Echo Simulation of Dual Polarization Doppler Weather Radar Based on the Physical Model

Using the scattering characteristics of particles to simulate the radar echo can supply the test signals close to the real precipitation echo for the weather radar, and save the time and cost of the research and development and maintenance of the weather radar. In this paper, the precipitation echo of weather radar is simulated based on the theoretical basis that the falling raindrops have a shape well approximated by an oblate spheroid in the atmosphere. The Marshal-Palmer distribution is applied to describe the raindrop spectrum distribution of precipitation particles. It is assumed that the raindrop particles of different sizes have the random distribution in the radar resolution volume, and then the spatial distribution of precipitation particles in the resolution volume is modeled. The echo signals of horizontal and vertical polarization channels of dual polarization weather radar are obtained by vector superposition of backscattering echoes of each particle. The experimental results show that this method can describe the microphysical characteristics of precipitation particles more completely and can be used to test the signal processing module of dual polarization Doppler weather radar.


Introduction
The signal receiving and processing part is one of the core parts of the Doppler weather radar system, and its technical indicators play an extremely significant role in the performance evaluation of the radar. In the process of production and research and development of the weather radar, the evaluation of the signal processing performance of radar receiver is an indispensable link. In this process, it is essential to provide an effective signal source for weather radar. If we can simulate the echo sequence received by the real weather radar, not only can it optimize the design parameters of the radar system, but also help to improve the performance and advance the research and development efficiency of the weather radar. Therefore, how to quickly obtain an effective test signal source has become an important research topic.
Capsoni C et al. [1] proposed a multi-parameter radar simulator based on the physical characteristics of rainfall, but did not generate radar echo signals with polarization information.
Capsoni C et al. [2] simulated the multi-parameter polarization radar by setting raindrop particle distribution units with different sizes and shapes. The contribution of each particle to the radar echo is accurately estimated, but the method is lack of real-time. Cheong B L et al. [3] used the Advanced Regional Prediction System (ARPS) to simulate and generate an atmospheric field for avoiding the shortcomings of the Capsoni C research methods, and designed and implemented a real-time radar simulator, which can generate three-dimensional time series. The atmospheric field consists of many small particles, each of which corresponds to a vector echo signal, and the radar echoes are presented by superposition of these vector signals. But this simulation method only provides the spectral moments.
Lupidi A et al. [4] simulated the echoes of X-band polarization weather radar under rain and hail, and calculated the electromagnetic scattering field by utilizing the transmission matrix T-Matrix. Lischi S et al. [5] designed a full polarized Doppler weather radar simulator that can generate the original time series, whose scattering and propagation effects of electromagnetic waves between particles are taken into account. The real weather radar observation data, combined with the power spectrum and propagation characteristics of the radar echo, were employed to produce the radar echo signals. Lischi S et al. [6] combined the Weather Research and Forecasting (WRF) model and the transfer matrix coding to generate radar echo signals, and the 3D scenes of graupel, rain and hail are modeled. The echo signals are gained through simulating the weather radar antenna scanning mode, and the consistent result is obtained by comparing with the real radar test data. However, in this framework, the microphysical model does not output the parameters such as oblateness and descent rate, and there are some differences between the observed results and the simulated moisture variables, especially for ice particles. Tang et al. [7] implemented a weather radar echo simulation system based on power spectrum model which takes 3 the real radar data as the input. The simulated echo map of weather radar is produced according to the generated I and Q signals, but the algorithm limits the range of input velocity. Based on the time autocorrelation between the transmitted pulses of weather radar, Zhi-Chao B U et al. [8] proposed a time-domain signal simulation algorithm of weather radar under the atmospheric stratification mode, and the accuracy of the simulation algorithm was verified by taking real weak precipitation echo data of C-band dual polarization weather radar as comparative data.
Providing weather radar echo signals with polarization information are particularly important for the research and development and testing on dual polarization Doppler weather radar. However, the existing simulation methods require the foundation of real weather radar data and a great amount of computation. This paper, in the light of related research, applies the scattering characteristics of precipitation particles to simulate the echo of dual polarization Doppler weather radar. This method realizes the modeling of meteorological scene based on the theoretical basis that the falling drops are not spherical but have an oblate spheroidal shape at the atmosphere. The electromagnetic scattering matrix of raindrop particles is calculated by T-matrix method, and the radar echo simulation of different polarization waves is achieved. In this paper, this method is applied to the simulation of dual polarization Doppler weather radar echo signals. It is shown that the simulation results can more completely describe the microphysical characteristics of precipitation particles and can be employed in the test of dual polarization Doppler weather radar signal processing module.

