330 likes | 492 Views
Future QPE: Dual-Pol and Gap-Filler Radars. Kevin Scharfenberg University of Oklahoma/CIMMS and NOAA National Severe Storms Laboratory. Quantitative Precipitation Estimation. Dual-polarization in one slide. Current state: linear horizontal E pulses: — — — …
E N D
Future QPE: Dual-Pol and Gap-Filler Radars Kevin Scharfenberg University of Oklahoma/CIMMS and NOAA National Severe Storms Laboratory
Dual-polarization in one slide • Current state: linear horizontal E pulses: — — — … • Original WSR-88D contract specified capability for later upgrade to dual-pol • After upgrade, WSR-88D will transmit simultaneous horizontal/vertical pulse (“slant 45º”): ∕ ∕ ∕ ∕ … • Separate receivers will listen for horizontal and vertical backscatter
Early dual-pol QPE results Areal (basin) Estimates Point Estimates
Early dual-pol QPE results Bias of radar areal rainfall estimates Spring hail cases Cold season stratiform rain
Differential reflectivity (Zdr) Reflectivity (Zh) Similar reflectivity – very different differential reflectivity! Northeast – Mostly large rain drops Southwest – Mostly hail
Differential reflectivity (Zdr) Reflectivity (Zh) Similar reflectivity – very different differential reflectivity! Northwest – relatively large number of relatively small drops Southeast – relatively small number of relatively large drops
Quantitative Precipitation Estimation Warm rain case – A very unusual DSD!
Quantitative Precipitation Estimation Hail case – Z-R relations break down!
Hydrometeor Classification Z ZDR RHI in stratiform rainfall rhv KDP
No Echo Lgt/mod rain Heavy rain Hail “Big drops” Graupel Ice crystals Dry snow Wet snow Unknown AP or Clutter Biological Hydrometeor classification algorithm
Dual-pol QPE Algorithm Operational strategy Where HCA detects Use R= Ground clutter / AP / biologicals 0 Rain R(Z, Zdr) Possible hail below melting layer R(KDP) Wet snow 0.6R(Z) Graupel/hail above melting layer 0.8R(Z) Dry snow / ice crystals 2.8R(Z) R(Z) is from standard WSR-88D R(Z) equations.
Dual-pol and partial attenuation NCAR SPOL radar ; From Vivekanandan et al. 1999, JTech 16, 837-845 Partial terrain blockage
“Gap-Filler” Boundary Layer Radars Courtesy CASA project
“Gap-Filler” Boundary Layer Radars CASA radars Nearest WSR-88D
Radar-based QPE: The Future • Dual-pol WSR-88D upgrade • Dual-pol, low-power “gap-filler” radars • Multiple-radar data mergers incorporating NWP • Corrections for dual-pol radar QPE using rain gages • Incorporation of dual-pol base data vertical profiles • Incorporate corrections for partial beam attenuation (including partial terrain blockage!)
Questions? Kevin.Scharfenberg@noaa.gov
Quantitative Precipitation Estimation R(Z) on a 2 km x 2 km grid
Quantitative Precipitation Estimation Dual-pol QPE on a 2 km x 2 km grid
Hydrometeor Classification * * Height * Increasing value
Quantitative Precipitation Estimation • Operational QPE algorithm • Significant improvement over R(Z), particularly inside 150 km and in heavy rain (and possible hail) • Measurable improvement 150-230 km • Measurable improvement over adjusted R(Z) using vertical Zh profiles/mean-field bias (MFB) corrections • Later work to incorporate multiple radars, corrections using MFB, vertical dual-pol profiles, beam attenuation
Differential Reflectivity (Zdr) • Indicates the presence of larger liquid drops • Hail shafts without a lot of liquid water
Differential Reflectivity (Zdr) Differential reflectivity Zdr = 10 log (Eh/Ev) = Zh - Zv [dB] The reflectivity-weighted mean axis ratio of scatterers in a sample volume Zdr > 0 Horizontally-oriented mean profile Zdr < 0 Vertically-oriented mean profile Zdr ~ 0 Near-spherical mean profile Ev Eh
Differential Phase Shift (fDP) Differential Phase Shift fDP = fh – fv (fh, fv≥ 0) [deg] The difference in phase between the horizontally- and vertically-polarized pulses at a given range along the propagation path. - Two-way process - Independent of partial beam blockage, attenuation - Independent of absolute radar calibration - Immune to propagation effects on calibration - Independent of system noise
Specific Differential Phase Shift (KDP) Specific Differential Phase Shift fDP(r2) – fDP(r1) KDP = [deg/km] 2 (r2 – r1) The range derivative of differential phase shift - Identify areas with significantly non-spherical scatterers (usually, rain) - Can estimate rain amount in rain/hail mixture
Specific Differential Phase Shift (KDP) Result: The KDP dilemma - Using a long-distance derivative for calculating KDP can oversmooth heavy rain features but reduces noise - Using a short-distance derivative for calculating KDP retains features in heavy rain but is also noisy
Specific Differential Phase Shift (KDP) Calculating KDP: current practice - If Z > 40 dBZ, use a KDP calculation range of 9 gates (2 km). - Otherwise, use a range derivative of 25 gates (6 km) - Filter the final KDP product at 0.9 rhv
Outline Differential phase shift (fDP)
Outline Specific differential phase shift (KDP)
Quantitative Precipitation Estimation Rainfall estimation using polarimetric variables R(Z, ZDR) = 0.0142 Z0.77 ZDR-1.67 [mm/h] R(KDP) = 44|KDP|0.822 sign(KDP) [mm/h]