Theoretical and Experimental Study on Echo Fluctuation Suppression of a Cirrus Cloud by Millimeter Wave MIMO Radar

Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Key Laboratory for Aerosol-Cloud-Precipitation of China Meteorological Administration, Nanjing University of Information Science and Technology, Nanjing 210044, China Jiangsu Key Laboratory of Meteorological Observation and Information Processing, Nanjing University of Information Science and Technology, Nanjing 210044, China National Demonstration Center for Experimental Atmospheric Science and Environmental Meteorology Education, Nanjing University of Information Science and Technology, Nanjing 210044, China Key Laboratory of Middle Atmosphere and Global Environment Observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China


Introduction
Cirrus clouds are globally distributed and play an important role in regulating the energy budget in the earth-atmosphere system, but they are still a source of major uncertainties in satellite and climate modeling studies [1][2][3][4][5].
To study macro-and microphysical properties such as cloud height, cloud thickness, and ice water content (IWC), active and passive remote sensing instruments that cover visible, infrared, and millimeter and submillimeter wavelengths have been adopted.These instruments include satellite instruments, such as TRMM (tropical rainfall measuring mission) [6] and CloudSat [7]; airborne-based instruments, such as WCR (U. of Wyoming) [8] and ACR (UMass and NASA) [9]; and ground-based instruments, such as MMCR (ARM millimeter wave cloud radar) [10] and CPRS (cloud profiling radar system) [11].Among them, millimeter wavelength radars have an advantage in detecting the dynamic and structural properties of cirrus clouds [12].To retrieve the physical properties of clouds, the interaction between electromagnetic waves and ice particles should be studied.
Before studying the scattering properties, the size and shape of cirrus clouds have to be confirmed.Cirrus clouds are composed of almost nonspherical ice crystals with various sizes and shapes.Some numerical methods have been widely used to calculate the scattering properties of cirrus cloud particles, such as the T-matrix [13], the finite difference time domain (FDTD) [14], the discrete dipole approximation (DDA) [15], the finite element method (FEM) [16], and the method of moments (MOM) [17].However, these ice particles are smaller compared to the wavelength of a millimeter wave and are usually randomly orientated, so the scattering properties of an ice cloud can be approximately calculated by the Rayleigh approximation [18].
At present, few researchers discuss the characteristics of a scattering electromagnetic wave when ice particles are in motion because modern weather radar has been equipped with video integrators, which can obtain the average value of received power by integral calculation [19].Video integrators can transform the video signal, which can amplify the video signal and conduct analog-to-digital conversion and integration processing.Integration processing consists of calculating the distance integral and azimuth integral.However, video integrators have the disadvantage of reducing the distance resolution and have high average time requirements [20].The radial resolution Δr of radar has the following relationship with the pulse width τ and velocity of light C: Δr = 1/2 * τ * C. If we use a video integrator, then the integrated times of distance K would be introduced, with which the radial resolution Δr of radar has the following relationship along with the pulse width τ and velocity of light C: Δr = 1/2 * K * τ * C. When the pulse width τ remains unchanged, the video integrator would reduce the resolution Δr of the radar range due to K.
These factors led to an important question: how can we suppress echo fluctuation (RCS fluctuation) without reducing the distance resolution?To solve this problem, first, the echo fluctuation phenomenon of the moving spherical particles was simulated by using the vector potential method under the Rayleigh approximation.Second, we studied the reason why video integrators would reduce the resolution of the radar range and easily lose changed information and we suggested that the best way to suppress echo fluctuation is MIMO radar.To validate the suppression effect of MIMO radar, the RCS values of horizontally orientated hexagonal columns, ellipsoids, and spherical particles were simulated theoretically by HFSS software and compared with the results of detecting moving particles with vertically pointing radar.To experimentally validate the echo fluctuation suppression ability of MIMO radar, the RCS values of cuboid candles based on MIMO spatial diversity were measured in a microwave anechoic chamber at Nanjing University of Information Science and Technology.

