# Regional Attenuation Correction of Weather Radar Using a Distributed Microwave-Links Network

Authors: Yang Xue; Xi-chuan Liu; Tai-chang Gao; Chang-ye Yang; Kun Song
Journal: Advances in Meteorology
 (2017)

---

## Abstract

The complex temporal-spatial variation of raindrop size distribution will affect the precision of precipitation quantitative estimates (QPE) produced from radar data, making it difficult to correct echo attenuation. Given the fact that microwave links can obtain the total path attenuation accurately, we introduce the concept of regional attenuation correction using a multiple-microwave-links network based on the tomographic reconstruction of attenuation coefficients. Derived from the radar-based equation, the effect of rainfall distribution on the propagation of radar and microwave link signals was analyzed. This article focuses on modeling of the tomographic reconstruction of attenuation coefficients and regional attenuation correction algorithms. Finally, a numerical simulation of regional attenuation correction was performed to verify the algorithms employed here. The results demonstrate that the correction coefficient (0.9175) falls between the corrected and initial field of radar reflectivity factor (root mean square error, 2.3476 dBz; average deviation, 0.0113 dBz). Compared with uncorrected data, the accuracy of the corrected radar reflectivity factor was improved by 26.12%, and the corrected rainfall intensity distribution was improved by 51.85% validating the region attenuation correction algorithm. This method can correct the regional attenuation of weather radar echo effectively and efficiently; it can be widely used for the radar attenuation correction and the promotion of quantitative precipitation estimation by weather radar.


---

## Body

## 1. Introduction

The attenuation of precipitation on the propagation of electromagnetic waves is one of the main factors that affect the quantitative precipitation estimation (QPE) by weather radar [1, 2]. The attenuation effect will cause a decrease in radar echo intensity and detection area. In particular, at relatively far distances, radar reflectivity factors are lower than actual values and cannot reflect the actual distribution of rainfall. For short-wavelength radar, such as X-band radar, this type of attenuation is particularly serious. The development of correction techniques that address the attenuation of radar echoes has been extensively analyzed.The most important aspect of attenuation correction is to obtain the characteristic of microwave attenuation. In some weather radar applications, researchers desire to measure attenuation along the propagation path to improve the effect of attenuation correction. Lin and Lv [3] proposed the use of a microwave radiometer to measure the total path attenuation. Serrar et al. [4] suggested the estimation of total attenuation of the echoes via the echo difference between sunny and rainy days. However, variations in terrain and fluctuations in the refraction in the path limit the usefulness of this method. Perez [5] introduced a method that can be used to obtain the attenuation of X-band radar by using dual-wavelength radar (S, X) based on the assumption that the S-band radar is with no attenuation. However, during relatively intense rainfall events, the attenuation difference between two bands is notably small and does not provide additional effective X-band attenuation information. The surface reference technique (SRT) is commonly used with space-borne precipitation radar to calculate the path attenuation by comparing the measured echo differences between areas of rainfall and nonrainfall. This method provides the advantage that the relative error decreases with an increase in integral path, but the reflection and attenuation of the electromagnetic waves are different because of the apparent fluctuation of the ground surface and humidity, which causes the SRT to have some limitations [6].A new approach for rainfall measurement was recently presented using the attenuation caused by rainfall that affects the microwave signals for wireless data exchange [7, 8]. Microwave propagation attenuation can be obtained by detecting the signal level of both the transmitter and receiver of a microwave link. During a rainfall event, the microwave is attenuated by falling raindrops when passing through the rain area. The degree of microwave attenuation can be obtained by measuring the difference between clear and rainy days, and these data can be used to retrieve rainfall intensity and its spatial distribution [9, 10] and correct the radar reflectivity factor as a constraint. Krämer et al. [11, 12] described forward and backward iterative algorithms that can be used to correct the X-band radar reflectivity factor using a 10.5/17.5 GHz dual-frequency microwave link. Krämer and Verworn [13] used a dual-frequency microwave link to correct the C-band radar reflectivity factor. Cummings et al. [14] explored the use of two single-frequency microwave links in weather radar correction. However, these studies only enable corrections of the radar signal along the ray under the link path with a narrow range of a radial beam. Therefore, expanding the correction range using microwave links is worth exploring.To improve the correction of radar attenuation by using microwave links, based on a microwave link for single-beam correction, this paper proposes extending the region of radar attenuation correction using a multiple-microwave links network. In a radar detection area with simultaneously multiple communication links, several microwave links form a network that satisfies a certain topology. Deriving the attenuation coefficient of the grids based on the attenuation information of the microwave links network is the key to obtaining the total attenuation of each radial path of the radar as a reference. Therefore, in view of the above analysis, this paper presents a new method that is similar to that used in medical computerized tomography (CT) imaging technology. The area of radar and microwave link monitoring is addressed using discrete grids.Section2 of this paper analyzes the effects of rainfall distribution on the propagation of radar signals and microwave link signals. Section 3 analyzes and establishes the models of tomographic reconstruction of the attenuation coefficients while Section 4 proposes the regional attenuation correction algorithm, where a numerical simulation experiment spanning 20 × 20 km was performed to verify the algorithm. Section 5 provides a summary and concluding remarks.


