The Chinese Meteorological Society
Article Information
- ZOU, Haibo, Shuwen ZHANG, Xudong LIANG, et al., 2018.
- Improved Algorithms for Removing Isolated Non-Meteorological Echoes and Ground Clutters in CINRAD. 2018.
- J. Meteor. Res., 32(4): 584-597
- http://dx.doi.org/10.1007/s13351-018-7176-9
Article History
- Received December 5, 2017
- in final form May 14, 2018
2. State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081;
3. Meteorological Disaster Emergency Warning Center of Jiangxi Province, Nanchang 330046
Doppler weather radar data have high temporal and spatial resolutions and have been widely applied in operational forecasts and research on mesoscale weather phenomena (Anagnostou and Krajewski, 1999; Alfieri et al., 2010; Pan et al., 2010; Qin et al., 2014; Liu et al., 2016), presenting significant advantages compared with other observational data (Serafin and Wilson, 2000; Zheng et al., 2004; Xiao et al., 2008). For instance, radar data can be used to generate high-resolution quantitative precipitation estimations (Fulton et al., 1998; Qi and Zhang, 2017) and are often assimilated into numerical weather models for improving the prediction of convective systems (Jones et al., 2015; Ridal and Dahlbom, 2017). In radar data, however, there are not only meteorological echoes, but also non-meteorological echoes, such as ground clutters and isolated echoes that may comprise biological echoes, electronic interference, and so on (Steiner and Smith, 2002; Zhang et al., 2005). The occurrence of isolated non-meteorological echoes and ground clutters can often degrade the performance of applications of radar data in forecasting and research. Therefore, it is important to remove isolated non-meteorological echoes and ground clutters from radar data before they can be used.
Various techniques have been applied to identify and remove isolated non-meteorological echoes and ground clutters from radar data. For identifying isolated non-meteorological echoes, Fulton et al. (1998) simply defined them as gates with reflectivity exceeding a certain threshold (18 dBZ as default) and having no more than 1 in 8 surrounding neighbors above that threshold. Bergen and Albers (1988) and Krishnapuram and Keller (1993) presented a fuzzy logic based classification scheme to identify isolated points in a picture. Their scheme has been used in meteorological fields to identify isolated non-meteorological echoes, and was incorporated into the Weather Surveillance Radar-1988 Doppler (WSR-88D) system (Zhang et al., 2004) in the USA and the Severe Weather Automatic Nowcast (SWAN) system in China (Wu et al., 2013). Specifically, the method involves computing the percentage of gates with non-missing values in a box (window) and comparing the results with a given default threshold. If the percentage is less than the threshold, the center point (gate) of the box is regarded as an isolated echo. Although the method performs well in eliminating isolated non-meteorological echoes, it simultaneously removes the edges of meteorological echoes (Wu et al., 2013).
Many studies have also focused on methods to iden-tify ground clutters. For example, using the notion that ground clutters are characterized by near-zero Doppler velocities, Chrisman et al. (1995) developed an automated procedure to reduce or eliminate ground clutters using the radar-observed radial velocity. Unfortunately, however, the procedure also eliminates meteorological echoes with near-zero Doppler velocities, because they exhibit characteristics similar to those of ground clutters. In a later work, Fabry and Gadoury (2009) indicated that although the peak contribution from ground clutters is centered around zero velocity, ground clutters can extend to other velocities and even exceed those of meteorological echoes if they are strong. Meanwhile, Hubbert et al. (2009a) pointed out that the actual measured velocity from ground clutters can vary somewhat from zero and the spectrum width from the clutter can be widened (theoretically, it should be narrow) due to rotational scanning. Hence, such an approach would bias the resulting reflectivity. In radar-based quantitative precipitation estimation, Fulton et al. (1998) simply defined ground clutters as gates whose reflectivity exceed a fixed threshold (65 dBZ in default). Pamment and Conway (1998) devised a scheme to detect ground clutters by combining multiple independent observational data, such as surface observations, satellite images, and so on. However, it is difficult to apply this scheme in practice because of the sparse spatial and temporal resolutions of independent observations. Steiner and Smith (2002) used a decision tree comprised of the vertical extent of echoes, the spatial variability of reflectivity, and the vertical gradient of reflectivity, to identify ground clutters. However, this method also removes boundary layer features (e.g., gust fronts) and cannot handle strong, clear-air echoes, and the ground clutters appearing at two adjacent tilts because of the small vertical gradient of reflectivity. Kessinger et al. (2003) developed a detection algorithm for ground clutters based on fuzzy logic classification (depending on the echo variation on the vertical and azimuth-radial plane), which has since been introduced in the WSR-88D system in the USA. The algorithm has worked well for anomalous propagation and most ground clutters (Liu L. P. et al., 2007; Liu S.-Y. et al., 2008). However, a disadvantage is that it cannot identify the ground clutters presented at two adjacent tilts because they can induce small vertical gradients of reflectivity. Zhang et al. (2004) improved this algorithm by introducing a reference tilt with 3–4.5-km height, instead of an upper tilt, so that it can identify ground clutters appearing at two adjacent tilts. The method has also been introduced into the WSR-88D system in the USA (Zhang et al., 2004) and the SWAN system in China (Wu et al., 2013). Nevertheless, this improved method can only be used for the gates within a certain range from the radar, as the reference tilt is limited by the radar beam height of 4.5 km.
Recently, Hubbert et al. (2009a) constructed a new parameter [called “clutter phase alignment” (CPA)], based on the received base-band complex signal [referred to as the in-phase and quadrature-phase (I and Q) components] at the signal-processing level (i.e., level-I data). It is a good indicator in the identification of ground clutters. For example, the CPA is greater than 0.90 for the vast majority of ground clutters, whereas it is less than 0.5 for most weather echoes. Subsequently, Hubbert et al. (2009b) proposed a clutter mitigation decision (CMD) algorithm (a fuzzy logic recognition algorithm) using a combination of the CPA and other measures (such as the texture of reflectivity, the spatial variability of reflectivity, and so on). The CMD algorithm is designed to run in real time on operational radars and has been implemented in the updated WSR-88D Open Radar Product Generator (Hubbert et al., 2009b). Torres and Warde (2014) introduced a clutter environment analysis using adaptive processing (CLEAN-AP) filter, which is also based on the received base-band complex signals (I and Q) and can adjust its suppression characteristics in real time to match dynamic atmospheric environments. With the operational application of more advanced dual-polarization radars, polarimetric quantities such as the cross-correlation coefficient, the standard deviation of differential propagation phase measurements, differential reflectivity, etc., have been used to classify ground clutters and shown better results (Friedrich et al., 2006; Gourley et al., 2007; Rico-Ramirez and Cluckie, 2008; Lakshmanan et al., 2014).
The China New Generation Weather Radar (CINRAD) network consists of more than 170 radars comprising S-band (SA, SB, and SC types) and C-band (CB, CC, and CD types) instruments manufactured by three companies (Xu et al., 2017). These radars are now single-polarization radars (Su et al., 2018) and their received base-band complex signals (I and Q) at the signal-processing level are not stored in operation. Therefore, although the CMD algorithm (Hubbert et al., 2009b), CLEAN-AP filter (Torres and Warde, 2014), and methods based on polarimetric quantities (Friedrich et al., 2006; Gourley et al., 2007; Rico-Ramirez and Cluckie, 2008; Lakshmanan et al., 2014) have performed well in identifying ground clutters, they cannot be used in CINRAD (i.e., level-II data of single polarization) to eliminate ground clutters. Thus, for CINRAD data, the algorithms of Kessinger et al. (2003) and Zhang et al. (2004) are still a good choice for removing isolated non-meteorological echoes and ground clutters. However, as mentioned above, these methods can also remove the edges of meteorological echoes and cannot identify the ground clutters appearing at two adjacent tilts. Thus, our main goal here is to improve these algorithms such that they can remove isolated non-meteorological echoes whilst preserving the edges of meteorological echoes, as well as identify ground clutters (especially those appearing at two adjacent tilts) on the overall azimuth–range plane.
Following this introduction, the data and processing procedures are described in Section 2. The improved algorithms for identifying isolated echoes and ground clutters are presented in Sections 3 and 4, respectively. A statistical analysis is presented in Section 5, followed by conclusions in Section 6.
2 Data and processingThree CINRADs located in Jiangxi Province (Fig. 1a) are used to test the methods of removing isolated non-meteorological echoes and ground clutters from radar base (i.e., level-II) data in this study. The locations of these three radars and altitude of their surroundings are shown in Fig. 1b. The radars located at Nanchang and Shangrao are S-band with type-A CINRADs (i.e., CINRAD/SA), which is a similar instrument to the WSR-88D radar. The one located at Ganzhou is an S-band with type-C CINRAD (i.e., CINRAD/SC), which is produced by 784 Factory in Chengdu, China. Up to 973 automatic rainfall observational stations around Nanchang are used to calculate the rainfall frequency from 0000 UTC 5 to 1600 UTC 14 March 2015. These data are obtained from the Meteorological Information Center of Jiangxi Province, China. The radar data are processed in polar coordinates (radar coordinates), but the corresponding displays are all in Cartesian coordinates. The method of horizontal interpolation with the nearest neighbor on the range–azimuth plane (Xiao et al., 2008) is selected to convert the radar data from polar to Cartesian coordinates with a resolution of 0.01° × 0.01°.
3 Isolated echoes 3.1 Original methodIn the raw data observed by Doppler weather radar, there are sometimes many isolated non-meteorological echoes (i.e., noise), derived from sources such as insects, an aircraft, and so on. These isolated non-meteorological echoes obviously affect the application of radar data in estimating rainfall and monitoring weather, as well as data assimilation. Bergen and Albers (1988) and Krishnapuram and Keller (1993) presented a scheme based on fuzzy logic classification to identify isolated points in a picture, and the scheme was then applied in the quality control of radar data to remove isolated non-meteorological echoes (Zhang et al., 2004). The expression of the scheme is
${P_{\rm X}} = \frac{N}{{{N_{\rm {total}}}}}, $ | (1) |
where N is the number of grid points with valid reflectivity in a window (default 5 × 5, i.e., 5 radials and 5 azimuths) whose center is located at point X, and
The reflectivity measurements at 0.5° elevation observed by Nanchang CINRAD/SA radar at 2203 UTC 18 June 2016 are selected to test the method (Fig. 2a), and the results show that the method can eliminate the isolated echoes remarkably well (Fig. 2b). However, the method also partially removes the margins of meteorological echoes, such as the encircled area in Fig. 2b. Wu et al. (2013) also pointed out this shortcoming of the me-thod. Intuitively, this deficiency is not particularly problematic in synoptic- and meso-α-scale meteorological echoes, because the area of the removed echoes is relatively small. However, in meso-β- and meso-γ-scale meteorological echoes (i.e., small-scale convective cells), the situation is quite different. Figure 3 depicts the primitive convective echoes observed by the Ganzhou CINRAD/SC radar and the convective echoes filtered by the current method. As shown in Fig. 3, the echo area is significantly reduced after the removal of isolated non-meteorological echoes based on the current method (Fig. 3b). The areas of all the convective cells (clouds A–F) are reduced by more than 30% (Table 1). Particularly, the area of cloud C decreases from 32 to 9 km2 (a reduction of 72%), such that it is difficult to find in Fig. 3b. Moreover, the method can also remarkably decrease the coverage of radar echoes in nascent squall lines. This is because their shape is often long and narrow, meaning that the edge of echoes can be removed more effectively. For example, the area of the nascent squall line that occurred in northeastern Nanchang at 1003 UTC 10 June 2017 (Fig. 4a) decreases by about 30%, from 2189 to 1543 km2. In particular, west of the squall line, the echo area decreases by about 50% (Fig. 4b). It is clear that the method has a shortcoming in that, when it eliminates isolated echoes, it partially removes the margins of meteorological echoes too. This situation is quite obvious in small-scale convective cells and squall lines with a long and narrow shape, because the edges of meteorological echoes account for a large proportion of the total echo area.
Cloud | ||||||
A | B | C | D | E | F | |
Area | 123 | 81 | 32 | 87 | 374 | 56 |
Area after removing isolated echoes | 81 | 51 | 9 | 36 | 223 | 26 |
Percentage reduction (%) | 34 | 37 | 72 | 58 | 40 | 53 |
What has caused the problem in the original method? In fact, when a window comprising 5 × 5 grid points is located at the edge of a storm, only half of all 25 grid points have valid radar reflectivity (Fig. 5a). The variable
${P_{\rm o}} = \frac{{{N_{\rm {pv}}}}}{{{N_{\rm {pt}}}}}, $ | (2) |
where
Theoretically, when echoes in the window are disconnected with external echoes (i.e., isolated echoes),
Figures 2a, 3a, 4a, and 6a show that the intensity of most isolated non-meteorological echoes is less than 20 dBZ, and the echoes with reflectivity larger than 20 dBZ are considered as meteorological echoes. Therefore, the areas of echoes over 20 dBZ can imply the range of meteorological echoes. The rate of reduction of echoes over 20 dBZ is defined as
${P_{\rm r}} = 1 - \frac{{{A_{\rm a}}}}{{{A_{\rm b}}}}, $ | (3) |
where
|
||||||||||
0.042 | 0.084 | 0.125 | 0.167 | 0.208 | 0.25 | 0.292 | 0.333 | 0.375 | 0.416 | |
Case 1 | 8.2 | 13.1 | 19.4 | 26.2 | 36.8 | 51.9 | 79.4 | 96.5 | 99.9 | 100.0 |
Case 2 | 4.1 | 8.2 | 12.5 | 16.5 | 24.3 | 38.1 | 65.2 | 89.5 | 99.8 | 100.0 |
Case 3 | 2.3 | 4.3 | 6.9 | 9.8 | 17.2 | 28.1 | 53.3 | 89.4 | 98.3 | 100.0 |
Case 4 | 0.5 | 1.4 | 2.7 | 4.2 | 7.9 | 11.2 | 20.3 | 44.2 | 80.2 | 98.2 |
Case 5 | 0.0 | 0.0 | 0.0 | 0.0 | 0.1 | 0.2 | 0.5 | 1.2 | 4.3 | 12.8 |
Note: “Bold” is used to divide small values (left) and big values (right). |
Although the above improved algorithm is very useful for removing small, isolated non-meteorological echoes, it is impossible for it to eliminate ground clutters because ground clutters often occupy a relatively large area. However, ground clutters usually show worse spatial continuity than meteorological echoes, in both horizon-tal and vertical directions (especially the latter). Therefore, the difference in the spatial discontinuity between ground clutters and meteorological echoes can be used to identify and remove the ground clutters.
4.1 Original methodMeteorological echoes usually extend from the lower troposphere to the middle or upper troposphere, having a relatively deep vertical extent. Thus, the variation in reflectivity in the vertical direction is relative small (i.e., the vertical gradient of echoes is relatively small). In contrast, ground clutters often appear at only one or two elevation angles, with an abrupt change in intensity in the vertical direction. Based on this, the vertical discontinuity of reflectivity is widely used to filter out ground clutters (Steiner and Smith, 2002; Kessinger et al., 2003; Zhang et al., 2004; Wu et al., 2013). In the WSR-88D system in the USA (Kessinger et al., 2003), the vertical gradient of reflectivity with range weighting (RGDZ) is selected as the main variable for the identification and removal of ground clutters from its single-polarization, level-II radar data. The expression of RGDZ is
${\rm {RGDZ}} = W\left({Z - {Z_{\rm {up}}}} \right), $ | (4) |
where Z is the reflectivity value at the selected elevation angle,
To address this problem, Zhang et al. (2004) designed a new parameter, VDZ, which can be expressed as
${\rm{VDZ}} = \frac{{Z - {Z_{\rm {ref}}}}}{{{H_{\rm {ref}}} - H}}, $ | (5) |
where H is the height of the radar beam at the selected elevation angle; Zref and Href are the reflectivity and height at the upper reference elevation angle, respectively. The first upper elevation with radar beam height between 3 and 4.5 km will be treated as the reference elevation (i.e., 3 km ≤ Href ≤ 4.5 km). Thus, VDZ can remove ground clutters better than RGDZ under the condition of ground clutters appearing at two adjacent elevation angles. This promising performance of VDZ is the result of Zref not being a ground clutter (because Href is more than 3 km and ground clutters are often below 3 km). Similar to the situation with RGDZ, a large VDZ indicates existence of ground clutters, while a small VDZ corresponds to meteorological echoes. The VDZ parameter is applied in the WSR-88D radar operational system in the USA (Zhang et al., 2004) and the SWAN system in China (Wu et al., 2013).
Although VDZ has the advantage over RGDZ in terms of removing ground clutters, it is confined by the Href being less than 4.5 km. There are 9 elevation angles for the VCP21 mode (common mode) in the CINRAD SA/SC radar, and these elevation are 0.5°, 1.5°, 2.4°, 3.4°, 4.3°, 6.0°, 9.9°, 14.6°, and 19.5°. The corresponding range (
${({R^2} + R_{_W}^2 + 2R \cdot {R_{_W}} \cdot \sin \theta)^{1/2}}\,{\rm{ - }}\,{R_{_W}} \,= \, 4.5, $ | (6) |
where
The parameters RGDZ and VDZ have their own respective advantages and disadvantages. Based on their advantages, a new parameter, NDZ, can be constructed as follows:
${\rm{NDZ}} = \left\{ \begin{gathered} W\left({Z - {Z_{\rm {ref}}}} \right)\quad R \leqslant {R_{\max }} \\ W\left({Z - {Z_{\rm {up}}}} \right)\quad R > {R_{\max }} \\ \end{gathered} \right..$ | (7) |
Here, W is the same as the weighting function in Eq. (4), and
The removal of isolated echoes and ground clutters is performed by following two steps. First, under the condition of NDZ ≥ 20 dBZ, the ground clutter is removed from the radar base (level-II) data. Second, combining the two conditions
To further evaluate the performance of the improved algorithms in removing ground clutters and isolated non-meteorological echoes, the frequency images of original reflectivity over 20 dBZ at 0.5° evaluation, echo top (ET) height more than 3 km, precipitation more than 0.1 mm h–1, and processed (using the improved algorithm) reflectivity over 20 dBZ at 0.5° elevation, are all generated from 0000 UTC 5 to 1600 UTC 14 March 2015 (232 scans with a 1-h temporal resolution; Fig. 10). This pe-riod encompasses a series of precipitation and clear-sky echoes. The frequency image has been used previously by Steenburgh et al. (2000) to describe the echoes induced or enhanced by the Great Salt Lake in the USA. In the frequency image of the raw reflectivity over 20 dBZ (Fig. 10a), most regions are filled by points with low frequency (less than 0.15). However, several scattered points with high frequencies (larger than 0.2, and even reaching 0.5) are apparent in the surroundings of Nanchang, as well as some flaky areas with high frequencies (larger than 0.2, and even reaching 0.8) near the mountains located to the west, north, and south of Nanchang. These regions with high frequency are evidently caused by isolated echoes (possibly from the tall buildings of Nanchang city) and ground clutters, because they are more stationary compared with meteorological echoes. The ET height can reflect the development of a storm and the intensity of the associated precipitation system (Bedka et al., 2015). Normally, the ET height of meteorological echoes is more than 3 km. Therefore, the frequency image of ET height more than 3 km (Fig. 10b) can to a large extent depict the distribution of meteorological echoes. As shown in Fig. 10b, the frequencies in the whole region are less than 0.2. Particularly, in the regions where the frequencies in Fig. 10a are more than 0.2, the frequencies of ET height are less than 0.15, or even 0.1. This further verifies that the regions with high frequency in Fig. 10a are created by non-meteorological echoes. The frequency distribution of precipitation with a 1-h temporal resolution (Fig. 10c) also illustrates this viewpoint, because there are no high frequencies in the regions shown in Fig. 10a. After removing the ground clutters and isolated non-meteorological echoes using the improved algorithms, in the frequency image of reflectivity over 20 dBZ at 0.5° elevation, there is no point with a frequency more than 0.2 (Fig. 10b) and the distribution of frequencies less than 0.15 coincides with that of unprocessed frequencies (Fig. 10a). It is clear that the improved algorithms perform well in removing the isolated non-meteorological echoes and ground clutters, whilst preserving the structure of storms at the same time. The areas with the highest frequency (> 0.15) are located over the mountain region in southern Nanchang inFig. 10d, but the areas with the highest frequency in Figs. 10b, c are not. This indicates that, although the improved algorithm for identifying ground clutters has performed well, it cannot completely remove all ground clutters. Note that in Figs. 10a, b, d, the areas with low frequency (< 0.05, and even 0.1) are mainly induced by the blocking of the radar beam at 0.5° elevation in the southeast of Nanchang.
The above analyses show that the high frequencies (> 0.2) inFig. 10a are mainly induced by isolated non-meteorological echoes and ground clutters (near mountain regions). Clearly, the ground clutters in the mountain regions characterized by high frequency are the important part of the observed echoes. Meanwhile, the isolated non-meteorological echoes in other regions with high frequency are also important. In order to quantitatively evaluate the performance of the improved algorithms, the ratio of removal of isolated non-meteorological echoes and ground clutters is defined as
${P_{\rm {rm}}} = 1 - \frac{{{A_{\rm {ha}}}}}{{{A_{\rm {hb}}}}}, $ | (8) |
where
The present study proposes improvements to two algorithms commonly used for removing isolated non-meteorological echoes and ground clutters from single-polarization, level-II radar data. In view of the shortcomings of the original methods, which can remove the edges of meteorological echoes but cannot identify ground clutters appearing at two adjacent elevation angles, the current algorithms offer improvements in these respects. To remove isolated non-meteorological echoes as much as possible whilst preserving the edges of meteorological echoes, the improved algorithm introduces a constraint parameter, Po, to determine whether a window of 5 × 5 grid points is in isolation. To identify the ground clutters appearing at two adjacent elevation angles without range limitation, the improved algorithm establishes a new vertical difference of reflectivity (NDZ) based on the parameters RGDZ and VDZ in the WSR-88D system. Correspondingly, a set of criteria are proposed for the improved algorithms to remove isolated non-meteorologi-cal echoes and ground clutters.
A statistical analysis of 5 cases (each with 30 scans, making a total of 150 scans) with different echo types (3 small-scale convection cases, 1 typhoon case, and 1 large-scale synoptic system case) demonstrates that a constraint parameter of Po ≤ 0.167 is suitable for the improved algorithm to remove isolated non-meteorological echoes as much as possible whilst at the same time preserving the edges of meteorological echoes. The improved algorithms with Po ≤ 0.167 and NDZ ≥ 20 dBZ are tested in 4 cases (3 for isolated non-meteorological echoes and 1 for ground clutter) and the statistics of 232 scans of radar data (temporal resolution of 1 h) measured at Nanchang from 0000 UTC 5 to 1600 UTC 14 March 2015. The results show that the improved algorithms not only remove most (> 95% under clear-sky conditions) of the isolated non-meteorological echoes and ground clutters (including that occurring at two adjacent elevation angles), but also preserve the structure of storms well.
Although the improved algorithms have shown promise in removing isolated non-meteorological echoes and ground clutters, certain residual non-meteorological echoes (especially for ground clutter) still exist, such as those in Figs. 10d, 11 (under clear-sky conditions, the ratios of removal of ground clutters are sometimes less than 100%). This indicates that there is still some room for the algorithms to be further improved. Particularly, other independent data, such as radar radial velocity, radar spectrum width, etc., should be introduced to jointly improve the algorithms even further.
Acknowledgments. We thank the two anonymous reviewers and the editor for their constructive comments, which have certainly helped improve the manuscript from its original version. Dr. Zhiqun Hu provided useful comments and suggestions regarding the received base-band complex signal at the signal-processing level.
Alfieri, L., P. Claps, and F. Laio, 2010: Time-dependent Z–R relationships for estimating rainfall fields from radar measurements . Nat. Hazards Earth Sys., 10, 149–158. DOI:10.5194/nhess-10-149-2010 |
Bedka, K. M., C. Wang, R. Rogers, et al., 2015: Examining deep convective cloud evolution using total lightning, WSR-88D, and GOES-14 super rapid scan datasets. Wea. Forecasting, 30, 271–590. DOI:10.1175/WAF-D-14-00062.1 |
Anagnostou, E. N., and W. F. Krajewski, 1999: Real-time radar rainfall estimation. Part I: Algorithm formulation. J. Atmos. Oceanic Technol., 16, 189–197. DOI:10.1175/1520-0426(1999)016<0189:RTRREP>2.0.CO;2 |
Bergen, W. R., and S. C. Albers, 1988: Two- and three-dimensional de-aliasing of Doppler radar velocities. J. Atmos. Oceanic Technol., 5, 305–319. DOI:10.1175/1520-0426(1988)005<0305:TATDDA>2.0.CO;2 |
Chrisman, J., D. Rinderknecht, and R. Hamilton, 1995: WSR-88D clutter suppression and its impact on meteorological data interpretation. Preprints, First WSR-88D User’s Conference, Norman, OK, WSR-88D Operational Support Facility, USA, 11–14 October, 9–20. |
Fabry, F., and J. Gadoury, 2009: Another method for ground target identification and filtering using spectral processing of dual-polarization returns. 34th Conference on Radar Meteorology, Williamsburg, VA, 6 October, Amer. Meteor. Soc. |
Friedrich, K., M. Hagen and T. Einfalt, 2006: A quality control concept for radar reflectivity, polarimetric parameters, and Doppler velocity. J. Atmos. Oceanic Technol., 23, 865–887. DOI:10.1175/jtech1920.1 |
Fulton, R. A., J. P. Breidenbach, D. J. Seo, et al., 1998: The WSR-88D rainfall algorithm. Wea. Forecasting, 13, 377–395. DOI:10.1175/1520-0434(1998)013<0377:TWRA>2.0.CO;2 |
Gourley, J. J., P. Tabary, and J. Parent 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. DOI:10.1175/jtech2035.1 |
Hubbert, J. C., M. Dixon, S. M. Ellis, et al., 2009a: Weather radar ground clutter. Part I: Identification, modeling, and simulation. J. Atmos. Oceanic Technol., 26, 1165–1180. DOI:10.1175/2009jtecha1159.1 |
Hubbert, J. C., M. Dixon, and S. M. Ellis, 2009b: Weather radar ground clutter. Part II: Real-time identification and filtering. J. Atmos. Oceanic Technol., 26, 1181–1197. DOI:10.1175/2009jtecha1160.1 |
Jones, T. A., D. Stensrud, L. Wicker, et al., 2015: Simultaneous radar and satellite data storm-scale assimilation using an ensemble Kalman filter approach for 24 May 2011. Mon. Wea. Rev., 143, 165–194. DOI:10.1175/mwr-d-14-00180.1 |
Kessinger, C. K., S. Ellis, J. Van Andel, et al., 2003: The AP clutter mitigation scheme for the WSR-88D. Preprints, 31st Conference on Radar Meteorology, Seattle WA, 6–12 August, Amer. Meteor. Soc., 526–529. |
Krishnapuram, R., and J. M. Keller, 1993: A possibilistic approach to clustering. IEEE Transactions on Fuzzy Systems, 1, 98–110. DOI:10.1109/91.227387 |
Lakshmanan, V., C. Karstens, J. Krause, et al., 2014: Quality control of weather radar data using polarimetric variables. J. Atmos. Oceanic Technol., 31, 1234–1249. DOI:10.1175/jtech-d-13-00073.1 |
Liu, L., M. X. Chen, and Y. C. Wang, 2016: Numerical nowcasting experiments for the simulation of a mesoscale convective system using a cloud model and radar data assimilation with 4DVar. Acta Meteor. Sinica, 74, 213–228. DOI:10.11676/qxxb2016.021 |
Liu, L. P., L. L. Wu, and Y. M. Yang, 2007: Development of fuzzy-logical two-step ground clutter detection algorithm. Acta Meteor. Sinica, 65, 252–260. DOI:10.3321/j.issn:0577-6619.2007.02.011 |
Liu, S.-Y., L.-F. Yan, and J. Sun, 2008: Retrieval of horizontal wind and quality controlling of single Doppler radar data on landing typhoon. J. Trop. Meteor., 24, 105–110. DOI:10.3969/j.issn.1004-4965.2008.02.001 |
Pamment, J. A., and B. J. Conway, 1998: Objective identification of echoes due to anomalous propagation in weather radar data. J. Atmos. Oceanic Technol., 15, 98–113. DOI:10.1175/1520-0426(1998)015<0098:oioedt>2.0.co;2 |
Pan, Y. J., K. Zhao, and Y. N. Pan, 2010: Single-Doppler radar observations of a high precipitation supercell accompanying the 12 April 2003 severe squall line in Fujian province. Acta Meteor. Sinica, 24, 50–65. |
Qi, Y. C., and J. Zhang, 2017: A physically based two-dimensional seamless reflectivity mosaic for radar QPE in the MRMS system. J. Hydrometeor., 18, 1327–1340. DOI:10.1175/jhm-d-16-0197.1 |
Qin, Y. Y., J. D. Gong, Z. C. Li, et al, 2014: Assimilation of Doppler radar observations with an ensemble square root filter: A squall line case study. J. Meteor. Res., 28, 230–251. DOI:10.1007/s13351-014-2046-6 |
Rico-Ramirez, M. A., and I. D. Cluckie, 2008: Classification of ground clutter and anomalous propagation using dual-polarization weather radar. IEEE Trans. Geosci. Remote Sens., 46, 1892–1904. DOI:10.1109/tgrs.2008.916979 |
Serafin, R. J., and J. W. Wilson, 2000: Operational weather radar in the United States: Progress and opportunity. Bull. Amer. Meteor. Soc., 81, 501–518. DOI:10.1175/1520-0477(2000)081<0501:owritu>2.3.CO;2 |
Ridal, M., and M. Dahlbom, 2017: Assimilation of multinational radar reflectivity data in a mesoscale model: A proof of concept. J. Appl. Meteor. Climatol., 56, 1739–1751. DOI:10.1175/jamc-d-16-0247.1 |
Steenburgh, W. J., S. F. Halvorson, and D. J. Onton, 2000: Climatology of lake-effect snowstorms of the Great Salt Lake. Mon. Wea. Rev., 128, 709–727. DOI:10.1175/1520-0493(2000)128<0709:coleso>2.0.co;2 |
Steiner, M., and J. A. Smith, 2002: Use of three-dimensional reflectivity structure for automated detection and removal of nonprecipitating echoes in radar data. J. Atmos. Oceanic Technol., 19, 673–686. DOI:10.1175/1520-0426(2002)019<0673:uotdrs>2.0.co;2 |
Su, T. J., J. X. Ge, and H. B. Zhang, 2018: Review of dual polarization weather radar system development in China. J. Mar. Meteor., 38, 62–68. DOI:10.19513/j.cnki.issn2096-3599.2018.01.008 |
Tang, F., 2011: Quality control of radar reflectivity data and its application to GRAPES 3DVAR. Master dissertation, Nanjing University of Information Science & Technology, China, 48 pp. (in Chinese) |
Torres, S. M., and D. A. Warde, 2014: Ground clutter mitigation for weather radars using the autocorrelation spectral density. J. Atmos. Oceanic Technol., 31, 2049–2066. DOI:10.1175/jtech-d-13-00117.1 |
Wang, Y. B., and Y. F. Wan, 2006: Automatic quality control for radar volume-scanning reflectivity fields. Meteor. Sci. Tech-nol., 34, 615–619. DOI:10.3969/j.issn.1671-6345.2006.05.021 |
Wu, T., Y. F. Wan, W. F. Wo, et al., 2013: Design and application of radar reflectivity quality control algorithm in SWAN. Meteor. Sci. Technol., 41, 809–817. DOI:10.3969/j.issn.1671-6345.2013.05.004 |
Xiao, Y. J., L. P. Liu, and Y. Shi, 2008: Study of methods for three-dimensional multiple-radar reflectivity mosaics. Acta Meteor. Sinica, 22, 351–361. |
Xu, M. Y., F. Li, Y. C. Xia, et al., 2017: Analysis of CINRAD radar operation status during 2009–2014. Meteor. Mon., 43, 365–372. DOI:10.7519/j.issn.1000-0526.2017.03.013 |
Yang, M. L., L. P. Liu, D. B. Su, et al., 2011: The application of an automated 2-D multi-pass Doppler radar velocity dealiasing and the research of its effect. Meteor. Mon., 37, 203–212. DOI:10.7519/j.issn.1000-0526.2011.2.010 |
Yu, X. D., X. P. Yao, T. N. Xiong, et al., 2006. Doppler Weather Radar Principles and Operational Applications. China Meteorological Press, Beijing, 25 pp. (in Chinese) |
Zhang, J., S. Wang, and B. Clarke, 2004: WSR-88D reflectivity quality control using horizontal and vertical reflectivity structure. Preprints, 11th Conference on Aviation, Range and Aerospace Meteorology, Hyannis, MA, 4–8 October, Amer. Meteor. Soc. |
Zhang, J., K. Howard, and J. J. Gourley, 2005: Constructing three-dimensional multiple-radar reflectivity mosaics: Examples of convective storms and stratiform rain echoes. J. Atmos. Oceanic Technol., 22, 30–42. DOI:10.1175/jtech-1689.1 |
Zheng, Y. Y., X. D. Yu, C. Fang, et al., 2004: Analysis of a strong classic supercell storm with Doppler weather radar data. Acta Meteor. Sinica, 62, 317–328. DOI:10.3321/j.issn:0577-6619.2004.03.006 |