Analysis of Echo Fluctuation
When the radar transmitter transmits an electromagnetic pulse, an accurate timer is started to measure the time between the starting transmitting pulse and the received echo signal, which is scattered by atmospheric particles.The distance between the radar and the atmospheric particle can be deduced from the consumed time and the speed of light.The atmospheric particles are randomly distributed in space.The magnitude and phase of the electromagnetic waves scattered by atmospheric particles are randomly changed with the enhancement and attenuation of the interference.Atmospheric particles are not fixed at a certain position: assuming the amplitude of the electromagnetic wave scattered to the radar is A, the fluctuating distribution follows the Rayleigh probability distribution, which leads to probability distributions of backscattered power in the forms of A 2 , and then, average processing would be used in order to retrieve the parameters of particles, which are proportional to the backscattering cross section of all particles in the unit volume [21].
To analyze the echo fluctuation of moving ice particles, the scattering properties of moving ice particles must first be computed.Mie theory [22] can be used to compute the scattering properties of ice particles because an ice cloud can be made of approximately spherical particles.Logan has reviewed the early history of light scattering [23], and the book written by Stratton has outlined this theory in terms of spherical vector wave functions.Professor Shanjie Zhang discussed the plane electromagnetic wave scattering of spherical particles using vector potential [24].Rayleigh approximation is a special case of Mie scattering that describes the elastic scattering of light by spherical particles, which is much smaller than the wavelength of incident electromagnetic waves [25].The monographs written by Van de Hulst and Bohren and Huffman have provided a concrete derivation [26,27].The vector potential technique was used, and the phenomenon of echo fluctuation produced by moving particles is studied in this paper.Because the ensembles of ice particles are sparse, the mutual interaction would not be considered.
According to Maxwell equations in a linear, isotropic, homogeneous medium, the charge density and current density are ρ = 0 and J = 0, respectively, and the time-harmonic electromagnetic fields ( E , H ) will satisfy the wave equation; thus, the scattering electric field intensity of the spherical particles in the spherical coordinate system can be obtained by the vector potential technique [24]: where E 0 is the amplitude of an incident electromagnetic wave, k 0 = ω ε 0 μ 0 is the wave number of the incident electromagnetic wave, ε 0 and μ 0 are the dielectric constant and magnetic permeability (usually air) outside the sphere, respectively, r is the distance between the center of sphere and the scattering observation point, a is the radius of the spherical particle, ϕ is the azimuth angle in the spherical coordinate system, and ε r and μ r are the relative permittivity and permeability of the spherical particles, respectively.

International Journal of Antennas and Propagation
Whether they are "E-plane" or "H-plane," the total scattering electromagnetic field can be expressed as follows: The received radar power P r is proportional to the scattering electromagnetic fields E sca total 2 or H sca total 2 , and P r can be obtained by E sca total 2 [19].To compute more conveniently, taking two rotating particles as an example and assuming an incident electromagnetic field is E 0 = 10 kV/m, the frequency of the millimeter wave radar is 94 GHz, and the equivalent dielectric constant of the particle is ε = 3 1307 − j0 0111 [28].The radii of spherical particles A and B are 10 μm and 8 μm, respectively, the distance between A and B is 0.8 m, and the distance between the radar antenna and the center of the ensembles of ice particles is 6000 m, which can be found in Figure 1.
Assuming the two particles are spatially rotated clockwise at angular velocity (Ω t = 1 rad/s) under the wind vector gradient, the starting position can be as shown in Figure 1, the change in distance between the radar antenna and the particles can be computed at 0 to 6.28 seconds, and the received electric field as a function of time can be found as shown in Figure 2 using equation (2).
Figure 2 shows the electric field received at the radar antenna along with changes in time in the situation of a cosine function when the particles move at an angular velocity of 1 rad/s; the changed received electric field makes it impossible to accurately obtain the received power, and the physical parameters of the cloud cannot be accurately retrieved based on this power.The reason cosine function fluctuates is that the distance between the spherical particles and the radar antenna depends on the cosine of the angle θ in Figure 1.To reduce the influence of RCS fluctuation on the retrieval of cloud physical parameters, the weather radar uses the video integrator to suppress the echo fluctuation.However, if the average time is not sufficient, the measured average power still fluctuates greatly; if the average time is too long, the actual change characteristic of the ice particles will be smoothed out and lose their original information.It should be noted that the video integrator can not only suppress echo fluctuation but also reduce the distance resolution [20].