## 2. Theory of Radar Echo Attenuation Correction

When electromagnetic waves propagate in the atmosphere, the physical effects of weather target scattering, absorption, and reflection cause the energy of the propagation path to become attenuated. Suppose thatP¯r0 is the average echo power without considering the attenuation of the meteorological target and P¯r is the average echo power after the attenuation is considered. The attenuation characteristic can be described by the attenuation factor Kr; then(1)P¯r=P¯r0×Kr.SupposedP¯r is the attenuation value of the received power, which is affected by rainfall and other factors between radar and the target; this attenuation value can be expressed as(2)dP¯r=-2KtP¯r0dr.Next, the integral from 0 toR is calculated, where R is the distance between the radar and the target, and the average echo power is(3)P¯r=P¯r0×e-2∫0RKtdr.Therefore, the dimension ofKt is 1/km. Because the attenuation of the received power is usually expressed in decibels (dB), the attenuation coefficient is converted into kt in dB/km. According to lgM=0.4343lnM, (3) can be written as(4)lnP¯rP¯r0=-2∫0RKtdr.Further derivation shows that(5)P¯r=P¯r0×10-0.2∫0Rktdr.Next, compare (1) with (5)(6)Kr=10-0.2∫0Rktdr.If we use the radar reflectivity observationZm (mm6/m3) and the real radar reflectivity Zr (mm6/m3) of the target instead of the echo power in the previous derivation process, then(7)Zmr=Zrr×Kr=Zrr×10-0.2∫0Rktdr,where R is the distance between the radar and the detection target, Kr is the attenuation factor, and kt is the attenuation coefficient in dB/km. The logarithm of (7) is expressed as(8)lgZm=lgZr-0.2∫0Rktdr.In the process, we use the radar reflectivity factor. The conversion formula between radar reflectivityZ (mm6/m3) and radar reflectivity factorz (dBz) is(9)z=10logZ.Further derivation shows that(10)zr=zm+2∫0Rktdr,where kt is the attenuation coefficient in dB/km; 2∫0Rktdr is the total attenuation value of the radar path after integration in units of dB; zr is the corrected radar reflectivity factor, and zm is the radar-measured reflectivity factor. The accurate determination of the radar path integral attenuation is the key to attenuation correction. Therefore, using the attenuation of microwave link along the path attrain, the attenuation correction process can be expressed as(11)zr=zm+2×attrain.


## 3. Tomographic Model of the Attenuation Coefficient

### 3.1. Reconstruction of the Grid Attenuation Using Microwave Links

