Next Article in Journal
Applications of Cell-Ratio Constant False-Alarm Rate Method in Coherent Doppler Wind Lidar
Next Article in Special Issue
Ensemble Classification for Anomalous Propagation Echo Detection with Clustering-Based Subset-Selection Method
Previous Article in Journal
Chemical Composition of PM10 at Urban Sites in Naples (Italy)
Previous Article in Special Issue
Detection of Ground Clutter from Weather Radar Using a Dual-Polarization and Dual-Scan Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Rain Attenuation Correction of Reflectivity for X-Band Dual-Polarization Radar

1
Key Laboratory of Cloud-Precipitation and Severe Storms, Center of Disaster Reduction, Institute of Atmospheric Physics, Chinese Academy of Sciences, Huayanli No.40, Qijiahuozi, Deshengmenwai Avenue, Beijing 100029, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Tianjin Meteorological Information Center, Tianjin 300074, China
4
Dalian Meteorological Equipment Support Center, Dalian 116001, China
*
Author to whom correspondence should be addressed.
Atmosphere 2016, 7(12), 164; https://doi.org/10.3390/atmos7120164
Submission received: 19 September 2016 / Revised: 6 December 2016 / Accepted: 8 December 2016 / Published: 17 December 2016
(This article belongs to the Special Issue Radar Meteorology)

Abstract

:
In order to improve the performance of X-band dual-polarization radars, it is necessary to conduct attenuation correction before using the X-band radar data. This study analyzes a variety of attenuation correction methods for the X-band radar reflectivity, and proposes a high-resolution slide self-consistency correction (SSCC) method, which is an improvement of Kim et al.’s method based on Bringi et al.’s original method. The new method is comprehensively evaluated with the observational data of convective cloud, stratiform cloud, and the stratiform cloud with embedded convection. Comparing with the intrinsic reflectivity at X-band calculated from the reflectivity at S-band, it is found that the new method can effectively reduce the correction errors when calculating differential propagation shift increments using the conventional self-consistency attenuation correction method. This method can efficiently correct the X-band dual-polarization radar reflectivity, in particular, for the echoes with reflectivity greater than 35 dBZ.

1. Introduction

