## 1. Introduction

Range oversampling followed by a decorrelation transformation is a recently suggested method for increasing the number of independent samples from which to estimate the Doppler spectrum and its moments, as well as several polarimetric variables on pulsed weather radars (Torres and Zrnić 2003a,b). Range-oversampling techniques rely on the precise knowledge of the range correlation of oversampled signals, which is a function of the transmitter pulse envelope, the receiver filter impulse response, and the distribution of scatterers illuminated by the radar. Theoretical and simulation studies demonstrating the advantages of these techniques have been successfully verified on weather data collected with a single-transmitter dual-polarization radar (Ivić et al. 2002; Torres and Ivić 2005). In contrast, recent experimental results on a dual-transmitter system have revealed some difficulties: if the amplitude and/or phase mismatch between transmission pulses is disregarded in the formulation of the decorrelation transformation, processing of range-oversampled dual-polarization signals with the standard whitening transformation can produce abnormally biased^{1} polarimetric variable estimates (Choudhury and Chandrasekar 2007). These authors concluded that matching the correlation of samples in range for the horizontal and vertical channels is critical to effectively use the whitening transformation introduced by Torres and Zrnić (2003b). Further, they recognized “the need to develop a variant of the whitening transformation algorithm” for dual-polarization systems with mismatched channels. Recently, Hefner and Chandrasekar (2008) proposed a Hermitian symmetric whitening transformation as a means to mitigate the biases observed when applying the originally proposed whitening transformation. In addition, a Hermitian transformation was shown to fix numerical inconsistencies that could arise in the construction of whitening transformations from mismatched polarimetric channels. Although these results are very encouraging, the proposed whitening transformation does not remove the biases completely and the authors conclude that more work is needed to determine unbiased whitening (UWTB) methods. This is the purpose of this work.

Although having a dual-polarization radar system with matched channels is ideal for proper measurement of the polarimetric variables, there are situations in which, despite the best design efforts, one must deal with a system with mismatched polarization channels. Systems with one transmitter and a power splitter, such as the National Severe Storms Laboratory KOUN radar (Zahrai and Zrnić 1993), are relatively immune to mismatches in the horizontal and vertical channels, as demonstrated by Torres and Ivić (2005). Still, differences may arise because of variations in the hardware paths specific to each channel. In contrast, systems with dual transmitters, such as the Colorado State University–University of Chicago–Illinois State Water Survey (CSU–CHILL) radar (Brunkow et al. 2000), are more susceptible to waveform mismatches, as reported by Choudhury and Chandrasekar (2007). This problem is aggravated with magnetron-based radars, because precise control of transmitted waveform phases and frequencies is not possible. Despite this limitation, dual-transmitter systems may be preferred as a way to increase the sensitivity of the radar or to exploit the ability to control each transmitted waveform individually, such as required by the method suggested by Chandrasekar et al. (2007) to simultaneously perform co- and cross-polarization measurements.

This work demonstrates that, by properly accounting for the amplitude and/or phase differences in the two polarization channels, it is always possible to obtain acceptable polarimetric variable estimates from transformed range-oversampled data. Nonetheless, the variance of these estimators increases as the degree of mismatch between the horizontally and the vertically polarized transmitted pulses increases. In such cases, estimators that achieve maximum variance reduction can be obtained by solving constrained minimization problems.

The paper is organized as follows: section 2 reviews the theory behind range-oversampling techniques on dual-polarimetric radars. Section 3 examines the bias in auto- and cross-correlation estimates for systems with matched channels. This is followed by a similar analysis for systems with mismatched channels in section 4, in which a formulation that leads to normally biased estimates of the polarimetric variables is presented. Section 5 discusses the construction of optimum transformations for each of the polarimetric variables. The final section demonstrates the performance of the different transformations using simulated data.

## 2. Range oversampling in dual-polarimetric radars

Traditional sampling of weather radar signals *V* occurs at a rate of *τ*^{−1}, where *τ* is the duration of the transmitted pulse. Oversampling in range entails acquiring polarimetric time series data at increased rates so that *L* complex samples are collected during the time *τ*. This is termed as oversampling by a factor of *L* and has become feasible with the advent of commercial single-board digital receivers (Ivić et al. 2003a) and digital signal processors (Zahrai et al. 2002).