CT (computerized tomography) imaging technology is mainly based on ray scanning, where the attenuation of a ray power pass is obtained through the wave field. Tomographic techniques include techniques aimed toward reconstructing the cross-sectional distribution of a parameter from a set of one-dimensional transmission or reflection data, measured and collected along many different paths crossing the spatial domain where the object field has to be reconstructed. The theoretical basis of a CT image reconstruction algorithm is based on a Radon transform and Fourier slicing theorem [15]. This is characterized as follows: (1) the projection data contain the characteristic information of the original image, which can be reconstructed with the information; (2) in order to realize the reconstruction of the image, theoretically, an infinite number of continuous projection data is needed.Radon transform is a linear integral projection transformation, and assumes the distribution function of a two-dimensional target isf(x,y); then Radon transformp for the function f(x,y) is integral along the straight linez:(12)p=∫-∞∞fx,ydz=∫-∞∞f^r,θdz=∫-∞∞fl2+z2,ϕ+arctanzldz.Therefore, the problem of tomographic image reconstruction involves the calculation of the image functionf(x,y) by the projection data p. The meanings of the other quantities in (12) are shown in Figure 1. Meanwhile, the Radon inverse transform satisfies(13)f^r,θdz=12π2∫0π∫-∞∞1rcosθ-ϕ-l∂p∂ldldϕ.Figure 1
Coordinate system used in the Radon formula, wherer,z,l refer to the length, and θ,ϕ refer to the angle.The corresponding tomography is the inverse of the Radon transform.f(x,y) is a continuous function, and the actual projection data are finite. Therefore, f(x,y) as a solution is not unique. In practical work, a limited number of projection data is used to get a satisfactory reconstruction effect. Therefore, in the tomography, the continuous image is usually discretized to the pixel level initially, and the average value of f(x,y) in each pixel is obtained.With reference to the principle of CT imaging technology, the monitoring area is discretized. For the general measurement of the two-dimensional distribution, the nonuniform characteristic measurement area can be divided intom×n grids, with several imaging units (Figure 2). Assuming that the distribution within each grid is uniform, each grid cell can be referred to as a pixel. The S site is a ray transmitter designed to emit with a fixed power, the R site is an energy receiver designed to receive the ray energy after attenuation. The energy attenuation of the ray that passes through each grid is measured via the difference between the transmitted and received power. The pixel distribution f(x,y) in the measurement area was reconstructed with the inversion technique.Figure 2
Reconstruction of the regional discrete diagram. This region was divided intoN equal square meshes Cj (j=1,…,N). A microwave link Li is taken as an example, where the transmitting and receiving terminals are S and R, respectively.The energy attenuation coefficient of each pixel isK(f(x,y)), where the pixel-energy attenuation coefficient conversion operator is denoted as K(·). Assuming that the propagation path of theith ray is Li and the total energy attenuation is Ai,(14)Ai=∫LiKfx,ydl.The key to discrete image reconstruction is related to the problem of calculating the area attenuation coefficient vectorK by measuring the total attenuation A of the area through a series of rays.


### 3.2. Establishment of the Attenuation Coefficient Model

In the weather radar monitoring area, a network of microwave links pass through the rain area. The transmitting end of the microwave link realizes the transmission of a microwave signal with the fixed power, frequency, and polarization, and the receiving end receives the power after attenuation. Figure3 provides a schematic representation of the attenuation coefficient tomography for a hypothetical region. This region was divided intoN equal square meshes Cj (j=1,…,N). A microwave link Li is taken as an example, where the transmitting and receiving terminals are S and R, respectively. The link passes through the grid at C5, C6, C3, and C4, and their lengths in the grid are denoted as li5, li6, li3, and li4, respectively. In practical applications, the monitoring area was divided into the appropriate size and shape based on a combination of local geographical conditions, link density, and other information.Figure 3
Schematic diagram of the chromatography. The link passes through the grid atC5, C6, C3, and C4, and their lengths in the grid are denoted as li5, li6, li3, and li4, respectively.If the number of grids after discretization isN, the rain attenuation of theith microwave link can be expressed as(15)Ai=∑j=1Nlijkj,where lij is the length of theith microwave link that passes through thejth grid; the unit of kj is (dB/km), and i=1⋯M (for M wireless microwave links).According to the power law relation of rainfall attenuation, the rain attenuation of the link can be related to the link length, the average rain rate of the link, and the power law parameter of rain attenuation(16)Ai=aiRibidi,where Ai is the total attenuation of the propagation of theith microwave link; Ri is the average observed rainfall intensity on theith link;di is the length of theith link; ai and bi are the conversion constants of theith link.Suppose that there areM microwave links in the monitoring area. By combining (15), the rain attenuation of the M microwave links can be expressed in vector form:(17)LK=A,where the total attenuation column vector is A=AiM×1; the distance matrix is L=lijM×N; the attenuation coefficient column vector is K=(kj)N×1; kj is the attenuation coefficient of the grid Cj (assuming that the attenuation coefficient in the same grid cell is identical). lij is the length of theith link in the grid Cj. Then, the chromatographic model can be further established as follows:(18)l11k1+l12k2+⋯+l1jkj+⋯+l1NkN=A1,l21k1+l22k2+⋯+l2jkj+⋯+l2NkN=A2,⋮li1k1+li2k2+⋯+lijkj+⋯+liNkN=Ai,⋮lM1k1+lM2k2+⋯+lijkj+⋯+lMNkN=AM.


