Identification and Removal of Non-meteorological Echoes in Dual-polarization Radar Data Based on a Fuzzy Logic Algorithm

Bo-Young YE,GyuWon LEE?,,and Hong-Mok PARK

1Department of Astronomy and Atmospheric Sciences,Research and Training Team for Future Creative Astrophysicists and Cosmologists,Kyungpook National University,Daegu 702–701,Korea

2Center for Atmospheric REmote sensing(CARE),Kyungpook National University,Daegu 702–701,Korea

Identification and Removal of Non-meteorological Echoes in Dual-polarization Radar Data Based on a Fuzzy Logic Algorithm

Bo-Young YE1,GyuWon LEE?1,2,and Hong-Mok PARK2

1Department of Astronomy and Atmospheric Sciences,Research and Training Team for Future Creative Astrophysicists and Cosmologists,Kyungpook National University,Daegu 702–701,Korea

2Center for Atmospheric REmote sensing(CARE),Kyungpook National University,Daegu 702–701,Korea

A major issue in radar quantitative precipitation estimation is the contamination of radar echoes by non-meteorological targets such as ground clutter,chaff,clear air echoes etc.In this study,a fuzzy logic algorithm for the identification of non-meteorological echoes is developed using optimized membership functions and weights for the dual-polarization radar located at Mount Sobaek.For selected precipitation and non-meteorological events,the characteristics of the precipitation and non-meteorological echo are derived by the probability density functions of five fuzzy parameters as functions of reflectivity values.The membership functions and weights are then determined by these density functions.Finally,the nonmeteorological echoes are identified by combining the membership functions and weights.The performance is qualitatively evaluated by long-term rain accumulation.The detection accuracy of the fuzzy logic algorithm is calculated using the probability of detection(POD),false alarm rate(FAR),and clutter–signal ratio(CSR).In addition,the issues in using filtered dual-polarization data are alleviated.

dual-polarization radar,non-meteorological echo,quality control,fuzzy logic algorithm

1.Introduction

The weather radar is commonly used to estimate rainfall at high spatial and temporal resolution over large regions.However,radarmomentdata are contaminated by nonmeteorological echo(NME)such as ground clutter,anomalous propagation(AP),chaff etc.Smith et al.(1996)showed that AP echoes lead to systematic overestimation of rainfall. In addition,chaff echo can lead to radar precipitation on clear days.Thus,quality control of radar data is essential for improving radar quantitative precipitation estimation.

Several methods have been used to identify and remove NME:(1)average ground echo mask;(2)time domain filter;(3)frequency domain filter;and(4)moment-based suppression using the temporal and spatial variability of moment data.Moment-based suppression is widely known and accepted for operational application with various statistical techniques(e.g.neural network,Bayes classifier,and fuzzy logic).

The neural network is one of the most commonly used techniques in identifying NME(Haykin and Deng,1991; Grecu and Krajewski,2000;Krajewski and Vignal,2001; Lakshmanan et al.,2007).Grecu and Krajewski(2000)proposed the use of neural networks to identify AP echoes by applying the spatial fluctuation of the reflectivity.Moszkowicz et al.(1994)suggested the Bayes classifier to detect NME by applying discriminating functions of Bayesian theorem. Rico-Ramirez and Cluckie(2008)compared the Bayes classifier with the fuzzy classifier for C-band dual-polarization radardata.TheirBayesclassifiergave a slightly betterperformance than the fuzzy classifier.Recently,Nicol et al.(2011) showed that ground clutter identification is likely to benefit from measurements of the power ratio or clutter phase alignment even when dual-polarization parameters are available in the Bayes classifier.

Fuzzy logic algorithmshave been widely used forthe mitigation of NME due to their simplicity.Cho et al.(2006)derived the normalized frequency distributions of three feature parameters(vertical gradient of reflectivity,standard deviation of reflectivity,and absolute value of radial velocity)to define the characteristics of ground clutter echo,AP echo, and precipitation echo(PRE).Membership functions(MFs) and their weights are derived from these characteristics for given parameter and reflectivity intervals.Their algorithm identifies NME using the total membership value calculatedfrom each MF and weight.The performance of their algorithm is comparable with that from a polarimetric approach. Berenguer et al.(2006)also used a fuzzy logic algorithm that employs the echo top,vertical gradient of reflectivity,spin change(Steiner and Smith,2002),and texture of reflectivity.