Echo Fluctuation Suppression by Single Millimeter Wave Radar
At present, the single radar uses the video integrator to suppress the echo fluctuation of the moving particles, using average technology to suppress the echo fluctuation as shown in Figure 2. Assume that the radar collects 128 points in 6.28 seconds and define three average numbers: (1) normal average, the average number of times is 128; (2) over average, the average number of times is 130; and (3) under average, the average number of times is 126.The simulation results of echo fluctuation under above three average methods can be found in Figure 3.
Figure 3 shows that the results of echo fluctuation suppression using the normal average are the best.The results of echo fluctuation suppression using the over average are lower, and those with the under average are higher than the actual values (RCS fluctuation), which are results consistent with the conclusions of the video integrator in radar meteorology [19].

International Journal of Antennas and Propagation
In summary, the conventional video integrator of meteorological radar can suppress the echo fluctuation with the situation of a normal average.However, a video integrator will reduce the radar range resolution, and under-and oversuppression phenomena easily appeared when the average number of times was lesser or greater than the normal average, respectively.In addition, the video integrator would reduce the radar range resolution.Therefore, it is of great significance to explore the millimeter wave MIMO radar to suppress the echo fluctuation of moving particles.

Echo Fluctuation Suppression Theory and
Experiment by MIMO 4.1.MIMO Simulation of Echo Fluctuation Suppression.To verify the echo fluctuation suppression ability of MIMO radar, it is assumed that a nonspherical ice particle is a hexagonal column, which is one of the most common particles in cirrus clouds [29].A theoretical simulation of echo fluctuation suppression using a 4 * 4 MIMO antenna array can be found in Figure 4.
In Figure 4, the diameter of the circular array is 4000 m, the distance between the center O of the circular array and the center K of the cloud particles is 6000 m, OC is the positive direction of x, OF is the positive direction of y, and OK is the positive direction of z.The parameters of the transmitting antennas and receiving antennas set in HFSS (high-frequency structure simulator) software can be found in Tables 1  and 2, respectively.
We assume that the targets detected by the MIMO radar array in Figure 4 are ensembles of hexagonal columns, long ellipsoid particles, and spherical particles, which have the same volume.The complex refraction indices of these particles   4 International Journal of Antennas and Propagation are m = 1 782 + j * 0 0027, and the corresponding complex permittivities are ε = 3 1755 + j * 0 0096 [30].The long side and short side of the hexagonal column are 400 μm and 69.6000 μm, respectively; the axis of the long half of the ellipsoid is 200 μm, and the axis of the short half is 77.5186 μm.
The radius of the spherical particle is 106.3197μm, and the volumes of the three particles are 5 0342 * 10 6 μm.The simulated models of ensembles of ice particles in the HFSS environment can be found in Figure 5.
The arrangements of ensembles of hexagonal columns, ellipsoids, and spherical particles are shown in Figures 5(a)-5(c), respectively.These models of particles consist of 147 identical particles.These particles are arranged in the upper and lower left and right spaces with the distance of 1 mm.The definitions of transmitting antenna and receiving antenna of MIMO radar are shown in Tables 1 and 2, respectively, and the incident electromagnetic wave frequency is 94 GHz.The positions of rotated ensembles of particles at 1 s are shown in Figure 5.The positions of rotated ensembles of particles at 2 s, 3 s, 4 s, 5 s, 6 s, 7 s, and 8 s are rotated counterclockwise by 20 °, 40 °, 60 °, 80 °, 100 °, 120 °, and 140 °, with Figure 5 chosen as a benchmark.The simulated RCS values of ensembles of hexagonal columns, ellipsoids, and spherical particles can be found in Figure 6.    3.
Table 3 shows that the maximum RCS fluctuation value occurs at the combination of the transmitting antenna 4 and receiving antenna 1 (24.55 dB) for horizontally orientated hexagonal columns; for horizontally orientated ellipsoids, the maximum RCS fluctuation value occurs at the T stands for transmitting antenna, R stands for receiving antenna, and + stands for combination.Figure 7 shows that the maximum RCS fluctuations are 0.5374 dB, 1.3126 dB, and 0.3856 dB after the coherent cumulative average processing for horizontally rotated hexagonal column particles, ellipsoidal particles, and spherical particles, respectively.The simulation results of ensembles of particles in Figure 5 vertically detected by single radar (the polarization direction of the electric field along the x-axis and the propagation direction of the electromagnetic wave along the z-axis) can be found in Figure 8.
Figure 8 shows that the maximum RCS fluctuation values of the horizontally rotated hexagonal column particles, ellipsoidal particles, and spherical particles are 1.2668 dB, 3.9765 dB, and 0.4342 dB, respectively.Comparing the results of Figure 7 with those of Figure 8, the coherent accumulation averaging results of horizontally rotated hexagonal column particles, ellipsoidal particles, and spherical particles are smaller than those of the vertically pointing single radar.Therefore, MIMO radar can reduce the RCS fluctuation of the moving particles.The following experiment was used to verify the ability of MIMO radar to suppress echo fluctuation.