### 3.3. Solving the Algorithm of the Attenuation Coefficient Model

The solution of the tomography model for attenuation coefficients is found by solving the linear system (18). Because the microwave links only pass through a small number ofN grids, their distribution is sparse with respect to the monitoring area. The number of unknown parameters is much larger than the number of equations. Therefore, it is an ill-posed linear inverse problem if one attempts to use the attenuation information of the microwave link to perform the attenuation coefficient tomographic inversion.According to the characteristics of the linear equation, the attenuation coefficient field is iteratively solved. The related basic theory was recently explained and described in detail by Giuli et al. [16, 17] who provide more detailed information. To minimize the error of (17), the Simultaneous Iterative Reconstruction Technique was used to construct the cost function using the least-square criterion. The objective function was minimized by repeated iterations, and the optimal solution was obtained:(19)JK=A-LKTA-LK=min!.To obtain the minimum cost function,(20)LTA=LTLK,the iterative formula needed to solve (18) can be written as(21)kjt=kjt-1+1M∑i=1NAit-∑p=1Mliptkpt-1·lijt∑p=1Mlipt2.Here is the process of the specific iterative algorithm. First, lett=0 and give the initial value of the vector. Then, let t=t+1 and estimate k(t) by (21). Finally, stop the iteration when k(t)-k(t-1)<ε or t>T, where ε, T is a predefined threshold. The process is repeated until the optimal solution is iterated.


## 3.1. Reconstruction of the Grid Attenuation Using Microwave Links

CT (computerized tomography) imaging technology is mainly based on ray scanning, where the attenuation of a ray power pass is obtained through the wave field. Tomographic techniques include techniques aimed toward reconstructing the cross-sectional distribution of a parameter from a set of one-dimensional transmission or reflection data, measured and collected along many different paths crossing the spatial domain where the object field has to be reconstructed. The theoretical basis of a CT image reconstruction algorithm is based on a Radon transform and Fourier slicing theorem [15]. This is characterized as follows: (1) the projection data contain the characteristic information of the original image, which can be reconstructed with the information; (2) in order to realize the reconstruction of the image, theoretically, an infinite number of continuous projection data is needed.Radon transform is a linear integral projection transformation, and assumes the distribution function of a two-dimensional target isf(x,y); then Radon transformp for the function f(x,y) is integral along the straight linez:(12)p=∫-∞∞fx,ydz=∫-∞∞f^r,θdz=∫-∞∞fl2+z2,ϕ+arctanzldz.Therefore, the problem of tomographic image reconstruction involves the calculation of the image functionf(x,y) by the projection data p. The meanings of the other quantities in (12) are shown in Figure 1. Meanwhile, the Radon inverse transform satisfies(13)f^r,θdz=12π2∫0π∫-∞∞1rcosθ-ϕ-l∂p∂ldldϕ.Figure 1
Coordinate system used in the Radon formula, wherer,z,l refer to the length, and θ,ϕ refer to the angle.The corresponding tomography is the inverse of the Radon transform.f(x,y) is a continuous function, and the actual projection data are finite. Therefore, f(x,y) as a solution is not unique. In practical work, a limited number of projection data is used to get a satisfactory reconstruction effect. Therefore, in the tomography, the continuous image is usually discretized to the pixel level initially, and the average value of f(x,y) in each pixel is obtained.With reference to the principle of CT imaging technology, the monitoring area is discretized. For the general measurement of the two-dimensional distribution, the nonuniform characteristic measurement area can be divided intom×n grids, with several imaging units (Figure 2). Assuming that the distribution within each grid is uniform, each grid cell can be referred to as a pixel. The S site is a ray transmitter designed to emit with a fixed power, the R site is an energy receiver designed to receive the ray energy after attenuation. The energy attenuation of the ray that passes through each grid is measured via the difference between the transmitted and received power. The pixel distribution f(x,y) in the measurement area was reconstructed with the inversion technique.Figure 2
Reconstruction of the regional discrete diagram. This region was divided intoN equal square meshes Cj (j=1,…,N). A microwave link Li is taken as an example, where the transmitting and receiving terminals are S and R, respectively.The energy attenuation coefficient of each pixel isK(f(x,y)), where the pixel-energy attenuation coefficient conversion operator is denoted as K(·). Assuming that the propagation path of theith ray is Li and the total energy attenuation is Ai,(14)Ai=∫LiKfx,ydl.The key to discrete image reconstruction is related to the problem of calculating the area attenuation coefficient vectorK by measuring the total attenuation A of the area through a series of rays.