Many ground clutter detection algorithms based on dualpolarimetric radar data have also been developed using fuzzy logic algorithms.Gourley et al.(2007)used an algorithm that derives the MFs of the correlation coefficient,standard deviation of differential reflectivity,and standard deviation of differential phase using Gaussian kernel density estimation proposed by Silverman(1981)for C-band dual-polarization radar data.They evaluated the performance using accumulation maps from PRE and NME.In addition,the clutter mitigation decision(CMD)combines three discriminants(clutter phase alignment,texture of reflectivity,and spin)using the fuzzy logic algorithm to identify NME(Hubbert et al., 2009a,2009b)and then applies an adaptive frequency domain filter.The hydrometeor classification algorithm(HCA), which discriminates hydrometeor type as well as ground clutter,also uses a fuzzy logic classifier(Liu and Chandrasekar, 2000;Park et al.,2009;Mahale et al.,2014).

Recently,several ground clutter filters using spectral parameters have been demonstrated in the radar research community.Warde and Torres(2014)developed CLEAN-AP (Clutter Environment Analysis using Adaptive Processing) using the phase of the auto-correlation spectral density.It shows better results than CMD(Torres et al.,2012).STEP (Spectrum-Time Estimation and Processing)can improve the quality of polarimetric radar data using three novel algorithms:Spectrum Clutter Identification(SCI),bi-Gaussian clutter filter,and multi-lag moment estimation(Cao et al., 2012).The SCI algorithm has four discriminants(spectral power distribution,spectral phase fluctuations,spatial texture of echo power,and spatial texture of spectrum width)based on the Bayesian method(Li et al.,2013).Li et al.(2014)also detected ground clutter using a Bayesian classifier for polarimetric radar.

The Ministry of Land,Infrastructure and Transport (MOLIT)has been operating the Mount Sobaek S-band polarimetric radar,situated in complex terrain,since November 2011.In this paper,a moment-based fuzzy logic algorithm is developed by optimizing the MFs and weights for Mount Sobaek dual-polarization radar measurements.Our approach is similar to Cho et al.(2006),except for the fuzzy parameters.However,their approach is adapted into dualpolarization data.In addition,this radar does not archive unfiltered polarimetric parameters and only unfiltered reflectivity is obtained.Thus,the algorithm is further adapted to resolve issues with an absence of unfiltered polarimetric measurements.The detection accuracy of the proposed method is evaluated using accumulated rainfall,POD,FAR,and CSR.

The data used are summarized in section 2 and the methodology is discussed in section 3.Section 4 presents the results of determining the MFs and weights,as well as validation of the proposed algorithm.

2.Data

Moment data from the Mount Sobaek S-band dual polarization radar are used.This radar has been operated by MOLIT since 15 November 2011.The peak power of the Klystron transmitter is 750 kW.The Mount Sobaek radar operates a full volume scan composed of six elevation angles (?0.5?,?0.1?,0.3?,0.7?,1.0?,1.3?)in a simultaneous transmitter and receiver(STAR)mode.The unambiguous range is 150 km in the rain mode.The gate size is 125 m and the pulse repetition frequency is 1000 Hz.