The shape of raindrops
The rainfall measurement of dual polarization weather radar is based on the basis that the raindrops in the process of falling are not spherical but have an oblate spheroidal shape, and the shape is determined by the equivalent diameter and axis ratio. The axis ratio b/a , the ratio of the major axis length to minor axis length, is related to the diameter e D of an equivalent-volume spherical raindrop as [9]  where s T is the surface tension of water, g is the acceleration due to gravity, and w  is the water density.

Raindrop spectrum
Raindrop spectrum represents the distribution of the number of raindrops in the unit scale interval and per unit volume. The microphysical structure and evolution properties of precipitation can be derived by analyzing the raindrop spectrum. Marshall-Palmer (M-P) distribution function is widely applied to describe the raindrop spectrum, which is given by [10]

Radar meteorological equation
Radar meteorological equation is a relationship between radar operating range, radar parameters, backscattering cross section and average echo power, which is defined as for horizontal and vertical incident fields are given by [11]   2 56 2 , 4 2 where m is the refractive index of water, and e is the eccentricity of the ellipsoid,      system. The  and  are the elevation of the particle position and the azimuth of the particle position, respectively. The G is the gain of radar during transmitting and receiving, that is, tr G G G  , and the a represents the direction of the electromagnetic scattering field.

(b) The horizontal and vertical reflectivity factor
Some polarization parameters can be gained directly from inversion of raindrop spectrum.
When the rotation axis of raindrop particles is vertically oriented, the horizontal and vertical reflectivity factors can be expressed by , hv  as: (c) The differential reflectivity factor The differential reflectivity factor DR Z is the logarithm of the ratio of horizontal reflectivity factor h Z to vertical reflectivity factor v Z , which reflects the intensity difference between the two polarization directions, and can be obtained by using the horizontal and vertical reflectivity The unit of DR Z is dB , which has nothing to do with the number distribution of precipitation particles in the radar resolution volume, but is related to the shape distribution of precipitation particles. It reflects the spatial orientation of precipitation particles deviated from the sphere and particle swarm.
(d) The specific differential phase Raindrop particles will cause phase shift to electromagnetic waves. Due to the resistance of air, and the long axis of oblate spheroidal raindrops is horizontal, the phase shift of horizontal electric field is relatively larger. The difference between the phase change of horizontal and vertical wave propagation is represented by the specific differential phase ( DP K ), which is not affected by beam blocking or ground clutter and is a little affected by attenuation, so the application of where n V is radial velocity, d f is Doppler shift, and N is the number of pulses. 7

Radar resolution volume
The pulse signal beam emitted by radar can be regarded as conical, and the radar resolution volume is a segment of the volume of the cone, which is drawn in Fig 1.

The matrix of raindrop particle size distribution
The T-matrix can be constructed by raindrop parameters, which are listed in Table 1.

The matrix raindrop position
The position of each raindrop particle in the resolution volume is randomly distributed. The position information is described by the spherical coordinate system and stored in a raindrop position matrix, which can be written as where r represents the radial distance from the radar station,  is the elevation of the particle position, and  is the azimuth of the particle position. The angle is the angle between the particle and the centerline of the beam, which is used to determine which raindrop particles are inside the beam.
where D is the particle diameter in m, and the unit of v is m/s. Using the pulse repetition time, we can determine the moving distance of the raindrop particles, and the raindrop particle position 9 corresponding to each pulse time is saved in the position matrix.

Polarization radar echo simulator
According to the particle scattering model and the detection principle of dual-polarization weather radar, a dual-polarization weather radar echo simulator can be constructed. The input parameters of the simulator mainly include two parts, namely, radar parameters and weather parameters, as shown in Table 2. The output of the simulator is the radar echo signals of vertical polarization and horizontal polarization.
The simulator uses different size resolutions to determine the equivalent diameter of particles in different size ranges according to the size range of raindrops. In the experiment, the axis ratio b/a needs to be computed from the size of the equivalent diameter, but it is difficult to solve Eq.
(1) directly. Hence the least square algorithm is introduced to find the best match of the data by calculating the square of the error minimized. The distribution curves of axis ratio and equivalent diameter are fitted by the least square algorithm, and the resolution of the equivalent diameter of the original data points is 0.01 cm, which can ensure that the corresponding axis ratio can be obtained by setting different precision and range of equivalent diameter (Table 1).
Firstly, the distribution of the original data points is drawn, and then the least square algorithm is used to fit the original discrete data points according to the distribution. The fitting effect is shown in Fig. 2, it can be seen clearly that the fitting result of the least square algorithm is consistent with the theoretical data.