## 3.2. Establishment of the Attenuation Coefficient Model

In the weather radar monitoring area, a network of microwave links pass through the rain area. The transmitting end of the microwave link realizes the transmission of a microwave signal with the fixed power, frequency, and polarization, and the receiving end receives the power after attenuation. Figure3 provides a schematic representation of the attenuation coefficient tomography for a hypothetical region. This region was divided intoN equal square meshes Cj (j=1,…,N). A microwave link Li is taken as an example, where the transmitting and receiving terminals are S and R, respectively. The link passes through the grid at C5, C6, C3, and C4, and their lengths in the grid are denoted as li5, li6, li3, and li4, respectively. In practical applications, the monitoring area was divided into the appropriate size and shape based on a combination of local geographical conditions, link density, and other information.Figure 3
Schematic diagram of the chromatography. The link passes through the grid atC5, C6, C3, and C4, and their lengths in the grid are denoted as li5, li6, li3, and li4, respectively.If the number of grids after discretization isN, the rain attenuation of theith microwave link can be expressed as(15)Ai=∑j=1Nlijkj,where lij is the length of theith microwave link that passes through thejth grid; the unit of kj is (dB/km), and i=1⋯M (for M wireless microwave links).According to the power law relation of rainfall attenuation, the rain attenuation of the link can be related to the link length, the average rain rate of the link, and the power law parameter of rain attenuation(16)Ai=aiRibidi,where Ai is the total attenuation of the propagation of theith microwave link; Ri is the average observed rainfall intensity on theith link;di is the length of theith link; ai and bi are the conversion constants of theith link.Suppose that there areM microwave links in the monitoring area. By combining (15), the rain attenuation of the M microwave links can be expressed in vector form:(17)LK=A,where the total attenuation column vector is A=AiM×1; the distance matrix is L=lijM×N; the attenuation coefficient column vector is K=(kj)N×1; kj is the attenuation coefficient of the grid Cj (assuming that the attenuation coefficient in the same grid cell is identical). lij is the length of theith link in the grid Cj. Then, the chromatographic model can be further established as follows:(18)l11k1+l12k2+⋯+l1jkj+⋯+l1NkN=A1,l21k1+l22k2+⋯+l2jkj+⋯+l2NkN=A2,⋮li1k1+li2k2+⋯+lijkj+⋯+liNkN=Ai,⋮lM1k1+lM2k2+⋯+lijkj+⋯+lMNkN=AM.


## 3.3. Solving the Algorithm of the Attenuation Coefficient Model

