An Improved Clutter Suppression Method for Weather Radars Using Multiple Pulse Repetition Time Technique

This paper describes the implementation of an improved clutter suppression method for the multiple pulse repetition time (PRT) technique based on simulated radar data.The suppression method is constructed using maximum likelihood methodology in time domain and is called parametric time domain method (PTDM). The procedure relies on the assumption that precipitation and clutter signal spectra follow a Gaussian functional form.The multiple interleaved pulse repetition frequencies (PRFs) that are used in this work are set to four PRFs (952, 833, 667, and 513Hz). Based on radar simulation, it is shown that the newmethod can provide accurate retrieval ofDoppler velocity even in the case of strong clutter contamination.Theobtained velocity is nearly unbiased for all the range of Nyquist velocity interval. Also, the performance of the method is illustrated on simulated radar data for plan position indicator (PPI) scan. Compared with staggered 2-PRT transmission schemes with PTDM, the proposed method presents better estimation accuracy under certain clutter situations.


Introduction
The velocity and range ambiguity [1] is still one of the major limitations for Doppler weather radar observation.Given a radar wavelength , the relationship between the maximum unambiguous range   and maximum unambiguous velocity V  can be presented as V    = /8, where  is the speed of light.With a uniform pulse repetition frequency (PRF)  PRF observation, V  can be shown as V  =  PRF /4.That means large maximum unambiguous velocity V  can be obtained by increasing the PRF.Along with the increase of maximum unambiguous velocity V  comes a decrease in maximum range   , and vice versa.
Various techniques have been proposed to overcome this problem [2][3][4][5], like dual-pulse repetition time (PRT) and staggered PRT [3,6].The dual-PRT and staggered dual-PRT techniques utilize two PRTs in an  pulses transmission sequence.The dual-PRT method transmits two blocks of uniform PRT pulses with each block at a different PRT ( 1 ,  1 , . . .,  1 ,  2 ,  2 , . . .,  2 ).The optimal velocity estimates can be recovered from the joint analysis of the two Doppler measurements.This method is suitable for the low antenna rotation rate, because the radial velocity remains the same during the  pulse intervals.The staggered PRT method is not restricted by this condition.The PRT alternates between two values ( 1 ,  2 ,  1 , . . .,  1 ,  2 ).The interval of the maximum unambiguous velocity V  was well extended and successively adopted for the Next Generation Weather Radar (NEXRAD) network [7].In order to improve the dealiasing success rate and obtain larger Nyquist velocity, the multiple PRT (-PRT) technique which utilizes multiple interleaved PRTs was proposed.Tabary et al. [8] presented a triple-PRT scheme implemented on the C-band French Radar Network.This technique works well when spectral width is smaller than 3 m/s and SNR is high.A simultaneous multiple pulse repetition frequency (SMPRT) code [9] was introduced and evaluated.This technique can provide enough information to produce a high-resolution measured spectrum for each range gate.Tahanout et al. [10] proposed an optimal 9-PRT scheme for the suitable shape reproduction of the power spectrum of the radar signal.These techniques which used nonuniform time series make filtering of clutter from the radar signal more complicated.Standard clutter filters cannot be applied directly 2 Advances in Meteorology to the nonuniform sampling scheme.Banjanin and Zrnic [11] developed a scheme which consists of two filters that operate sequentially; the overall filter is time-varying with periodically changing coefficients.Cho and Chornoboy [12] introduced a finite impulse response time-varying filter that is applied to -PRT transmission.Siggia and Passarelli [13] proposed a Gaussian model adaptive processing (GMAP) algorithm that filters clutter and recovers weather spectrum after notching the spectrum around zero Doppler.Nguyen et al. [14] and Moisseev et al. [15] presented a parametric time domain method (PTDM) for staggered PRT observations, which can accurately estimate the spectral moments of precipitation echoes even with strong clutter.
In this work, we extend PTDM to -PRT observations.This paper is organized as follows.In Section 2, staggered 2-PRT and -PRT observation techniques are described.The PTDM for clutter suppression and spectral moments estimation is then presented in Section 3. Based on radar data simulations, Section 4 focuses on performance analysis of the proposed method.The error analyses of the proposed method and the staggered 2-PRT technique with PTDM for different noise and clutter scenarios are presented.In order to test the performance of the proposed method, different sample schemes techniques are used for plan position indicator (PPI) observations for comparison.In Section 5, discussions and conclusions are provided.

𝑀-PRT Techniques
2.1.Uniform PRT.Assuming that  is the PRT of the traditional uniform PRT technique, the autocorrelation calculated from return signal time series at lag  is R.The maximum unambiguous velocity is V  = /4.The estimated velocity V can be deduced from R [1] as where arg( ) is argument function.