Comparing with the conventional Doppler weather radar, dual-polarization radars can measure more valuable polarized information precipitation systems. The polarized information allows to improve the accuracy of radar-based quantitative precipitation estimation, raindrop size distribution (DSD) retrieval and precipitation particle identification [1,2,3,4,5,6,7,8,9,10]. Previous studies on the application of dual-polarization radars are mostly designed for S, C-band radars [11,12,13,14,15]. Research on X-band dual-polarization radars is limited, since X-band radars experience severe attenuation compared to S, C-band radars. However, due to their low cost, small antennas, easy mobility and high temporal and spatial resolution, X-band dual-polarization radars have become an important detection equipment in the areas of cloud and precipitation physics and weather modification. In order to improve the performance of X-band dual-polarization radars, the attenuation needs to be corrected before application.
Atlas and Banks [16] showed there were two main factors resulting in attenuation. One is detection range, whereby the echo power received by the radar will decrease with increasing range, and this applies to all radar wavelengths; the other is rain attenuation. With the exception of intense storms, rain attenuation for electromagnetic waves with a wavelength greater than approximately 7 cm is negligible. Zhang et al. [17] noted that the first attenuation factor is mainly because gas molecules absorb electromagnetic waves and the influence of scattering can be neglected. The absorption attenuation of electromagnetic waves with a wavelength greater than 2 cm is generally small, but when the wavelength is approximately 1 cm or the detection range is large, the attenuation caused by the first factor still needs to be considered. Therefore, the attenuation at X-band (about 3 cm in wavelength) is mainly due to the second factor.
The conventional method of attenuation correction for single-polarization radars is mostly based on the empirical relationship between horizontal reflectivity factor ZH and rainfall intensity R(ZH = aRb; a, b are empirical constants). This method retrieves Zret using the measured rainfall intensity R and then calculates horizontal specific attenuation coefficient AH (AH = ZretZH) [18]. However, the relationship between ZH and R is not stable, depending on not only different locations, seasons and precipitation patterns, but also precipitation process and time. It is mainly because of the variability of drop size distributions (DSDs). Meanwhile, the empirical relationship is also affected by radar calibration and beam blockage. Therefore, it is not accurate to correct rain attenuation by using this method.
Dual-polarization radars can avoid the shortcomings of the single-polarization radar attenuation correction method, because they can provide differential propagation phases (φDP) and specific differential phases (KDP). The two parameters are independent of radar calibration, rain attenuation and partial beam blockage. Therefore, dual-polarization radars can provide a stable rain attenuation correction relationship using φDP and KDP. Bringi et al. [19] found that there was almost a linear relationship (AH = αHKDP) between the attenuation (AH) and specific differential phase (KDP) by scattering simulation. Zrnic and Ryzhkov [20] pointed out that KDP was unaffected by attenuation and relatively immune to the beam blockage. Based on this fact, Ryzhkov and Zrnic [21] proposed an empirical correction method, where the coefficients were determined as a mean slope between φDP and ZH or differential reflectivity (ZDR) in a sampling area. Their method was evaluated with the S-band dual-polarization radar data and improved for the C-band dual-polarization radar data by Carey et al. [22]. He et al. [23] adopted this correction method and introduced Kalman filter for filtering the measured φDP, then obtained the relation coefficient α’H between AH and φDP, finally corrected the stratiform case, which was detected by an X-band dual-polarization radar. Although Carey et al. [22] and He et al. [23] have greatly improved the method of Ryzhkov and Zrnic [21], the method is only applied to stratiform precipitation. Hu et al. [24] compared the correction method by KDP with the convectional correction method and found that the correction by KDP was better than by ZH. However, the KDP correction method would cause errors since the KDP may contain errors when rainfall intensity is small. Thus he proposed a comprehensive ZH-KDP correction method to overcome shortcomings of the correction by ZH or by KDP. However, ZH-KDP correction method still uses a fixed coefficient to correct rain attenuation.
Smyth and Illingworth [25] introduced a constraint to correct the differential reflectivity (ZDR) for S-band dual-polarization radars. In this method, the coefficient (αH) of the relationship between KDP and ADP is not fixed, but determined by the constraint that ZDR at the edge of a rain cell should be 0 dB (assuming that edge of the rain cell is drizzle). However, this is not applicable in some cases, particularly for the shorter wavelengths, high-resolution X-band dual-polarization radar observations. Due to rain attenuation, the rain edge of the radar display is not necessarily the actual edge of the rain cell. The rain edge may be drizzle, moderate or even heavy rain. Thus, it is inappropriate to set ZDR as 0 dB at the farther edge of a rain zone. It is necessary to create a new ZDR constraint according to the actual situation. Testud et al. [26] proposed an attenuation correction method, called ZPHI method. The core idea assumes that the differential propagation phase calculated by AH should be equal to the increments of the measured radial differential propagation phase. This method achieves a better performance, but still needs to set the coefficient αH of the relationship between AH and KDP.
Bringi et al. [3] proposed an algorithm referred to as “the self-consistent method with constraints”, which can resolve the limitations of Smyth and Illingworth [25] and Testud et al. [26]. The algorithm improved the method of Testud et al. [26] for ZH correction and the method of Smyth and Illingworth [25] for ZDR correction. One of the advantages of the algorithm is that the coefficient αH of the relationship between AH and KDP is estimated from the radar data rather than scattering simulation. Park et al. [27,28] extended the algorithm to the X-band dual polarization radar, and calculated the range of αH by scattering simulation using drop size distributions. Kim et al. [29] corrected ZDR by the horizontal reflectivity ZH and the vertical reflectivity ZV using the method of Bringi et al. [3]. The method turned the range resolution into 1.5 km. Later, Kim et al. [4] improved the resolution to 0.5 km further. For stratiform cloud and the stratiform cloud with embedded convection, a resolution of 0.5 km may be appropriate, because the KDP is not large in the two kinds of cloud for X-band dual polarization radars. However, it is large for convective cloud, for example, the KDP can reach 10°/km or more in convective cores. Such a resolution may result in errors when correcting convective cloud. In addition, φDP would be used to correct ZH and ZV in the method of Bringi et al. [3]. This method needs to seek an initial phase and a terminal phase for every radial in the corrected process. This may result in phase errors due to radar system noise and finally result in correction errors.
In this paper, we propose a high-resolution slide self-consistency correction method to improve the method of Bringi et al. [3]. The new algorithm applies a slide window to avoid seeking the initial and terminal phases. The accuracy of the correction results is evaluated with convective cloud, stratiform cloud and the stratiform cloud with embedded convection by comparing with the intrinsic reflectivity at X-band, which is calculated from the reflectivity at S-band.

2. Radar Feature

The IAP-714XDP-A mobile dual-polarization weather radar has been operated since 2006 by the Key Laboratory of Cloud-Precipitation Physics and Severe storms (LACS), Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences (CAS). The signal processor of the radar is RVP8. The scanning strategy includes plane position indicator (PPI), radar height indicator (RHI) and volume coverage pattern (VCP). The main specifications of the IAP-714XDP-A mobile radar system are listed in Table 1.

3. The Slide Self-Consistency Correction Method

