## Abstract

A new methodology for estimating the depolarization ratio (DR) by dual-polarization radars with simultaneous transmission/reception of orthogonally polarized waves together with traditionally measured differential reflectivity *Z*_{DR}, correlation coefficient *ρ*_{hυ}, and differential phase Φ_{DP} in a single mode of operation is suggested. This depolarization ratio can serve as a proxy for circular depolarization ratio measured by radars with circular polarization. The suggested methodology implies the use of a high-power phase shifter to control the system differential phase on transmission and a special signal processing to eliminate the detrimental impact of differential phase on the estimate of DR. The feasibility of the suggested approach has been demonstrated by retrieving DR from the standard polarimetric variables and the raw in-phase *I* and quadrature *Q* components of radar signals and by implementing the scheme on a C-band radar with simultaneous transmission/reception of horizontally and vertically polarized waves. Possible practical implications of using DR include the detection of hail and the determination of its size above the melting layer, the discrimination between various habits of ice aloft, and the possible identification and quantification of riming, which is associated with the presence of supercooled cloud water. Some examples of these applications are presented.

## 1. Introduction

Historically, the first dual-polarization weather radars transmitted radio waves with either linear or circular polarizations and received both the copolar and cross-polar components of the reflected signals. The ratio of powers of these components represents the depolarization ratio. By definition, the linear depolarization ratio *L*_{dr} is the ratio of the power of a cross-polar radar return to the power of a copolar radar return if a wave with horizontal (or vertical) polarization is transmitted (e.g., Bringi and Chandrasekar 2001). Circular depolarization ratio *C*_{dr} is the ratio of the power of the copolar component of the reflected signal (e.g., left-hand circular) to its cross-polar component (e.g., right-hand circular) (McCormick and Hendry 1975). Depolarization ratios are usually represented in a logarithmic scale [LDR = 10 log_{10}(*L*_{dr}) and CDR = 10 log_{10}(*C*_{dr})]. In previous studies, depolarization ratios have been used primarily for the identification of hail and the melting layer wherein both *L*_{dr} and *C*_{dr} are much higher than in pure rain and dry snow (e.g., Minervin and Shupyatsky 1963; Barge 1974; Bringi et al. 1986; Holler et al. 1994; Hubbert et al. 1998; Kennedy et al. 2001). It was also shown by Matrosov et al. (2001, 2012) that the CDR dependency on antenna elevation angle can be used to distinguish between planar and columnar types of ice cloud hydrometeors and to estimate the hydrometeor degree of nonsphericity (i.e., the deviation of their axis ratios from 1).

There are several important differences between LDR and CDR. First, CDR depends primarily on the scatterers’ shape and only slightly on their orientations (in fact, CDR does not depend on the scatterer orientation in the incident wave propagation plane), whereas LDR depends both on the shape and (very strongly) on the orientation of hydrometeors. Second, CDR is generally larger than LDR (if the hydrometeors are nonspherical) and can be more reliably measured at low signal-to-noise ratio (SNR). For randomly oriented hydrometeors, CDR is 3 dB higher than LDR, and the difference between CDR and LDR increases for oriented nonspherical hydrometeors (Holt et al. 1999; Ryzhkov 2001). In rain, CDR can exceed LDR by more than 10 dB. On the negative side, CDR is much more sensitive to propagation effects than LDR. In other words, CDR strongly depends on differential phase Φ_{DP}; therefore, separating the contributions of backscattering and forward propagation effects to Φ_{DP} in interpreting CDR is a notoriously difficult task (e.g., Al-Jumily et al. 1991; Torlaschi and Holt 1993, 1998). This may limit the utilization of CDR close to the radar where Φ_{DP} is small.

Radar variables measured in linear and circular polarization bases are interrelated, and any single variable in a circular polarization basis, such as CDR, can be expressed via a combination of the variables measured in a linear horizontal–vertical (HV) polarization basis provided that all elements of the covariance scattering matrix in the HV basis are known along with the elements of the propagation matrix describing effects of differential attenuation and differential phase (e.g., Jameson 1987; Tragl 1990; Zrnić 1991; Krehbiel et al. 1996; Holt et al. 1999; Bringi and Chandrasekar 2001). Holt et al. (1999) attempted to estimate CDR using simultaneous and switched transmission of vertical and horizontal polarization implemented on the research CSU–CHILL polarimetric radar, which utilizes two transmitters and two receivers for orthogonal polarizations to obtain a full covariance scattering matrix (Brunkow et al. 2000).

Although circular polarized meteorological radars have been used in the past (e.g., McCormick and Hendry 1979; Kropfli et al. 1990; Krehbiel et al. 1996), modern operational dual-polarization radars increasingly use simultaneous transmission and reception of the waves of horizontal and vertical polarizations (SHV) with a single transmitter (Doviak et al. 2000), which allows measuring differential reflectivity *Z*_{DR}, differential phase Φ_{DP}, and copolar correlation coefficient *ρ*_{hυ} but not LDR (in the SHV mode of operation). Some polarimetric weather radars optionally measure LDR in a special LDR mode of operation; then only a horizontally polarized wave is transmitted, and both the horizontally and vertically polarized radar returns are received. This mode of operation, however, requires a polarization switch on transmission, and simultaneous measurements of *Z*_{DR} and LDR are not possible. Most users do not utilize the LDR mode because switching back and forth between the SHV and LDR modes slows down the radar data update and causes premature wear of the high-power switch. Chandrasekar and Bharadwaj (2009) suggested a technique to estimate LDR in the SHV mode using phase coding. According to this technique, the horizontally and vertically polarized transmit waveforms are coded with orthogonal phase sequences. This phase coding is not implemented on operational polarimetric weather radars, and such an opportunity has been explored only on the research CSU–CHILL radar.