### a. Characterization of range-oversampled dual-polarimetric signals

**v**

*and*

_{H}**v**

*be the sets of*

_{V}*L*oversampled signals in range for the horizontal (H) and vertical (V) polarization channels for a given sample time

*m*. In vector notation,

*H*,

*V*(read as

*H*or

*V*) denote signals corresponding to either the horizontal or vertical channels. The first index in the time series corresponds to range time; the second corresponds to sample time. The two-dimensional correlation of range-oversampled signals considering both range- and sample-time lags is defined as a separable function:

*k*and

*n*are range- and sample-time lags,

*E*[.] is the expected value operation, the superscript * denotes complex conjugation, superscripts (

*R*) and (

*T*) designate range- and sample-time correlations, and subscripts

*Y*and

*Z*can be either

*H*or

*V*denoting signals from the horizontal or vertical channels (e.g.,

*p*and

_{H}*p*be the normalized “modified” pulse envelopes for the horizontal and vertical channels (i.e., the transmitted pulses after each channel’s receiver filter). For a calibrated system, the amplitudes of the pulses are such that they do not bias the power in each channel; that is,

_{V}_{i,j}denotes the element in the

*i*th row and

*j*th column of the Hermitian matrix 𝗖.

### b. Estimation of auto- and cross correlations

*L*oversampled range gates are suitably combined. With this technique, auto- and cross correlations are estimated at each of the

*L*oversampled range gates. These

*L*correlation estimates are averaged to produce one correlation estimate with reduced variance. As with traditional sampling, averaged auto- and cross correlations for the first few lags are used to compute the spectral moments and the polarimetric variables. The focus of this paper is on the estimation of the polarimetric variables: differential reflectivity

*Z*

_{DR}, differential phase Φ

_{DP}, and magnitude of the cross-correlation coefficient

*ρ*. Classical estimators are given by

_{HV}*M*is the number of samples in the dwell time and

*Y*and

*Z*can again be either

*H*or

*V*. Equation (8) can be rewritten as

### c. Transformation of range-oversampled signals

A whitening transformation on range-oversampled time series data can be used to decorrelate these signals before averaging; that is, through a linear transformation, a set of *L* correlated complex samples is transformed into a set of *L* decorrelated (or whitened) complex samples. Because data are uncorrelated, averaging covariances after whitening oversampled signals reduces the variance of estimates by a factor of *L* (Torres and Zrnić 2003a).

^{T}. It is important to note that this decomposition is not unique. As argued by Torres et al. (2004), a family of whitening transformations can be obtained by premultiplying the inverse of a given matrix square root of 𝗖 with any unitary matrix. This fact was exploited by Hefner and Chandrasekar (2008) to construct Hermitian symmetric whitening transformations based on the eigenvalue decomposition of 𝗖.

**x**of

*L*transformed oversampled data at a given sample time is obtained as

*Y*and

*Z*can again be either

*H*or

*V*.

## 3. Range oversampling on systems with matched polarimetric channels

For a radar system with perfectly matched channels, the normalized modified pulses for the H and V channels are the same; that is, *p _{H}* =

*p*. This is not an unrealistic assumption for dual-polarization radars with one transmitter. In this situation, the normalized auto- and cross-correlation matrices are the same [cf. (3)]; that is, 𝗖

_{V}*= 𝗖*

_{VHVH}*= 𝗖*

_{VVVV}*.*

_{VHVV}### a. Autocorrelation estimation

*as*

_{VHVH}*= 𝗖*

_{VHVH}*and 𝗪 also whitens the V channel data.*

_{VVVV}### b. Cross-correlation estimation

*= 𝗖*

_{VHVH}*.*

_{VHVV}## 4. Range oversampling on systems with mismatched polarimetric channels

*p*≠

_{H}*p*. This is more likely to occur in dual-polarization radars with dual transmitters. In this case, auto- and cross-correlation matrices are generally different. Therefore, a whitening matrix that works for the H channel may not work for the V channel and vice versa. Then, it makes sense to consider two independent whitening transformations, 𝗪

_{V}*and 𝗪*

_{H}*—one for each channel. Transformed data*

_{V}**x**are obtained as [cf. (14)]

### a. Autocorrelation estimation

*defined by (13). The same is true for the V channel autocorrelation estimator, with 𝗪*

_{H}*derived from an analogous decomposition of 𝗖*

_{V}*.*

_{VVVV}### b. Biased cross-correlation estimation

_{H}*𝗖

*𝗪*

_{VHVV}*) ≠*