Staggered 2-PRT.
The staggered PRT method applies a pulse transmission sequence that changes intervals between  1 and  2 .For  2 >  1 , the maximum unambiguous velocity V  will be extended as The autocorrelations calculated from return signal time series at lags  1 and  2 are R1 and R2 .The Doppler velocities V 1 and V 2 can be deduced from R1 and R2 , respectively, with (1).
Two approaches have been introduced to deal with the velocity aliasing.One algorithm uses the ratio of autocorrelations [7], given by If the argument of R1 / R2 exceeds the range of (−, ), the estimated velocity will be ambiguous.The other technique utilizes the maximum unambiguous velocity to revise the aliasing velocity.Suppose V 1 and V 2 are the maximum unambiguous velocity for uniform PRTs  1 and  2 .The dealiased velocity can be shown as where  and  are the correct Nyquist interval number and  1 and  2 are estimate errors.Similar to the staggered 2-PRT method, the maximum unambiguous velocity can be shown as The relationship between the measured true velocity and the estimated velocities V 1 , V 2 , . . ., V  is given by Finding the appropriate integers ( 1 ,  2 , . . .,   ) which minimize the error  1 + 2 +⋅ ⋅ ⋅+  is the way for the measured true velocity retrieval.

Parametric Time Domain Method (PTDM)
Bringi and Chandrasekar [16] have shown that the real and imaginary parts of radar signals follow zero mean normal distribution.Suppose V = [ 1 ,  2 , . . .,   ]  is  samples of received radar signals in one radar resolution volume with sample time  1 ,  2 , . . .,   .Hence, the complex vector can be expressed as the sum of two real Gaussian vectors representing in-phase and quadrature components, respectively, like The multivariate density function of V [15] can be represented as where R = [V(V * )  ] is the covariance matrix and R = V(V * )  is the sample covariance matrix.tr( ) is the trace function.det(R) is the determinant of R.
Assuming that Doppler spectra of clutter and precipitation obey Gaussian distribution, the observed spectrum of Doppler velocity V can be expressed as the summation of independent spectra coming from precipitation   (V), clutter  cl (V), and white noise  N (V), as where   and  cl are the precipitation signal and clutter power,   and  cl are the precipitation and clutter spectrum width, respectively,  2 N is the noise power, V is mean velocity, and  s is sample time.Then, the covariance function can be obtained by fast Fourier transform (FFT) as where  is the vector of unknown parameters, .  is the temporal variable.For -PRT transmission, the element in the covariance matrix R can be written as where Here, int is integral function.mod is the remainder of a division.
Using the multivariate density function (8), the loglikelihood estimation of  is as follows: The parameter vector  can be estimated by minimizing min[− ln ()].1.

