## 1. Introduction

*r*

_{max}and maximum unambiguous velocity

*u*

_{max}of a conventional pulsed Doppler radar are related by (Doviak and Zrnić 1984; Miller and Rochwarger 1972)

*r*

_{max}

*u*

_{max}

*cλ*

*c*is the speed of light and

*λ*is the radar wavelength. While this ambiguity relation applies to all weather radars, millimeter-wavelength radars face a particularly stringent limit due to the increase in Doppler frequency shift with increasing operating frequency. This high Doppler shift often forces the operation of millimeter-wave radars to sacrifice range and accept multiple folding in the velocity data, even for moderate winds. For example, at 3 mm (95 GHz) wavelength the product of the maximum velocity and range is so small that for a

*u*

_{max}of ±8 m s

^{−1},

*r*

_{max}is reduced to 14 km. These limitations may be acceptable during zenith observations of some clouds and precipitation; longer-range capability is often necessary during horizontal measurements, however, and weather phenomena such as convective clouds and storms require the ability to measure much higher wind speeds.

Conventional estimators of Doppler velocity, such as the Fourier transform and pulse-pair techniques, can measure beyond *u*_{max} by using multiple (staggered) pulse-repetition frequencies (PRF) (Doviak and Zrnić 1984); however, the performance of these estimators degrades when the standard deviation of the Doppler spectrum is large compared to the folding velocity. An increase in PRF increases the folding velocity but also decreases the maximum range, which can make interference from scatterers located past *r*_{max} more likely. Golestani et al. (1995) successfully reduced this interference by combining multiple polarizations with staggered PRFs, but this method did not remove the limitation imposed by (1).

Doviak and Sirmans (1973) proposed a polarization diversity pulse-pair (PDPP) technique that can decouple *r*_{max} from *u*_{max}. The technique takes advantage of 1) the isolation between orthogonally polarized signals to prevent ambiguity and 2) the high degree of correlation between the orthogonal copolarized backscatter coefficients (*S*_{υυ} and *S*_{hh}) of atmospheric particles to measure velocity. Figure 1 compares a conventional pulse-pair transmit polarization sequence with the PDPP sequence proposed by Doviak and Sirmans. The technique allows a radar system to measure the first two moments of the Doppler spectrum by simultaneously using two receiver channels to receive the scattered field from two closely spaced, orthogonally polarized pulses. The received signal from the two pulses then is processed like conventional pulse-pair signals. The resulting *u*_{max} measurement is limited only by the receiver bandwidth, and *r*_{max} has no theoretical limit because the pairs can be spaced as far apart as desired. In practice, however, operating with an excessively high *u*_{max} will unnecessarily increase the standard deviation of the Doppler measurements (as shown in section 3), and a very long *r*_{max} implies operation at low average power, which degrades sensitivity.

To the best of our knowledge, PDPP has not been implemented with a weather radar until this project, in spite of the potential benefit of the technique. There are two possible reasons for this. First, PDPP requires a complex radar system with two receiver channels that can simultaneously measure the orthogonal polarization components of the received electric field and can switch the transmit polarization from pulse to pulse. Second, complex, high-frequency weather radars with these capabilities only recently have been developed and low-frequency (S, C, X band) radars almost always can eliminate problems associated with the coupling of *r*_{max} with *u*_{max} using conventional staggered PRF techniques.

The PDPP technique was applied using the University of Massachusetts (UMass) 95-GHz polarimetric cloud radar (Pazmany et al. 1994a, 1994b; Bluestein et al. 1997) to the study of tornadoes and of substructures in supercell storms as part of a joint UMass and University of Oklahoma (OU) experiment. The PDPP technique is necessary because the combination of high operating frequency of the radar, the extreme wind speeds inside a tornado, and the difficulty and risks associated with getting close to a tornado requires high values of both *r*_{max} and *u*_{max}.

In the next section, the conventional pulse-pair technique and PDPP estimation of the first two moments of the reflectivity-weighted Doppler spectrum are described. In section 3, an expression for the standard deviation of the PDPP mean velocity estimate is derived and is compared to that of conventional pulse pair. In section 4, an interleaved PDPP and conventional pulse-pair sequence is described that minimizes Doppler measurement errors while extending both *u*_{max} and *r*_{max}. Images of the reflectivity and Doppler velocity of a severe storm that contained a cyclonic and an anticyclonic hook also are presented. Furthermore, section 4 presents comparisons of the conventional pulse-pair velocity images with the PDPP images and the standard deviation of the PDPP velocity measurements with the standard deviation predicted by the expression derived in section 3. The paper is concluded in section 5.