_{V}^{T}*L*. A simple example suffices to demonstrate this problem. Consider the case of mismatched H and V channels such that

*p*=

_{V}*e*, where

^{jα}p_{H}*α*is a real number. Then, according to (3), 𝗖

*= 𝗖*

_{VVVV}*= 𝗛*𝗛*

_{VHVH}^{T}and 𝗖

*=*

_{VHVV}*e*𝗖

^{jα}*. If 𝗪*

_{VHVH}*= 𝗪*

_{H}*= 𝗛*

_{V}^{−1}are the whitening transformations for each channel, it is easy to see that

*α*is an integer multiple of 2

*π*.

### c. Unbiased cross-correlation estimation

*H*and

*V*transformation matrices; that is, let a new set of scaled transformation matrices for the cross-correlation estimator be

_{H,V}=

*γ̃*

_{H,V}𝗪

_{H,V}, where

*γ̃*

_{H,V}are complex constants (the tilde is used throughout to discriminate the cross-correlation estimator transformations from the ones corresponding to the autocorrelation estimator). With these new transformations, (25) becomes

### d. General unbiased correlation estimation

As will be shown in section 6, scaled transformations based on whitening transformations for each channel may not result in polarimetric variable estimators with the lowest possible variance. Thus, it may be advantageous to explore other transformation structures. Other transformations such as pseudowhitening have been proposed as a way to increase the effective number of independent samples while minimizing the noise enhancement effects inherent in the whitening transformation (Torres et al. 2004). In general, we need a procedure to determine the best transformation matrices that result in unbiased correlation estimates for any given situation without being constrained to choosing whitening transformations. The general formulation below produces unbiased auto- and cross-correlation estimates for any transformation matrix structure.

_{H,V}be the set of transformations for the autocorrelation estimator, and

_{H,V}be the set for the cross-correlation estimator, where the basic structure of each matrix can be determined using different criteria (e.g., 𝗔

_{H,V}can be chosen as the whitening matrices for the

*H*and

*V*channel data, whereas

_{H,V}can be chosen as the matrices that diagonalize the range cross-correlation matrix 𝗖

*). To produce unbiased auto- and cross-correlation estimates, these four matrices require scaling given by respective factors, which could be determined as*

_{VVVH}_{H,V}=

*γ*

_{H,V}𝗔

_{H,V}and

_{H,V}=

*γ̃*

_{H,V}

_{H,V}. Note that with this scaling, auto- and cross-correlation estimates are always unbiased. However, the variance of these estimators depends on the basic structure chosen for each transformation matrix.

## 5. Optimum unbiased correlation estimation

*θ̂*, where

*θ*is a function of one or more correlation estimates, we must find the solution to the following constrained minimization problem:

*θ̂*] is a function of one or more transformation matrices 𝗪

_{H}, 𝗪

_{V},

_{H}, and

_{V}. Note that the constraint in the minimization problem guarantees polarimetric variable estimates with acceptable biases (see the appendix). As presented in the previous section, this constraint is easily satisfied by the scaling in Eqs. (29) and (30) for generic (not necessarily whitening) transformation matrices 𝗪

_{H,V}and

_{H,V}corresponding to the auto- and cross-correlation estimators, respectively.

*Y*

_{1},

*Y*

_{2},

*Z*

_{1}, and

*Z*

_{2}can be either

*H*or

*V*to denote transformed signals from either polarimetric channel). This generic quantity can be expanded as

*Z*

_{DR}, Φ

_{DP}, and

*ρ*as functions of the

_{HV}*H*and

*V*transformation matrices.

*Z*

_{DR}estimator, we need a set of transformations that leads to unbiased estimates of autocorrelations (see the appendix); that is, transformed signals are obtained as

**x**

_{H,V}= 𝗪

_{H,V}

**v**

_{H,V}, where 𝗪

_{H,V}satisfies

_{DP}estimator, we need a set of transformations that lead to unbiased estimates of the cross correlation (see the appendix). In this case, transformed signals are obtained as

*X̃*

_{H,V}=

_{H,V}

**v**

_{H,V}, where

_{H,V}satisfy tr(

_{H}*𝗖

_{VHVV}_{V}

^{T}) =

*L*. The variance of the estimator in (40) was given by Ryzhkov and Zrnić (1998) and is adapted here as

*X*s and 𝗪s carry a tilde.

*ρ*, we need two sets of transformations—one for unbiased autocorrelation estimates and another for unbiased cross-correlation estimates—that is, two sets of transformed signals are needed. As shown in (43), the tilde is again used to distinguish the set of transformed data used in the estimation of the cross correlation. The variance of

_{HV}*ρ*estimates for simultaneous transmission and reception of horizontally and vertically polarized signals was given by Torlaschi and Gingras (2003) and is adapted here as

_{HV}*=*

_{X̃HXV}_{H}*𝗖

*𝗪*

_{VHVV}_{V}

^{T}).

It is important to recall that the use of a single optimum transformation matrix applies only if polarimetric channels are perfectly matched. If this is not the case, each polarimetric variable requires its own transformation matrix set for optimum estimates (i.e., two transformation matrices for *Z*_{DR}, another two for Φ_{DP}, and yet another four for *ρ _{HV}*). These variable-specific transformation sets are obtained by solving (31) with the functions in (36), (42), and (45). As in the case of matched channels, these transformation sets can be precomputed to reduce the computational complexity of the optimum unbiased estimators. However, because the transformations arising from the constrained minimization of (36), (42), and (45) depend on

*ρ*, a transformation set should be precomputed for every single value of

_{HV}*ρ*between 0 and 1—a daunting task. In practice, it suffices to precompute transformation sets for a finite set of

_{HV}*ρ*values (e.g., 20 sets would cover the range of

_{HV}*ρ*values with a resolution of 0.05) and select (in real time) the set that best matches the situation at hand to get an “almost optimum” estimator performance. Additionally,

_{HV}*ρ*is not known a priori but is one of the variables to be estimated from oversampled data. To handle this apparent paradox, an initial estimate of

_{HV}*ρ*could be obtained from nonoversampled (or decimated) data to select the proper transformation set using a lookup table. Even with precomputed sets of transformations, a real-time implementation of optimum unbiased transformation (OUTB) estimators on polarimetric radars with mismatched channels would be computationally more expensive than biased whitening (WTB) estimators on systems with matched channels; that is, oversampled time series data for the H and V channels have to be transformed 4 times each, instead of just once as in (14). Also, estimates of auto- and cross correlations cannot be “shared” among the three polarimetric variable estimators. Each variable requires its own correlation estimators. All in all, a real-time implementation of OUTB estimators should be feasible with modern digital signal processing technology.

_{HV}## 6. Simulation results

*α*

_{0}and

*β*

_{0}control constant amplitude and phase mismatches and unitless factors

*α*

_{1}and

*β*

_{1}control time-varying amplitude and phase mismatches (

*p*and

_{α}*p*are unitless functions of range time

_{β}*l*). Note that perfect matching is obtained if

*α*

_{0}= 1 and

*α*

_{1}=

*β*

_{0}=

*β*

_{1}= 0. For this simulation study, we consider two types of functions for

*p*and

_{α}*p*: a linearly increasing function

_{β}^{2}

Polarimetric range-oversampled weather-like data are simulated for varying degrees of channel mismatch to illustrate the performance of the transformations developed in the previous sections. For each mismatch case, 1000 realizations of time series data are generated with the following parameters: the oversampling factor *L* = 5, the number of samples in the dwell time *M* = 64, the signal-to-noise ratio is very large, and the Nyquist velocity is 32 m s^{−1}. The spectral moments are *S _{H}* = 0 dB,

*υ*= 0 m s

^{−1}, and

*σ*= 4 m s

_{v}^{−1}; the polarimetric variables are

*Z*

_{DR}= 1 dB, Φ

_{DP}= 30°, and

*ρ*= 0.985. Range-oversampled data are processed using four matrix transformation sets:

_{HV}oversampling and averaging (OAB), in which oversampled signals are not transformed; that is, 𝗪

_{H,V}=_{H,V}= 𝗜;WTB, in which the same whitening transformation (the one for the H channel) is used disregarding channel mismatch (section 4b); that is, 𝗪

_{H,V}=_{H,V}= 𝗛^{−1}, where 𝗖= 𝗛*𝗛_{VHVH}^{T};UWTB, in which the proper scaling factors are applied (section 4d), 𝗪

_{H,V}are the whitening transformations for each channel, and_{H,V}diagonalize the normalized range cross-correlation matrix 𝗖; and_{VHVV}OUTB, in which three sets of transformations are obtained by solving the constrained minimization problems presented in section 5. These problems are solved using the sequential quadratic programming (SQP) method implemented in MATLAB function “fmincon.” SQP methods are best suited for solving problems with significant nonlinearities in the constraints (Nocedal and Wright 2006), which is the case in our formulation.

Table 1 shows a list of abnormally biased polarimetric variable estimators on range-oversampled data for different channel mismatch cases. Case 1 is the ideal case of matched channels in which all the estimators are unbiased for all the transformations under consideration. Cases 2, 3, and 4 correspond to a constant mismatch in amplitude, phase, and both, respectively. Unless mismatches are properly accounted for through unbiasing factors (UWTB and OUTB), amplitude mismatches lead to an abnormal bias in *Z*_{DR} and phase mismatches lead to an abnormal bias in Φ_{DP}. Note that these constant mismatches can be corrected via *Z*_{DR} and Φ_{DP} calibration constants; thus, biases in these polarimetric variables are typically of little concern. Conversely, biases in *ρ _{HV}* are more difficult to correct through calibration; nevertheless, in these cases

*ρ*estimates remain normally biased. Cases 5, 6, and 7 correspond to time-varying mismatches of the ramp form [cf. (48)] in amplitude, phase, and both, respectively. Finally, case 8 is similar to case 7 but uses (49) for the time-varying amplitude mismatch. Analogously to the previous cases, amplitude (phase) mismatches lead to abnormally biased OAB and WTB estimates of

_{HV}*Z*

_{DR}(Φ

_{DP}), which can be removed through calibration. However, in all these cases, OAB and WTB estimates of

*ρ*are abnormally biased and removal of this bias is not straightforward using calibration techniques.

_{HV}Next, varying degrees of time-dependent mismatches are evaluated. Figures 1 and 2 show the bias and standard deviation of the polarimetric variables corresponding to varying degrees of phase and amplitude mismatch, respectively. These figures show that although WTB has lower standard deviation than OAB (except for *ρ _{HV}* estimates with phase mismatches with

*β*

_{1}larger than about

*π*/9), it produces abnormally biased estimates of polarimetric variables as the degree of mismatch between the H and V channels increases (note that whereas significant

*Z*

_{DR}biases correspond to amplitude mismatches and Φ

_{DP}biases only occur for phase mismatches,

*ρ*biases are of practical concern in both types of mismatches). On the other hand, UWTB has almost the same or better variance reduction as WTB but produces normally biased estimates irrespective of the degree of mismatch between the H and V channels. Notice, however, that the standard deviations of estimates with UWTB (and WTB) increase relative to those of OAB as the degree of mismatch between the H and V channels increases. Indeed, for large degrees of mismatch, UWTB performs worse than OAB (in this specific case, this is mainly observed for

_{HV}*ρ*with phase mismatches with

_{HV}*β*

_{1}larger than about

*π*/9). Conversely, in all cases, OUTB achieves maximum variance reduction of polarimetric variable estimates with minimum bias as predicted theoretically.

## 7. Conclusions

This paper demonstrates that, by properly accounting for the amplitude and/or phase differences in the transmission channels (i.e., by properly scaling the transformation matrices), it is always possible to obtain unbiased polarimetric variable estimates from range-oversampled signals. Nevertheless, as shown by the simulations, the variance of these estimators degrades as the degree of mismatch between the horizontally and the vertically polarized transmitted pulses increases. Still, by solving a constrained minimization problem, it is possible to find transformation structures that result in unbiased auto- and cross-correlation estimates and at the same time achieve maximum variance reduction.

Although it was shown that polarimetric channel mismatches can be properly accounted for, the implementation of optimum unbiased estimators on range-oversampled signals comes at a price because of the additional complexity and computational requirements. However, depending on the maximum acceptable polarimetric variable biases, the conventional use of whitening transformations may still be possible on systems with small channel mismatches; that is, assuming that *Z*_{DR} and Φ_{DP} biases can be effectively removed through calibration procedures, the maximum acceptable *ρ _{HV}* bias can be used to determine the worst channel mismatch conditions for which a whitening transformation is still viable (e.g., for the cases depicted in Figs. 1, 2, a maximum acceptable

*ρ*bias of 0.01 would result in whitening being applicable for amplitude mismatches less than about 0.175 or phase mismatches less than about

_{HV}*π*/12). However, on systems with significant channel mismatches, one might be forced to implement the more-involved optimum unbiased estimators. Fortunately, on polarization-diverse single-transmitter radar systems, such as the planned upgrades of the NEXRAD network (Doviak et al. 2000), significant channel mismatches are not likely to occur.

## Acknowledgments

The author wishes to thank three anonymous reviewers for excellent comments that substantially improved the quality of this paper. Also, I would like to thank Chris Curtis and Dusan Zrnić for engaging in fruitful discussions and providing many constructive comments. Funding was provided by NOAA/Office of Oceanic and Atmospheric Research under NOAA–University of Oklahoma Cooperative Agreement NA17RJ1227, U.S. Department of Commerce.

## REFERENCES

Brunkow, D., Bringi V. N. , Kennedy P. C. , Rutledge S. A. , Chandrasekar V. , Mueller E. A. , and Bowie R. K. , 2000: A description of the CSU–CHILL National Radar Facility.

,*J. Atmos. Oceanic Technol.***17****,**1596–1608.Chandrasekar, V., Bharadwaj N. , and George J. , 2007: Orthogonal channel coding for simultaneous co- and cross-polarization measurements. Preprints,

*33rd Int. Conf. on Radar Meteorology,*Cairns, Queensland, Australia, Amer. Meteor. Soc., P7.10. [Available online at http://ams.confex.com/ams/pdfpapers/123516.pdf].Choudhury, S., and Chandrasekar V. , 2007: Wideband reception and processing for dual-polarization radars with dual transmitters.

,*J. Atmos. Oceanic Technol.***24****,**95–101.Doviak, R. J., Bringi V. N. , Ryzhkov A. V. , Zahrai A. , and Zrnić D. S. , 2000: Considerations for polarimetric upgrades to operational WSR-88D radars.

,*J. Atmos. Oceanic Technol.***17****,**257–278.Hefner, E., and Chandrasekar V. , 2008: Whitening dual-polarized weather radar signals with a Hermitian transformation.

,*IEEE Trans. Geosci. Remote Sens.***46****,**2357–2364.Ivić, I., Zrnić D. S. , and Torres S. M. , 2002: Demonstration of optimum processing of oversampled signals in range to improve Doppler spectral moment estimates. Preprints,

*18th Int. Conf. on IIPS for Meteorology, Oceanography, and Hydrology,*Orlando, FL, Amer. Meteor. Soc., 5.22. [Available online at http://ams.confex.com/ams/pdfpapers/30085.pdf].Ivić, I., Zahrai A. , and Zrnić D. , 2003a: Digital IF receiver—Capabilities, tests and evaluation. Preprints,

*31st Conf. on Radar Meteorology,*Seattle, WA, Amer. Meteor. Soc., 9.3. [Available online at http://ams.confex.com/ams/pdfpapers/64211.pdf].Ivić, I., Zrnić D. , and Torres S. , 2003b: Whitening in range to improve weather radar spectral moment estimates. Part II: Experimental evaluation.

,*J. Atmos. Oceanic Technol.***20****,**1449–1459.Melnikov, V., and Zrnić D. S. , 2004: Simultaneous transmission mode for the polarimetric WSR-88D: Statistical biases and standard deviations of polarimetric variables. NOAA/NSSL Rep., 84 pp. [Available online at http://cimms.ou.edu/~schuur/jpole/SHV_statistics.pdf].

Nocedal, J., and Wright S. , 2006:

*Numerical Optimization*. 2nd ed. Springer-Verlag, 664 pp.Papoulis, A., 1984:

*Probability, Random Variables, and Stochastic Processes*. 2nd ed. McGraw-Hill, 576 pp.Reed, I. S., 1962: On a moment theorem for complex Gaussian processes.

,*IRE Trans. Inf. Theory***8****,**194–195.Ryzhkov, A., and Zrnić D. , 1998: Polarimetric rainfall estimation in the presence of anomalous propagation.

,*J. Atmos. Oceanic Technol.***15****,**1320–1330.Sachidananda, M., and Zrnić D. , 1985:

*Z*measurement considerations for a fast scan capability._{DR},*Radio Sci.***20****,**907–922.Torlaschi, E., and Gingras Y. , 2003: Standard deviation of the copular correlation coefficient for simultaneous transmission and reception of vertical and horizontal polarized weather radar signals.

,*J. Atmos. Oceanic Technol.***20****,**760–766.Torres, S., and Zrnić D. , 2003a: Whitening in range to improve weather radar spectral moment estimates. Part I: Formulation and simulation.

,*J. Atmos. Oceanic Technol.***20****,**1433–1448.Torres, S., and Zrnić D. , 2003b: Whitening of signals in range to improve estimates of polarimetric variables.

,*J. Atmos. Oceanic Technol.***20****,**1776–1789.Torres, S., and Ivić I. , 2005: Demonstration of range oversampling techniques on the WSR-88D. Preprints,

*32nd Int. Conf. on Radar Meteorology,*Albuquerque, NM, Amer. Meteor. Soc., 4R.5. [Available online at http://ams.confex.com/ams/pdfpapers/96151.pdf].Torres, S., Curtis C. , and Cruz J. R. , 2004: Pseudowhitening of weather radar signals to improve spectral moment and polarimetric variable estimates at low signal-to-noise ratios.

,*IEEE Trans. Geosci. Remote Sens.***42****,**941–949.Zahrai, A., and Zrnić D. S. , 1993: The 10-cm-wavelength polarimetric weather radar at NOAA’s National Severe Storms Laboratory.

,*J. Atmos. Oceanic Technol.***10****,**649–662.Zahrai, A., Torres S. , Ivić I. , and Curtis C. , 2002: The open radar data acquisition (ORDA) design for the WSR-88D. Preprints,

*18th Int. Conf. on IIPS for Meteorology, Oceanography, and Hydrology,*Orlando, FL, Amer. Meteor. Soc., 5.10. [Available online at http://ams.confex.com/ams/pdfpapers/29242.pdf].

## APPENDIX

### Bias of Polarimetric Variable Estimators

Throughout this paper, it is assumed that unbiased correlation estimates are required to preserve the integrity of the polarimetric variables; this condition motivates the focus on correlation estimates in section 2 and is used as a constraint to the minimization problem in section 5. However, polarimetric variable estimators are nonlinear functions of auto- and cross-correlation estimates; thus, unbiased correlation estimates do not necessarily lead to unbiased *Z*_{DR}, *ρ _{HV}*, and Φ

_{DP}estimates. In fact, polarimetric variable estimators are inherently biased (or normally biased), as shown by Melnikov and Zrnić (2004). Nevertheless, a condition to ensure that biases of polarimetric variables obtained from range-oversampled signals are not any larger than those obtained with standard processing (no oversampling) is that correlation estimators based on range-oversampled signals be unbiased. This is shown next.

*ŷ*derived as a function of

*n*primary estimators:

*x̂*

_{1},

*x̂*

_{2}, … ,

*x̂*; that is,

_{n}*ŷ*=

*g*(

*x̂*

_{1},

*x̂*

_{2}, … ,

*x̂*), where

_{n}*g*is any function. According to Papoulis (1984), if

*g*is sufficiently smooth near the point (

*E*[

*x̂*

_{1}],

*E*[

*x̂*

_{2}], … ,

*E*[

*x̂*]), the expected value of the derived estimator can be approximated as

_{n}*g*and its derivatives are evaluated at (

*E*[

*x̂*

_{1}],

*E*[

*x̂*

_{2}], … ,

*E*[

*x̂*]) and Cov is the covariance operator. Note that if

_{n}*g*is a linear function, the terms from the double sum vanish because the second derivatives of a linear function are zero.

*ŷ*] =

*E*[

*ŷ*] −

*y*can be approximated as

*E*[

*x̂*] =

_{i}*x*(

_{i}*i*= 1, … ,

*n*)], the bias of the derived estimator depends only on the covariance terms; that is,

*x̂*,

_{i}*x̂*]| ≤

_{j}*x̂*]Var[

_{i}*x̂*]

_{j}*g*

_{1},

*g*

_{2}, and

*g*

_{3}are nonlinear complex functions. For these estimators, based on transformed oversampled signals, the variances in (A5) can be computed using (33) as

*Y*and

*Z*can be either

*H*or

*V*.

*𝗖*

_{XYXY}*)/*

_{XZXZ}*L*

^{2}that originate from using (A6) in (A5). On one hand, for estimators based on standard processing, these factors are always equal to 1 (because

*L*= 1 and the normalized range correlation matrices reduce to scalars; i.e., 𝗖

*= 𝗖*

_{XYXY}*= 1). On the other hand, for estimators based on transformed oversampled signals, these factors are at most 1; that is,*

_{XZXZ}In summary, provided that correlation estimates are unbiased, the biases of polarimetric variable estimators based on transformed oversampled signals are never larger than those of estimators based on standard processing.

Bias (left) and standard deviation (right) of polarimetric variable estimates for oversampled data processed using: (1) oversampling and averaging (OAB), (2) biased whitening (WTB), (3) unbiased whitening (UWTB), and (4) optimum unbiased transformation (OUTB). In addition, the biases of estimates for standard processing (non-oversampled data) without calibration are included as a reference. The H and V channels exhibit a progressively increasing amplitude mismatch with *α*_{1} varying from 0 to 0.2 and *α*_{0} + *α*_{1} = 1.

Citation: Journal of Atmospheric and Oceanic Technology 26, 7; 10.1175/2009JTECHA1209.1

Bias (left) and standard deviation (right) of polarimetric variable estimates for oversampled data processed using: (1) oversampling and averaging (OAB), (2) biased whitening (WTB), (3) unbiased whitening (UWTB), and (4) optimum unbiased transformation (OUTB). In addition, the biases of estimates for standard processing (non-oversampled data) without calibration are included as a reference. The H and V channels exhibit a progressively increasing amplitude mismatch with *α*_{1} varying from 0 to 0.2 and *α*_{0} + *α*_{1} = 1.

Citation: Journal of Atmospheric and Oceanic Technology 26, 7; 10.1175/2009JTECHA1209.1

Bias (left) and standard deviation (right) of polarimetric variable estimates for oversampled data processed using: (1) oversampling and averaging (OAB), (2) biased whitening (WTB), (3) unbiased whitening (UWTB), and (4) optimum unbiased transformation (OUTB). In addition, the biases of estimates for standard processing (non-oversampled data) without calibration are included as a reference. The H and V channels exhibit a progressively increasing amplitude mismatch with *α*_{1} varying from 0 to 0.2 and *α*_{0} + *α*_{1} = 1.

Citation: Journal of Atmospheric and Oceanic Technology 26, 7; 10.1175/2009JTECHA1209.1

Summary of abnormally biased range-oversampling polarimetric variable estimators for different amplitude and phase mismatch cases for oversampling and averaging (OAB), biased whitening (WTB), unbiased whitening (UWTB), and optimum unbiased transformation (OUTB). The amplitude and phase mismatch parameters are those in (47), and ramp and triang. correspond to Eqs. (48) and (49), respectively.

^{1}

Throughout this work, the term abnormally biased is used to denote range-oversampling polarimetric variable estimators having biases larger than their standard (no oversampling) counterparts. Range-oversampling estimators with the same or smaller biases than their standard counterparts are termed normally biased.

^{2}

This formula is for odd *L*.