The radar produces unfiltered reflectivity(DZ),filtered reflectivity(CZ),Doppler velocity(Vr,spectrum width,differential reflectivity(ZDR),correlation coefficient(ρHV),differential phase(ΦDP),specific differential phase(KDP),and signal quality index(SQI).Figure 1 shows plan position indicator(PPI)images of all radar parameters for the precipitation case at the elevation angle of?0.5?at 0802 LST(local standard time)6 July 2012.The signal processor of the radar provides one unfiltered parameter,DZ only,and all other parameters are filtered by a Gaussian model adaptive processing(GMAP;Siggia and Passarelli Jr.,2004)clutter filter in the signal processor.Thus,other parameters show areas of“no data”when compared with DZ data.

Two different data sets are used in this study:one for building MFs and weights,and the other for validating the proposed algorithm.The MFs and weights are obtained from seven NME events(377 volume scans)and four PRE events (375 volume scans),listed in Table 1.The NME events include ground clutter,clear air,and chaff echoes.The precipitation events are composed of widespread rain and monsoon fronts.The validation of the fuzzy logic algorithm is performed for three dominant echo types(chaff,clear air,and precipitation),shown in Table 2.The selection of the various events is due to the increase in diversity in the feature parameters.

3.Methodology

3.1.Determination of MFs and weights

The MFs and weights of the Mount Sobaek radar are necessary to detect NME.The MFs and weights are determined by the same procedure described in Cho et al.(2006).Radar moment data are divided into two areas(the NME and PRE areas)based on the average ground echo map.Then,the probability density functions(PDFs)of the five fuzzy parameters are derived for the two areas to define the characteristics of the NME and PRE.Finally,the MFs and weights are derived from these PDFs.

The average ground echo map is produced by the mean reflectivity of ground clutter echoes during 0345–0445 LST 8 July 2012.This period is composed of 25 volume scans and is representative of the normal propagation of the radar beam. No precipitation and clouds are present.The derived average ground echo map at the elevation angle of?0.5?shows significant ground echoes in northwest,south,and southeastareas,where several mountain ranges are located(Fig.1). The NME and PRE areas are classified using this average ground echo map as reference.For the PRE events in Table 1,the PRE area is limited to pixels where DZ＜0 dB Z in the average ground echo map.On the other hand,the NME area is determined for the NME events(Table 1)with the condition of DZ≥5 dB Z in the average ground echo map.

For the selected NME and PRE areas,five fuzzy parameters(FPs)are calculated.They are:the standard deviation of reflectivity[SD(DZ)];the standard deviation of differential reflectivity[SD(ZDR)];the standard deviation of the correlation coefficient[SD(ρHV)];the standard deviation of differentialphase[SD(ΦDP)];and the correlation coefficient(ρHV). The standard deviation of a given variable(V:DZ,ZDR,ΦDP, andρHV)is defined as follows:

Table 1.Events for deriving the MFs and weights of the fuzzy logic algorithm.

Table 2.Events used for evaluation of the fuzzy logic algorithm.

Here,the c is the gate number at which the SD is calculated. The i indicates the gate number.

The PDFs of the five FPs are then calculated with five reflectivity intervals for both the NME and PRE areas.The four intervals of reflectivity(DZ)are:0≤DZ＜10;10≤DZ＜20;20≤DZ＜30;30≤DZ.

For given fuzzy parameters(FPi)and reflectivity intervals (DZi),the MF isderived using the PDFsofthe NME and PRE with the following equation:

Here,the FPstands forSD(DZ),SD(ZDR),SD(ρHV),SD(ΦDP), andρHV.

The MFs of the fuzzy parameters are combined with weights to derive the total membership value as Eq.(3).The total membership value is the same concept as in another study using a fuzzy logic algorithm(Cho et al.,2006):

Here,W indicates the weight,which is derived below:

where the A indicates the overlapping area of PDFNME[FPi, DZj]and PDFPRE[FPi,DZj].A point of intersection between PDFNME[FPi,DZj]and PDFPRE[FPi,DZj]is determined.The overlapping area is then calculated by adding smaller values between normalized PDFNME[FPi,DZj]and PDFNME[FPi,DZj]on both sides of the point.Thus,the FPiwith the large overlapping area is considered as a less significant parameter in the detection of NME.

3.2.Detection of NME using MFs and weights

The five fuzzy parameters are calculated from radar moment data at each gate.The fuzzy membership values of fuzzy parameters are calculated by predetermined MFs and are combined with weights as in Eq.(3)to lead to the total fuzzy membership value(MFtotal).When MFtotalis larger than the threshold(=0.5),the pixel is identified as the NME, and vice versa for the PRE.

Since all parameters of Mount Sobaek radar except DZ are filtered,the derived FP does not fully represent the characteristics of NME.Thus,the performance of the algorithm can be significantly undermined.To avoid this problem, the difference between DZ and CZ is additionally applied to eliminate residual NME.The difference between DZ and CZ reflects the degree of contamination from ground clutter. The larger the difference,the more contaminations PRE will receive.The threshold of 5 dB Z is used to further remove echoes from PRE in this study.

4.Results

4.1.Determination of MFs and weights

For the selected NME and PRE areas from pure NME and pure PRE events in Table 1,the PDFs of NME and PRE are calculated for five fuzzy parameters with reflectivity intervals(Fig.2).Figure 2a shows that the SD(DZ)of the PRE is within 5 dB.The higher the values of the reflectivity interval,the broader and the more skewed the distribution of NME.The PDFs of SD(ZDR)show a distinctive discrepancy between PRE and NME.The PDFPRE[SD(ZDR)]appears mostly at SD(ZDR)≤2 dB,except for 0≤DZ＜10 dB Z(Fig.2b).The higher value of SD(ZDR)for 0≤DZ＜10 dB Z is due to the edge of precipitation echoes.The PDFPRE[SD(ρHV)]shows a smaller value[SD(ρHV)＜0.2] (Fig.2c).The most distinctive separation is shown in the PDFs of SD(ΦDP)(Fig.2d).The PDFPRE[SD(ΦDP)]mostly appears at SD(ΦDP)＜10?.When DZ≥30 dB Z,the overlap area between the PDFs of NME and PRE is almostnegligible. The higher the value of the reflectivity interval,the larger the value ofρHV.(Fig.2e).The PDFs clearly show that the SD of the three dual-polarimetric variables will play a key role in identifying NME.In particular,the SD(ΦDP)is the most important indicator at higher DZ.

Figure 3 shows the MFs of each fuzzy parameter that are calculated from the PDFs of NME and PRE.Each line represents the MF of different reflectivity intervals.The MFs become broader and less steep as the reflectivity decreases, except for SD(ZDR).For example,when SD(ρHV)=0.1,MF =0.9 in DZ≥30 dB Z,whereas MF=0.65 in 0≤DZ＜10 dB Z.This indicates the MFs are good indicators at higher DZ rather than smaller DZ.

Table 3 shows the weights of the five fuzzy parameters derived from the overlapping area of the PDFs of NME and PRE.When DZ≥30 dB Z,SD(ZDR)and SD(ΦDP)have higher weights,whereas the weight ofρHVis the lowest.The SD(ΦDP)shows higher values of weights at overall ranges of DZ intervals.The SD(ZDR)is also an important feature parameter for the classification of NME.

4.2.Determination of NME using MFs and weights

To detect NME for a selected case,the five fuzzy parameters are first calculated from the radar moment data.The NME and PRE are distinctive in the PPI images of the fuzzy parameters,especially SD(DZ)and SD(ZDR)(Fig.4).In particular,the high value of SD(DZ)and SD(ZDR)is distinctive in the areas of normal ground echoes,while these areas are not present in the other parameters.This is due to the GMAPfiltering of power spectra near zero Doppler velocity.In addition,the values of SD(ρHV)and SD(ΦDP)are larger in the residual areas of normal ground echoes.This indicates that the GMAP filtering is not suf ficient to completely eliminate the contamination of ground echoes.Furthermore,the partial beam blocking areas are clearly identi fied with the SD(ρHV) and SD(ΦDP).

Table 3.The weights of the five fuzzy parameters for Mount Sobaek radar.

The value of MFs for each parameter are shown in Fig.5. When MF=1.0(0.0),a favorable condition is indicated for NME(PRE).The membership value of SD(DZ)in the areas of normal ground clutter is close to 1.0,while the membership value for PRE is below 0.6(Fig.5a).The NME and PRE areas are clearly discernible in the membership value of SD(ZDR)(Fig.5b).The membership values of SD(ΦDP) show smaller areas of high MF(Fig.5d).In addition,the value of MF in the areas of partial beam blocking is small, indicating less of an effect of beam blocking in ΦDP.The values of SD(ρHV)andρHVbetter represent the partial beam blockage areas and ground clutter compared with SD(ΦDP). In particular,the high values are shown in the rays along the southwest direction when the value ofρHVis low.The total membership value is higher in the areas of partial beam blocking,residual clutter,and edges.The MFtotalin most PRE areas has a value as small as 0.1.

4.3.Issues with filtered dual-polarization parameters

In the dual-polarization fuzzy logic algorithm,the unfiltered parameter(DZ)and filtered parameters(ZDR,ρHV,and ΦDP)are used.Figure 6a is the PPI image of hourly average reflectivity at the elevation angle of?0.5?from 0802 to 0859 LST 6 July 2012.After applying the fuzzy logic algorithm, ground clutter echoes above 40 dB Z still remain in areas of residual clutter due to filtered dual-polarimetric parameters (Fig.6a).This is more prominent for areas of mixed ground echoes and precipitation.After applying GMAP filtering,the NME characteristics of dual-polarimetric parameters become weaker and less distinctive.

Figure 6c shows the fraction of removed pixels to overall pixels as a function of the difference between DZ and CZ. The minimum fraction exists at DZ?CZ=0 dB.As the values of DZ?CZ increase,the fraction also increases,up to DZ?CZ=25 dB.This indicates that mixed areas of clutter and precipitation echoes are not properly removed with such high values of DZ?CZ.The fraction is then less than 1 at DZ?CZ=25–50 dB.This may be due to residual clutter and precipitation echoes at near zero Doppler velocity.The threshold of 5 dB is further applied to eliminate residual clutters.If the difference between DZ and CZ in a pixel is higher than the threshold,the pixel is removed as NME.After applying the threshold,all residual ground clutter echoes shown in Fig.6a are removed(see Fig.6b).This removed area is similar to the ground echo map above 40 dB Z in Fig.1.In addition,the average reflectivities in the zero isodop area,indicated by the red dashed line in Fig.6b,are also decreased by about 10 dB Z.This is unavoidable due to the use of the DZ?CZ threshold.

4.4.Verification of the dual-polarization fuzzy logic algorithm

4.4.1.Verification with rainfall accumulation

The performance of the dual-polarization fuzzy logic algorithm is validated using rainfall accumulation.The follow-ing rainfall estimator is used for each PPI:

Then,total rainfall accumulations are obtained for the entire period.The total rainfall amount from the original DZ and CZ,and the DZ after applying the fuzzy algorithm,are compared.

Mount Sobaek radar produces filtered dual-polarization data using GMAP.This frequency domain filter has a weakness in that all echoes are filtered regardless of their echo types.These filtered data are significantly biased and can lead to systematic bias in rainfall accumulation.Moreover,the broad clutter spectrum cannot be completely filtered.These drawbacks are shown in the verification using rainfall accu-mulation.

The performance of the dual-polarization fuzzy classifier is validated using rainfall accumulation from 0100 to 1300 LST 22 October 2012(Fig.7).The first column shows accumulated rainfall derived from unfiltered reflectivity(DZ); the second shows that calculated from filtered reflectivity (CZ);and the last column from DZ after applying the dualpolarization fuzzy classifier(FZ).High rainfall above 200 mm associated with ground clutter echoes exists in the rainfall accumulation derived from DZ at all elevation angles. Most of these contaminated pixels are removed in both the rainfall accumulations derived from CZ and FZ.However, R(CZ)and R(FZ)show several different features in terms of residual ground clutter,rainfall amount,bright band,andsecond trip echo.The features are easily discernable in Fig. 7 and Fig 8.High rainfall areas(＞100 km),shown in red dashed circles(Fig.7),are residual ground clutter echoes that cannot be removed by GMAP,but they are properly removed by FZ(black dashed circles).

In addition,the rainfall accumulation of R(FZ)shows higher values than that of R(CZ).Figure 8 shows differences among different accumulations.The difference in rain accumulation is larger over precipitation areas in R(DZ)–R(CZ) than in R(DZ)–R(FZ).This is due to the over-filtering of precipitation echoes by GMAP,whereas the precipitation areas are properly presented in FZ.The relatively larger values of the difference in the northwest and southeast directions(Fig. 8)are due to filtering of precipitation echoes near zero velocity.The dashed circle areas in the right panel show overdeduction of precipitation by FZ.Detailed examination reveals that these areas are due to the second trip echoes(not shown).The SQI values in the second trip echo area are low. The dual-polarization parameters of second trip echo are already removed by the SQI threshold before applying our algorithm in several cases.For this reason,if the measurements of dual-polarization variables are void in these areas,then FZ is removed.Thus,FZ can be used to remove second trip echoes.In addition,the proposed algorithm can reduce the overestimated rainfall due to bright band at＞100 km for the reason thatρHVis low and SD(ρHV)is high in the melting layer.The rain accumulation in Figs.7g and h clearly shows overestimation at far ranges,as shown by the rings close to the 150 km.This overestimation is not shown in Fig 7i and the large value of R(DZ)–R(FZ)in the same area supports this argument.

4.4.2.Verification with skill score

The detection accuracy of the fuzzy logic algorithm is evaluated using the POD and FAR:

The PODand FARare derived based on the contingency table shown in Table 4.For the PRE events,the hit(H)is defined as the PRE echoes with SQI＞0.7 in actual measurements and DZ＜10 dB Z in the average ground echo map.The miss(M)is the eliminated pixel by the fuzzy logic algorithm with the same condition as in H.The false(F)is the remaining pixel in the NME area after applying the fuzzy algorithm.For NME events,the M is the area of DZ＜10 dB Z in the average ground echo map that is not detected by the algorithm.The H is the detected area in the algorithm.Thus,the POD in NMEevaluates the removal efficiency of ground clutter,chaff,and AP echoes.For the NME events,the F cannot be calculated because the precipitation pixels do not exist.

Table 4.Contingency table for the POD and FAR.

Validation is also performed as functions of CSR for the PRE events only.The CSR is defined as

where T0is the unfiltered signal power and C0is the filtered signal power.Since the real-time moment data do not provide T0and C0,the CSR is approximated with DZ and CZ.For a given CSR,the fraction(F)of removed pixels to overall pixels is calculated by:

The term Pallrepresents the overall number of pixels with radar echoes and Preis the number of removed pixels among them.

The validation of the dual-polarization fuzzy classifier is performed for the events in Table 2.Figure 9a shows the PODs as a function of reflectivity thresholds.The average POD for the clear echo and chaff echo events is about 1.00.The average POD for the PRE events is 0.7–0.85.This slightly lower POD is due to the exclusion of mixed areas of ground clutter and precipitation and the removal of edge regions.The average FAR for the PRE events decreases from 0.06 to 0.02 with reflectivity thresholds.In the case of ground clutter,chaff and precipitation echo,NMEs are properly removed by our algorithm(Fig.10).The spatial variations of ZDRand ΦDPin the chaff echo are large andρHVis less than 0.8.

The fraction of removed pixels to overall pixels is shown in Fig.9b as a function of CSR.The fraction is about 0.5 at CSR=0 dB.As the CSR increases,more pixels are removed.When the echoes are more contaminated,they are easily eliminated.

5.Summary and conclusions

A moment-based fuzzy logic algorithm is developed to identify and remove NME by optimizing membership functions and weights for Mount Sobaek dual-polarization radar data.We use five fuzzy parameters[SD(DZ),SD(ZDR), SD(ρHV),SD(ΦDP),andρHV]and reflectivity intervals from selected NME and PRE cases for determination of MFs and weights.For the selected NME and PRE areas,the PDFs of NME and PRE clearly show that the five fuzzy parameters show discernible features to identify NME,especially ΦDP.The MFs and weights are derived using the PDFs of the NME and PRE.The MFs become broader and less steep as the reflectivity interval decreases.The SD(ΦDP)shows higher values of weights at overall ranges of DZ intervals. The SD(ΦDP)is an important feature parameter for classif ication of NME.The fuzzy membership values at each pixel are calculated using predetermined MFs and are combined with weights to lead to the total membership value.When the total membership value is higher than 0.5,this pixel is identified and removed as NME.

After applying the fuzzy logic algorithm,ground clutter echoes still remain in areas of mixed ground and precipitation echoes.This is due to the absence of unfiltered dualpolarimetric measurements provided from Mount Sobaek radar.The signal processor of Mount Sobaek radar provides one unfiltered parameter only.All the others,including dualpolarization parameters are filtered.After GMAP filtering, the NME characteristics of dual-polarimetric measurements become weaker or are removed.Thus,the dual-polarimetric fuzzy parameters cannot be identified as NME after applying the fuzzy logic algorithm.To avoid this problem,our algorithm further applies the threshold of difference between DZ and CZ.The DZ–CZ indicates the degree of contamination from ground clutter.Average reflectivity after applying the threshold(5 dB)shows that residual ground clutter echoes.

The performance of the fuzzy logic algorithm is evaluated using accumulated rainfall,POD,FAR,and CSR.In the comparison of rainfall accumulation,the proposed algorithm has several different results to GMAP filtering of Mount Sobaek radar.The GMAP filters all echoes regardless of echo type. Thus,significant underestimation appears in high rainfall ar-eas compared with the proposed algorithm.Moreover,residual ground clutter echoes still remain and cause high rainfall that cannot be removed by GMAP due to the broad clutter spectrum.These areas are properly removed by R(FZ).In addition,the proposed algorithm can reduce the overestimated rainfall due to bright band and second trip echoes.Values of high SD(ρHV)and lowρHVin the bright band can be detected in the proposed algorithm.Moreover,removal of second trip echo is possible using the SQI threshold.The average POD for NME(PRE)events is about 1.00(0.76).The average FAR for PRE events is 0.05,and decreases as a function of the reflectivity threshold.The fraction of removed pixels to overallpixels increases as a function of CSR and reaches F=1 at CSR=4 dB.

Acknowledgements.This research was supported by a grant (14AWMP-B079364-01)from Water Management Research Program funded by Ministry of Land,Infrastructure and Transport of Korean government.

REFERENCES

Berenguer,M.,D.Sempere-Torres,C.Corral,and R.S′anchez-Diezma,2006:A fuzzy logic technique for identifying nonprecipitating echoes in radar scans.J.Atmos.Oceanic Technol.,23,1157–1180.

Cao,Q.,G.F.Zhang,R.D.Palmer,M.Knight,R.May,and R. J.Stafford,2012:Spectrum-Time Estimation and Processing (STEP)for imprving weather radar data quality.IEEE Trans. Geosci.Remote Sens.,50,4670–4683.

Cho,Y.H.,G.W.Lee,K.E.Kim,and I.Zawadzki,2006:Identification and removal of ground echoes and anomalous propagation using the characteristics of radar echoes.J.Atmos. Oceanic Technol.,23,1206–1222.

Gourley,J.J.,P.Tabary,and J.P.du Chatelet,2007:A fuzzy logic algorithm for the separation of precipitating from nonprecipitating echoes using polarimetric radar observations.J.Atmos. Oceanic Technol.,24,1439–1451.

Grecu,M.,and W.F.Krajewski,2000:An efficient methodology fordetection of anomalous propagation echoes in radarreflectivity data using neural networks.J.Atmos.Oceanic Technol.,17,121–129.

Haykin,S.,and C.Deng,1991:Classification of radar clutter using neural networks.IEEE Transactions on Neural Networks,2,589–600.

Hubbert,J.C.,M.Dixon,S.Ellis,and G.Meymaris,2009a: Weather radar ground clutter.Part I:Identification,modeling, and simulation.J.Atmos.Oceanic Technol.,26,1165–1180.

Hubbert,J.C.,M.Dixon,and S.Ellis,2009b:Weather radar ground clutter.Part II:Real-time identification and filtering. J.Atmos.Oceanic Technol.,26,1181–1197.

Krajewski,W.F.,and B.Vignal,2001:Evaluation of anomalous propagation echo detection in WSR-88Ddata:Alarge sample case study.J.Atmos.Oceanic Technol.,18,807–814.

Lakshmanan,V.,A.Fritz,T.Smith,K.Hondl,and G.Stumpf, 2007:An automated technique to quality control radar reflectivity data.J.Appl.Meteor.Climatol.,46,288–305.

Li,Y.,G.Zhang,R.J.Doviak,L.Lei,and Q.Cao,2013:A new approach to detect ground clutter mixed with weather signals. IEEE Trans.Geosci.Remote Sens.,51,2373–2387.

Li,Y.G.,G.F.Zhang,and R.J.Doviak,2014:Ground clutter detection using the statistical properties of signals received with a polarimetric radar.IEEE Transactions on SignalProcessing,62,597–606.

Liu,H.P.,and V.Chandrasekar,2000:Classification of hydrometeors based on polarimetric radar measurements:Development of fuzzy logic and neuro-fuzzy systems,and in situ verification.J.Atmos.Oceanic Technol.,17,140–164.

Mahale,V.N.,G.F.Zhang,and M.Xue,2014:Fuzzy logic classification of S-band polarimetric radar echoes to identify threebody scattering and improve data quality.J.Appl.Meteor.Climatol.,53,2017–2033.

Moszkowicz,S.,G.J.Ciach,and W.F.Krajewski,1994:Statistical detection of anomalous propagation in radar reflectivity patterns.J.Atmos.Oceanic Technol.,11,1026–1034.

Nicol,J.C.,A.J.Illingworth,T.Darlington,and J.Sugier,2011: Techniques for improving ground clutter identification.Proc. Symp.Weather Radar Hydrol.,IAHS Press,351 pp.

Park,H.S.,A.V.Ryzhkov,D.S.Zrni′c and K.-E.Kim,2009: The hydrometeor classification algorithm for the polarimetric WSR-88D:Description and application to an MCS.Wea. Forecasting,24,730–748.

Rico-Ramirez,M.A.,and I.D.Cluckie,2008:Classification of ground clutter and anomalous propagation using dualpolarization weather radar.IEEE Trans.Geosci.Remote Sens.,46,1892–1904.

Siggia,A.D.,and R.E.Passarelli Jr.,2004:Gaussian model adaptive processing(GMAP)for improved ground clutter cancellation and moment calculation.Proc.3rd European Conf.on Radar in Meteorol.and Hydrology,Visby,Sweden,67–73.

Silverman,B.W.,1981:Using kernel density estimates to investigate multimodality.Journal of the Royal Statistical Society: Series B,43,97–99.

Smith,J.A.,D.J.Seo,M.L.Baeck,and M.D.Hudlow,1996:An intercomparison study of NEXRAD precipitation estimates. Water Resour.Res.,32,2035–2045.

Steiner,M.and J.A.Smith,2002:Use of three-dimensional reflectivity structure for automated detection and removal of nonprecipitating echoes in radardata.J.Atmos.Oceanic Technol.,19,673–686.

Torres,S.,D.Warde,and D.Zrnic,2012:Signal Design and Processing Techniques for WSR-88D Ambiguity Resolution: Part 15.The CLEAN-AP Filter,National Severe Storms Lab., Norman,OK,65 pp.

Warde,D.A.,and S.M.Torres,2014:The autocorrelation spectral density for Doppler-weather-radar signal analysis.IEEE Trans.Geosci.Remote Sens.,52,508–518.

:Ye,B.-Y.,G.Lee,and H.-M.Park,2015:Identification and removal of non-meteorological echoes in dualpolarization radar data based on a fuzzy logic algorithm.Adv.Atmos.Sci.,32(9),1217–1230,

10.1007/s00376-015-4092-0.

27 September 2014;revised 31 December 2014;accepted 19 January 2015)

?Corresponding author:GyuWon LEE

Email:gyuwon@knu.ac.kr