How much information do we lose by abandoning depolarization measurements, or is it still possible to infer depolarization characteristics of atmospheric particles in the standard SHV mode of operations? Matrosov (2004) was the first to propose the idea of measuring a proxy for the circular depolarization ratio by an X-band radar operating in the SHV mode. This idea was echoed in more recent investigations by Melnikov and Matrosov (2013) and Ryzhkov et al. (2014). In this study, we further explore such an approach and demonstrate how CDR—or its proxy, depolarization ratio (DR)—can be estimated from operational dual-polarization weather radar observations and thus complement *Z*_{DR}, Φ_{DP}, and *ρ*_{hυ} data without slowing down or compromising the standard mode of operation. Moreover, we show that our method automatically eliminates the impact of propagation differential phase on DR (and thus CDR) estimates at the signal processor level. The best results are achieved if the system differential phase on transmission is controlled using a high-power phase shifter to ensure that the polarization state of the transmitted wave is close to circular.

## 2. Theoretical background

Typical signal processors of polarimetric radars estimate the powers of the received signals in the orthogonal channels *P*_{h} and *P*_{υ} and the complex covariance *R*_{hυ} from which the radar reflectivity, the differential reflectivity, the magnitude of correlation coefficient, and its phase (i.e., differential phase) are computed. If *V*_{h} and *V*_{υ} are complex voltages of radar returns at the output of the two receivers for backscattered waves with horizontal and vertical polarizations, then

where overbars indicate averaging over a number of radar samples, and the superscript * denotes the complex conjugate. In the SHV mode of operation whereby the horizontal (*h*) and vertical (*υ*) waves are transmitted and received simultaneously, the complex voltages *V*_{h} and *V*_{υ} are expressed as (Doviak et al. 2000; Matrosov 2004; Ryzhkov and Zrnić 2007)

where *S*_{ij} are the elements of the scattering matrix of hydrometeors in the horizontal–vertical linear polarization basis, *T*_{i} are the elements of the propagation matrix accounting for signal attenuation and propagation differential phase , and and are the system differential phases on transmission and reception, respectively. A zero-mean canting angle of the hydrometeors along a propagation path is assumed so that the transmission matrix has a diagonal form. The propagation terms are ensemble averaged along the propagation path as opposed to the elements of the scattering matrix , which are ensemble averaged in the radar resolution volume. According to Eq. (2), the voltages in the two orthogonal channels are