The Slide Self-Consistency Correction (SSCC) method is mainly based on the self-consistent method with constraints proposed by Bringi et al. [3]. The radar reflectivity factor Zh (mm6m−3) in linear scale and ZH (dBZ) in logarithm scale have the following relationship:
Z H = 10 lg Z h
The corrected reflectivity Z H A (dBZ) at a range r is related to the attenuated (measured) reflectivity Z H as follows:
Z H A ( r ) = Z H ( r ) + 2 0 r A H ( s ) d s
where AH is specific attenuation in decibels per kilometer. The change in differential propagation phase is as follows:
Δ ϕ D P = ϕ D P ( r 1 ) ϕ D P ( r 0 )
where r0 and r1 are the beginning and ending range gate, respectively.
In Bringi et al. [3], specific attenuation AH is determined with a constraint that the cumulative attenuation from range r0 to r1 should be consistent with the total change in differential propagation phase Δ ϕ D P . Under the assumption that there is a linear relationship between AH and KDP, the final formula of AH is given by
A H ( r ) = [ Z h ( r ) ] b × [ 10 0.1 ( b α ) Δ ϕ D P 1 ] I ( r 0 , r 1 ) + [ 10 0.1 ( b α ) Δ ϕ D P 1 ] × I ( r , r 1 )
where,
I ( r 0 , r 1 ) = 0.46 b r 0 r 1 [ Z h ( s ) ] b d s
I ( r , r 1 ) = 0.46 b r r 1 [ Z h ( s ) ] b d s
In the above equations, α and b are the empirical parameters of the following Equations (8) and (7) that can be obtained by scattering simulation by raindrop size distribution. Bringi et al. [3,19] found an exponent relationship between AH and Zh and a linear relationship between specific attenuation AH and specific differential phase KDP at frequencies from 2.8 to 9.3 GHz; that is, the exponent c in Equation (8) is close to 1.
A H = a Z h b
A H = α K D P c
where KDP is in °·km−1 and c is set as a constant 1.
Therefore, if AH(r) is calculated by Equation (4) and substituted into Equation (2), the corrected reflectivity ZHA(r) at a range r is obtained. However, α and b need to be set to a fixed value before calculating AH(r). Carey et al. [22] noted that the coefficient α can vary widely with temperature and drop shape. Park et al. [27] found that it changes from 0.139 to 0.335 dB(°)−1 at X-band. Comparing with α, the exponent b is less influenced. Delriu et al. [30] found that b varies from 0.76 to 0.84 at X-band. Thus, in this paper, b is set as a constant 0.8.
When calculating AH(r) using a fixed α value, the correction errors are introduced in the process of attenuation correction. To eliminate the impact of the α variability, Bringi et al. [3] proposed a self-consistent method with constraints. This method does not set α as a fixed constant, but seek an optimal α within a predetermined scope (αmin, αmax), which is obtained from scattering simulation under various temperatures and raindrop size distributions.
For each α, A H ( r ; α ) at each range is calculated by Equation (4), and then ϕ D P c a l ( r ; α ) is calculated as follows:
ϕ D P c a l ( r ; α ) = 2 r 0 r 1 A H ( s ; α ) α d s
The optimal α is the value that leads Equation (10) to the minimum.
ϕ D P e r r o r ( α ) = i = 1 N | ϕ D P c a l ( r i ; α ) ϕ D P ( r i ) |
where i denotes the range gate from r0 to r1. The main advantage of the self-consistent method with constraints is estimating an optimal α rather than setting a fixed value.
According to the scattering simulation results at X-band by Park et al. [27], α is set between 0.1 and 0.5, in a step of 0.03. Kim et al. [4] sets the distance between r0 and r1 as 1.0 km with an overlap of 0.5 km (ultimately, α has a resolution of 0.5 km), referring to Figure 1a. In the paper, the SSCC method employs a slide window processing shown in Figure 1b by setting the distance between r0 and r1 as 1.5 km (10 gates), thus α has a high-resolution of 0.15 km, improving the resolution of α estimation. After developing the method of Bringi et al. [3], Equations (3) and (10) become as follows:
Δ ϕ D P _ 10 g a t e s = ϕ D P ( G i + 10 ) ϕ D P ( G i )
ϕ D P _ 10 g a t e s e r r o r ( α ) = i = 0 9 | ϕ D P c a l ( r i ; α ) ϕ D P ( r i ) |
In addition, the new method does not seek the initial and terminal phase of each radial as the method of Bringi et al. [3] does, which would finally cause correction errors.

4. Result