Experimental Measurements of Echo Fluctuation
Suppression.The echo fluctuation suppression of moving particles by MIMO at 10 GHz using the PNA vector network analyzer (model: E8363C, test frequency: 10 MHz-40 GHz) from Agilent Technologies, USA, can be transformed into results at 94 GHz by electromagnetic similarity theory [31].A set of MIMO echo fluctuation suppression systems based on a vector network analyzer in the microwave anechoic chamber of Nanjing University of Information Science and Technology can be found in Figure 9.
The vector network analyzer in Figure 9 has a two-port device, which cannot achieve multi-input and multioutput characteristics of MIMO radar at the same time; however, the time-based transmitting mode of MIMO spatial diversity can be achieved by the different combinations of transmitting ports and receiving ports at different times.The moving particles are divided into eight different states (8-second motion states), and the ensembles of particles consist of 19 cuboid candles (length and width are 23 mm, height is 49 mm, and the dielectric constant is 2.1); we should notice that the cuboid candles have a different dielectric constant as compared to ice.  Figure 11 shows the positions of the 3 * 3 MIMO antenna array, and the distances between the transmitting antenna and the center of the target support structure, the receiving antenna, and the center of the target support structure are all 0.5 m.The bottom of the target support structure is a rectangular box with a wave-absorbing characteristic (both the length and width are 0.5 m, and the height is 0.3 m).The overall height of the target support structure is 2.5 m.The horn antenna is placed on a 0.15 m thick foam pad, and the elevation angle of all horn antennas is 78 °.The specific positions of the transmitting horns and receiving horns in the MIMO testing experiment can be found in Figure 12.
The location parameters of the transmitting and receiving horns of the 3 * 3 MIMO antenna array in the spherical coordinate system in Figure 12 are shown in Table 4.
The S21 values of the cuboid candles in Figure 10 are detected by the MIMO antenna array in Figure 12 and can be found in Figure 13.
Figure 13 showed the S21 values of moving cuboid candles within 8 s using the vector network analyzer at Nanjing University of Information Science and Technology; these values will be used to compute the RCS values of cuboid candles.
An aluminum plate with dimensions of 50 cm * 50 cm and a thickness of 1.35 mm was used as a calibrator when conducting the MIMO testing experiment.The testing environment can be seen in Figure 14.
The bistatic RCS results of the calibration plate in Figure 14 simulated by HFSS software are shown in Figure 15.
According to the RCS transformation formula [30] and the results of Figures 13 and 15, the RCS fluctuation of the cuboid candles in 8 states can be found in Figure 16.
Figure 16 shows the MIMO detection results of cuboid candles.In order to study the RCS fluctuation values of cuboid candles at different detection methods, we introduced the maximum RCS fluctuation detected by each combination of the MIMO transmitting and receiving antennas, which can be shown in Table 5.
Table 5 shows that the RCS fluctuation is the largest (28.96 dB) for cuboid candles in the combination of transmitting antenna 2 and receiving antenna 2. The results in Figure 16 also showed that RCS fluctuation in the combination of transmitting antenna 2 and receiving antenna 2 is higher than those of other combinations of transmitting antenna and receiving antenna.
The coherent cumulative averaging results of cuboid candles detected by MIMO radar in Figure 16 can be found in Figure 17.
The maximum RCS fluctuation of cuboid candles is 21.24 dB by the coherent cumulative average process, which can be computed in Figure 17.Compared with the maximum RCS fluctuation (transmitting antenna 2, receiving

Conclusions
The echo (RCS) fluctuation phenomenon of the moving spherical particles is calculated by using the vector potential method under the Rayleigh approximation.The RCS fluctuation can be suppressed by the video integrator of single millimeter radar, but the average time is very important to the suppression effect, which can not only easily lose changed information but also reduce the distance resolution.In addition, a video integrator would reduce the resolution of the

Figure 1 :
Figure 1: The geometric relationship between two particles and radar antenna.

Figure 2 :
Figure 2: Variations in received electric field as a function of time when spherical particles move at the angular velocity of 1 rad/s.

Figure 3 :
Figure 3: Simulation results of echo fluctuation under three average methods.

Figure 4 :
Figure 4: A schematic of a cirrus cloud detected by round MIMO array.

Figure 6 :
Figure 6: RCS simulation results of horizontally orientated (a) hexagonal columns, (b) ellipsoids, and (c) spheres by means of MIMO.

Figure 7 :
Figure 7: Coherent cumulative averaging results of (a) horizontally rotated hexagonal column particles, (b) horizontally rotated ellipsoidal particles, and (c) horizontally rotated spherical particles by means of MIMO.

Figure 8 :
Figure 8: RCS fluctuations of three particles detected by vertically pointing single radar.

Figure 11 :
Figure 11: Overhead view of the MIMO antenna positions.

Figure 10 :
Figure 10: Eight different states of cuboid candles corresponding to different times.

Table 4 :Figure 13 :
Figure 13: S21 values of moving cuboid candles within 8 s by means of MIMO array.

9
International Journal of Antennas and Propagation radar range.To solve these problems, the RCS values of horizontally orientated hexagonal columns, ellipsoids, and spherical particles were simulated by HFSS software.Compared with the results of detecting moving particles with vertically pointing radar, RCS fluctuation can be effectively reduced by MIMO radar after conducting coherent accumulative averaging processes for the above results.To experimentally validate the RCS fluctuation suppression ability of MIMO radar, RCS values of cuboid candles based on the time-based transmitting mode of MIMO spatial diversity have been measured in a microwave anechoic chamber at Nanjing University of Information Science and Technology; the results showed that the RCS fluctuation of moving particles can be effectively reduced by MIMO radar.

Figure 14 :Figure 15 :Figure 16 :
Figure 14: The calibrator for testing the cuboid candles by means of the MIMO array.

Figure 17 :
Figure 17: Coherent cumulative averaging values of cuboid candles by means of MIMO array.

Table 2 :
Parameter settings of receiving antennas in HFSS software.

Table 3 :
The difference between the maximum and minimum RCS fluctuation by means of MIMO (units: dB).

Table 5 :
The difference between the maximum and minimum values of RCS fluctuation by means of MIMO (units: dB).stands for transmitting antenna, R stands for receiving antenna, and + stands for combination. T