## 2. Pulse-pair techniques

*τ.*Since a Gaussian-shaped power spectrum is a good approximation of the Doppler spectrum of various turbulent media, the autocorrelation function may be expressed as (Zrnić 1977)

*R*

*τ*

*R*(0)

*βe*

^{j2πfoτ}

*R*(0) is the signal power,

*f*

_{o}is the mean Doppler shift, and

*β*is the normalized signal correlation. Here

*β*is a function of the power spectrum standard deviation

*σ*

_{f}and the pulse-pair spacing

*τ*according to

*β*

*e*

^{−2π2σ2fτ2}

*τ*≪ 1/

*σ*

_{f}, the Doppler velocity mean and standard deviation can be approximated by (Doviak and Zrnić 1984; Miller and Rochwarger 1972)

*μ̂*

_{u}

*σ̂*

_{u}

*υ*) and horizontal (

*h*) polarization basis, the cross-correlation function of the scattered field for a

*υ–h*pair

*R*

_{υh}can be expressed as

*R*

_{υh}

*τ*

*R*

_{υh}(0)|

*βe*

^{j(2πfoτ+ϕ)}

*h–υ*pair,

*R*

_{hυ}can be expressed as

*R*

_{hυ}

*τ*

*R*

_{hυ}(0)|

*βe*

^{j(2πfoτ−ϕ)}

*ϕ*accounts for differential phase shift during propagation between the radar and the scatterers and for any difference in the transmission line lengths between the two polarization channels in the radar. Since this bias has an opposite sign when the polarization sequence is reversed, it can be eliminated by interleaving

*υ–h*pulse pairs with

*h–υ*pairs. After combining

*R*

_{hυ}(

*τ*) with

*R*

_{hυ}(

*τ*), the unbiased estimate of the mean velocity can be expressed as

*R*

_{υh}(

*τ*)| = |

*R*

_{hυ}(

*τ*)| is assumed. The practical use of (9) is limited, however, because |

*R*

_{υh}(0)| is difficult to measure, especially in conditions when conventional pulse-pair technique breaks down. A technique outlined by Doviak and Zrnić (1984) for the estimation of |

*R*

_{υh}(0)|, for example, requires the measurement of the lag 2 autocorrelation

*R*

_{hh}(2

*τ*). However, when this function can be measured, PDPP technique is not needed. The technique by Doviak and Zrnić (1984) could be modified to replace |

*R*

_{hh}(2

*τ*)| with |

*R*

_{υh}(2

*τ*)|, which can be obtained from lag 2 PDPP measurements, but an a priori model for the power spectrum shape also is needed, which model is difficult to estimate in high-shear regions such as a tornado. Noncoherent polarimetric measurement technique (Mead et al. 1993; Pazmany and McIntosh 1994) may be used to estimate |

*R*

_{υh}(0)|, but this technique requires a complex radar that can transmit three or more linearly independent polarizations and can measure the polarization of the scattered field. Without this capability, |

*R*

_{υh}(0)| may be approximated from the copolarized received power at vertical and horizontal polarizations and by assuming that

*S*

_{υυ}and

*S*

_{hh}are completely correlated. This approximation, however, leads to an imprecise estimate of

*σ*

_{u}; therefore, it is useful only to locate regions that have very high spectral width, such as tornadoes.

The implementation of PDPP requires a polarimetric Doppler radar system with some unique capabilities. As seen in Fig. 1, the close pulse spacing and varying transmit polarization during PDPP measurements require a radar with fast transmit/receive and polarization switches, and the radar must have two receivers to measure simultaneously the scattering from the two orthogonally polarized transmit pulses. Such a radar system is the UMass W-band cloud radar, documented in Pazmany et al. (1994a), which was constructed with two receiver channels to speed the measurement of the polarimetric scattering properties of clouds and precipitation from an aircraft.

## 3. Polarization diversity pulse-pair estimate errors