The SSCC method is evaluated using the data collected by the X-band dual-polarization radar (IAP-714XDP-A), which is located in Shunyi District of Beijing City (BJ) from June to September 2015. The data contain observations of convective cloud, stratiform cloud and the stratiform cloud with embedded convection. The corrected reflectivity is compared with the intrinsic reflectivity at X-band calculated from the reflectivity at S-band, which is obtained by the CINRAD/SA S-band single-polarization weather radar located in Daxing District of Beijing City. Figure 2 shows the locations of the two radars. The X-band radar (at ShunyiSY, 116.68° E, 40.19° N) is located at the northeast of the S-band radar (at DaxingDX, 116.47° E, 39.81° N). The straight-line distance between the two radars is about 46 km.
Hubbert and Bringi [31] pointed out that the measured differential propagation phase ΨC consists of true differential propagation phase φDP and backward scattering differential phase shift δ. Since δ could result in errors of KDP estimation, δ needs to be eliminated before using ΨC. In this paper, the method of Hubbet and Bringi [31] is applied to filter δ out.
A strong convective weather event swept Beijing City from north to south on 19 June 2015. Figure 3 shows the plane position indicator (PPI) of the X-band radar at elevation 3° at 14:45 Beijing Time (BJT), 19 June 2015. Three strong echo cores of convective cloud are situated between the two radars and the largest reflectivity observed by the S-band radar is about 60 dBZ.
The intrinsic reflectivity at X-band is not equal to the reflectivity at S-band. Chandrasekar et al. [32] proposed three different methodologies for simulating X-band radar observations from the S-band radar data and the empirical conversion method is used in the paper. Figure 4 shows a plot of the intrinsic reflectivity at X and S bands for a monodispersed drop size distribution using the shape mode proposed by Beard and Chuang [33]. The relationship between the intrinsic reflectivity at X-band and the reflectivity at S-band is obtained by curve fitting, which is divided into three parts as shown in Equation (13) where subscripts X and S indicate simulated radar variables at X-band and measured radar measurements at S-band. Note that the reflectivity at S-band is assumed to be non-attenuated.
Z H , X = { 0.9696 Z H , S 0.0145 Z H , S 25 d B Z 1.1982 Z H , S 5.7726 25 d B Z < Z H , S < 45 d B Z 0.8206 Z H , S + 11.7934 Z H , S 45 d B Z
In order to analyze the accuracy of the SSCC method, the corrected X-band radar reflectivity is compared with the intrinsic reflectivity at X-band, which is calculated from the reflectivity at S-band. The S-band radar reflectivity from the volume scan data is interpolated into the coordinate of the X-band radar. The X-band radar PPI is shown in Figure 5a. Figure 5b is the intrinsic reflectivity at X-band. As shown in Figure 5a, there is a strong convective cloud band with three strong echo cores in the southwest of the X-band radar. Due to severe attenuation, the X-band radar cannot observe the echo after the intensive rain region, which is shown in circle A in Figure 5b. The echo of the circle B is also not be detected by the X-band radar, resulting from partial beam blockage.
The shapes of the two radar echoes are similar to each other, but the X-band radar reflectivity is seriously attenuated. The maximum reflectivity of convective cores at X-band is about 50 dBZ, while the corresponding intrinsic reflectivity is about 60 dBZ, indicating that the X-band radar echo has a serious distortion due to rain attenuation.
The reflectivity at X-band is corrected by the SSCC method, which is shown in Figure 6. Compared with the uncorrected X-band radar reflectivity in Figure 5a, the corrected reflectivity is effectively compensated and the scope of strong echois extended. For further analysis, the reflectivity is mapped into a 1000 × 1000 matrix grids with a resolution of 150 m.
Figure 7 shows the scatter diagrams of the uncorrected and corrected X-band radar reflectivity versus the intrinsic reflectivity. Figure 7a shows that the uncorrected reflectivity significantly deviates from the intrinsic reflectivity, especially when the echoes are strong. The fitting line between the uncorrected reflectivity and the intrinsic reflectivity (the green line) is y = 0.5824x + 9.2234, while the fitting line between the corrected reflectivity and the intrinsic reflectivity becomes y = 0.8036x + 3.8382, referring to Figure 7b. After attenuation correction, the slope of fitting line turns 0.5824 into 0.8036, indicating that the corrected reflectivity is much closer to the intrinsic reflectivity.
In order to show the validity of the correction for any path, the reflectivity at X-band with an azimuth angle of 228° is analyzed herein. As shown in Figure 3, the electromagnetic wave successively passes from A to D. Due to the impact of distance, antenna elevation and earth curvature, there are not echoes of the S-band radar in the area A, referring to Figure 8a. Compared with the intrinsic reflectivity (Intrinsic), the uncorrected reflectivity (UnC) nearly has no attenuation in areas A, B and C while with serious attenuation in area D. Figure 8b also illustrates this phenomenon, the φDP increases by nearly 50°, corresponding to intensive rain region, while no increase in areas A, B and C. After attenuation correction using the SSCC method, the corrected X-band radar reflectivity at 33 km has compensated about 10 dBZ (SSCC). The corrected reflectivity at X-band is consistent with the intrinsic reflectivity. However, the corrected reflectivity using the method of Kim et al. [4] (Kim) is lower than the intrinsic reflectivity, indicating the correction is not enough.
The corrected X-band radar reflectivity at 37 km in Figure 8a is 15 dBZ larger than the intrinsic reflectivity. This results from the rapid development and fast moving speed of the convective cloud. Because the intrinsic reflectivity shown in Figure 5b is interpolated by the 6-min volume scan data of the S-band radar, the two factors causes a slight deviation from the intrinsic reflectivity. This influence is significant at the edge of convective cloud but negligible for the stratiform cloud and the stratiform cloud with embedded convection. As shown in Figure 8a, the corrected X-band radar reflectivity both the SSCC and the Kim is close to the intrinsic reflectivity. However, the SSCC method has an advantage over Kim et al. [4] in correcting the reflectivity of the convective cloud. In order to analyze the two methods comprehensively, all the radials are used.
Figure 9 shows four cumulative distributions of the radar reflectivity. Comparing with the cumulative distribution of the uncorrected reflectivity (UnC), the method of Kim et al. [4] (Kim), and the SSCC method (SSCC) both shift to the right, indicating that the low cumulative value of reflectivity decreases, while the high cumulative value increases. Both the cumulative distribution of the Kim and SSCC are closer to that of the intrinsic reflectivity than the UnC. The average biases (AB) of the reflectivity are calculated for the two methods.
The average bias (AB) shown in Figure 10 is defined as below:
A B = R R s | x
where | x is average value above parameter x, R is the reflectivity at X-band and RS is the intrinsic reflectivity. As shown in Figure 10a, the AB between the uncorrected reflectivity and the intrinsic reflectivity (line UnC) is decreasing with increasing reflectivity. The AB of the UnC is greater than 10 dB, indicating that the attenuation is significantin the rain area. The AB of the Kim and SSCC significantly reduces the difference from the intrinsic reflectivity. To accurately retrieve meteorological products, a resolution of 1 dB for the reflectivity is necessary. Figure 10b shows that there are more than 1 dB differences between the Kim and SSCC from 35 dBZ, illustrating that the SSCC method has a better performance than Kim et al.’s method at correcting convective cloud, especially with reflectivity greater than 35 dBZ.
In order to analyze the impact of different sampling resolutions for the SSCC method, the range resolution is set at 0.45 km (SSCC_450) and 0.75 km (SSCC_750) as shown in Figure 11, respectively. The SSCC_750 is further away from the Intrinsic than the SSCC_450, which is closer to the method by Kim et al. [4] (the range resolution is 0.5 km). Figure 11 shows that the decreasing range resolution would lead to reduced correction effect and the SSCC method performs better than the method by Kim et al. [4] for convective cloud.
Compared with the method by Bringi et al. [3], the SSCC method does not require an initial and a terminal differential propagation phase of each radial, which could avoid correction errors. To illustrate this problem, we assume the terminal phase is true and examine the errors due to the wrong initial phase. Figure 12 shows correction errors with various initial phases using the method by Bringi et al. [3], whereby 275° (Bringi_275) is the actual radar initial phase and 255° (Bringi_255) and 265° (Bringi_265) are not. The SSCC is consistent with the Bringi_275. In contrast, the Bringi_255 and the Bringi_265 are far away from the Bringi_275, indicating the SSCC method does not require seeking the initial phase and terminal phase but the cumulative distribution is also consistent with that of the true correction results.
To verify the applicability of the SSCC method under various precipitation conditions, the stratiform cloud with embedded convection on 26 June 2015 (Figure 13 and Figure 14), and the stratiform cloud on 16 June 2015 (Figure 15 and Figure 16) are analyzed. Both Figure 17 and Figure 18 show that the reflectivity at X-band is corrected effectively and the corrected cumulative distribution closer to that of the intrinsic reflectivity. Figure 19 and Figure 20 show that the SSCC method is consistent with the method by Kim et al. [4]. The change of the resolution does not lead to correction biases, because the KDP of the stratiform cloud with embedded convection and the stratiform cloud is lower than that of the convective cloud. Figure 21 and Figure 22 illustrate that the corrected cumulative distribution using the SSCC method are consistent with that of the 275°, which is the true initial phase of the radar. The correction verification of the three different precipitation cases indicates that the SSCC method is also applicable for the stratiform cloud and the stratiform cloud with embedded convection. Note that both the SSCC method and the method by Kim et al. [4] may have no significant effect or lead to slightly worse attenuation correction due to the error of the integration resolution when correcting stratiform rain if the rainfall is very small.

5. Conclusions and Discussion

Based on Bringi et al. [3], the paper proposed a high-resolution slide self-consistency correction (SSCC) method for the X-band dual-polarization radar reflectivity, which is an improvement from Kim et al.’s method. The proposed method improved the correction resolution and effect, adapting a slide window consisting of 10 gates.
In the paper, the SSCC method is evaluated with the reflectivity of the convective cloud, the stratiform cloud with embedded convection and the stratiform cloud, comparing with the correction results from the methods by Bringi et al. [3] and Kim et al. [4], as well as the intrinsic reflectivity at X-band calculated from the reflectivity at the S-band. It is found that the reflectivity at X-band can be corrected effectively by the SSCC method. The corrected reflectivity is closer to the intrinsic reflectivity and has a better performance than the method by Kim et al. [4] in correcting the convective cloud. However, the correction results of the two methods are very similar for the stratiform cloud with embedded convection and the stratiform cloud. This may be because the KDP of the two kinds of precipitation cloud is much less than that of the convective cloud. For this reason, the SSCC method and the method by Kim et al. [4] may have no significant effect or may lead to slightly worse attenuation correction when correcting stratiform rain if the rainfall is very small. In addition, the SSCC method has better results than Bringi et al. [3] for the three cases due to the reduced correction errors when computing differential propagation shift increments.
In summary, the SSCC method has three advantages as follows:
  • Improving the correction resolution;
  • Having no need for seeking the initial and terminal differential phases;
  • Good performance in correcting convective cloud.
Meanwhile, it must be noted that the accuracy of the attenuation correction is restricted in the SSCC method by the length (1.5 km) of the sliding window and this is a more significant effect than the resolution.

Acknowledgments

This work was partially supported by the National Key Basic Research 973 Program of China (Grant Nos. 2014CB441403, 2013CB430105), the National Natural Science Foundation of China (Grant Nos. 41575037, 41605019), Guizhou Province Scientific Research Joint Project (Grant No. G [2013]4001), and the Special Scientific Research Project of Meteorological Public Welfare Profession of China (Grant No. GYHY201006031).

Author Contributions

The authors contribute equally to this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, L.P.; Qian, Y.F. The study of spacial distribution of phase and size of hydrometeors in cloud by dual linear polarization radar. Acta Meteorol. Sin. 1996, 54, 590–599. [Google Scholar]
  2. Zrnic, D.S.; Ryzhkov, A.V. Polarimetry for weather surveillance radars. Bull. Am. Meteorol. Soc. 1999, 80, 389–406. [Google Scholar] [CrossRef]
  3. Bringi, V.N.; Keenan, T.D.; Chandrasekar, V. Correcting C-band radar reflectivity and differential reflectivity data for rain attenuation: A self-consistent method with constraints. IEEE Trans. Geosci. Remote Sens. 2001, 39, 1906–1915. [Google Scholar] [CrossRef]
  4. Kim, D.-S.; Maki, M.; Lee, D.-I. Retrieval of three-dimensional raindrop size distribution using X-band polarimetric radar data. J. Atmos. Ocean. Technol. 2010, 27, 1265–1285. [Google Scholar] [CrossRef]
  5. Snyder, J.C.; Bluestein, H.B.; Zhang, G.F. Attenuation correction and Hydrometeor classification of high-resolution, X-band, Dual-polarized mobile radar measurements in severe convective storms. J. Atmos. Ocean. Technol. 2010, 27, 1979–2001. [Google Scholar] [CrossRef]
  6. Thompson, E.J.; Rutledge, S.A.; Dolan, B. A dual-polarization radar hydrometeor classification algorithm for winter precipitation. J. Atmos. Ocean. Technol. 2014, 31, 1457–1481. [Google Scholar] [CrossRef]
  7. Bechini, R.; Chandrasekar, V. A semisupervised robust hydrometeor classification method for dual-polarization radar applications. J. Atmos. Ocean. Technol. 2015, 32, 22–47. [Google Scholar] [CrossRef]
  8. Wen, G.; Protat, A.; May, P.T.; Wang, X.; Moran, W. A Cluster-based method for hydrometeor classification using polarimetric variables. Part I: Interpretation and analysis. J. Atmos. Ocean. Technol. 2015, 32, 1320–1340. [Google Scholar] [CrossRef]
  9. Wen, G.; Protat, A.; May, P.T.; Moran, W.; Dixon, M. A Cluster-based method for hydrometeor classification using polarimetric variables. Part II: Classification. J. Atmos. Ocean. Technol. 2016, 33, 45–60. [Google Scholar] [CrossRef]
  10. Wen, G.; Oue, M.; Protat, A.; Verlinde, J.; Xiao, H. Ice particle type identification for shallow Arctic mixed-phase cloud using X-band polarimetric radar. Atmos. Res. 2016, 182, 114–131. [Google Scholar] [CrossRef]
  11. Liu, H.P.; Chandrasekar, V. Classification of hydrometeors based on polarimetric radar measurements: Development of fuzzy logic and neuro-fuzzy systems, and in situ verification. J. Atmos. Ocean. Technol. 2000, 17, 140–164. [Google Scholar] [CrossRef]
  12. Keenan, T. Hydrometeor classification with a C-band polarimetric radar. Aust. Met. Mag. 2003, 52, 23–31. [Google Scholar]
  13. Paultsch, H.; Teschl, F.; Randeu, W.L. Dual-polarization C-band weather radar algorithms for rain rate estimation and hydrometeor classification in an Alpine region. Adv. Geosci. 2009, 20, 3–8. [Google Scholar] [CrossRef]
  14. Kumjian, M.R.; Ryzhkov, A.V.; Reeves, H.D.; Schuur, T.J. A dual-polarization radar signature of hydrometeor refreezing in winter storms. J. Appl. Meteorol. Climatol. 2013, 52, 2549–2566. [Google Scholar] [CrossRef]
  15. Junyent, F.; Chandrasekar, V. An examination of precipitation using CSU-CHILL dual-wavelength, dual-polarization radar observations. J. Atmos. Ocean. Technol. 2016, 33, 313–329. [Google Scholar] [CrossRef]
  16. Atlas, D.; Banks, H.C. The interpretation of microwave reflections from rainfall. J. Meteorol. 1951, 8, 271–282. [Google Scholar] [CrossRef]
  17. Zhang, P.C.; Du, B.Y.; Dai, T.P. Radar Meteorology; China Meteorological Press: Beijing, China, 2001; p. 511. (In Chinese) [Google Scholar]
  18. Hitschfeld, W.; Bordan, J. Errors inherent in the radar measurement of rainfall at attenuating wavelengths. J. Meteorol. 1954, 11, 58–67. [Google Scholar] [CrossRef]
  19. Bringi, V.N.; Balakrishnan, N.; Zrnic, D.S. An examination of propagation effects in rainfall on polarimetric variables at microwave frequencies. J. Atmos. Ocean. Technol. 1990, 7, 829–840. [Google Scholar] [CrossRef]
  20. Zrnic, D.S.; Ryzhkov, A.V. Advantages of rain measurements using specific differential phase. J. Atmos. Ocean. Technol. 1996, 13, 464–476. [Google Scholar] [CrossRef]
  21. Ryzhkov, A.V.; Zrnic, D.S. Precipitation and attenuation measurements at a 10-cm wavelength. J. Appl. Meteorol. 1995, 34, 2121–2134. [Google Scholar] [CrossRef]
  22. Carey, L.D.; Rutledge, S.A.; Ahijevych, D.A.; Keenan, T.D. Correcting propagation effects in C-band polarimetric radar observations of tropical convection using differential propagation phase. J. Appl. Meteorol. 2000, 39, 1405–1433. [Google Scholar] [CrossRef]
  23. He, Y.X.; Lv, D.R.; Xiao, H.; Lei, H.C.; Liu, D.M.; Duan, S. Attenuation correction of reflectivity for X-band dual polarization radar. Chin. J. Atmos. Sci. 2009, 33, 1027–1037. (In Chinese) [Google Scholar]
  24. Hu, Z.Q.; Liu, L.P.; Chu, R.Z. Comparison of different attenuation correction methods and their effects on estimated rainfall using X-band dual linear polarimetric radar. ActaMeteorol. Sin. 2008, 66, 251–261. (In Chinese) [Google Scholar]
  25. Smyth, T.J.; Illingworth, A.J. Correction for attenuation of radar reflectivity using polarization data. Quart. J. R. Meteorol. Soc. 1998, 124, 2393–2415. [Google Scholar] [CrossRef]
  26. Testud, J.; Bouar, E.L.; Obligis, E.; Ali-Mehenni, M. The rain profiling algorithm applied to polarimetric weather radar. J. Atmos. Ocean. Technol. 2000, 17, 332–356. [Google Scholar] [CrossRef]
  27. Park, S.-G.; Bringi, V.N.; Chandrasekar, V.; Maki, M.; Iwanami, K. Correction of radar reflectivity and differential reflectivity for rain attenuation at X-band. Part I: Theoretical and Empirical Basis. J. Atmos. Ocean. Technol. 2005, 22, 1621–1632. [Google Scholar] [CrossRef]
  28. Park, S.-G.; Maki, M.; Iwanami, K.; Bringi, V.N.; Chandrasekar, V. Correction of radar reflectivity and differential reflectivity for rain attenuation at X-band. Part II: Evaluation and Application. J. Atmos. Ocean. Technol. 2005, 22, 1633–1655. [Google Scholar] [CrossRef]
  29. Kim, D.-S.; Maki, M.; Lee, D.-I. Correction of X-band radar reflectivity and differential reflectivity for rain attenuation using differential phase. Atmos. Res. 2008, 90, 1–9. [Google Scholar] [CrossRef]
  30. Delrieu, G.; Caoudal, S.; Creutin, J.D. Feasibility of using mountain return for the correction ofground-based X-band weather radar data. J. Atmos. Ocean Technol. 1997, 14, 368–385. [Google Scholar] [CrossRef]
  31. Hubbert, J.; Bringi, V.N. An iterative filtering technique for the analysis of copolar differential phase and dual-frequency polarimetric variables. J. Atmos. Ocean. Technol. 1995, 12, 643–648. [Google Scholar] [CrossRef]
  32. Chandrasekar, V.; Lim, S.; Gorgucci, E. Simulation of X-band rainfall observations from S-band radar data. J. Atmos. Ocean. Technol. 2006, 23, 1195–1205. [Google Scholar] [CrossRef]
  33. Beard, K.V.; Chuang, C. A new model for the equilibrium shape of raindrops. J. Atmos. Sci. 1987, 44, 1509–1524. [Google Scholar] [CrossRef]
Figure 1. (a) The Kim et al.’s method; (b) the SSCC method, assuming αi = 0 (i = 0, 4; 995,999).
Figure 1. (a) The Kim et al.’s method; (b) the SSCC method, assuming αi = 0 (i = 0, 4; 995,999).
Atmosphere 07 00164 g001
Figure 2. The locations of the X-band radar and the S-band radar,the X-band radar is at Shunyi (SY, 116.68° E, 40.19° N), the S-band radar is at Daxing (DX, 116.47° E, 39.81° N). The symbols BJ, HB and TJ in the Figure are the abbreviation of Beijing City, Hebei Province and Tianjin City, respectively.
Figure 2. The locations of the X-band radar and the S-band radar,the X-band radar is at Shunyi (SY, 116.68° E, 40.19° N), the S-band radar is at Daxing (DX, 116.47° E, 39.81° N). The symbols BJ, HB and TJ in the Figure are the abbreviation of Beijing City, Hebei Province and Tianjin City, respectively.
Atmosphere 07 00164 g002
Figure 3. The reflectivity at X-band (elevation angle: 3°; range: 75 km) at 14:45 BJT, 19 June 2015.
Figure 3. The reflectivity at X-band (elevation angle: 3°; range: 75 km) at 14:45 BJT, 19 June 2015.
Atmosphere 07 00164 g003
Figure 4. Scattering simulation of intrinsic reflectivity at X and S bands. The blue line stands for the intrinsic reflectivity at X-band is equal to the reflectivity at S-band. The red curve is obtained by scattering simulation.
Figure 4. Scattering simulation of intrinsic reflectivity at X and S bands. The blue line stands for the intrinsic reflectivity at X-band is equal to the reflectivity at S-band. The red curve is obtained by scattering simulation.
Atmosphere 07 00164 g004
Figure 5. (a) Original reflectivity at X-band for convective cloud on 19 June 2015; (b) the intrinsic reflectivity at X-band. The color shade guide is the reflectivity factor in unit of dBZ and the X-Y axis is the range in unit of km in the figure and the following figures.
Figure 5. (a) Original reflectivity at X-band for convective cloud on 19 June 2015; (b) the intrinsic reflectivity at X-band. The color shade guide is the reflectivity factor in unit of dBZ and the X-Y axis is the range in unit of km in the figure and the following figures.
Atmosphere 07 00164 g005
Figure 6. Corrected reflectivity at X-band.
Figure 6. Corrected reflectivity at X-band.
Atmosphere 07 00164 g006
Figure 7. Scatter diagrams of the uncorrected (a) and corrected (b) X-band radar reflectivity versus the intrinsic reflectivity. The red line is an ideal line, indicating that the corrected X-band radar reflectivity is equal to the intrinsic reflectivity. The green line is the fitting curve.
Figure 7. Scatter diagrams of the uncorrected (a) and corrected (b) X-band radar reflectivity versus the intrinsic reflectivity. The red line is an ideal line, indicating that the corrected X-band radar reflectivity is equal to the intrinsic reflectivity. The green line is the fitting curve.
Atmosphere 07 00164 g007
Figure 8. (a) Range profile of different reflectivity along the azimuth of 228°, uncorrected reflectivity (UnC), reflectivity corrected by Kim et al.’s method (Kim), reflectivity corrected by the SSCC method (SSCC) and the intrinsic reflectivity (Intrinsic); (b) the measured ΨC and the filtered φDP.
Figure 8. (a) Range profile of different reflectivity along the azimuth of 228°, uncorrected reflectivity (UnC), reflectivity corrected by Kim et al.’s method (Kim), reflectivity corrected by the SSCC method (SSCC) and the intrinsic reflectivity (Intrinsic); (b) the measured ΨC and the filtered φDP.
Atmosphere 07 00164 g008
Figure 9. Cumulative distribution of reflectivity for the convective cloud. In the figure, line Unc is uncorrected reflectivity, line Kim stands for reflectivity corrected by Kim et al.’s method, line SSCC represents reflectivity corrected by the SSCC method, and line Intrinsic stands for intrinsic reflectivity.
Figure 9. Cumulative distribution of reflectivity for the convective cloud. In the figure, line Unc is uncorrected reflectivity, line Kim stands for reflectivity corrected by Kim et al.’s method, line SSCC represents reflectivity corrected by the SSCC method, and line Intrinsic stands for intrinsic reflectivity.
Atmosphere 07 00164 g009
Figure 10. Average biases of the reflectivity for the convective cloud. (a) The AB of the uncorrected reflectivity, the reflectivity corrected by Kim et al.’s method and the reflectivity corrected by the SSCC; (b) the AB of the difference between the Kim and the SSCC.
Figure 10. Average biases of the reflectivity for the convective cloud. (a) The AB of the uncorrected reflectivity, the reflectivity corrected by Kim et al.’s method and the reflectivity corrected by the SSCC; (b) the AB of the difference between the Kim and the SSCC.
Atmosphere 07 00164 g010
Figure 11. Cumulative distribution of reflectivity with different resolutions for the convective cloud.
Figure 11. Cumulative distribution of reflectivity with different resolutions for the convective cloud.
Atmosphere 07 00164 g011
Figure 12. Cumulative distribution of reflectivity with different initial phases for convective cloud.
Figure 12. Cumulative distribution of reflectivity with different initial phases for convective cloud.
Atmosphere 07 00164 g012
Figure 13. Original reflectivity at X-band for the stratiform cloud with embedded convection on 26 June 2015.
Figure 13. Original reflectivity at X-band for the stratiform cloud with embedded convection on 26 June 2015.
Atmosphere 07 00164 g013
Figure 14. (a) Corrected reflectivity at X-band; (b) the intrinsic reflectivity at X-band.
Figure 14. (a) Corrected reflectivity at X-band; (b) the intrinsic reflectivity at X-band.
Atmosphere 07 00164 g014
Figure 15. Original reflectivity at X-band for the stratiform cloud on 16 June 2015.
Figure 15. Original reflectivity at X-band for the stratiform cloud on 16 June 2015.
Atmosphere 07 00164 g015
Figure 16. (a) Corrected reflectivity at X-band; (b) the intrinsic reflectivity at X-band.
Figure 16. (a) Corrected reflectivity at X-band; (b) the intrinsic reflectivity at X-band.
Atmosphere 07 00164 g016
Figure 17. Cumulative distribution of reflectivity for the stratiform rain with embedded convection.
Figure 17. Cumulative distribution of reflectivity for the stratiform rain with embedded convection.
Atmosphere 07 00164 g017
Figure 18. Cumulative distribution of reflectivity for the stratiform rain.
Figure 18. Cumulative distribution of reflectivity for the stratiform rain.
Atmosphere 07 00164 g018
Figure 19. Average bias of reflectivity for the stratiform rain with embedded convection.
Figure 19. Average bias of reflectivity for the stratiform rain with embedded convection.
Atmosphere 07 00164 g019
Figure 20. Average bias of reflectivity for the stratiform rain.
Figure 20. Average bias of reflectivity for the stratiform rain.
Atmosphere 07 00164 g020
Figure 21. Cumulative distribution of reflectivity with different initial phases for the stratiform rain with embedded convection.
Figure 21. Cumulative distribution of reflectivity with different initial phases for the stratiform rain with embedded convection.
Atmosphere 07 00164 g021
Figure 22. Cumulative distribution of reflectivity with different initial phases for the stratiform rain.
Figure 22. Cumulative distribution of reflectivity with different initial phases for the stratiform rain.
Atmosphere 07 00164 g022
Table 1. System Characteristics of the IAP-714XDP-A radar.
Table 1. System Characteristics of the IAP-714XDP-A radar.
ItemIAP-714XDP-A Radar
Frequency9.370 GHz
Antenna type2.4 m diameter parabolic antenna
Antenna gain44.78 dB
Beam width
Pulse width0.5/1/2 μs
Pulse repetition frequency500~2000 Hz
PolarizationHorizontal/Vertical
Observation range75/150/300 km
Observation parametersZH, ZDR, φDP, KDP, ρHV, V, W
Doppler processingPPP/FFT

Share and Cite

MDPI and ACS Style

Feng, L.; Xiao, H.; Wen, G.; Li, Z.; Sun, Y.; Tang, Q.; Liu, Y. Rain Attenuation Correction of Reflectivity for X-Band Dual-Polarization Radar. Atmosphere 2016, 7, 164. https://doi.org/10.3390/atmos7120164

AMA Style

Feng L, Xiao H, Wen G, Li Z, Sun Y, Tang Q, Liu Y. Rain Attenuation Correction of Reflectivity for X-Band Dual-Polarization Radar. Atmosphere. 2016; 7(12):164. https://doi.org/10.3390/atmos7120164

Chicago/Turabian Style

Feng, Liang, Hui Xiao, Guang Wen, Zongfei Li, Yue Sun, Qi Tang, and Yanan Liu. 2016. "Rain Attenuation Correction of Reflectivity for X-Band Dual-Polarization Radar" Atmosphere 7, no. 12: 164. https://doi.org/10.3390/atmos7120164

APA Style

Feng, L., Xiao, H., Wen, G., Li, Z., Sun, Y., Tang, Q., & Liu, Y. (2016). Rain Attenuation Correction of Reflectivity for X-Band Dual-Polarization Radar. Atmosphere, 7(12), 164. https://doi.org/10.3390/atmos7120164

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop