remote

: It is commonly believed that clear-air echoes detected by weather radars are caused by atmobios migration. However, clear-air echoes are sometimes inconsistently related to the activity of living creatures. In some cases, the characteristics of radar products seem to conform to biological scattering, but the movement of echoes cannot be observed. For these reasons, we sought to expand the cause of clear-air echoes from a Chinese Doppler S-band Weather Radar (CINRAD/SA) in Beijing. Some contradictions were discovered in a case which diverged from previous conclusions. It was found that the progression and movement of clear-air echoes do not conform to the rules of biological activities. The frequency distribution of dual-wavelength ratio peaks is 21.5 dB, which is in accordance with Villars–Weisskopf’s turbulence theory. From 1 May to 20 May, the 58% dual-wavelength ratio between the S-band and the X-band was distributed between 18 dB and 24 dB. These results show that more than half of the clear-air echoes of CINRAD/SA at night were caused by turbulence in Beijing. A new model of troposcatter propagation, the reﬂecting-layers model, was then introduced to explain the radar observations. According to the reﬂecting-layers model, the echoes’ diurnal variation and reﬂectivity characteristics are attributed to the effects of turbulent mixing. Excessive turbulent mixing affects the generation of the reﬂective layer, thereby weakening the echo signal. It is necessary to re-examine the position of turbulence in clear-air echoes.


Introduction
Billions of atmobios cross the sky each year in search of food, partners, and habitats.With the imminent threat to the ecological state driven by human activity, broad-scale biological monitoring may prove crucial to successful conservation efforts.However, efforts to monitor atmobios are hampered by the unpredictability of their movements [1].Traditional methods, such as visual and auditory observations, laboratory research, trapping, and ringing studies, have left a wide gap caused by limitations of space, time, and labor [2].Since weather radars, which were originally used to observe clouds and rain, have been found to be able to observe the bio-scattering of atmobios after World War II, many researchers deemed that atmobios are responsible for clear-air echoes.For instance, Wilson et al. compared the reflectivity at different wavelengths to conclude that particulate scattering dominates in the boundary layer [3].Martin et al. used data from the Weather Surveillance Radar-1988 Doppler (WSR88D) and X-and W-band research radars and deemed that the targets of nocturnal clear-air echoes are almost insects [4].Further, Broeke found that biological scatterers, consisting of birds and insects, may become trapped near the circulation center of tropical cyclones [5].Westbrook et al. used a WSR-88D radar to detect corn earworm moth migration [6].Now, radars are believed to be a practical tool and an important data source for monitoring atmobios activity.
Although two mechanisms, turbulence-scattering mechanisms caused by turbulent inhomogeneities and biological scatterers exemplified by insects and birds, can dominate the scattering process of clear-air echoes [7], it is recognized that most of the echoes at centimeter wavelengths are primarily caused by insects and birds.This viewpoint is based on the theory of locally homogeneous isotropic turbulence developed by Kolmogorov [8,9].
The Kolmogorov-Obukhov theory holds that turbulent motion is homogeneous and isotropic, and its average properties are uniquely determined by the average rate of dissipation of the turbulent kinetic energy per unit mass of fluid within a subrange or regime of turbulent eddy sizes [9].According to this theory, Ottersten clarified the relationship of turbulence scattering (frequently referred to as Bragg scattering) from refractive index irregularities to the atmospheric structure [8].In light of this relationship, Wilson compared the reflectivity at different wavelengths and examined the differential reflectivity at the S-band, concluding that biological scatterers dominate clear-air echoes because a smaller reflectivity difference and a nonzero value of the differential reflectivity are not consistent with Bragg scattering [3].
However, Landau's query (1957) caused Kolmogorov and Obukhov (1961) to introduce important modifications to Kolmogorov's theory.They took into account spatial fluctuations in the turbulent energy dissipation and chose a specific form (lognormal) for the probability density as a third hypothesis [10].However, Mandelbrot (1976) pointed out that the lognormal assumption is only a special, probably physically unrealistic, case of weighted curdling [11].Moreover, many experimental cases for the fine-structure intermittency of turbulence showed the inhomogeneity of physical quantities in space and time [12][13][14][15][16][17][18][19][20][21].Batchelor and Towsend firstly observed that the turbulence and its energy dissipation were very unevenly distributed over space, and the velocity gradients became increasingly intermittent as the Reynolds number increased [13].Siggia made a numerical simulation and found that 95% of the energy dissipation is concentrated in a tiny region of space [15].Turbulent mixing of a passive scalar is an extremely intermittent phenomenon [21].Intermittency has been shown to be one of the fundamental mechanisms of turbulence.Conclusions regarding clear-air echoes are overshadowed by theoretical defects.
Although it is essentially appropriate to use Kolmogorov's theory in atmospheric science, it is a fact that non-Kolmogorov turbulence is widely present in the boundary layer [22][23][24][25][26][27].Experiments show that the Kolmogorov theory is sometimes incomplete to describe atmospheric statistics properly, in particular, in portions of the troposphere and stratosphere [25].The power-law exponent for the inverse spatial frequency dependence has been observed experimentally to be both larger and smaller than the value of 11/3 that derives from Kolmogorov's model [24].Consortini, Ronchi, and Stefanutti illustrated in the laboratory that the statistics of laser beam wander differed for horizontal and vertical orientations, implying the presence of turbulent anisotropy [23].Other observations also have shown that the atmospheric structures of refractive index irregularities often differ between the horizontal and vertical planes at the same height and distance [27].This means that a nonzero value of differential reflectivity is also consistent with turbulent scattering.Additionally, skepticism towards scattering mechanisms has been exacerbated by these contradictions.
With the continued decline of global species biodiversity, China wants to shoulder more environmental responsibilities as a global economic power and implements the strictest possible systems for environmental protection.Ecological monitoring is an important part of environmental protection.Similar to the Next-Generation Weather Radar (NEXRAD) network in the US, the China Meteorological Administration (CMA) deployed the China Next-Generation Weather Radar (CINRAD) network.The CMA wants to use the radar network for monitoring the movement and abundance of animals in the airspace.However, a few irrationalities were found in the monitoring.The characteristics of clear-air echoes are not completely consistent with the law of seasonal biological activities in China.The source of the contradictions seems to be the influence of turbulence, which can be disregarded in clear-air echoes.Therefore, a rethinking of scattering mechanisms regarding clear-air echoes is needed for an expanded set of causes on the cause of clear-air echoes.
If Bragg scattering was the cause of echoes, the dependence of the echo strength on the radar wavelength would be expected.Since 2015, the Beijing Meteorological Service has built multiple X-POL radars (wavelength λ~3 cm) in the observation coverage area of the CINRAD/SA radar (λ~10 cm) in the Daxing district.Moreover, the CINRAD/SA radar has completed a polarization upgrade in April 2021.The feasible conditions for studying the cause of clear-air echoes in Beijing have been met.Thus, this study is focused on determining the cause of clear-air echoes and their scattering mechanism in Beijing to clarify the mechanism of clear-air echo causes and help the quantitative observation of biology.Observations of the multi-time and dual-wavelength characteristics of clear-air echoes were used to confirm the diagnosis of clear-air echoes.A troposcatter mechanism was introduced to provide some explanations of the observed phenomena.
Section 2 introduces some basic concepts and theories.The data and methods are described in Section 3. The clear-air echo characteristics are analyzed in Section 4. In Section 5, a theory in communication engineering is used to explain the phenomenon of clear-air echoes.Section 6 presents the conclusion.