The PDPP technique provides an unbiased estimate of the Doppler velocity and can be used to measure very high wind speeds effectively in the presence of high shear, but it sacrifices accuracy to achieve this. The high unambiguous velocity of PDPP is one of the main contributors to the relatively high estimate error because, with increasing *u*_{max}, the error in the measured phase of the correlation function maps into a larger velocity error. These errors in the measured phase are generated by thermal and phase noise and by interference between the orthogonally polarized signals resulting from the finite polarization isolation of the antenna and orthomode transducer (OMT). These interfering “noiselike” signals add to thermal noise to decrease the effective signal-to-noise ratio, and therefore must be considered when estimating the standard deviation of PDPP velocity estimates. This analysis is presented next.

**V**

_{m}that corresponds to scattering from range

*r*at time

*t*can be written in terms of the signal and noise components according to

*V*

_{ij}is the signal component of the voltage at the output of the

*i*-polarized receiver when

*j*polarization was transmitted, and

*N*

_{υ}and

*N*

_{h}represent system noise in the vertical and horizontal receiver channels, respectively. The signal component is a function of the target scattering matrix and matrices that characterize the polarimetric transformation of the signal during propagation according to Wood (1986)

**D**

_{R}is the normalized [

**D**

_{R}(1, 1) = 1] complex receiver distortion matrix,

**D**

_{P}is the normalized propagation distortion matrix that accounts for differential attenuation

*a*and phase shift

*ϕ*

_{p}according to

**D**

_{T}is the normalized complex transmitter distortion matrix,

**S**

*r, t*) is the complex scattering matrix corresponding to the range cell located at a distance

*r*from the radar at time

*t,*and the complex constant

*k*accounts for system gain and propagation effects. The estimate of the cross-correlation function

*R̂*of the vertical and horizontal signals from

*M*independent, closely spaced (small

*τ*relative to the decorrelation time of the scattered signal)

*υ–h*pairs can be expressed as

*r*=

*cτ*/2,

*t*

_{i}is the sample time of the

*i*th pair’s first sample, and the superscript (*) represents a complex conjugate. Since the cross-polarized signals are scattered from cells located at ranges

*r*± Δ

*r,*only

*V*

_{hh}(

*r*) and

*V*

_{υυ}(

*r*) are correlated, and therefore the expectation (〈 〉) of the estimated cross-correlation function reduces to

*R̂*

_{υh}

*r, τ*

*V*

_{υυ}

*r, t*

*V*

_{hh}

*r, t*

*τ*

*h–υ*pairs, the expression for 〈

*R̂*

_{hυ}(

*r, τ*)〉 takes a similar form except that the

*υ*’s and

*h*’s are interchanged in (14).

*R̂*〉 and therefore do not bias the estimated mean velocity; however, they do increase the standard deviation of the velocity estimate. From (13) it can be shown (see the appendix) that, if the number of independent pulse-pair samples

*M*is sufficiently large, the variance of the unbiased mean velocity estimate,

*P*

_{co}(

*r*) = 〈|

*V*

_{υυ}(

*r*)|

^{2}〉 = 〈|

*V*

_{hh}(

*r*)|

^{2}〉,

*P*

_{cross}(

*r*) = 〈|

*V*

_{υh}(

*r*)|

^{2}〉 = 〈|

*V*

_{hυ}(

*r*)|

^{2}〉,

*P*

_{n}= 〈|

*N*

_{υ}|

^{2}〉 = 〈|

*N*

_{h}|

^{2}〉, and

*M*represents

*M*/2 independent

*υ–h*pairs combined with

*M*/2 independent

*h–υ*pairs. To simplify (15), the following assumptions have been made: no differential reflectivity, reciprocity (

*S*

_{υh}=

*S*

_{hυ}) (Zrnić 1991), and equal gain in the two orthogonally polarized channels of the radar.

*σ*

_{μu}

*P*

_{co}and

*P*

_{cross}can be measured directly, or (11) can be evaluated if

**D**

_{T},

**D**

_{R},

**S**

**D**

_{P}are known or can be approximated; matrices

**D**

_{T}and

**D**

_{R}can be obtained as a result of a polarimetric calibration of the radar (Wood 1986) and

**S**

**D**

_{P}can be approximated based on the shape and composition of the particles that are located in the radar beam. Note that, if

*P*

_{cross}is zero and the cross-correlation function