where *A*_{hυ} = |*T*_{h}/*T*_{υ}|^{2}, Φ_{DP} =1/2[arg(*T*_{υ}) *−* arg(*T*_{h})], and depends on the radar parameters and the distance from the radar.

In addition to the standard polarimetric variables routinely measured in the SHV mode, it is suggested that the ratio

can serve as a proxy for CDR (Melnikov and Matrosov 2013; Ryzhkov et al. 2014). Equation (5) can be rewritten as

where

is the measured differential reflectivity expressed in linear scale and

is the measured copolar correlation coefficient. It is important to emphasize that the values of *Z*_{dr} and *ρ*_{hυ} measured in the SHV mode are generally affected by cross coupling between the H and V polarizations and are not identical to the intrinsic *Z*_{dr} and *ρ*_{hυ} ( and ) defined through the elements *S*_{hh}, *S*_{υυ}, and *S*_{hυ} of the covariance scattering matrix as

In Eqs. (9) and (10), angular brackets mean averaging over an ensemble of hydrometeors. Additionally, the measured value of *Z*_{dr} is generally biased because of differential attenuation [factor *A*_{hυ} in Eqs. (3) and (4)]. As can be deduced from Eqs. (3) and (4) containing cross-polarization term *S*_{hυ}, the apparent values of *Z*_{dr} and *ρ*_{hυ} in Eqs. (7) and (8) depend to some extent on the linear depolarization ratio *L*_{dr}, and the magnitude of *ρ*_{hυ} is also dependent on the system differential phase upon transmission , that is, the polarization state of the transmitted wave. None of the existing operational dual-polarization radars have the capability to control the polarization state of the transmitted wave.

Using Eqs. (1), (3), and (4), the ratio *D*_{r} in Eq. (5) can be expressed via the intrinsic values of and as (Ryzhkov et al. 2014)

where

Equation (11) is obtained by assuming that the mean canting angle of hydrometeors is equal to zero (i.e., the moments 〈*S*_{hh} × *S*_{hυ}〉 and 〈*S*_{υυ} × *S*_{hυ}〉 vanish) and ≪ 1. This assumption, however, may not be valid in the areas of crystals canted in the presence of strong electric fields (Ryzhkov and Zrnić 2007).

Jameson (1987) showed that the circular depolarization ratio can be written in terms of linear polarization parameters as

Accounting for differential attenuation effects, we get

It is obvious from Eq. (15) that *C*_{dr} in Eq. (15) is strongly affected by differential phase (i.e., propagation effects), which is a major obstacle to the practical utilization of circular depolarization ratio. Comparing Eq. (11) for *D*_{r} and Eq. (15) for *C*_{dr}, it follows that *D*_{r} is very close to *C*_{dr} for small values of if = ±90°, that is, when the radar transmits a circularly polarized wave. In this case,

Moreover, Eqs. (15) and (16) become identical if = 0. If, however, the transmitted wave has a slant 45° linear polarization [i.e., = 0°], then

and the terms containing *L*_{dr} have a much smaller impact on *D*_{r} than in Eq. (16) provided that Φ_{DP} is relatively small so that sin(Φ_{DP}/2) ≈ 0. In that case, *D*_{r} might be quite different from *C*_{dr}.

The difference resulting from the transmitted wave being circularly or linearly polarized can be illustrated for randomly oriented hydrometeors. Then

and = 1 (Bringi and Chandrasekar 2001). It can be easily shown that *D*_{r} ≈ 2*L*_{dr} and *D*_{r} = *L*_{dr} for circularly and linearly polarized transmitted waves, respectively, provided that Φ_{DP} = 0. In the former case, *D*_{r} is equal to *C*_{dr}, and in the latter case, it is about 3 dB lower. A more detailed analysis for different types of hydrometeor orientation and various polarizations of the transmitted wave (generally elliptical) indicates that the highest values of DR are achieved for circular polarization of the transmitted wave although the variability of DR with does not exceed 3 dB (if the mean canting angle of hydrometeors is equal to 0°). One has to keep in mind that the impact of antenna cross coupling on the quality of the measurements of polarimetric radar variables is maximized if the polarization of the transmitted wave is circular (Hubbert and Bringi 2003; Wang and Chandrasekar 2006; Wang et al. 2006; Hubbert et al. 2010). The biases in *Z*_{DR}, *ρ*_{hυ}, and Φ_{DP} are also highest for circular polarization on transmission if the mean canting angle is different from zero (Ryzhkov and Zrnić 2007). However, if the antenna cross talk is small and the mean canting angle of hydrometeors is close to zero (which is generally a valid assumption except for low-inertia ice crystals in electrically active zones of clouds), then the use of circular polarization upon transmission (i.e., = ±90°) has an advantage in maximizing depolarization ratio DR. Anyway, the depolarization ratio *D*_{r} defined by Eq. (11) contains information useful for the discrimination between different types of ice and mixed-phase hydrometeors even if the polarization of the transmitted wave is not circular, as will be shown in the next sections.

It is important to note that *D*_{r} is much less affected by the propagation differential phase than *C*_{dr}, which is a great advantage of *D*_{r}. Both *D*_{r} and *C*_{dr} are, however, affected by another propagation effect, differential attenuation, which is described by the factor *A*_{hυ} in Eqs. (12) and (13). Therefore, the advantage of measuring *D*_{r} instead of *C*_{dr} is greater at longer radar wavelengths where attenuation/differential attenuation in rain is smaller.

## 3. Simulations based on experimental data

Our theoretical analysis demonstrates that the depolarization ratio, DR = 10 log_{10}(*D*_{r}), which serves as a good proxy for the circular depolarization ratio CDR, can be estimated by radars with simultaneous transmission/reception along with the standard set of polarimetric variables *Z*_{DR}, *ρ*_{hυ}, and Φ_{DP}. The best correspondence between DR and the “true” CDR is achieved if the radar transmits the wave with circular polarization, which means that the differential phase upon transmission is equal to ±90°. The implementation of DR measurements requires the computation following Eq. (5) and the control of the polarization of the transmitted wave using a high-power phase shifter in one of the orthogonal transmission channels to maximize DR and make it close to CDR if needed. According to Eq. (15), the circular depolarization ratio CDR can be computed from the measurements of *Z*_{DR}, *ρ*_{hυ}, and Φ_{DP} only if LDR is also available (which is not the case for operational polarimetric weather radars). For this purpose, we resort to the data collected by the National Center for Atmospheric Research (NCAR) research S-band dual-polarization Doppler radar (S-Pol), which measured all needed variables. A thunderstorm in Florida on 14 August 1998 is selected for the following analysis, a case already examined in the paper by Ryzhkov et al. (2002).

An example of vertical cross sections of the measured *Z*, , LDR, Φ_{DP}, and the estimates of two depolarization ratios is displayed in Fig. 1 for the azimuth angle 172.2°. The first depolarization ratio (tagged as CDR) is obtained from Eq. (15) and represents what would be measured by a true circularly polarized radar. The impact of propagation (or Φ_{DP}) on CDR is obvious: the intrinsic CDR is grossly overestimated in areas with even modest Φ_{DP} values (up to 10°–20°). Such an overestimation is expected to be much more dramatic at shorter radar wavelengths where Φ_{DP} is larger. Different approaches for correcting CDR for differential phase shift were discussed in the literature (e.g., Torlaschi and Holt 1993, 1998), but none proved to be efficient.

The second estimate of DR is computed from Eq. (16) for *D*_{r}, which is obtained by assuming that = π/2 in Eq. (11). This depolarization ratio is not affected by propagation while being very consistent with intrinsic CDR (in the absence of the propagation impact) and results of theoretical simulations, which can be found in literature (Jameson 1987; Al-Jumily et al. 1991; Torlaschi and Holt 1993, 1998; Holt et al. 1999; Matrosov et al. 2001, 2012). The DR seems to be more informative for hydrometeors above the melting layer than the traditionally utilized *Z*_{DR} and *ρ*_{hυ} and thus may complement other polarimetric measurements. The microphysical interpretation of DR is less ambiguous than deciphering LDR because the latter depends on particle shape, density, and orientation, whereas DR serves as a proxy for true CDR, which is only weakly affected by orientation. For example, DR allows distinguishing between the updraft associated with the *Z*_{DR} column centered at *R* = 20 km and the nearby downdraft filled with nearly spherical dry graupel particles centered at *R* = 23 km above the freezing level, whereas LDR does not. This is mainly because raindrops lofted in a convective updraft are relatively well oriented, which reduces LDR despite their quite nonspherical shape. The DR exhibits a very pronounced signature in the melting layer and offers a valuable resource for polarimetric microphysical retrievals.

Vertical cross sections of *Z*, LDR, DR, and Φ_{DP} for the same case but for a different azimuth (184.2°) in a more stratiform part of the storm are displayed in Fig. 2. Again, as in Fig. 1, and above the melting layer do not reveal pronounced patterns or signatures, whereas both LDR and DR exhibit a coherent structure in the ice part of the storm. Obviously, DR correlates well with LDR while being a few dB higher.

The value of DR displayed in Figs. 1 and 2 is computed from Eq. (11) assuming that the polarization of the transmitted wave is circular [ = 90°]. We examined the influence of on the estimated DR and compare the DR cross sections with changing from 90° to 20° (Fig. 3) with a 10° increment. It is evident that the highest DR corresponds to = 90° (i.e., the transmitted wave is circularly polarized) and gradually decreases with approaching 0°.

Vertical profiles of DR, LDR, and *ρ*_{hυ} through the melting layer at *Az* = 184.2° and a distance of 18 km from the radar are displayed in Fig. 4. The quantities LDR and DR are well correlated above the melting layer. Within the melting layer, the DR’s pronounced maximum coincides with the *ρ*_{hυ} minimum, whereas LDR increases monotonically with height. As expected, DR in rain is much higher than LDR. The difference between ρ_{hυ} in pure rain below the melting layer and frozen hydrometeors above it is vanishingly small, whereas DR in rain is very different from DR above the melting layer.

Ryzhkov (2001) showed that the difference between CDR and LDR depends primarily on particle orientations, so joint measurements of different depolarization ratios can be used to infer hydrometeor fall attitudes (e.g., Matrosov et al. 2005). If the hydrometeors are oriented randomly, then the difference CDR − LDR is roughly 3 dB and does not depend on the radar elevation angle as mentioned previously and is visible in many regions above the melting layer in the figures shown. For relatively narrow distributions of the canting angles,

where *s*_{h} and *s*_{υ} are the scattering amplitudes in orthogonal directions, 〈*α*〉 is the mean canting angle, and *σ* is the width of the canting angle distribution. In rain, 〈*α*〉 ≈ 0, and the factor containing *s*_{h} and *s*_{υ} varies between 0.25 and 0.35; hence, the ratio *C*_{dr}/*L*_{dr} (or the difference CDR − LDR expressed in dB) is a measure of *σ*. In the examples shown in Figs. 1 and 2, DR − LDR varies between 10 and 15 dB in rain below the melting layer. This corresponds to *σ* varying between 5° and 10° in full agreement with the disdrometer studies (Huang et al. 2008). This is another evidence indicating that our procedure yields realistic estimates of CDR using DR as its proxy.

In general, DR positively correlates with *Z*_{DR}, and such correlation is much stronger in rain than in snow. This is illustrated in Fig. 5, where the scatterplots of DR versus *Z*_{DR} in rain below the melting layer (blue dots) and ice/snow above the melting layer (red dots) estimated for two values of (π/2 and 0) are displayed for the S-Pol case. The DR in snow/ice is noticeably higher than in rain for a given *Z*_{DR}. Another important conclusion from Fig. 5 is that using circular polarization on transmission [ = π/2] instead of slant linear polarization [ = 0] results in a higher DR and a better discrimination between rain and ice in the *Z*_{DR}–DR plane. This emphasizes the value of the phase shifter to control the polarization of the transmitted wave in order to ensure maximal values of DR.

In the previous examples, CDR and its proxy DR were computed from the radar variables (i.e., , , LDR, and Φ_{DP}). In the next example, CDR and DR are estimated directly from the in-phase and quadrature components *I* and *Q* of the radar signals (Doviak and Zrnić 1993) measured by the polarimetric S-band KOUN WSR-88D in Norman, Oklahoma, following Eq. (5). Figure 6 shows composite range–height indicators (RHIs) of *Z*, *Z*_{DR}, Φ_{DP}, LDR, CDR, and DR for the storm observed on 3 July 2007. The quantities *Z*, *Z*_{DR}, Φ_{DP}, CDR, and DR are estimated in the SHV mode of operation, while LDR is measured in the LDR mode whereby only the H-polarized wave is transmitted. The system differential phase on transmission was about 82° at the time of measurements. This means that the polarization of the transmitted wave in the SHV mode was very close to circular. Again, DR estimated from Eq. (11) does not exhibit a bias attributed to propagation differential phase.

## 4. Direct measurements of DR with a prototype of an operational radar

The suggested methodology for DR measurements prescribing the use of a high-power phase shifter, and a signal processing according to Eq. (5) has been implemented by the Enterprise Electronics Corporation on its C-band dual-polarization radar with simultaneous transmission/reception in Enterprise, Alabama. The phase shifter controlling the differential phase upon transmission was calibrated using alternate measurements of differential phase in the SHV and LDR modes in rain at ranges close to the radar. The estimated system differential phase in the SHV mode of operation is equal to the sum of and , whereas it is equal to in the LDR mode. Subtracting the two yields the estimate of . The shifter was tuned to achieve = 90°.

The quality of DR measurements can be evaluated using the expected consistency between DR and *Z*_{DR} in pure rain. Figure 7 (left panel) shows the scatterplot of DR versus in rain simulated from a large dataset of disdrometer measurements in central Oklahoma (Schuur et al. 2005). The DR was computed using Eq. (16) with Φ_{DP} = 0 from simulated , , and LDR at C band, assuming that the width of the canting angle distribution is 10°. The scatterplot of DR versus *Z*_{DR} obtained from measured C-band data collected in stratiform rain in Alabama on 25 February 2015 is presented in the top right panel in Fig. 7. It is consistent with the scatterplot of the simulated values in the left panel (at least for *Z*_{DR} > 0.5 dB). Observed DR for *Z*_{DR} < 0.5 dB corresponding to the areas dominated by small quasi-spherical raindrops are higher than the theoretical ones. This can be attributed to antenna cross coupling in the SHV mode. A possible impact of noise can also not be excluded, although the noises in the two orthogonal channels were subtracted from the measured powers before using Eq. (5). The scatterplot of DR versus *Z*_{DR} in ice above the melting layer for the same event exhibits generally higher values of DR than in the rain below the melting layer for a given *Z* in full agreement with the estimates from the S-Pol measurements displayed in Fig. 5. The fact that raindrops are well oriented with the width of the canting angle distribution within a 5°–10° range makes DR and *Z*_{DR} tightly correlated in rain. This is not the case in snow and ice, where the two variables are less correlated because of a more random orientation of snowflakes and ice crystals. Indeed, dry snow aggregates have nonspherical shape, but their *Z*_{DR} is close to zero because of their low density.

Average vertical profiles of *Z*, *Z*_{DR}, *ρ*_{hυ}, and DR obtained via azimuthal averaging of these variables at elevation 15° for the above case are shown in Fig. 8. The depolarization ratio has a pronounced maximum within the melting layer with DR approaching −16 dB. Additionally, DR exhibits a small bump at about 6-km height within the dendritic growth layer in the temperature interval between −10° and −20°C, which is also marked by an increase in *Z*_{DR} (e.g., Kennedy and Rutledge 2011; Williams et al. 2015). Overall, the depolarization ratio for this winter stratiform event is lower than DR estimated from the S-Pol for a convective summer storm at all heights, as the comparison of Figs. 4 and 8 demonstrates. Further scrutiny is required to explain these differences.

In deep convective storms, the linear depolarization ratio LDR has been traditionally used for hail detection above the freezing level (e.g., Barge 1974; Bringi et al. 1986; Holler et al. 1994; Carey and Rutledge 1998; Hubbert et al. 1998; Kennedy et al. 2001). Kennedy et al. (2001) suggest an S-band LDR threshold of −25 dB at mid- to upper levels of a thunderstorm for the detection of growing hail. Hail growth in the convective updrafts is often seen as an “LDR cap” on the top of the *Z*_{DR} column, which is aligned with the base of such updrafts (Jameson et al. 1996; Bringi et al. 1997; Hubbert et al. 1998; Carey and Rutledge 1998; Kumjian et al. 2014).

The depolarization ratio DR was estimated with the Enterprise Electronics Corporation (EEC) C-band radar data using Eq. (5) in a hail-bearing storm in Alabama on 7 August 2015. Figure 9 displays composite RHI and plan position indicator (PPI) plots of *Z*, *Z*_{DR}, *ρ*_{hυ}, and DR for a convective cell containing hail aloft. The RHI scans were reconstructed from the series of successive PPI sweeps; thus, their vertical resolution is somewhat compromised. Nevertheless, the convective cell with maximal reflectivity of 58 dB*Z* is marked by a well-pronounced *Z*_{DR} column and a much taller column of DR with a maximal value as high as −10 dB, which is indicative of wet, large hail growing in the convective updraft. Obviously, the DR column can better indicate a strong convective updraft than the *Z*_{DR} column. Indeed, the DR column seems to combine the base of updraft where *Z*_{DR} is high because of the presence of large liquid or partially frozen raindrops and the LDR cap in the middle or at the top of updraft where water-coated or spongy hail is growing in a wet growth regime. Therefore, DR can potentially be used for hail detection aloft and for the estimation of its size (e.g., Mirkovic 2016). The composite PPIs of *Z*, *Z*_{DR}, and DR at different antenna elevations (and thus heights) in Fig. 10 demonstrate that the updraft signatures in terms of *Z*_{DR} disappear at 7 km, whereas the DR column is still clearly detectable up to 9 km.

## 5. Potential practical utilization of DR

Three major challenging practical tasks can be addressed using DR measurements: 1) detection of hail and determination of its size above the melting layer, 2) differentiating between various ice habits aloft, and 3) quantification of riming associated with the presence of supercooled cloud water hinting at potential icing hazard for airplanes. These are addressed next.

### a. Detection of hail and estimation of its size

The hail size discrimination algorithm (HSDA) recently developed at NSSL for implementation on the WSR-88D network distinguishes between small (*D* < 2.5 cm), large (2.5 < *D* < 5.0 cm), and giant (*D* > 5.0 cm) hail (Ryzhkov et al. 2013b; Ortega et al. 2016). The algorithm capitalizes primarily on polarimetric signatures below the melting layer, where melting hail is mixed with rain. This implies that hail has already been formed aloft and falls to the ground. For hail nowcasting and suppression, however, it is important to detect large hail earlier, when it is just formed in the upper part of the storm where radar reflectivity *Z* is a primary discrimination parameter in HSDA and where differential reflectivity *Z*_{DR} is usually close to zero. In contrast, both LDR and CDR (or DR) above the freezing level vary significantly depending on hail size and shape. Simulations of CDR based on the microphysical model of hail described in Ryzhkov et al. (2013a) and Mirkovic (2016) show that CDR can increase by about 10 dB if the maximal diameter of dry hail changes from 8 to 50 mm. In computations, it was assumed that the width of the mean canting angle distribution is equal to 40°, and the aspect ratio of hailstones is 0.8 if their size exceeds 1 cm and varies between 0.8 and 1.0 for smaller sizes as specified in Ryzhkov et al. (2011).

The steady increase of CDR with the maximal hail size is caused by progressively stronger resonance scattering effects for larger hailstones and by the decrease of the slope of the exponential size distribution in case of larger hail (Aydin and Zhao 1990; Ryzhkov et al. 2013a). Very big hail with diameters larger than 5 cm usually grows in a wet growth regime and may have quite significant water fraction even at negative temperatures while residing in a convective updraft. The magnitude of dielectric constant for wet hail is higher than for dry hail, which results in higher *Z*_{DR}, LDR, CDR, and lower *ρ*_{hυ}.

The utilization of the DR column discussed in the previous section offers an additional opportunity to identify large or giant hail growing in a wet growth regime within the convective updraft. As was mentioned in section 2, the depolarization ratio DR best approximates the true circular depolarization ratio CDR (unbiased by the impact of differential phase) if the polarization of the transmitted wave is circular. Then the value of DR and its discriminatory power are highest (see Fig. 3). However, we believe that the depolarization ratio derived from SHV radar measurements can be also useful for arbitrary polarization of the transmitted wave, which is the case for the majority of existing operational and research dual-polarization radars. Moreover, the depolarization ratio computed from the measured *Z*_{DR} and *ρ*_{hυ} using Eq. (6) may provide new insights into the microphysical properties of frozen and mixed-phase hydrometeors because it combines the information contained in *Z*_{DR} and *ρ*_{hυ} in a single signature or feature. To demonstrate this, we examined several weather events observed with dual-polarization radars for which the differential phase upon transmission was not known.

A hailstorm producing large hail in central Oklahoma on 17 June 2005 was observed by the KOUN S-band radar, which served as a prototype of the operational dual-polarization WSR-88D sets. The radar antenna was scanning in the vertical plane at a number of azimuthal directions within a narrow azimuthal sector so that high-resolution RHIs have been produced. Examples of the composite RHIs of *Z*, *Z*_{DR}, *ρ*_{hυ}, and DR computed from *Z*_{DR} and *ρ*_{hυ} are shown in Figs. 11 and 12.

A classical signature of melting hail falling toward the ground is displayed in Fig. 11. It is centered at range 75 km from the radar and marked by high *Z* (exceeding 60 dB*Z* next to the surface), low *Z*_{DR}, and depressed *ρ*_{hυ}. A “*Z*_{DR} hole” is clearly visible at this range. Of primary interest, however, is the signature of a new convective updraft indicated by a white oval, centered 68 km from the radar. The updraft can be identified by a relatively short *Z*_{DR} column stretching 1.5 km above the environmental freezing level at 4 km. This updraft is also characterized by a much taller area of *ρ*_{hυ} depression, which is a proxy for an LDR cap. These *Z*_{DR} and *ρ*_{hυ} signatures are combined in a strong DR column, which is very similar to the one illustrated in Fig. 9. Note that radar reflectivity is very modest in this area, indicating that hailstones growing in the upper part of this updraft are still small.

Figure 12 shows the signature of very large hail observed in the same storm but 40 min later. The hail core centered 35 km from the radar is characterized by very high values of *Z* exceeding 65 dB*Z*, very low *ρ*_{hυ} (<0.8), and negative *Z*_{DR} (as low as −1.3 dB) aloft at the height exceeding 3 km AGL. Negative values of *Z*_{DR} definitely signal the presence of giant hailstones with diameters larger than 5 cm as a result of resonance scattering, as shown by Ryzhkov et al. (2013a). The combination of anomalously low *Z*_{DR} and *ρ*_{hυ} yields very high values of DR exceeding −10 dB. The examples in Figs. 11 and 12 demonstrate that combining the *Z*_{DR} and *ρ*_{hυ} signatures into a single feature, that is, depolarization ratio, may help to better identify large hail aloft and complement traditional reflectivity-based criteria of the vertically integrated liquid (VIL) or maximum expected size of hail (MESH) (Witt et al. 1998).

### b. Identification of different ice particle habits

Matrosov et al. (2001) and Matrosov (2015) showed that the dependence of the absolute values of CDR [(and/or closely related to it slant 45° linear depolarization ratio (SLDR)] on the antenna elevation angle can be used for the discrimination between different ice hydrometeor habits and the estimation of particle shape parameters (e.g., aspect ratios). The CDR increases with decreasing antenna elevation for planar (oblate) types of snow crystals (hexagonal plates, thick plates, dendrites), whereas the elevation dependency of CDR is relatively “flat” for columnar (prolate) type crystals (columns, needles, etc.). At low elevation angles, CDR varies in a large range (from the radar system depolarization limit to about −8 dB) depending on the ice habit. According to cloud radar observations by Matrosov et al. (2001, 2012), pristine dendrites and hexagonal plates have highest values of CDR up to −8 dB (at slant viewing) followed by columns and aggregates of dendrites (from −23 to −16 dB) with graupel indicating the lowest CDR below −25 dB. Note that the specific values of actually measured depolarization ratios are affected (especially for lower depolarizations) by the particular radar antenna cross-coupling characteristics. Since the suggested DR parameter represents a CDR proxy, it can be useful for future studies of ice hydrometeor habits.

Another potential application may involve the discrimination between aggregated and rimed snow and possible quantification of the degree of riming. Riming denotes the freezing of small supercooled liquid droplets (with sizes from microns to tens of microns) on faster falling ice crystals and snowflakes. Supercooled liquid water cannot be observed directly by weather surveillance radars because of the small size of liquid droplets, but its presence can be detected indirectly by estimating the degree of riming of ice crystals. Riming tends to increase the density and aspect ratio of ice particles, that is, it makes them denser and more spherical. The density effect increases CDR while the shape effect changes CDR in the opposite direction. Both theoretical and experimental studies backed by in situ particle observations indicate that the shape effect usually dominates so that riming reduces CDR at low elevation angles (e.g., Matrosov et al. 2001, 2012).

### c. Quasi-vertical profiles of depolarization ratio

Another possible way to discriminate between aggregated and rimed snow implies the use of the concept of quasi-vertical profiles (QVPs) of radar variables including the depolarization ratio DR. Such a concept has been recently introduced by Trömel et al. (2013) and Ryzhkov et al. (2016). It suggests azimuthal averaging of the radar data on conical scans at high antenna elevations (between 10° and 30°) to reduce the statistical uncertainty of the radar variables estimates. Trömel et al. (2013) demonstrated the QVP strategy as a pathway to reliably estimate the backscatter differential phase in the melting layer. The QVP profiles are represented in a height versus time format, and the evolution of the microphysical processes of snow formation above the freezing level can be captured in relative detail. Since the statistical errors of *Z*_{DR} and *ρ*_{hυ} in the QVP profiles are very small, Eq. (6) can be effectively used for the computation of QVPs of DR with high accuracy for most dual-polarization radars.

An example of the composite QVP of *Z*, *Z*_{DR}, *ρ*_{hυ}, and DR estimated from the KMOB WSR-88D data using Eq. (6) for a winter stratiform event near Mobile, Alabama, is presented in Fig. 13. The QVP methodology allows delineating the melting layer and the dendritic growth layer usually found in the temperature interval between −10° and −15°C (or 5–6 km in height) and their temporal evolution with unprecedented accuracy and vertical resolution. The depolarization ratio reaches −15 dB within the melting layer and −19 dB in the dendritic growth layer. The latter is also denoted by strong *Z*_{DR} enhancement and a reduction of *ρ*_{hυ} (Williams et al. 2015).

In Fig. 13, DR follows the pattern of Z_{DR}, which is not surprising from the physical standpoint and because DR is computed from *Z*_{DR} and *ρ*_{hυ}. Nevertheless, there are certain differences due to the dependence of DR on *ρ*_{hυ} as well. How close is DR computed from Eq. (6) to the real circular depolarization ratio CDR? This depends on the system differential phase upon transmission , which changes from radar to radar on the WSR-88D fleet and generally is not known. Regardless of the difference between DR and true CDR, DR turns out to be a useful radar parameter characterizing microphysical properties of snow. In general, the range of values of DR in Fig. 13 is consistent with what was reported by Matrosov et al. (2001, 2012), who measured SLDR, which is practically equivalent to CDR for axisymmetric particles predominantly oriented with their major dimensions in the horizontal plane.

Similar to *Z*_{DR}, DR tends to decrease in rimed snow relative to aggregated snow. Vogel et al. (2015), Ryzhkov et al. (2016), Kumjian et al. (2016), and Giangrande et al. (2016) claim that *Z*_{DR} of rimed snow is slightly lower than *Z*_{DR} of aggregated snow, but this difference is really small (typically 0.2–0.4 dB), which can be clearly visible in the QVP profiles. The corresponding difference in DR is significantly larger (at least 2–4 dB) as can be seen in Fig. 13. Both Z_{DR} and DR exhibit their lowest values in rimed snow just above the melting layer, which is apparently “sagging” during intense riming periods between 1530 and 1630 UTC and between 2100 and 2300 UTC [see also Kumjian et al. (2016) and Xie et al. (2016)]. It is likely that the relative differences of DR between various types of snow/ice may not be much affected by the uncertainty in , at least in a qualitative sense. We showed that these differences are maximized if the radar transmits the wave with truly circular polarization, which motivates the utilization of a phase shifter for controlling . This, however, does not prevent any dual-polarization radar with simultaneous transmission/reception and arbitrary to use DR for snow classification.

Similar computations of DR from the quasi-vertical profiles of *Z*_{DR} and *ρ*_{hυ} using Eq. (6) have been made for a number of storms observed with the EEC X-band radar owned by the University of Bonn, Germany (e.g., Diederich et al. 2015). Typical results are illustrated in Fig. 14. Again, the patterns of DR in its QVP profiles at X band are qualitatively and quantitatively similar to the ones obtained from the WSR-88D measurements although the polarization state of the wave transmitted by the X-band radar is not known. The advantage of using DR instead of *Z*_{DR} for discrimination between rimed and aggregated snow is clearly shown in Fig. 14. Indeed, DR exhibits large differences in snow between the first 30 volume scans (green shades of DR) and the subsequent 40 scans (mainly blue shades of DR). This means that snow was likely more rimed during the later period of observations, which implies a higher supercooled liquid water content favoring more extensive riming of snow.

## 6. Conclusions

A new methodology for estimating depolarization ratio by dual-polarization radars with simultaneous transmission/reception along with the traditionally measured *Z*_{DR}, *ρ*_{hυ}, and Φ_{DP} in a single mode of operation has been suggested. This depolarization ratio can serve as a proxy for circular depolarization ratio measured by radars with circular polarization. The suggested methodology implies the use of a high-power phase shifter to control the system differential phase on transmission and a special signal processing to eliminate the detrimental impact of propagation differential phase on the estimate of DR. In the past, the bias related to the differential phase on propagation was a principal obstacle to the use of the directly measured CDR and its earlier estimates from SHV mode measurements (Matrosov 2004). An important advantage of DR (i.e., a CDR proxy) suggested here is also that they are available in the same regions of the radar echo where *Z*, *Z*_{DR}, *ρ*_{hυ}, and Φ_{DP} are measured, whereas true CDR measurements can be reliably made only in the areas where the signal-to-noise ratio is significant (i.e., >30 dB) since the echo power in the “weak” polarization channel is significantly smaller than the one in the “strong” polarization channel.

The feasibility of the recommended approach has been demonstrated by retrieving DR from other polarimetric variables, the raw *I*–*Q* data, as well as by implementing the scheme on a C-band radar with simultaneous transmission/reception of the horizontally and vertically polarized waves. It is also shown that the proxy of DR (or CDR) can be roughly estimated from the combination of *Z*_{DR} and *ρ*_{hυ} measured by radars with arbitrary polarization of the transmitted wave. As opposed to linear depolarization ratio (LDR) (which requires a special mode of operation), DR depends on hydrometeor orientation rather weakly and is less affected by noise.

The suggested approach for estimating depolarization ratios can be implemented with operational weather radars (e.g., the WSR-88D network) and numerous research radars, for example, cloud and precipitation radars operated by the U.S. Department of Energy’s Atmospheric Radiation Measurement (DOE ARM) program. The addition of depolarization ratio estimates will improve the utility of the radar measurements to infer microphysical information.

Among potential applications, three major challenging practical tasks can be addressed using DR measurements: 1) detection of hail and determination of its size above the melting layer, 2) discrimination between various habits of ice aloft, and 3) identification and quantification of riming processes associated with the presence of supercooled cloud water and possible icing hazard to aviation.

It is demonstrated that strong convective updrafts where large hail is growing can possibly be better identified as DR columns encompassing most of the updraft as opposed to *Z*_{DR} columns indicative only of its bottom part. The DR columns are taller because they combine Z_{DR} column and “LDR cap” or *ρ*_{hυ} reduction in a single feature, which can be better detected at higher antenna elevations.

The depolarization ratio proxy estimates can potentially be utilized for the discrimination between different snow types similar to directly measured depolarization ratios. These ratios seem to be noticeably lower in rimed snow than in aggregated snow, and there is a possibility for utilizing DR for the detection of supercooled cloud water, which is an essential ingredient for snow riming.

Meaningful DR values can be estimated from the combination of *Z*_{DR} and *ρ*_{hυ} measured by polarimetric radars with arbitrary polarization of the transmitted wave for which the differential phase upon transmission is not controlled. It is likely that relative differences of DR between various types of snow/ice may not be much affected by the uncertainty in , at least in a qualitative sense. We assume that these differences are maximized if the radar transmits the wave with truly circular polarization, which suggests the utilization of a phase shifter for controlling but does not prevent any dual-polarization radar with simultaneous transmission/reception and arbitrary to use DR for classification of snow.

## Acknowledgments

Funding for some of the authors was provided by NOAA/Office of Oceanic and Atmospheric Research under NOAA–University of Oklahoma Cooperative Agreement NA11OAR4320072, U.S. Department of Commerce. Sergey Matrosov was sponsored in part by the U.S. Department of Energy Atmospheric System Research Program under award number DE-SC0013306. Research based on the Bonn X-band radar was carried out in the framework of the Hans-Ertel-Centre for Weather Research (http://www.herz-tb1.uni-bonn.de/) funded by the BMVBS (Federal Ministry of Transport, Building and Urban Development) and the SFB/TR 32 (Transregional Collaborative Research Centre 32, http://www.tr32.de/) funded by the DFG (German Research Foundation).

## REFERENCES

*Polarimetric Doppler Weather Radar: Principles and Applications.*Cambridge University Press, 636 pp.

*Doppler Radar and Weather Observations.*2nd ed. Dover Publications, Inc., 562 pp.

*29th Int. Conf. on Radar Meteorology*, Montreal, QC, Canada, 214–217.

*Conf. on Cloud Physics*, San Francisco, CA, Amer. Meteor. Soc., 341–345.

_{DR}columns: Investigating a polarimetric radar signature with a spectral bin microphysical model

*36th Conf. on Radar Meteorology*, Breckenridge, CO, Amer. Meteor. Soc., 245. [Available online at https://ams.confex.com/ams/36Radar/webprogram/Paper228291.html.]

*Computational Electromagnetics for Polarimetric Radar Scatterers: An Approach to Polarimetric Variable Modeling of Hydrometeors and Biota.*Lambert Academic Publishing, 272 pp.

*Eighth European Conf. on Radar in Meteorology and Hydrology*, Garmisch-Partenkirchen, Germany, ERAD, 11.3.

*32nd Conf. on Radar Meteorology*, Albequerque, NM, Amer. Meteor. Soc., 15R.4. [Available online at https://ams.confex.com/ams/32Rad11Meso/webprogram/Paper95995.html.]

*37th Conf. on Radar Meteorology*, Norman, OK, Amer. Meteor. Soc., 6B.6. [Available online at https://ams.confex.com/ams/37RADAR/webprogram/Paper275246.html.]

## Footnotes

© 2017 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).