Dual-Polarization Radar Products
A dual-polarization weather radar is an advanced sensor with a high observation accuracy and many products.It can provide multiple products at the same time.The reflectivity factor (denoted by Z in dBZ) is the most conventional product in radar meteorology, which is a more meteorologically meaningful way of expressing the radar reflectivity η.The radar's property is independent of the radar's wavelength in the case of small scatterers such as precipitation.The reflectivity factor Z is expressed as where λ is the wavelength; K 2 is a dielectric constant and is often taken to be 0.9 for water and 0.2 for ice, respectively.Another conventional product of the weather radar is the Doppler velocity (or called the radial velocity) which is the velocity vector of a scattering object along the radial direction as observed by the Doppler radar.In general, the Doppler velocity is positive for object motion away from the radar.
Prior to the dual-polarization upgrade, the weather radar transmitted and received in a purely horizontal polarization and could not measure any polarization-dependent attributes of a target.When the upgrade was finished, the radar could transmit and receiving two orthogonal signals.The ratio of the received wave powers between horizontal and vertical polarization is called the differential reflectivity (Z DR in units of dB), which helps to infer the shape of the scatterer.Z DR can be defined as the difference between the measured radar reflectivity factor in the horizontal polarization (Z H in dBZ) and the vertical polarization (Z V in dBZ): The correlation coefficient is another useful polarimetric product which depends on the similarity of the received signals at the horizontal and vertical polarizations across multiple pulses.It is a measure of the variety of hydrometeor shapes in a pulse volume.
There are also some other radar products that have not been introduced, and their definitions can be found in the Glossary of Meteorology, American Meteorological Society (http://glossary.ametsoc.org/wiki/climatology,accessed on 27 January 2023).

Turbulence
During the 1940s, Kolmogorov developed a model to illustrate how energy is transported from large-scale turbulent eddies to small-scale turbulent eddies [9,10].Kolmogorov's model provides a spatial power spectrum for index of refraction fluctuations.
Kolmogorov proposed that turbulence in the inertial subrange would reach a statistical equilibrium which is called the "homogeneous isotropic turbulence".Kolmogorov introduced a structure function defined as the squared mean of the difference between the meteorological elements from two independent points in space.For Kolmogorov's turbulence, the structure function of the index of refraction is where n(r 1 ) is the index of refraction at point r 1 ; r = |r|; C 2 n is the refractive index structure constant; and l 0 and L 0 are the inner and outer scales, respectively.
Tatarskii applied Kolmogorov's model and concluded that the structure constant of the refractive index C 2 n is related to meteorological parameters, as shown in Equation ( 4) [28]: where a 2 is a constant with the laboratory value being 2.8; L 0 is the outer scale of turbulence; T is the air temperature; P is the atmosphere pressure; θ is the potential temperature; and h is the altitude.The conversion relationship between θ and T is expressed in the form of Equation ( 5): In Kolmogorov's model, there is no preferred direction through the turbulence, and the turbulent fluctuations of the refractive index are homogeneous, statistically stationary, and isotropic.However, some experimental results have shown that atmospheric turbulence may not always obey Kolmogorov's law [24][25][26].The power spectrum exponent has been observed experimentally to be both larger and smaller than the value from Kolmogorov's model.Since Batchelor and Townsend's first observations of intermittency in 1949 [13], which found that turbulence and its energy dissipation are not space-filling but are intermittent in space, various turbulence theories have attempted to account for and reduce the intermittency geometry of turbulent dissipation [11,12,15,16,[19][20][21].

Bragg Scattering
Bragg scattering is caused by turbulent inhomogeneities with sizes around one-half of the transmitted wavelength of a radar [29].Technically, the atmospheric turbulence with refractivity gradients is essentially a dipole and causes scattering.Ottersten provides the radar reflectivity's relationship with the atmospheric structure constant of the refractive index C 2 n and the radar wavelength λ [8].The radar reflectivity η (or radar cross-section per unit volume) is given in Equation ( 6): According to Equation (4), the reflectivity factor Z is given in Equation ( 7): Therefore, the differential reflectivity Z DR can be equated to C 2 n as Based on isotropic turbulence, C 2 n is equal in the horizontal and vertical directions, and the value of Z DR will be zero in conventional dB form.Additionally, the ratio of the Z values for two radar wavelengths (also called dual-wavelength ratio, DWR) is Wilson used Equation ( 9) to study clear-air echoes over Florida and Colorado and concluded that the clear-air echoes over these areas were caused by living creatures [3].
However, Equation ( 6) is not the only relation between the radar reflectivity and the radar wavelength.The value of the radar reflectivity is variable based on the different applied turbulence theories, and Equation ( 4) is based on the Kolmogorov-Obukhov theory.According to the Villars-Weisskopf theory, the radar reflectivity's relation is given in Equation (10) [30]: where C is constant.Researchers hold different views regarding the relationship, which changes from η ∝ λ −1/3 to η ∝ λ, and the ratio of the Z values for the two radar wavelengths changes according to the researcher's view.

Biological Scattering
Unlike raindrops, biological scatterers have complex shapes that result in highly aspect-dependent scattering characteristics.As a result, radar cross-sections (RCS) are normally modeled by prolate spheroids of equivalent mass [7,31].It is apparent that insects and birds, which are highly non-spherical, would exhibit large Z DR signals and low copular correlation coefficient values.
However, there are polarimetric differences between birds and insects.Insects often have a high Z DR (up to 10 dB) and a relatively low differential phase, while birds may have a lower Z DR (1 to 3 dB) and a much larger differential phase [5].Moreover, for both types of echoes, the cross-correlation coefficient is between 0.3 and 0.5, which is lower than the hydrometeorological signal.Polarimetry becomes a technical standard of the application of a dual-polarization radar to delineate meteorological and non-meteorological signals in the areas of aeroecology.
Based on fuzzy logic, some researchers have differentiated bio-scatterers by using typical values of polarimetric variables [32].Kilambi proposed an estimate of the depolarization ratio for separate types of echoes [33].Overall, the polarimetric characteristics of bio-scatterers are the primary means by which to solve this problem.
For radars that cannot measure the polarization characteristics, the difference derived from the Doppler velocity is another indicator of the scatterer type.Insects reasonably represent actual winds, whereas birds have a much larger independent velocity component [34].Therefore, the radial velocity standard deviation obtained from the velocity-azimuth display retrieval is an indicator of migration.For high migration densities, the raw weather radar wind vectors could be up to 15 m/s for birds and 6 m/s for strong insects [7].However, for micro-insects, which are weak flyers, their motions are primarily wind-driven and behave as quasi-passive wind tracers.Thus, migrating and wandering birds are the main targets of Doppler wind measurement.
Biological scattering and Bragg scattering are always rivals.Compared with Bragg scattering discussed above, the ratio of the Z values for two radar wavelengths is more complicated.When the physical diameter of the spherical particles is considerably smaller than the radar wavelength (approximately 6.25 mm for the S-band), Rayleigh scattering can be used.On the other hand, for larger particles, Mie scattering occurs.The different scattering makes the relation between their size and Z values far from straightforward.Provided that Mie scattering is occurring at one or both of the wavelengths, the ratio of the Z values with a spherical diameter is nonlinear [3].
For example, using prolate spheroids of a spinal cord dielectric, a new model that is closer to real insects [35], the RCS could be simulated by the method of moments, as shown in Table 1, according to biological datasets provided by the Chinese Academy of Agricultural Sciences.Additionally, based on the simulation, the DWRs were 1.9 dB, 13.8 dB, and 17.0 dB between the wavelengths of the S-band and X-band.Further, all birds and many common insects are above the 2 mm size threshold for Mie scattering at the X-band.Measurements of volumes containing multiple scatterers are likely often biased toward the characteristics of the largest scatterers [31].Thus, comparing the value of the DWR is a valid way to be certain of the dominant mechanism of scattering.

Instruments
The data of clear-air echoes used in this study were collected by China's New Generation Weather Radar (CINRAD/SA radar) and three X-band dual-polarization (X-POL) radars with the same technical parameters [36].CINRAD/SA was developed from the American WSR-88D (NEXRAD) through a joint agreement between the two countries [37].The CINRAD/SA radar of Beijing is located in the Daxing district of Beijing and was fully upgraded with polarimetric capabilities in April 2021.The products of CINRAD/SA have a radial distance resolution of 250 m and an azimuthal resolution of 1 degree.The volume coverage pattern 21 (VCP21) scan mode was selected which sweeps 9 elevation angles of 0.5, 1.5, 2.4, 3.4, 4.3, 6.0, 9.9, 14.6, and 19.5 degrees in 6 min.
The X-POL radars were built by the Beijing Meteorologic Service to improve radar usage in weather monitoring.The stations of the X-POL radars used in this study are in the Fangshan, Shunyi, and Tongzhou districts of Beijing (the BJXFS, BJXSY, and BJXTZ sites, respectively), located around the CINRAD/SA radar sites.The positions of the four radars and terrain are shown in Figure 1.The scan strategy for the X-Pol radars is the same as the CINRAD/SA radar.The detailed system characteristics of the CINRAD/SA radar and the X-POL radars are shown in Table 2.
The operations of all these radars are under the supervision of the CMA.The CIN-RAD/SA and the X-POLs are also calibrated weekly and monthly according to the technical standard of the CMA, which includes the system internal calibration, the receiving link calibration, the rotary joint calibration, and others.Some meteorological parameters were used in this study.The profiles of the air temperature, the relative humidity, and the absolute humidity were collected using a microwave radiometer (RPG-HATPRO-G5, Meckenheim, Germany).The range resolution of the microwave radiometer was 50 m for heights below 1 km, 100 m for heights between 1 km and 2 km, and 250 m for heights from 2 km to 10 km, and the temporal resolution was 1 s.The wind vertical velocity and the wind shear were retrieved using a 3D Doppler wind lidar (Windcube 100 s, Leosphere, Saclay, France).The wind measurement products had a spatial resolution of 25 m with a temporal resolution of 20 s, and the products were not used when the carrier-to-noise ratio was less than −30 dB.A radar wind profiler (RWP) CFL-06, which was manufactured by the 23rd Institute of the China Aerospace Science and Industry Corporation, was also used to detect and process the profiles of the horizontal wind speed and direction.The radar wind profiler operates in the L-band (1290 MHz) and provides data with a vertical resolution of 120 m for heights between 150 m and 2.9 km.The microwave radiometer, the wind lidar, and the wind profile radar were located several tens of meters from the S-band weather radar in the same observation field.The operations of all these radars are under the supervision of the CMA.The CINRAD/SA and the X-POLs are also calibrated weekly and monthly according to the technical standard of the CMA, which includes the system internal calibration, the receiving link calibration, the rotary joint calibration, and others.Some meteorological parameters were used in this study.The profiles of the air temperature, the relative humidity, and the absolute humidity were collected using a microwave radiometer (RPG-HATPRO-G5, Meckenheim, Germany).The range resolution of the microwave radiometer was 50 m for heights below 1 km, 100 m for heights between 1 km and 2 km, and 250 m for heights from 2 km to 10 km, and the temporal resolution was 1 s.The wind vertical velocity and the wind shear were retrieved using a 3D Doppler wind lidar (Windcube 100 s, Leosphere, Saclay, France).The wind measurement products had a spatial resolution of 25 m with a temporal resolution of 20 s, and the products were not

Preprocessing
Preprocessing was needed before data comparison between multiple radars because of the differences temporally and spatially.It was necessary to preprocess and use multiple actual sounding data to create a time-height cross-section for the comparison about the DWR.The time-height cross-section is composed of vertical profiles of continuous time.Thus, the calculation of the vertical profile is introduced in the following sections.

Threshold
It is commonly believed that the signal of clear-air echoes is generally weak.A distorted signal would hamper the estimation of echoes.Hence, first, the signal-to-noise ratio (SNR) is set as the standard of thresholding.The minimum value of the SNR thresholding is 6 dB.
Secondly, the setting of the threshold needs to consider the limits of the radar's minimum detectable reflectivity.If the echoes are out of the range of the radar system's designed criteria, weak echoes may be distorted.The distortion values of the echoes then introduce errors into the statistics and the comparison.Meanwhile, the different minimums of the reflectivity values on different radar systems may affect the statistics as well.On the other hand, an excessive detection range would reduce the effectiveness of the comparison with the vertical observation system.Consequently, the radial distance from the radar to the objective is limited.
As shown in Table 2, the minimum detectable reflectivity of the CINRAD/SA and X-POLs was less than −7.5 dBZ at 50 km and 5 dBZ at 60 km.Therefore, balancing the minimum detectable reflectivity and the volume of the data, the range threshold of the data was 59 km for the S-band and 17 km for the X-band.The range threshold was set according to the radar equation to ensure that the radars with different wavelengths had the same minimum detectable reflectivity, which was −6 dBZ.

Vertical Profiles
The vertical profile of the reflectivity factor is useful to estimate the rainfall intensity because of the complexity of the vertical structure of radar echoes [38][39][40][41].From radar data recorded at multiple elevation angles, the mean value of the reflectivity at each altitude can be calculated.Thus, the mean vertical profile of reflectivity (MVPR) can be easily extracted from volume-scan data.Unlike the MVPR, which focuses on the precision of the radar precipitation estimation, a method for determining the state of clear-air echoes needs to be presented.
For weak clear-air echoes, an extreme value can cause fluctuations in the mean value.The limited number of antenna elevation angles also introduces a discretization of the sampling of echoes and lessens the accuracy of the profile.Therefore, a minor alteration is being made to the MVPR.
Assuming the value of the reflectivity factor in a certain altitude range obeys the Gaussian distribution, the expected value of the reflectivity factor can be used to accurately estimate the state of the clear-air echoes at the sampling altitude.By applying Gauss curve fitting to the frequency distribution function of the reflectivity factor, the expected value can be obtained.The fitting uses the bi-square method for robustness, and the adjusted R-square of the effective fitting needs to be larger than 0.95.
Of special note is that a multiplicative adjustment factor is applied to the frequency distribution function because of the characteristics of the radar beam that make the volume of echoes different.The adjustment factor (called the weight factor) is the ratio of the single-sampling volume to the sum of the volume of the scan in the altitude range, or simply denoted by the sampling volume of the single echo.
It is also noticed that the fit coefficients are affected by the number of samples, so the sample size of the reflectivity factor needs to be checked at the sampling altitude.The minimum ratio of the sample size is set to 10%, which means the ratio of the valid sample to the total of the scan volume at the sampling altitude needs to be greater than one-tenth experientially in this study.

Dual-Wavelength Ratio
Unlike the differential reflectivity measured by means of two orthogonal signals that differ in polarization, the DWR describes the ratio of the radar reflectivity measured with two signals of different wavelengths.In Section 2, it has been found that biological scattering and turbulent scattering have different characteristics of the DWR.The DWR of Bragg scattering is dependent on the radar wavelength, as shown in Equation (9).Moreover, the DWR change in biological scattering is unpredictable because of the biologically complex shape.Hence, the predominant scattering mechanism can be determined by the DWR.
However, it is also known that the value of Z from biological scatterers has an asymmetric pattern which depends on the angle between the radar beam and the main orientation of the biological scatterers.Moreover, values of the DWR may be more complicated and confusing because of the asymmetric pattern.Thus, it needs to examine the effect of the asymmetric pattern.
For the examination, the RCS of biological scatterers is simulated by the computer utilizing FEKO simulation software in this study [35].The simulative value of RCS is dependent on many factors.One of the factors is the incident angle of the radar beam.Although biological scatterers appear randomly at each azimuth angle of the radar, the incident angle of the simulation is only decided by the angle between the incident beam and the biological scatterer.By changing the coordinate origin of the simulation from the location of the scatterer to the location of the radar, the incident angle can be converted to the azimuth angle when the directions of the biological scatterers are same.
Other parameters of the simulation are obtained from Table 1 based on the model of prolate spheroids of the spinal cord dielectric.Further, supposing that insects fly horizontally in the east-west direction and the density is one per cubic meter, the Z value can be calculated by the simulated RCS shown in Figure 2. Figure 2 exemplifies the asymmetric pattern of the Z values and is almost consistent with previous studies [42].It is found that the asymmetric pattern may lead to confusion about the conclusions of Section 2.

Dual-Wavelength Ratio
Unlike the differential reflectivity measured by means of two orthogonal signals that differ in polarization, the DWR describes the ratio of the radar reflectivity measured with two signals of different wavelengths.In Section 2, it has been found that biological scattering and turbulent scattering have different characteristics of the DWR.The DWR of Bragg scattering is dependent on the radar wavelength, as shown in Equation ( 9).Moreover, the DWR change in biological scattering is unpredictable because of the biologically complex shape.Hence, the predominant scattering mechanism can be determined by the DWR.
However, it is also known that the value of  from biological scatterers has an asymmetric pattern which depends on the angle between the radar beam and the main orientation of the biological scatterers.Moreover, values of the DWR may be more complicated and confusing because of the asymmetric pattern.Thus, it needs to examine the effect of the asymmetric pattern.
For the examination, the RCS of biological scatterers is simulated by the computer utilizing FEKO simulation software in this study [35].The simulative value of RCS is dependent on many factors.One of the factors is the incident angle of the radar beam.Although biological scatterers appear randomly at each azimuth angle of the radar, the incident angle of the simulation is only decided by the angle between the incident beam and the biological scatterer.By changing the coordinate origin of the simulation from the location of the scatterer to the location of the radar, the incident angle can be converted to the azimuth angle when the directions of the biological scatterers are same.
Other parameters of the simulation are obtained from Table 1 based on the model of prolate spheroids of the spinal cord dielectric.Further, supposing that insects fly horizontally in the east-west direction and the density is one per cubic meter, the  value can be calculated by the simulated RCS shown in Figure 2. Figure 2 exemplifies the asymmetric pattern of the  values and is almost consistent with previous studies [42].It is found that the asymmetric pattern may lead to confusion about the conclusions of Section 2. However, the asymmetric pattern does not affect the results of the method in Section 3.2.3.The vertical profiles of the DWR between the S-and X-band are calculated and shown in Figure 3.The computation is deduced in terms of the scan mode of VCP21. Figure 3 shows that the vertical profiles of the DWR are same on each level of height, and the values are close to the simulation values of the insects' body sizes.This is because the RCS of the insect size changes more slowly with angle.However, it is also found that the Z difference values of large insects and the values of turbulence are close.Thus, the echo of turbulence is likely confused with that of large insects in the DWR.
difference values of large insects and the values of turbulence are close.Thus, the echo of turbulence is likely confused with that of large insects in the DWR.
Fortunately, large insects have a much stronger flight ability, and flight speed is correlated with body size in animals [43].Research on the migratory behavior of armyworms expresses that the flight speed of armyworms during migration is greater than 4 m/s [44].Thus, large insects can be easily distinguished by the velocity azimuth display (VAD).Analysis of the VAD is necessary to exclude confusion with large insects.

Plan Position Indicator
For the sake of discussion, the radar data from May 2 were taken as the case for further study.The focus of this case was from 11:30 (UTC), 2 May 2021, when the clear-air echo first appeared in the radar display, until it began to disappear at 20:30, when the echo became less distinct.The products in Figure 4 show that the echo was characterized by a low , a  higher than that typically observed in meteorological echoes, and a correlation coefficient  lower than that observed in meteorological echoes.Generally, similar echoes have been observed in bird scattering comprising many species flying over Southern Kansas from Wichita, KS, USA [5].Fortunately, large insects have a much stronger flight ability, and flight speed is correlated with body size in animals [43].Research on the migratory behavior of armyworms expresses that the flight speed of armyworms during migration is greater than 4 m/s [44].Thus, large insects can be easily distinguished by the velocity azimuth display (VAD).Analysis of the VAD is necessary to exclude confusion with large insects.

Plan Position Indicator
For the sake of discussion, the radar data from May 2 were taken as the case for further study.The focus of this case was from 11:30 (UTC), 2 May 2021, when the clear-air echo first appeared in the radar display, until it began to disappear at 20:30, when the echo became less distinct.The products in Figure 4 show that the echo was characterized by a low Z, a Z DR higher than that typically observed in meteorological echoes, and a correlation coefficient ρ HV lower than that observed in meteorological echoes.Generally, similar echoes have been observed in bird scattering comprising many species flying over Southern Kansas from Wichita, KS, USA [5].Some researchers have indicated that the ρ HV varies with scatterer orientation, with birds flying away producing low ρ HV values and birds flying toward producing higher ρ HV values.However, the velocity-azimuth display (VAD) confirmed that there was no intense bird migration.Analysis of the VAD is used to determine spatially averaged kinematic properties of the velocity field.According to the function of the radar azimuth angle, the VAD of the wind field resembles a sine function.Moreover, during intense bird migrations, a huge difference in scattering of the observed radial velocity around the modeled radial velocities will be evident.
Figure 5 shows a high-quality wind VAD, which was the same passage of a cold front as observed by an operational weather radar in De Bilt [45], Netherlands.Only a small part of the residual error exhibited signs of bird activity.Since the nonzero value of the Doppler velocity due to bio-scatterers is representative of biological target movement, the temporal change in Z is shown in Figure 6 to check for movement.
Figure 6 shows a continuous change in Z; the echoes began to increase rapidly at 11:30 and then remained basically unchanged after 12:06.It is strange that the echoes changed with the radar station as the center but not the "habitat", and the echoes with a larger Z also maintained their appearance.Whereas the value of the Doppler velocity was rhythmical, shown in Figure 7, the characteristics of the spatial distribution were also unchanged after 48 min, indicating that the scatterers did not move.
To confirm the state of the scatterers' movement, in Figure 8, the range-height crosssection of Z is displayed in the azimuth of the wind direction.As shown in Figure 8, the signal was enhanced along the height at each range and did not change following wind direction.The phenomena discovered by Adriaan were not observed [46].Paradoxically, the scatterer remained in the air, but the Doppler velocity was a nonzero value.Some researchers have indicated that the  varies with scatterer orientation, with birds flying away producing low  values and birds flying toward producing highe  values.However, the velocity-azimuth display (VAD) confirmed that there was no intense bird migration.Analysis of the VAD is used to determine spatially averaged kin ematic properties of the velocity field.According to the function of the radar azimuth an gle, the VAD of the wind field resembles a sine function.Moreover, during intense bird migrations, a huge difference in scattering of the observed radial velocity around the mod eled radial velocities will be evident.
Figure 5 shows a high-quality wind VAD, which was the same passage of a cold fron as observed by an operational weather radar in De Bilt [45], Netherlands.Only a smal part of the residual error exhibited signs of bird activity.Since the nonzero value of th Doppler velocity due to bio-scatterers is representative of biological target movement, th temporal change in  is shown in Figure 6 to check for movement.Figure 6 shows a continuous change in ; the echoes began to increase rapidly at 11:30 and then remained basically unchanged after 12:06.It is strange that the echoes changed with the radar station as the center but not the "habitat", and the echoes with a larger  also maintained their appearance.Whereas the value of the Doppler velocity was rhythmical, shown in Figure 7, the characteristics of the spatial distribution were also unchanged after 48 min, indicating that the scatterers did not move.x FOR PEER REVIEW 13 of 25 To confirm the state of the scatterers' movement, in Figure 8, the range-height crosssection of  is displayed in the azimuth of the wind direction.As shown in Figure 8, the signal was enhanced along the height at each range and did not change following wind direction.The phenomena discovered by Adriaan were not observed [46].Paradoxically, the scatterer remained in the air, but the Doppler velocity was a nonzero value.

Time-Height Cross-Section
Movement is the most important feature of bio-scatterers.To reveal the causes of clear-air echoes, the radar data were displayed as a time-height cross-section.The signal of the clear-air echoes showed significant diurnal variations, as shown in Figure 9.It was thought that the activities of nocturnal creatures caused greater echoes.However, the rapid growth in the signal in the time-height cross-section caught our attention.Surprisingly, the value of  grew fast.The 23 dBZ echo only took 13 min to increase from 500 m to 1000 m, and the increasing velocity was 0.64 m/s.When using the wind lidar, as Figure 10 shows, the vertical velocity of the wind was below 0.5 m/s and was smaller than the velocity of the echoes.This indicates that a small bio-scatterer, which does not fly well, could not be the scatterer.Despite the impractical state of motion, perhaps these clear-air echoes were still caused by large bio-scatterers.However, the vertical profiles of  shown in Figure 9 deny this view.

Time-Height Cross-Section
Movement is the most important feature of bio-scatterers.To reveal the causes of clear-air echoes, the radar data were displayed as a time-height cross-section.The signal of the clear-air echoes showed significant diurnal variations, as shown in Figure 9.It was thought that the activities of nocturnal creatures caused greater echoes.However, the rapid growth in the signal in the time-height cross-section caught our attention.Surprisingly, the value of Z grew fast.The 23 dBZ echo only took 13 min to increase from 500 m to 1000 m, and the increasing velocity was 0.64 m/s.When using the wind lidar, as Figure 10 shows, the vertical velocity of the wind was below 0.5 m/s and was smaller than the velocity of the echoes.This indicates that a small bio-scatterer, which does not fly well, could not be the scatterer.Despite the impractical state of motion, perhaps these clear-air echoes were still caused by large bio-scatterers.However, the vertical profiles of Z DR shown in Figure 9 deny this view.It is known that  is the radar reflectivity difference between the horizontal and vertical polarization and represents the dimension of the scatterer.If the measurements of the volumes contain multiple scatterers,  will be biased toward the characteristics of the largest scatterers.In Figure 9, the profiles of  increased with height and were temporarily greater than 2.5 dB above an altitude of 1200 m.Scatterers with a 2.5 dB  do not climb from the ground and do not land.This implies that the echoes did not float from the ground into the air, but they suddenly appeared in the middle of the air without immigration.

Velocity Analysis
Although we have preliminarily analyzed the VAD in Section 4.1, to avoid confusion with large insects and birds, the Doppler velocity and the results of the VAD analysis are displayed in Figure 11 again.In Figure 11, the velocity fields, which were obtained from  It is known that  is the radar reflectivity difference between the horizontal and vertical polarization and represents the dimension of the scatterer.If the measurements o the volumes contain multiple scatterers,  will be biased toward the characteristics o the largest scatterers.In Figure 9, the profiles of  increased with height and were tem porarily greater than 2.5 dB above an altitude of 1200 m.Scatterers with a 2.5 dB  do not climb from the ground and do not land.This implies that the echoes did not float from the ground into the air, but they suddenly appeared in the middle of the air without im migration.

Velocity Analysis
Although we have preliminarily analyzed the VAD in Section 4.1, to avoid confusion with large insects and birds, the Doppler velocity and the results of the VAD analysis ar displayed in Figure 11 again.In Figure 11, the velocity fields, which were obtained from It is known that Z DR is the radar reflectivity difference between the horizontal and vertical polarization and represents the dimension of the scatterer.If the measurements of the volumes contain multiple scatterers, Z DR will be biased toward the characteristics of the largest scatterers.In Figure 9, the profiles of Z DR increased with height and were temporarily greater than 2.5 dB above an altitude of 1200 m.Scatterers with a 2.5 dB Z DR do not climb from the ground and do not land.This implies that the echoes did not float from the ground into the air, but they suddenly appeared in the middle of the air without immigration.

Velocity Analysis
Although we have preliminarily analyzed the VAD in Section 4.1, to avoid confusion with large insects and birds, the Doppler velocity and the results of the VAD analysis are displayed in Figure 11 again.In Figure 11, the velocity fields, which were obtained from the VAD analysis, are similar to the radar Doppler velocity.The residual values of the velocity fields are shown in Figure 12a. in total, 86.4% of the residual values are less than 4 m/s, and 65.3% of the values are less than 2 m/s.The echoes were unlikely caused by large insects and birds mainly because large insects and birds have quite a great flight speed and would make the speed field messy and produce great deviations.12b, which is common in the results of the VAD wind profile [47], the changes in the speed profiles are consistent, generally.The profiles of the CINRAD/SA indicate that clear-air echoes are moving with the wind speed.However, the speed of biological autonomous movement is not detected by the weather radar.This means that these echoes are unlikely to be caused by creatures with a strong flight ability such as birds and large insects.

Comparison of the S-Band and X-Band
In Sections 4.1 and 4.2, some contradictions of bio-scatterers were revealed, and a quantitative analysis of the echoes was thus needed to convincingly demonstrate the cause    12b, which is common in the results of the VAD wind profile [47], the changes in the speed profiles are consistent, generally.The profiles of the CINRAD/SA indicate that clear-air echoes are moving with the wind speed.However, the speed of biological autonomous movement is not detected by the weather radar.This means that these echoes are unlikely to be caused by creatures with a strong flight ability such as birds and large insects.

Comparison of the S-Band and X-Band
In Sections 4.1 and 4.2, some contradictions of bio-scatterers were revealed, and a quantitative analysis of the echoes was thus needed to convincingly demonstrate the cause of the clear-air echoes.A comparison of the reflectivity factors of the two bands is another The comparison of the speed profiles with the RWP and the CINRAD/SA confirms the conclusion.The profiles of the horizontal wind speed which were sensed by the RWP are not different from the speed profiles calculated by the CINRAD/SA, as shown in Figure 12b.Despite the system deviation in Figure 12b, which is common in the results of the VAD wind profile [47], the changes in the speed profiles are consistent, generally.The profiles of the CINRAD/SA indicate that clear-air echoes are moving with the wind speed.However, the speed of biological autonomous movement is not detected by the weather radar.This means that these echoes are unlikely to be caused by creatures with a strong flight ability such as birds and large insects.

Comparison of the S-Band and X-Band
In Sections 4.1 and 4.2, some contradictions of bio-scatterers were revealed, and a quantitative analysis of the echoes was thus needed to convincingly demonstrate the cause of the clear-air echoes.A comparison of the reflectivity factors of the two bands is another method that can be used to determine the cause of clear-air echoes.Since the DWR of the turbulence echoes in different bands was more regular than biological, the time-height cross-section of Z using the data from BJXFS is shown in Figure 13, and the DWR between the S-band and the X-band is displayed in Figure 14.
turbulence echoes in different bands was more regular than biological, the time-height cross-section of  using the data from BJXFS is shown in Figure 13, and the DWR between the S-band and the X-band is displayed in Figure 14.
In Figure 13, intuitively, the signal of the X-band differed from the S-band signal in the value of ; the value of the X-band was much smaller than that of the S-band.The signal of the X-band also did not have an as pronounced diurnal variation as that of the Sband echoes.Although the signal of the X-band echoes appeared as fast as those of the Sband signal at dusk, the signal did not last long and slowly vanished from the top to the bottom after 13:00.Meanwhile, the value of  in Figure 13 is low at lower altitudes and increases with height, which is different in Figure 9.This is due to the minimum scale of turbulence which increases with height.When the minimum scale exceeds the half-wavelength of the X-band, which is still much smaller than the wavelength of the S-band, the scattering vanishes before the refractivity gradient weakens.Hence, the  value of the Xband seems to increase progressively without reducing.Figure 14a shows the time-height cross-section of the DWR between the  values of the CINRAD/SA and X-POLs, and Figure 14b shows the normalized frequency histograms for the radar fields in Figure 14a.The histogram of the DWR shows the characteristic values near 21.5 dB.Moreover, the  value distributions of the S-band and X-band also show that the peaks of the  values are at intervals of about 22 dB (in Figure 15), turbulence echoes in different bands was more regular than biological, the time-height cross-section of  using the data from BJXFS is shown in Figure 13, and the DWR between the S-band and the X-band is displayed in Figure 14.
In Figure 13, intuitively, the signal of the X-band differed from the S-band signal in the value of ; the value of the X-band was much smaller than that of the S-band.The signal of the X-band also did not have an as pronounced diurnal variation as that of the Sband echoes.Although the signal of the X-band echoes appeared as fast as those of the Sband signal at dusk, the signal did not last long and slowly vanished from the top to the bottom after 13:00.Meanwhile, the value of  in Figure 13 is low at lower altitudes and increases with height, which is different in Figure 9.This is due to the minimum scale of turbulence which increases with height.When the minimum scale exceeds the half-wavelength of the X-band, which is still much smaller than the wavelength of the S-band, the scattering vanishes before the refractivity gradient weakens.Hence, the  value of the Xband seems to increase progressively without reducing.Figure 14a shows the time-height cross-section of the DWR between the  values of the CINRAD/SA and X-POLs, and Figure 14b shows the normalized frequency histograms for the radar fields in Figure 14a.The histogram of the DWR shows the characteristic values near 21.5 dB.Moreover, the  value distributions of the S-band and X-band also show that the peaks of the  values are at intervals of about 22 dB (in Figure 15),  In Figure 13, intuitively, the signal of the X-band differed from the S-band signal in the value of Z; the value of the X-band was much smaller than that of the S-band.The signal of the X-band also did not have an as pronounced diurnal variation as that of the S-band echoes.Although the signal of the X-band echoes appeared as fast as those of the S-band signal at dusk, the signal did not last long and slowly vanished from the top to the bottom after 13:00.Meanwhile, the value of Z in Figure 13 is low at lower altitudes and increases with height, which is different in Figure 9.This is due to the minimum scale of turbulence which increases with height.When the minimum scale exceeds the half-wavelength of the X-band, which is still much smaller than the wavelength of the S-band, the scattering vanishes before the refractivity gradient weakens.Hence, the Z value of the X-band seems to increase progressively without reducing.
Figure 14a shows the time-height cross-section of the DWR between the Z values of the CINRAD/SA and X-POLs, and Figure 14b shows the normalized frequency histograms for the radar fields in Figure 14a.The histogram of the DWR shows the characteristic values near 21.5 dB.Moreover, the Z value distributions of the S-band and X-band also show that the peaks of the Z values are at intervals of about 22 dB (in Figure 15), whereas the Z DR distribution does not peak at 0 dB.In Figure 15, the frequency distribution of the Z value resembles the Gaussian distribution and slightly changes with height, which is in line with the condition of the hypothesis in Section 3.2.However, it is not expected that the correlation or other products should follow the Gaussian distribution.The distributions of turbulence variables deviate from normality because of the intermittence [48][49][50].
The large widths of these distributions match with the characteristics of bird echoes seemingly.However, the echo is unlikely to be mainly caused by birds for the following reasons.First, the analysis of the VAD has already denied the dominance of birds.Second, for north China, May is almost the end of the spring migration and is close to the period of preincubation, yet the scatterers of the clear-air echo still travel to the north until late July [51,52].Third, there are two international airports located at the azimuth angles of 22° and 188° and the distances of 31 km and 34 km, respectively.The large numbers of birds migrating would obviously be a threat to the safety of the flights if the clear-air echo was caused by birds.Therefore, the suspicion of birds and large insects has been ruled out.
It is worth noting that the CINRAD/SA and X-POL are located along the same migration path.It means that the two radars should observe the same group of scatterers with the same scattering characteristic, moving from upstream (XFS) to downstream (SDX) the whole night, continuously.Although Z values from atmospheric insects may have asymmetric patterns which depend on the angle between the radar beam and the main orientation of the biological scatterers, the time-height cross-section of the Z values will not be affected because the cross-sections are counted by the whole coverage volume of the radar observation and are not dependent on the azimuth.Thus, the asymmetric patterns cannot stain the ratio between the Z values of the CINRAD/SA and X-POL.It is already known that, according to Equation (6) and Equation (10), the  difference in the turbulence of the S-band (2870 MHz, 10.45 cm) and the X-band (9455 MHz, 3.17 cm) is about 19.0 dB and 21.4 dB, respectively.Since many creatures are within the The large widths of these distributions match with the characteristics of bird echoes seemingly.However, the echo is unlikely to be mainly caused by birds for the following reasons.First, the analysis of the VAD has already denied the dominance of birds.Second, for north China, May is almost the end of the spring migration and is close to the period of preincubation, yet the scatterers of the clear-air echo still travel to the north until late July [51,52].Third, there are two international airports located at the azimuth angles of 22 • and 188 • and the distances of 31 km and 34 km, respectively.The large numbers of birds migrating would obviously be a threat to the safety of the flights if the clear-air echo was caused by birds.Therefore, the suspicion of birds and large insects has been ruled out.
It is worth noting that the CINRAD/SA and X-POL are located along the same migration path.It means that the two radars should observe the same group of scatterers with the same scattering characteristic, moving from upstream (XFS) to downstream (SDX) the whole night, continuously.Although Z values from atmospheric insects may have asymmetric patterns which depend on the angle between the radar beam and the main orientation of the biological scatterers, the time-height cross-section of the Z values will not be affected because the cross-sections are counted by the whole coverage volume of the radar observation and are not dependent on the azimuth.Thus, the asymmetric patterns cannot stain the ratio between the Z values of the CINRAD/SA and X-POL.
It is already known that, according to Equation (6) and Equation (10), the Z difference in the turbulence of the S-band (2870 MHz, 10.45 cm) and the X-band (9455 MHz, 3.17 cm) is about 19.0 dB and 21.4 dB, respectively.Since many creatures are within the size threshold for Mie scattering and have different RCSs at different wavelengths, the reflectivity relationships to the wavelength are far from straightforward.Meanwhile, wide discrepancies are found between the characteristics of clear-air echoes and the law of biological activities.Therefore, this proves that the Beijing CINRAD/SA echoes are mainly caused by turbulence.
The DWR statistics for the region of Beijing reveal that, from 1 May 2021 to 20 May 2021, 58% of the DWR was distributed between 18dB and 24dB, which is the characteristic interval of turbulence.Moreover, the clear-air echoes maintain the features of the variations in the Z values and the VAD analysis.It means that nearly 58% of the echoes were probably caused by turbulence and proves that the influence of turbulence is erroneously ignored in clear-air echoes.Almost all these echoes did not show the classic features of Bragg scattering, such as a low Z and zero Z DR .These echoes could wrongly be identified as bioscatterers and cause a misestimation of biomass.The number of flying creatures observed by the weather radar was much smaller than previously thought.

Discussion
Many scholars deny or reject the domination of turbulence in clear-air echoes.First, they propose that a nonzero Z DR indicates that the scattering is due to creatures rather than turbulence because, according to the statement of Kolmogorov's theory that turbulence is homogeneous and isotropic, the Z DR of turbulent echoes should be zero.However, a problem with this viewpoint is that non-Kolmogorov turbulence widely exists in the atmospheric boundary layer.Some observations have already shown that the atmospheric structure constant of the refractive index differs between the horizontal direction and the vertical direction.Therefore, the Z DR cannot support their rejections.We agree that biological scatterers cause large values of Z DR , but the contribution from the effect of turbulence still needs to be investigated.
Second, some scholars believe that because the C 2 n calculated from echoes is larger than the high end of the observed values, clear-air echoes should not be attributed to turbulence.For example, clear-air echoes of a 10 to 20 dB Z require C 2 n to be greater than 10 −11 m −2/3 , which is far above the observed C 2 n .The phenomenon of clear-air echoes cannot be explained by Bragg scattering.Therefore, it is assumed that the domination of clear-air echoes is not turbulence.However, scholars have ignored the effect of another scattering mechanism.We found another reasonable scattering mechanism to explain clear-air echoes in communication and wireless areas.
In communication, tropospheric scattering (also known as troposcatter) is admitted as an efficient propagation method by the Radiocommunication Sector of the International Telecommunication Union, which describes the mechanism by which microwave radio systems inadvertently achieve beyond-the-horizon communications [53].One of the three models of the troposcatter is the reflecting-layers model.The other two models are scattering from turbulence and reflections from an exponential atmosphere.Tropospheric scattering is the result of the combination of these three models.Zhang points out that the reflectivity of the troposcatter can be written as follows [54]: where B, n, and m are constants and are measured by experiments; ε r is the dielectric constant; h is the height; ϕ is the glancing angle; the bar is the sign of the mean; and the relationship between ε r and h can refer to the assumption by K. Bullington [55].Interestingly, the scattering from turbulence is essentially the same as Bragg scattering, but few studies have highlighted the role of reflecting layers in clear-air echoes.
To demonstrate how reflecting layers affect clear-air echoes, we calculated the gradient Richardson number and the intensity of the turbulence (expressed by the ratio of the standard deviation of the wind speed to the mean wind speed) using a microwave radiometer and a wind lidar [56,57].Moreover, the atmospheric profiles were collected by a microwave radiometer.Although the geometrical structure of the reflecting layers remains obscure, it can be inferred that turbulent mixing is detrimental to the reflecting layers because a reflecting layer in the atmosphere is formed by relatively sharp gradients of the refractive index, but turbulent mixing makes the temperature and humidity homogeneous and reduces the gradient of the refractive index.Thus, the Richardson number and the turbulence intensity, which indicate turbulent mixing, are shown in Figure 16 to give the relationship between clear-air echoes and turbulent mixing.In Figure 16, large values of  correspond to large Richardson numbers and turbulent mixing.When the sun sets in the west, the turbulent mixing rapidly weakens.Additionally, the unmixed air masses with relatively sharp gradients become the reflecting layer and generate the scattering signal of the clear-air echoes.In contrast, in the daytime, the strong turbulent mixing breaks the structure of the reflecting layer, and only the turbulent scattering remains.Thus, turbulence has a greater impact on the echoes at night than in the daytime.
The same diurnal variation has been observed in the field-strength variation in the short radio wave propagation in Arizona [58].Other studies have found that troposcatter In Figure 16, large values of Z correspond to large Richardson numbers and turbulent mixing.When the sun sets in the west, the turbulent mixing rapidly weakens.Additionally, the unmixed air masses with relatively sharp gradients become the reflecting layer and generate the scattering signal of the clear-air echoes.In contrast, in the daytime, the strong turbulent mixing breaks the structure of the reflecting layer, and only the turbulent scattering remains.Thus, turbulence has a greater impact on the echoes at night than in the daytime.
The same diurnal variation has been observed in the field-strength variation in the short radio wave propagation in Arizona [58].Other studies have found that troposcatter propagation only occurs when the turbulence scale is larger than the wavelength of the radio signal, and the signal level with a wavelength of 9 cm is much larger than that with a wavelength of 3 cm [58,59].These phenomena have the same characteristics as clear-air echoes, with the value of the reflectivity factor at the S-band being bigger than that at the X-band.Thus, we propose that the reflecting layers cause a diurnal variation in clear-air echoes and enhance the signal of clear-air echoes at night.Therefore, the coexistence of the reflecting layers and Bragg scattering is the reason why the value of C 2 n calculated from the echoes is larger than the value from the theoretical calculations.Further, the reflecting layers explain some of the characteristics of Z DR .For example, the vertical signal of clear-air echoes is weaker than the horizontal signal because turbulent mixing is stronger in the vertical direction, and the vertical reflecting layers cannot easily survive.Thus, the values of Z DR are generally greater than zero.
An interesting note is that the high reflectivity area in the upper-left region of Figure 4 analyzes the characteristics of a low correlation coefficient, and in this area, Z DR presents a mixture of high and low values.This characteristic is ascribed to the atmospheric response to the underlying surface, which belongs to a mountainous region and is quite different from the others in Figure 1.
The rough terrain of mountains makes the turbulence more chaotic and intense, which produces a stronger refractive index and scattering.A thin layer of shear turbulence excited on the shear plane makes the Z DR larger, but the layer is not stable, bringing a mixture of high and low values to the Z DR .However, it is noted that, by upthrust, too high mountains may affect the turbulent scales, which are too large to match up to the radar wavelength for scattering.The relationship between the underlying surface and clear-air echoes is worth further exploration in the future.
Another interesting characteristic of the Z DR is that its value grows weak and becomes close to zero in the day-night shift scenes.We speculate that the turbulent mixing reaches a quasi-equilibrium state, and this will be further investigated in future studies.
Yet despite all this, some scholars still express a slightly different point of view.They claim that the reflection from turbulent air at the S-band has been studied by some studies [29,[60][61][62][63], and it has been well established with polarization radars that turbulent air has Z DR values close to 0 dB at the S-band.However, these echoes were commonly observed at the top of the convective boundary layer (also called the entrainment layer).The entrainment layer is essentially static in stability.Due to the effects of penetrative convection and entrainment, thermals reciprocate in the layer [64].The reciprocation brings turbulence generation on the one hand; on the other hand, the stable air brings the turbulence closer to a locally homogeneous isotropic state.Thus, the observed echoes can be regarded as special cases and cannot represent the whole characteristic of clear-air turbulent echoes.Figure 17 shows the vertical profiles of the Z values and rawinsonde data.It is observed that there is a layer of enhanced reflectivity at the entrainment layer, which is consistent with the known results [63].
Some scholars also conject that dust and other particulate matter can be the cause of the clear-air echoes.At the beginning of this study, we also inferred that nocturnal pollutant accumulation is the reason for the echo diurnal variation because of the uncleanliness of the atmosphere of the megacities.Yet, it is hard to explain the variation in the echo signal.
One piece of evidence is based on seasonal variation.In winter, the air quality is the worst because of heating, which uses fossil fuels, but the clear-air echo is hardly observed by radars during winter.The annual and monthly variation in the sand and dust also showed that the sand and dust weather is most frequent in spring, whereas the signal of the spring clear-air echo was generally weaker than summer and autumn.
lence generation on the one hand; on the other hand, the stable air brings the turbulence closer to a locally homogeneous isotropic state.Thus, the observed echoes can be regarded as special cases and cannot represent the whole characteristic of clear-air turbulent echoes.Figure 17 shows the vertical profiles of the  values and rawinsonde data.It is observed that there is a layer of enhanced reflectivity at the entrainment layer, which is consistent with the known results [63].Some scholars also conject that dust and other particulate matter can be the cause of the clear-air echoes.At the beginning of this study, we also inferred that nocturnal pollutant accumulation is the reason for the echo diurnal variation because of the uncleanliness of the atmosphere of the megacities.Yet, it is hard to explain the variation in the echo signal.
One piece of evidence is based on seasonal variation.In winter, the air quality is the worst because of heating, which uses fossil fuels, but the clear-air echo is hardly observed by radars during winter.The annual and monthly variation in the sand and dust also showed that the sand and dust weather is most frequent in spring, whereas the signal of the spring clear-air echo was generally weaker than summer and autumn.The other piece of evidence comes from the daily change in pollutants.The concentration of pollutants has two peaks because of traffic congestion in the metropolis, which is different from the Z value.The wind, which is related to pollutant diffusion and dust, also does not show a significant correlation with the clear-air echoes.Although it is known that fire plumes can cause clear-air echoes [65,66], it is unrealistic to detect plumes of wild fires in an urban area.Thus, dust and the pollutants are unlikely to be the main causes of clear-air echoes.
The issue of the influence of meteorological factors on troposcatter communication remains unknown and requires future examination.In troposcatter, the value of the signal level depends on the refractive index and its gradient, which are affected by the intensity of the turbulence fluctuation and atmospheric stratification.Thus, Gaoming Zhang proposed that the inversion of temperature leads to the diurnal variation in the signal [54].However, few studies focus on the effect of turbulent mixing on signals level.Thus, we plan to provide a more in-depth explanation of reflecting layers forming and the structure of the nocturnal boundary layer in future studies.Meanwhile, Guifu Zhang pointed that bistatic radars have an advantage in the sensitivity of clear echoes [67].Technically, experiments using bistatic radars are closer to the principle of troposcatter propagation.Moreover, experiments with bistatic radar data would provide more information for the analysis in future research.

Conclusions
In this study, clear-air echoes detected by CINRAD were analyzed to find their causes.Some observations diverge from the previous conclusion that bio-scatterers are the main reason for clear-air echoes.
Echoes with a larger Z are not deformed in the air, even if the scatterers are moving.The change in echoes in the vertical direction is also closer to the switching of the physical mode rather than biological flying.The analysis results of the DWR and the VAD support that turbulence plays an important role in clear-air echoes.In the case of May 2, the frequency distribution of the DWR peaks at 21.5 dB, which is consistent with the theory of turbulence.From 1 May to 20 May, 58% of the DWR between the S-band and the X-band is distributed between 18 dB and 24 dB, which means that more than half of the echoes at night were caused by turbulence.It was confirmed that the influence of turbulence is erroneously ignored in clear-air echoes.The reflecting-layers model of troposcatter propagation is the cause of the clear-air echoes, and this model can explain the main phenomena of the radar observations.

25 Figure 1 .
Figure 1.Distribution of radars (square signs and diamonds) and topography (coloring) of Beijing and its vicinity.The three square signs indicate the locations of the three X-POL radar sites (XFS, XSY, and XTZ).The diamond shows the location of the SDX site.The distance of each X-POL radar relative to the SDX site is labeled.The black dotted circle is the distance circle of the S-band with a 58 km radius, and the radius of the white dotted circles of the X-POL radar sites are 34 km, respectively.These circles show the detection zones of the CINRAD/SA and X-POLs where their minimum detectable reflectivity is less than −5 dBZ.

Figure 1 .
Figure 1.Distribution of radars (square signs and diamonds) and topography (coloring) of Beijing and its vicinity.The three square signs indicate the locations of the three X-POL radar sites (XFS, XSY, and XTZ).The diamond shows the location of the SDX site.The distance of each X-POL radar relative to the SDX site is labeled.The black dotted circle is the distance circle of the S-band with a 58 km radius, and the radius of the white dotted circles of the X-POL radar sites are 34 km, respectively.These circles show the detection zones of the CINRAD/SA and X-POLs where their minimum detectable reflectivity is less than −5 dBZ.

Figure 2 .
Figure 2. Azimuthal dependencies of the reflectivity factor for three species at the S-band (a,c) and X-band (b,d).As in (a,b), the elevation of the radar beam is 0.5 • , which is the minimum elevation of VCP 21, and in (c,d), the elevation is 19.5 • , which is the maximum elevation of VCP 21.

Figure 3 .
Figure 3.The simulated vertical profiles of the DWR between the S-band and X-band with the vertical height.The vertical profiles of turbulence fall in the grey-shaded area because of the different turbulence theories.The parameters of the three species are referenced fromTable 1.All species are major agricultural pests in North China.

Figure 3 .
Figure 3.The simulated vertical profiles of the DWR between the S-band and X-band with the vertical height.The vertical profiles of turbulence fall in the grey-shaded area because of the different turbulence theories.The parameters of the three species are referenced fromTable 1.All species are major agricultural pests in North China.

2 Figure 4 .
Figure 4. Radar products from the CINRAD/SA radar (SDX), Daxing, Beijing, on 2 May 2021 at 13:0 UTC for an elevation angle of 2.4°.The mapped domain is 75 km by 75 km.Some researchers have indicated that the  varies with scatterer orientation, with birds flying away producing low  values and birds flying toward producing highe  values.However, the velocity-azimuth display (VAD) confirmed that there was no intense bird migration.Analysis of the VAD is used to determine spatially averaged kin

Figure 5 .
Figure 5. (a) An example of VAD data from the CINRAD/SA in Beijing.The line is the modeled radial velocities as a sine function of the azimuth, and the dots are the data of the Doppler velocity The data samples are in the range of 30 km and the elevation of 2.4 deg.(b) The residual error (th dot) is essentially less than 2 m/s.The root mean squared error is 1.278, and the adjusted R-squar is 0.9924.

Figure 5 . 25 Figure 6 .
Figure 5. (a) An example of VAD data from the CINRAD/SA in Beijing.The line is the modeled radial velocities as a sine function of the azimuth, and the dots are the data of the Doppler velocity.The data samples are in the range of 30 km and the elevation of 2.4 deg.(b) The residual error (the dot) is essentially less than 2 m/s.The root mean squared error is 1.278, and the adjusted R-square is 0.9924.

Figure 6 .
Figure 6.Combined reflectivity factor for 2.4 h from the CINRAD/SA, 2 May 2021.(a-i) Continuous observations from the CINRAD/SA at 12 min intervals from 11:12 to 13:36.The horizontal and vertical coordinates are, respectively, the ranges (km) in the west-east direction and the south-north direction.

Figure 7 .
Figure 7. Mean Doppler velocity for 2.4 h from the CINRAD/SA, 2 May 2021.(a-i) Continuous observations from the CINRAD/SA at 12 min intervals from 11:12 to 13:36.The horizontal and vertical coordinates are, respectively, the ranges (km) in the west-east direction and the south-north direction.

Figure 7 .
Figure 7. Mean Doppler velocity for 2.4 h from the CINRAD/SA, 2 May 2021.(a-i) Continuous observations from the CINRAD/SA at 12 min intervals from 11:12 to 13:36.The horizontal and vertical coordinates are, respectively, the ranges (km) in the west-east direction and the south-north direction.

Figure 8 .
Figure 8. Range-height cross-section of  following the azimuthal direction (azimuth: 225°).(a-f) Continuous observations from the CINRAD/SA at 6 min intervals from 11:48 to 12:28, 2 May 2021.The horizontal axis is the range (km) along the wind direction, and the vertical axis is the height (km).

Figure 8 .
Figure 8. Range-height cross-section of Z following the azimuthal direction (azimuth: 225 • ).(a-f) Continuous observations from the CINRAD/SA at 6 min intervals from 11:48 to 12:28, 2 May 2021.The horizontal axis is the range (km) along the wind direction, and the vertical axis is the height (km).

Figure 9 .
Figure 9. Time-height cross-section of radar products for 24 h from the CINRAD/SA at Daxing, Beijing, 2 May 2021.The fill color is  (unit: dBZ); the black isopleths are valid data proportions; and the white isopleths are  (unit: dB).The times of sunset and sunrise were 11:08 and 21:13, UTC, respectively.

Figure 10 .
Figure 10.Vertical wind profiles measured using the Windcube 100 s on 2 May 2021.The positive and negative values represent the vertical upward and downward wind speeds (unit: m/s), respectively.

Figure 9 . 2 Figure 9 .
Figure 9. Time-height cross-section of radar products for 24 h from the CINRAD/SA at Daxing, Beijing, 2 May 2021.The fill color is Z (unit: dBZ); the black isopleths are valid data proportions; and the white isopleths are Z DR (unit: dB).The times of sunset and sunrise were 11:08 and 21:13, UTC, respectively.

Figure 10 .
Figure 10.Vertical wind profiles measured using the Windcube 100 s on 2 May 2021.The positiv and negative values represent the vertical upward and downward wind speeds (unit: m/s), respec tively.

Figure 10 .
Figure 10.Vertical wind profiles measured using the Windcube 100 s on 2 May 2021.The positive and negative values represent the vertical upward and downward wind speeds (unit: m/s), respectively.
Remote Sens. 2023, 15, x FOR PEER REVIEW 16 of 25 4 m/s, and 65.3% of the values are less than 2 m/s.The echoes were unlikely caused by large insects and birds mainly because large insects and birds have quite a great flight speed and would make the speed field messy and produce great deviations.

Figure 11 .
Figure 11.Doppler velocity fields (a-c) and the fields which were obtained from the VAD analysis(d-f).The elevation angles of (a-c) are 0.5°, 1.5°, and 2.4°, respectively, the same as (d-f).Except for some error points and point targets, the velocity fields (d-f) which were obtained from the VAD analysis are similar to the Doppler velocity (a-c).

Figure 12 .
Figure 12.(a) Histogram of the residual error of the VAD analysis at different altitudes.The volume of the radar beam is used to calculate the ordinate.(b) The profiles of the mean wind speed during 12:30 to 13:30 (UTC) on 2 May 2021.

Figure 11 .
Figure 11.Doppler velocity fields (a-c) and the fields which were obtained from the VAD analysis (d-f).The elevation angles of (a-c) are 0.5 • , 1.5 • , and 2.4 • , respectively, the same as (d-f).Except for some error points and point targets, the velocity fields (d-f) which were obtained from the VAD analysis are similar to the Doppler velocity (a-c).

Figure 11 .
Figure 11.Doppler velocity fields (a-c) and the fields which were obtained from the VAD analysis(d-f).The elevation angles of (a-c) are 0.5°, 1.5°, and 2.4°, respectively, the same as (d-f).Except for some error points and point targets, the velocity fields (d-f) which were obtained from the VAD analysis are similar to the Doppler velocity (a-c).

Figure 12 .
Figure 12.(a) Histogram of the residual error of the VAD analysis at different altitudes.The volume of the radar beam is used to calculate the ordinate.(b) The profiles of the mean wind speed during 12:30 to 13:30 (UTC) on 2 May 2021.The comparison of the speed profiles with the RWP and the CINRAD/SA confirms the conclusion.The profiles of the horizontal wind speed which were sensed by the RWP are not different from the speed profiles calculated by the CINRAD/SA, as shown in Figure 12b.Despite the system deviation in Figure12b, which is common in the results of the VAD wind profile[47], the changes in the speed profiles are consistent, generally.The profiles of the CINRAD/SA indicate that clear-air echoes are moving with the wind speed.However, the speed of biological autonomous movement is not detected by the weather radar.This means that these echoes are unlikely to be caused by creatures with a strong flight ability such as birds and large insects.

Figure 12 .
Figure 12.(a) Histogram of the residual error of the VAD analysis at different altitudes.The volume of the radar beam is used to calculate the ordinate.(b) The profiles of the mean wind speed during 12:30 to 13:30 (UTC) on 2 May 2021.

Figure 13 .
Figure 13.Time-height cross-section of  (the fill color) and the valid data proportions (the black isopleths) for 24 h from the X-band radar at FS, Beijing, 2 May 2021.

Figure 14 .
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.

Figure 13 .
Figure 13.Time-height cross-section of Z (the fill color) and the valid data proportions (the black isopleths) for 24 h from the X-band radar at FS, Beijing, 2 May 2021.

Figure 13 .
Figure 13.Time-height cross-section of  (the fill color) and the valid data proportions (the black isopleths) for 24 h from the X-band radar at FS, Beijing, 2 May 2021.

Figure 14 .
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.

Figure 14 .
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.
Figure 14.The DWR between the values of Figures 8 and 12.The time-height cross-section and the histogram of the DWR at nighttime are exhibited in (a) and (b), respectively.The black isopleth in (a) is the proportion of the X-band valid data.The bar chart (b) represents the normalized frequency distribution of the DWR from sunset to sunrise, and the data in the histogram from 19 dB to 25 dB occupy 80%.

Figure 15 .
Figure 15.Histograms of the radar products of the CINRAD/SA (S-band) and X-POL (X-band), respectively, on 2 May 2021 at 13:00 UTC for a height ranging between 300 m and 1.2 km.

Figure 15 .
Figure 15.Histograms of the radar products of the CINRAD/SA (S-band) and X-POL (X-band), respectively, on 2 May 2021 at 13:00 UTC for a height ranging between 300 m and 1.2 km.

, 15 ,
x FOR PEER REVIEW 20 of 25 microwave radiometer.Although the geometrical structure of the reflecting layers remains obscure, it can be inferred that turbulent mixing is detrimental to the reflecting layers because a reflecting layer in the atmosphere is formed by relatively sharp gradients of the refractive index, but turbulent mixing makes the temperature and humidity homogeneous and reduces the gradient of the refractive index.Thus, the Richardson number and the turbulence intensity, which indicate turbulent mixing, are shown in Figure16to give the relationship between clear-air echoes and turbulent mixing.

Figure 16 .
Figure 16.Turbulence intensities (a) and Richardson number (b) for 20 h at Daxing, Beijing, 2 May 2021.The black isopleth is the time-height cross-section of  (unit: dBZ), which is the same as the fill color in Figure 6.

Figure 16 .
Figure 16.Turbulence intensities (a) and Richardson number (b) for 20 h at Daxing, Beijing, 2 May 2021.The black isopleth is the time-height cross-section of Z (unit: dBZ), which is the same as the fill color in Figure 6.

Figure 17 .
Figure 17.Time-height cross-section of the Z values (a) from 22:30 to 24:00, and the vertical profiles of the radar products (b,c) and rawinsonde products (d,e) at 23:15 UTC on 2 May 2021.The profiles of (d,e) are measured by rawinsondes.The estimated entrainment layer is based on the maximum vertical gradient in each variable (orange line).

Figure 17 .
Figure 17.Time-height cross-section of the Z values (a) from 22:30 to 24:00, and the vertical profiles of the radar products (b,c) and rawinsonde products (d,e) at 23:15 UTC on 2 May 2021.The profiles of (d,e) are measured by rawinsondes.The estimated entrainment layer is based on the maximum vertical gradient in each variable (orange line).

Table 1 .
Parameters of several insects and their RCSs at different wavelengths.The biometric data were provided by the Chinese Academy of Agricultural Sciences from captured insects in North China.The RCSs were simulated by FEKO simulation software using the prolate spheroid model of a spinal cord dielectric.

Table 2 .
System characteristics of the CINRAD/SA radar and the X-POL radars.

Table 2 .
System characteristics of the CINRAD/SA radar and the X-POL radars.

Table 1 .
All species are major agricultural pests in North China.