The solution of the tomography model for attenuation coefficients is found by solving the linear system (18). Because the microwave links only pass through a small number ofN grids, their distribution is sparse with respect to the monitoring area. The number of unknown parameters is much larger than the number of equations. Therefore, it is an ill-posed linear inverse problem if one attempts to use the attenuation information of the microwave link to perform the attenuation coefficient tomographic inversion.According to the characteristics of the linear equation, the attenuation coefficient field is iteratively solved. The related basic theory was recently explained and described in detail by Giuli et al. [16, 17] who provide more detailed information. To minimize the error of (17), the Simultaneous Iterative Reconstruction Technique was used to construct the cost function using the least-square criterion. The objective function was minimized by repeated iterations, and the optimal solution was obtained:(19)JK=A-LKTA-LK=min!.To obtain the minimum cost function,(20)LTA=LTLK,the iterative formula needed to solve (18) can be written as(21)kjt=kjt-1+1M∑i=1NAit-∑p=1Mliptkpt-1·lijt∑p=1Mlipt2.Here is the process of the specific iterative algorithm. First, lett=0 and give the initial value of the vector. Then, let t=t+1 and estimate k(t) by (21). Finally, stop the iteration when k(t)-k(t-1)<ε or t>T, where ε, T is a predefined threshold. The process is repeated until the optimal solution is iterated.


## 4. Region Attenuation Correction Model and Numerical Simulation

### 4.1. Region Attenuation Correction Model

Based on the establishment of the tomographic model for attenuation coefficients, the processing of the regional attenuation correction model (with the microwave link uniform distribution as an example) is described as follows:(1)
Construct a joint observation network of microwave links and weather radar data.(2)
The observational radar data and microwave attenuation data are matched with space and time, and the area is discretized.(3)
Using the rain attenuation of the microwave link network, a joint iterative reconstruction algorithm is used to obtain the rain attenuation coefficientkj of each grid cell.(4)
Assuming that there areX lines between each grid point and the radar, we can substitute the attenuation coefficients obtained from (18) into (22). Then, we can get the path attenuation of the qth link between the radar and each grid point:(22)Attqj=∑q=1X∑j=1Nkj×dqj,where dqj is the length of theqth link between the radar and the grid point through thejth grid.(5)
The corrected radar reflectivity factorzqj can be calculated by summing the radar-measured reflectivity factor zmqj and total attenuation of the two-way path on theqth link between the radar and the grid point:(23)zqj=zmqj+2×Attqj.


### 4.2. Establishment of the Initial Simulation Field

First, the 2D simulation region that corresponded to the microwave link and the radar data was established. Figure4(a) shows the monitoring area of the microwave link network, which has a uniform distribution. The monitoring area was discretized, and the total distribution was 20 × 20 km, which was divided into many 2 × 2 km square grids (Figure 4(b)). There were 34 microwave links in the simulation area. Among them, the microwave transmitting terminal is labeled as a point, the microwave receiving terminal is marked as pentagonal, the distance between adjacent transmitting and receiving terminals is 2 km, and the density of the microwave link network is 0.17/km2. The microwave end transmitted the signal with a frequency of 9.47 GHz and vertical polarization.Figure 4
(a) The monitoring area of the microwave link network; (b) discretized monitor area. There were 34 microwave links in the simulation area. The monitoring area was divided into many 2 × 2 km square grids.
(a)
(b)The initial field of rainfall had two weak-strong rainfall centers (Figure5). Using the radar empirical relationship Z=300R1.4 for simulation and calculation, the radar reflectivity factor initial field that corresponds to each grid was obtained (Figure 6).Figure 5
Initial field of rainfall intensity in the simulation area.Figure 6
Initial field of radar reflectivity factor in the simulated area.The attenuation coefficient that corresponds to each grid was obtained through the 9.47 GHz rain attenuation relationshipko=0.0085R1.2645 in the International Telecommunication Union (ITU) rainfall attenuation model [18]. The attenuation of each grid was multiplied by the distance between the grid and the radar. The radar reflectivity factor observation field was gained by subtracting the initial field of the reflectivity factor from the attenuation that corresponds to each grid. Figure 7 shows the simulated radar reflectivity factor observation field.Figure 7
Observation field of radar reflectivity factor.


### 4.3. Numerical Simulation Results and Analysis