*R*

_{υh}is replaced by the autocorrelation function

*R,*(15) reduces to the expression for the variance of conventional pulse-pair mean velocity estimate, which is given by Doviak and Zrnić (1984) as

Figure 2 presents a comparison between PDPP and conventional pulse-pair estimate standard deviation as a function of range and spectral width using (15) and (16). The following conditions were assumed for these calculations:

3-mm (95 GHz) radar wavelength,

uniform scatterer distribution,

20-dB signal–to–thermal-noise ratio at a range of 1 km (14 dB at 2 km, 8 dB at 4 km, . . . ),

no attenuation,

20-dB polarization isolation between the co- and cross-polarized signals; this isolation includes the effect of radar system isolation and scatterer depolarization,

conventional pulse-pair spacing of 66

*μ*s (*u*_{max}= ±12 m s^{−1}at 3-mm wavelength),PDPP spacing of 10

*μ*s (*u*_{max}= ±79 m s^{−1}at 3-mm wavelength), and100 independent samples used for each estimate.

^{−1}because of the high unambiguous velocity. The effect of cross-polarization interference on the PDPP estimate is evident at close range (1.5–2 km), where cross-polarized signal scattered from

*r*− Δ

*r*interferes with the copolarized signal scattered from

*r.*This interference becomes worse with decreasing

*r*since it is proportional to the square of the ratio of

*r*to

*r*− Δ

*r.*

An interleaved pulse sequence of PDPP and conventional pulse pairs takes advantage of the benefits of both techniques. Moreover, conventional measurements may be extended several times beyond their maximum unambiguous velocity because, like standard staggered PRF techniques, PDPP measurements can be used to correct folded Doppler data. The next section presents experimental data obtained with such an interleaved pulse sequence.

## 4. Measurements

Since 1993, the UMass 95-GHz polarimetric cloud radar (Pazmany et al. 1994b, 1994a) has been involved in a joint UMass–OU experiment to study mesocyclones (2–5 km wide in thunderstorms) and tornadoes in the southern plains of the United States (Bluestein et al. 1997, 1995). The 95-GHz radar system, a computer-controlled positioner, and a manual hydraulic lift were installed in an OU van so the radar could be elevated through an opening in the roof to collect data. The radar system also was equipped with a video monitor and camera, which was aligned with the radar beam to help the radar operator scan the radar and record its pointing direction during data collection. Figure 3 illustrates the chase van and the radar as it is deployed during a measurement.

For the 1995 spring tornado season, the radar was configured to interleave both alternating *υ–h* and *h–υ* PDPP pulse pairs and conventional pairs. This pulse sequence, shown in Fig. 4, allowed the radar to map the reflectivity field from 1.5 to 11 km and to measure mean Doppler velocity up to ±80 m s^{−1} without ambiguity. Data from an observation of a rotating cloud base (Fig. 5) from the Verification of the Origins of Rotation in Tornadoes Experiment (VORTEX) on 17 May 1995 (Bluestein et al. 1997, 1996) illustrates PDPP measurement quality and the ability of a 95-GHz radar to map the reflectivity and Doppler velocity field in a severe storm (Figs. 6–8). The sampled storm had formed south of Tulsa, Oklahoma, and moved northeastward as a supercell (Doswell and Burgess 1993). When the radar crew intercepted the storm, a gust front extended from the southwest to the north, where a wall cloud (Bluestein 1985) had been observed as a lowered cloud base near the edge of the gust front. The gust front was wrapping up cyclonically as it moved by to the north and northeast (Fig. 5).

The radar-observed reflectivity image (not corrected for attenuation) in Fig. 6a shows two mirror-image hook echoes, each 1 km or less in diameter, along the leading edge of the storm. The conventional pulse-pair velocity image in Fig. 6b, although folded, confirms that the southern member of the pair (0.9, 2 km) corresponds to the cyclonically rotating wall cloud that was visible from the van; the northern member (0.1, 2.6 km), not visible from the van, was rotating in an anticyclonic direction. The hook echoes apparently were associated with counterrotating vortices.