PTDM Implementation.
Using the log-likelihood function (13) to retrieve the precipitation and clutter spectral moments, the convergence of a local minimum will lead to wrong estimations in some cases.In order to reduce the estimate errors, we should make sure that the results of minimum (13) come from the convergence of global minimum.According to Moisseev et al. [15], this problem can be solved by selecting the appropriate seed value, especially for mean velocity seed value selection which has the biggest effect on the convergence of the optimization procedure.Similar to Moisseev et al. [15], the search range of unknown parameters is given in Table 2.  tot is the total signal power, and  sys N represents the system noise floor.Using the Nelder-Mead simplex method [17,18], the estimated values of unknown parameters are retrieved.Figure 1 shows a scatter diagram of estimated velocities for 4-PRT technique using PTDM, when CSR = 20 dB.It can be revealed that the estimated velocities are accurately retrieved using this method.The mean velocities almost coincide with true velocities.The standard deviations and biases are presented in Figure 2. The results show that the standard deviations are less than 0.5 m/s and the biases are near zero without increase.The noise and clutter environment is complex and variable in radar detection.The simulation can be carried out In the first type of simulation, the CSR changes from −20 dB to 40 dB.SNR is fixed to +40 dB during the computation.Figure 3 shows the curves of standard deviation and bias of velocity estimates versus CSR, when applying 2-PRT and -PRT schemes.The clutter suppression method used in both transmission schemes is PTDM.The biases increase with CSR at a similar rate for the two PRT schemes.The absolute values of biases do not exceed 0.5 m/s.For the range from −20 dB to 20 dB, the standard deviation are smaller than 1 m/s and obviously the 4-PRT schemes provide better estimation accuracy.The estimation error increases with CSR, especially from 20 dB to 30 dB.When CSR > 30 dB, the values tend to be almost steady near 1.8 m/s.
The second simulation is implemented for 0 dB ≤ SNR ≤ 40 dB and CSR = −20 dB.In Figure 4, the biases and standard deviations are shown for the PTDM velocity estimates.In the [15 dB, 40 dB] interval, the measured velocities for both PRT schemes are sufficiently accurate.Then, the bias and standard deviation increase with the decrease of SNR and reach the maximum near SNR = 0 dB.For the range from 5 dB to 40 dB, the measured velocities are acceptable for the biases close to zero and the standard deviations are less than 0.6 m/s.The performance of 4-PRT schemes has better precision especially when the noise is strong (from 0 dB to 15 dB).
In order to test the effectiveness of the proposed method, the  and  data for plan position indicator (PPI) scan were simulated.Choosing a range-azimuth space of 30 km × 60 ∘ , the resolution is 150 m and 1 ∘ .The radial velocity of each cell in the background is set to one constant.The fan-shape range from the center to 12 km is defined as ground clutter with zero Doppler velocity and spectral width   = 0.25 m/s.At the range from 8.5 km to 25 km, a circle region is set to simulate the wind shear in which the velocity field is unstable and changes rapidly in velocity magnitude and direction.The maximum velocity in wind shear field is 10 m/s and spectral width is equal to 4 m/s.From Figure 4, we have already known that the velocities estimates are desirable over a large range of SNR.This simulation is focused on comparing the observation results between different techniques under different CSR scenarios.SNR is fixed to 40 dB.Figures 5(a Figures 5(a), 5(d), and 5(g) illustrate the estimates of radial velocity sampled at single PRT 1050  using the pulse pair processing (PPP) method.The corresponding maximum unambiguous velocity is 7.14 m/s.The velocity aliasing can be seen obviously in the region of wind shear field.Because of the ground clutter, the estimation of velocity is highly influenced.

Conclusions
In this paper, a technique applying the parametric time domain method to the -PRT transmission schemes is proposed.Based on the assumption that the Doppler spectra of clutter and precipitation follow Gaussian shape, simulation studies are carried out using a case of -PRT schemes.It is shown that the proposed method can effectively estimate Doppler velocity.Two cases are simulated to analyze the standard deviation and bias of measured velocities, which can help to compare the proposed method with 2-PRT transmission schemes with PTDM.The first one describes a situation where SNR is fixed to 40 dB and CSR varies.The biases are nearly the same for both of the methods.As long as CSR is lower than 20 dB, the standard deviation of measured velocity using the proposed method are almost three times less than the 2-PRT transmission schemes with PTDM.The second type of simulation is performed for CSR = −20 dB and 0 dB ≤ SNR ≤ 40 dB.The estimation error of measured velocities using the 4-PRT schemes is much better than that given by the 2-PRT schemes.
The PPI scan of measured velocity is built on simulated radar data with SNR = 40 dB and different CSRs.The results confirm that the proposed method not only solves Doppler velocity aliasing but also performs better in clutter suppression than 2-PRT schemes with PTDM especially when CSR < 30 dB.
The future works will focus on finding the optimal -PRT schemes which minimize the estimation error and its applications to real air observation.The method to improve the computing efficiency is also the next emphasis of our research.

Figure 1 :Figure 2 :
Figure 1: The true velocity (solid line) versus estimated velocities (hollow dot) from 30 simulations when CSR = 20 dB and SNR = 40 dB using PTDM.The mean of estimated velocities of each input velocity is shown as a red asterisk character.

Figure 3 :
Figure 3: Bias (a) and standard deviation (b) of velocity estimates for the 4-PRT (solid line) and 2-PRT (dashed line) with SNR = 40 dB using the PTDM.
[9].-PRT.Mostly, -PRT technique[9]can be regarded as an extension of the staggered 2-PRT technique.This technique chooses a transmission sequence which has several blocks, and every block is formed of  different PRTs.

Table 1 :
Parameters of the simulated radar signal. , −0.9V  , . . ., 0.9V  , V 4.1.Simulation of Radar Signal.The main task of this work is to discuss the accuracy of estimated velocity under different

Table 2 :
Range of unknown parameters.
sys + 10)clutter conditions.So, different values of mean velocity, signal-to-noise ratio (SNR), and clutter-to-signal ratio (CSR) are chosen to generate the observed Doppler spectra.We select a periodic scheme with 4 different PRTs: 1050, 1200, 1500, and 1950 s.Their greatest common factor is   = 150 s.The ratio of PRTs is 7/8/10/13.The corresponding PRFs can be written as 952, 833, 667, and 513 Hz.The parameters of the simulated X-band radar signal are shown in Table