To evaluate the correction algorithm, we used the correlation coefficient, mean square error, and average deviation to describe the correction effect:(24)ρ=∑j=1Nzrj-zr¯zj-z¯∑j=1Nzrj-zr¯2∑j=1Nzj-z¯2,RMSE=1N∑jNzrj-zj2,e=∑jNzj-zrj∑jNzrj,where zr and z are the initial and corrected radar reflectivity factor fields, respectively; zr¯ and z¯ are the average initial and corrected radar reflectivity factors, respectively.Figures8(a) and 9(a) show the uncorrected radar reflectivity factor and uncorrected rainfall intensity, while Figures 8(b) and 9(b) show the corrected radar reflectivity factor and corrected rainfall intensity, respectively. Figures 8 and 9 show that the attenuation of the microwave link network significantly corrects the radar reflectivity factor and rainfall intensity. Compared with the uncorrected radar reflectivity factor, the correlation coefficient between the corrected reflectivity factor and the initial field of reflectivity factor improved from 0.4842 to 0.9175. The mean square error decreased from 14.5655 to 2.3476 dBz, and the mean deviation decreased from 0.3011 to 0.0113 dBz. Compared with the uncorrected rainfall intensity, the correlation coefficient between the corrected rainfall intensity and the initial field of rainfall increased from 0.7315 to 0.9271, the mean square error declined from 21.3271 to 8.2453 mm/h, and the mean deviation decreased from 0.7564 to 0.0359 mm/h. The corrected rainfall field can well reflect the position of the peak center, and the correction of the peak center of the heavy rainfall was slightly larger. The results show that the tomographic reconstruction model improved upon the radar echo attenuation correction quite well.Figure 8
Radar reflectivity factor: (a) uncorrected; (b) corrected.
(a)
(b)Figure 9
Rainfall intensity: (a) uncorrected; (b) corrected.
(a)
(b)


## 4.1. Region Attenuation Correction Model

Based on the establishment of the tomographic model for attenuation coefficients, the processing of the regional attenuation correction model (with the microwave link uniform distribution as an example) is described as follows:(1)
Construct a joint observation network of microwave links and weather radar data.(2)
The observational radar data and microwave attenuation data are matched with space and time, and the area is discretized.(3)
Using the rain attenuation of the microwave link network, a joint iterative reconstruction algorithm is used to obtain the rain attenuation coefficientkj of each grid cell.(4)
Assuming that there areX lines between each grid point and the radar, we can substitute the attenuation coefficients obtained from (18) into (22). Then, we can get the path attenuation of the qth link between the radar and each grid point:(22)Attqj=∑q=1X∑j=1Nkj×dqj,where dqj is the length of theqth link between the radar and the grid point through thejth grid.(5)
The corrected radar reflectivity factorzqj can be calculated by summing the radar-measured reflectivity factor zmqj and total attenuation of the two-way path on theqth link between the radar and the grid point:(23)zqj=zmqj+2×Attqj.


## 4.2. Establishment of the Initial Simulation Field

First, the 2D simulation region that corresponded to the microwave link and the radar data was established. Figure4(a) shows the monitoring area of the microwave link network, which has a uniform distribution. The monitoring area was discretized, and the total distribution was 20 × 20 km, which was divided into many 2 × 2 km square grids (Figure 4(b)). There were 34 microwave links in the simulation area. Among them, the microwave transmitting terminal is labeled as a point, the microwave receiving terminal is marked as pentagonal, the distance between adjacent transmitting and receiving terminals is 2 km, and the density of the microwave link network is 0.17/km2. The microwave end transmitted the signal with a frequency of 9.47 GHz and vertical polarization.Figure 4
(a) The monitoring area of the microwave link network; (b) discretized monitor area. There were 34 microwave links in the simulation area. The monitoring area was divided into many 2 × 2 km square grids.
(a)
(b)The initial field of rainfall had two weak-strong rainfall centers (Figure5). Using the radar empirical relationship Z=300R1.4 for simulation and calculation, the radar reflectivity factor initial field that corresponds to each grid was obtained (Figure 6).Figure 5
Initial field of rainfall intensity in the simulation area.Figure 6
Initial field of radar reflectivity factor in the simulated area.The attenuation coefficient that corresponds to each grid was obtained through the 9.47 GHz rain attenuation relationshipko=0.0085R1.2645 in the International Telecommunication Union (ITU) rainfall attenuation model [18]. The attenuation of each grid was multiplied by the distance between the grid and the radar. The radar reflectivity factor observation field was gained by subtracting the initial field of the reflectivity factor from the attenuation that corresponds to each grid. Figure 7 shows the simulated radar reflectivity factor observation field.Figure 7
Observation field of radar reflectivity factor.