The corresponding PDPP Doppler velocity image in Fig. 7a, depicts a similar windfield with inferior precision but better accuracy in regions where conventional pulse pairs folded in Fig. 6b. The PDPP measurements were sufficiently precise, however, to locate these folded regions and to reveal the number of folds (wind speeds did not exceed 33 m s^{−1} so double- or higher-order folding did not occur in this dataset). The unfolded conventional pulse-pair image is shown in Fig. 7b. Since the storm feature did not contain regions with high spectral width, such as those associated with a tornado, the conventional pulse-pair measurements were more precise everywhere in the image than were the PDPP measurements.

To test the validity of (15), *σ*_{μu}*σ*_{μu}*P*_{cross}(*r* ± Δ*r*)]. The standard deviation of the radar-measured Doppler velocity was estimated spatially from Fig. 7a using a 3 × 3 pixel window (about 90 m × 90 m area at a range of 3 km). The window was centered on each point in the PDPP image and the standard deviation of the enclosed nine points was calculated. The same procedure then was implemented on the conventional pulse-pair image to estimate and to correct for the spatial variation in the actual wind field over the 3 × 3 pixel window. The resulting image of this “measured” standard deviation is shown in Fig. 8a, and, although noisier, in magnitude it agrees well with the calculated *σ*_{μu}

## 5. Conclusions

The data presented in this paper demonstrate the feasibility and the benefit of PDPP technique and how PDPP enables millimeter-wave radars to be used to image severe storm dynamics, including tornadoes and hurricane-force winds. The implementation of the PDPP technique requires a radar system with two orthogonally polarized transmitters and receivers, the capability to switch the transmit polarization rapidly from pulse to pulse, and the capability to receive simultaneously the scattered field from both pulses. High polarization isolation between the two channels is necessary to minimize measurement errors.

The benefit of the PDPP technique is the ability to measure Doppler velocity from long range up to extreme speeds even when the standard deviation of the Doppler spectrum is large; its weakness is the lack of accuracy compared to that of conventional pulse pairs in moderate to low wind conditions and at short ranges. To take advantage of the benefits of both techniques, PDPP may be interleaved with conventional pulse pairs. This pulse sequence may be the ideal “staggered” PRF because it can be effective in almost any weather condition and measurement configuration.

## Acknowledgments

The authors wish to thank Bruce Williams for preparing the radar for the experiment, and David Dowell, Herb Stein, Todd Crawford, and Todd Hutchinson for their help during the storm chases. This work was supported by NSF Grant ATM-9616730 and a supplement to NSF Grant ATM-9320672 at UMass and by supplements to NSF Grants ATM-9019821 and ATM-9302379 at OU.

## REFERENCES

Bluestein, H. B., 1985: Wall clouds with eyes.

*Mon. Wea. Rev.,***113,**1081–1085.——, A. L. Pazmany, J. C. Galloway, and R. E. McIntosh, 1995: Studies of the substructure of severe convective storms using a mobile 3-mm-wavelength Doppler radar.

*Bull. Amer. Meteor. Soc.,***76,**2155–2169.——, ——, D. C. Dowell, J. C. Galloway, R. E. McIntosh, H. Stein, and S. Gaddy, 1996: Observations of sub-storm scale vortices in supercells using a mobile, 3-mm wavelength, pulsed Doppler radar. Preprints,

*18th Conf. on Severe Local Storms Conf.,*San Francisco, CA, Amer. Meteor. Soc., 23–26.——, S. G. Gaddy, D. C. Dowell, A. L. Pazmany, J. C. Galloway, and R. E. McIntosh, 1997: Doppler radar observations of substorm-scale vortices in a supercell.

*Mon. Wea. Rev.***125,**1046–1059.Doswell, C. A., and D. W. Burgess, 1993: Tornadoes and tornadic storms: A review of conceptual models.

*The Tornado: Its Structure, Dynamics, Prediction, and Hazards, Geophys. Monogr.,*No. 79, Amer. Geophys. Union, 161–172.Doviak, R. J., and D. Sirmans, 1973: Doppler radar with polarization diversity.

*J. Atmos. Sci.,***30,**737–738.——, and D. S. Zrnić, 1984:

*Doppler Radar and Weather Observations.*Academic Press, 458 pp.Golestani, Y., V. Chanderasekar, and R. J. Keeler, 1995: Dual polarization staggered PRT scheme for weather radars: Analysis and application.

