Polarization I: Radar concepts and ZDR part I

1. Polarization I:Radar concepts and ZDR part I

2. Dual polarization radars can estimate several return signal properties beyond those available from conventional, single polarization Doppler systems. Consider a dual linear polarization coherent radar in which both transmission and reception are possible in the horizontal (H) and vertical (V) polarization states. (Polarization defined by plane in which electric field lies). Some useful quantities that such a radar can measure are: Ratio of the H and V signal powers (ZDR) Phase difference between the H and V returns (fdp) Degree of correlation between the H and V returns (rhv) Ratio of orthogonal to “on channel” signal power (LDR)

3. Seliga and Bringi, 1976 J. Appl. Meteor.

5. First measurements of Zdr made by the CHILL radar in Oklahoma, 1977. Alternating HHH, VVV, HHH, VVV……….

6. Polarization of the electric field transmitted by the radar (and incident on the scatterers) is imposed by internal microwave signal paths in the antenna feed horn.

7. Backscattered electric field from an individual scatterer is described by the scattering matrix. “S” values are complex numbers that depend on the scatterer shape, orientation and dielectric constant Here, subscripts are transmit, receive from the particle viewpoint Largest terms are “co-polar” (repeated subscript) matrix elements Incident field due to transmitted radar pulse Backscattered electric field; contains both H and V components Matrix scanned from Bringi et al. 1986 part I

8. Scattering matrix multiplication result • Ebv = SvvEincv + SvhEinch • Ebh = ShvEincv + ShhEinch • Svhmeans particle scatters in v due to illumination in h

9. Computation of scattering matrix elements is simplified by considering particle shapes to be spheroids From appendix 1 of Bringi and Chandra (2001) text

10. Polarimetric variable #1: Ratio of the co-polar H and V return signal powers: Differential reflectivity (Zdr) Differential reflectivity ratio: linear scale power ratio as defined in Bringi and Chandra (2001); dependence on axis ratio (r) and diameter (D) explicitly shown

11. Co-polar scattering matrix terms under Rayleigh-Gans conditions. V = particle volume; er = relative permittivity; spheroid axis lengths = a,b; lz = depolarizing factor Recall that n=√εr Depends on composition and shape! Bringi and Chandra (2001), eq 7.5b,c and appendix 1

12. Basic Shh/Svv ratios (white text) for specified oblate scatterers (er from BC2001; water = 81; ice = 3.5) Shh and Svv calculated from previous equations. Zdr calculated accordingly. Main point: Zdr sensitivity to aspect ratio decreases as particle bulk density (as expressed by er relative permittivity) decreases

13. Single particle Zdr expressed as dB, curves labeled by er Zdr (dB)=10 log10 ((Shh / Svv)2) Plot from Herzegh and Jameson (1992)

14. Raindrop Characteristics and Zdr Equilibrium drop shape due to balance of surface tension and aerodynamic pressure distribution around the drop

15. Beard and Chuang drop shapes From 2 mm to 6 mm in 0.5 mm steps

16. Size-shape relationships 1 < D < 9 mm; Pruppacher and Beard (1970); wind tunnel measurementsmm b/a = 1.0048 + 5.7x10-4(D) – 2.628x10-2(D)2 + 3.682x10-3(D)3 - 1.677x10-4 (D)4 0 < D < 7 mm Beard and Chuang (1987); polynomial fit to numerical simulations Equations agree for D> 4 mm PB relation gives slightly more oblate drops compared to poly fit for D < 4 mm

17. Rayleigh-Gans Zdr from equilibrium (single) drop shape Plot from Herzegh and Jameson (1992)

18. Drop oscillations occur as diameter exceeds ~1 mm Image of modeled drop oscillations from K. Beard UIUC

19. Vortices shed in drop wake flow can help induce / sustain drop oscillations Saylor and Jones, Physics of Fluids (2005)

20. On average, observed raindrop shapes are somewhat less oblate than equilibrium force balance would dictate. Turbulent air motions and drop collisions also broaden the observed axis ratio range at a given diameter Lab study Range of size/shape relationships b/a = 1.012 – 0.01445D – 0.01028 (D)2 1 < D < 4 mm Fit to lab data Andsager et al., (1999; laboratory study using 25m fall column)

21. Basic exponential DSD: N(D)=N0*e-LD (D=diameter; L is the slope parameter). Here, N0’s have been adjusted to give the same reflectivity. Zdr is the reflectivity factor-weighted mean axis ratio of the drop size distribution. Drop population shift towards smaller diameters on the right is revealed by lower Zdr. (Note: rain rate estimation based on Z alone would be the same for these two DSD’s.) (See Andsager et al, JAS, 1999 equations 3 and 4)

22. Integral quantities from the rain DSD: D0 = median drop diameter; divides total water content in half Dm = mass weighted mean drop diameter Mean fit equations for the above two quantities as f(Zdr) based on thunderstorm-type DSD’s (Bringi and Chandra 2001, 7.14):

23. Verification of Zdr-estimated and observed Dm (Bringi and Chandra (2001)

24. Courtesy: Kristen George

25. 8 Sept 2002 (tropical) 19 Sept 2002 (trailing squall line) Rainfall Measurement with Polarimetric -88D (JPOL 2003)

26. Zdr typically used to adjust basic reflectivity-based rain rate estimator for variations in D0, etc. Low Zdr in rain -> small D0 -> reflectivity estimator low; Use of Zdr in the denominator will adjust Z rain rate up Ryzhkov et al. (2005)

27. B. Amer. Met. Soc., February, 1999

28. Ft. Collins flood: 29 July 1997 0208 UTC CSU-CHILL 2 km MSL CAPPI Using the NEXRAD relationship R = 0.017Z 0.714 This relationship is normally truncated at 55 dBZ to avoid hail contamination Using the relationship R = cZa 10 0.1bZdr Where a = 0.93 b= -3.43 c= 6.7 x 10-3

29. R(Z) rain rates low relative to R(Z, Zdr) in tropical (small D0) rain A sample from the Ft. Collins flood……….

30. Hail signal, Z and Zdr. V hail rain H Zdr = 10 log10 (Zhh/Zvv)

31. Insects are typically more oblate than raindrops, giving highly positive Zdr values. Zdr magnitudes and spatial textures are useful in identifying non-meteorological targets.

32. Difference Reflectivity, Z dp For a rain-ice mixture, we can write: Ice is isotropic so there is no polarization dependence. Therefore we can form the difference reflectivity as: Difference reflectivity only sensitive to oblate drops in the mixture.

33. When ice is present, the measured Zh will exceed Zh for rain. Whereas Zdp will be approximately the same since it is insensitive to ice. So the data point formed by Zh and Zdp will lie to the right of the rain line. Bringi and Chandrasekar (2001) - The distance between the vertical red lines denotes ΔZ, that is, Zh measured minus Zh,rain. Here f is the so called ice fraction, the ratio of the ice reflectivity factor to the total reflectivity factor. For delta Z = 3 dB, f is 50%. For delta Z = 10 db, f is 90%. That is, 90% of Ztotal is due to ice. So the difference reflectivity is very useful for determining mixed phase microphysics.

34. Not all frozen hydrometeors have quasi-spherical shapes and ~0 dB Zdr. Two cold season examples: Positive Zdr’s often occur in pristine, un-aggregated crystals near echo top Positive Zdr layers also noted near -15oC level in active crystal growth regimes