Fig. 2 Fitting effect of the least square method
Thus the effective scattering area of each raindrop particle can be calculated according to Eq. (5).
When the rain rate and particle size are determined, the number concentration () ND corresponding to particles of different sizes can be gained by raindrop spectrum distribution function Eq. (2).
The radar effective scattering cross section can be obtained by the vector superposition of the scattering echo of each raindrop particle, which is given by where k t is the sampling time, and s is the diameter of raindrop particles.
The echo simulator finally outputs the weighted electric field echo, which can be applied to calculate the precipitation parameters such as echo power, reflectivity, differential reflectivity, velocity and so on.

Results and discussions
In this study, the technical indicators of dual polarization Doppler weather radar can be set according to the need. In the precipitation radar echo simulation experiment, it is assumed that the precipitation raindrops are not electrified, and uncharged water drops commonly do not unite on collide but rebound, so there is no uniting and loss of raindrops. And the main technical indicators of the dual polarization Doppler weather radar used are presented in Table 3. In order to prove the accuracy of the simulation signal, the time-frequency characteristics of the signal are analyzed, and the results are as follows. The form of echo voltage is shown in Fig. 3.
It is obvious that the variation rules of horizontal and vertical echo voltage are basically coincident. The meteorological echo is composed of a large number of discrete scattered particle echoes, the position of these scattered particles will change due to the movement of the scatterer, and the amplitude of the echo will fluctuate in that the change of the position. Therefore, the statistical characteristics of echo amplitude can reflect part of the particle variation information. Fig. 4 shows the statistical characteristics of echo amplitude in time-domain. It can be seen that the echo amplitude spectrum has the Rayleigh distribution, which is determined by the microphysical characteristics of precipitation itself. The statistical characteristics of the phase distribution are shown in Fig. 5, from where we can see that the phase has the uniform distribution, which is essentially in agreement with the statistical characteristics of precipitation echoes described in 12 reference [13]. The meteorological echo signal is a random fluctuation signal, and its power spectrum can be represented by power spectral density, as shown in Fig. 6. Fig. 6 (a) shows the normalized power spectral density for echo of horizontally polarized waves, which is calculated by FFT(fast Fourier transforms) algorithm and Fourier transform of autocorrelation algorithm respectively. The corresponding power spectral density of vertically polarized waves is shown in Fig. 6 (b), from which we can conclude that the distribution of power spectral density calculated by these two algorithms is basically the same for both vertical polarization and horizontal polarization. And the distribution regulation is approximate to Gaussian distribution, which conforms to the conclusion of reference [14].  larger than that of vertical polarization, which is identical with that described in reference [15].
And it can be concluded that the reflectivity calculated by the droplet size distribution function is slightly larger than the estimated reflectivity, which proves that the simulation results can accord with the theoretical values.

Fig. 8 Reflectivity estimation
In order to verify the accuracy of the precipitation model, the value of DP K is calculated by Eq. (14) in the light of the distribution of raindrop size, as shown in Fig. 9. It is obvious from Fig. 9 that the value of DP K increases with the increase of rain rate, the reason is that when the rain rate is large, the raindrop distribution function assumed by the model corresponds to more heavy raindrop particles, while DP K is related to raindrop size distribution [16].
15 Fig. 9 Relationship between specific differential phase and rain rate By setting different rain rates, the relationship with the reflectivity can be derived, as shown in Fig. 10. The red dots represent the reflectivity factor calculated by the horizontal echo data, and the green dots indicate the reflectivity factor calculated by the vertical echo data. The least square algorithm is applied to fit the horizontal and vertical reflectivity factor data to acquire the relationship between reflectivity factor and rain rate, and the fitting curve effect as shown in Fig.   10. The theoretical curves h Z , R and v Z , R relationship in the figure come from Crozier et al. [15]. It can be seen that the horizontal reflectivity factor is slightly larger than the vertical reflectivity factor, and its distribution law basically conforms to the theoretical research results of Crozier et al.

Fig. 10
Relationship between rain rate and reflectivity 16

Conclusions
In this paper, the simulation of precipitation echo of dual polarization weather radar is realized by using the theoretical basis that the falling drops are not spherical but have an oblate spheroidal shape at the atmosphere. The Marshall-Palmer distribution is used to describe the raindrop spectrum distribution of precipitation particles, and it is assumed that the raindrop particles have random distribution in the radar resolution volume, thus the spatial distribution of precipitation particles in the radar resolution volume is modeled. The least square algorithm is applied to solve the problem that it is difficult to calculate the axis ratio of precipitation particles.
The backscattering echo of each precipitation particle under horizontal and vertical polarization

Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.