## 4.3. Numerical Simulation Results and Analysis

To evaluate the correction algorithm, we used the correlation coefficient, mean square error, and average deviation to describe the correction effect:(24)ρ=∑j=1Nzrj-zr¯zj-z¯∑j=1Nzrj-zr¯2∑j=1Nzj-z¯2,RMSE=1N∑jNzrj-zj2,e=∑jNzj-zrj∑jNzrj,where zr and z are the initial and corrected radar reflectivity factor fields, respectively; zr¯ and z¯ are the average initial and corrected radar reflectivity factors, respectively.Figures8(a) and 9(a) show the uncorrected radar reflectivity factor and uncorrected rainfall intensity, while Figures 8(b) and 9(b) show the corrected radar reflectivity factor and corrected rainfall intensity, respectively. Figures 8 and 9 show that the attenuation of the microwave link network significantly corrects the radar reflectivity factor and rainfall intensity. Compared with the uncorrected radar reflectivity factor, the correlation coefficient between the corrected reflectivity factor and the initial field of reflectivity factor improved from 0.4842 to 0.9175. The mean square error decreased from 14.5655 to 2.3476 dBz, and the mean deviation decreased from 0.3011 to 0.0113 dBz. Compared with the uncorrected rainfall intensity, the correlation coefficient between the corrected rainfall intensity and the initial field of rainfall increased from 0.7315 to 0.9271, the mean square error declined from 21.3271 to 8.2453 mm/h, and the mean deviation decreased from 0.7564 to 0.0359 mm/h. The corrected rainfall field can well reflect the position of the peak center, and the correction of the peak center of the heavy rainfall was slightly larger. The results show that the tomographic reconstruction model improved upon the radar echo attenuation correction quite well.Figure 8
Radar reflectivity factor: (a) uncorrected; (b) corrected.
(a)
(b)Figure 9
Rainfall intensity: (a) uncorrected; (b) corrected.
(a)
(b)


## 5. Conclusions

A regional attenuation correction method for weather radar is proposed based on using a multiple-microwave-links network using a tomographic reconstruction for attenuation coefficients. Derived from the radar equation, the effect of rainfall distribution on the propagation of the radar signal and microwave link signal were analyzed. After the attenuation of microwave links was matched with the weather radar data in the area, which was discretized, the model of the tomographic reconstruction was established for the attenuation coefficient of each grid cell, based on the multiple-microwave-links network and using Simultaneous Iterative Reconstruction Technique. The established model of regional attenuation correction achieved the attenuation correction for the radar reflectivity factor, where the total attenuation value of the path was obtained by adding the distance between the radar sites and each grid.A numerical simulation of regional attenuation correction (20 × 20 km) was performed to verify the algorithm. There were 34 microwave links in the simulation area, and the microwave transmitting end performed the signal transmission at the frequency of 9.47 GHz with vertical polarization. The results show that, compared with the uncorrected radar reflectivity and the uncorrected rainfall intensity, the accuracy of corrected radar reflectivity factor improved by 26.12%, and the corrected rainfall intensity distribution improved by 51.85%. Therefore, the region attenuation correction algorithm was validated.The advantage of using regional attenuation correction for radar with a multiple-microwave-links network is that this technique enlarges the scope of correction when compared with the method of single-path correction. The method has important application value for radar attenuation correction effect and improves the precision of quantitative precipitation estimation by weather radar.At present, commercial microwave links in cellular communication networks, electric power communication, maritime telecommunication, and microwave relay communication system have widely covered the urban and rural areas; on the basis of above equipment, there is no need to construct additional hardware; there is only a need to establish a network for signal acquisition from the existing microwave link devices; the economic cost is relatively low; therefore it has a unique advantage in the financial and engineering aspects. On the other hand, the regional attenuation correction method has a wider coverage and a higher spatial resolution than the single-point correction method, and the spatial representativeness is also much better. This paper has verified the feasibility and validity of above method. Aiming at the irregular distribution structures of the actual microwave links, we will investigate the attenuation coefficients reconstruction and regional attenuation correction under the condition of irregular topology network in the future.



---

*Source: 8621239-2017-05-22.xml*