*IEEE Trans. Geosci. Remote Sens.,***33,**239–246.Mead, J. B., P. S. Chang, S. P. Lohmeier, P. M. Langlois, and R. E. McIntosh, 1993: Polarimetric observations and theory of millimeter-wave backscatter from snow cover.

*IEEE Trans. Antennas Propag.,***41,**38–46.Miller, K. S., and M. M. Rochwarger, 1972: A covariance approach to spectral moment estimation.

*IEEE Trans. Inf. Theory,***18,**588–597.Pazmany, A. L., and R. E. McIntosh, 1994: The use of estimation techniques to reduce noncoherent polarimetric measurement errors.

*IEEE Trans. Antennas Propag.,***42,**1325–1328.——, ——, R. D. Kelly, and G. Vali, 1994a: An airborne 95 GHz dual-polarized radar for cloud studies.

*IEEE Trans. Geosci. Remote Sens.,***32,**731–739.——, J. B. Mead, R. E. McIntosh, M. Hervig, R. Kelly, and G. Vali, 1994b: 95-GHz polarimetric radar measurements of orographic cap clouds.

*J. Atmos. Oceanic Technol.,***11,**140–153.Wood, A. M., 1986: A theoretical study of calibration procedures for coherent and noncoherent polarimetric radars. Tech. Rep. 86011, Royal Signals and Radar Establishment, Malverne, 43 pp.

Zrnić, D. S., 1997: Spectral moment estimates from correlated pulse pairs.

*IEEE Trans. Aerosp. Electron. Syst.,***13,**344–345.——, 1991: Complete polarimetric and Doppler measurements with a single receiver radar.

*J. Atmos. Oceanic Technol.,***8,**159–165.

## APPENDIX

### Standard Deviation of the Polarization Diversity Pulse Pair Estimate of the Mean Velocity

*R̂*

_{υh}from (13) corresponding to range and pulse spacing can be expressed in terms of copolarization, cross-polarization and thermal noise voltages at the output of a radar receiver as

*M*is the number of independent pulse-pair samples,

*r*is range, the subscripts

*υh*of

*R*

_{υh}(

*r*) denote vertical–horizontal pulse pairs,

*τ*is the pulse spacing within the pairs, Δ

*r*=

*cτ*/2,

*V*

_{ij}is the received voltage at

*i*polarization when

*j*polarization was transmitted, and

*N*

_{i}is thermal noise voltage in the

*i*-polarized receiver. Since only

*V*

_{υυ}and

*V*

_{hh}are correlated in (A1), the expected value (〈〉) of

*R̂*

_{υh}reduces to the cross-correlation function

*R*

_{υh}:

*R*

_{υh}

*r*

*R̂*

_{υh}

*r*

*V*

_{υυ}

*r, t*

*V*

^{*}

_{hh}

*r, t*

*τ*

*R*

_{υh}and,

*R̃*

_{υh}, a complex, random component with uniformly distributed phase (see Fig. A1). The variance of the magnitude of this random component

*σ*

^{2}

_{|R̃υh|}

*P*

_{ij}is the average power at the output of the

*i*-polarized receiver when

*j*polarization was transmitted, and

*P*

_{nυ}and

*P*

_{nh}are the average noise power in the vertical and horizontal receivers, respectively.

*ϕ*of

*R̂*

_{υh}can be expressed as

*θ*is the difference between the phase of

*R̃*

_{υh}and the phase of

*R*

_{υh}, and is a uniformly distributed random variable from −

*π*to

*π.*If sufficient numbers of independent pairs are averaged so 〈|

*R̃*

_{υh}|〉 ≪ |

*R*

_{υh}|, (A4) simplifies to

*θ*is distributed uniformly and is independent of |

*R̃*

_{υh}| and |

*R*

_{υh}|, the expectation of the measurement phase reduces to

*ϕ*

*R*

_{υh}

*σ*

^{2}

_{ϕ}

*f*

_{|R̃υh(r)|}

*R̃*

_{υh}|. The resulting expression for the variance of the measurement phase is

*u*is related to the measurement phase according to

*u*

*ϕ*

*λ*

*π*

*τ*

*σ*

_{μ̂υh}

*h–υ*pairs

*σ*

_{μ̂hυ}

*υ–h*and

*h–υ*pulse-pair phase measurements, the standard deviation of the unbiased velocity estimate,

*σ*

_{μ̂u}