1. Introduction
Unwanted echoes represent a significant problem in many radar applications, such as the surveillance of airports, harbors, and inhabited areas. These echoes, globally referred to as “clutter” [
1], can be generated by both continuous (i.e., distributed) or discrete sets of scatterers. Discrete sources of clutter are, for example, road/rail infrastructures, buildings, water towers, and large lampposts. On the other hand, grass, vegetation, terrain, etc. represent distributed clutter sources, whose echoes are measured as the radar cross section per unit area. Clutter may be orders of magnitude more intense than the useful targets of interest, strongly affecting their visibility. Moreover, the dynamic range of the clutter echoes and their probability distribution widely vary both in time (at a given resolution cell) and in space (i.e., from adjacent radar resolution cells at the same time). These effects strongly depend on the radar installation in the operational area.
Ground clutter models are widely studied, see for instance [
1,
2,
3], but, in spite of the available literature and databases, dedicated measurements are needed in order to obtain the clutter levels in a given environment, especially when both discrete and distributed sources are present.
In order to generate valuable clutter maps, it is necessary to operate with a well calibrated radar sensor, that is, setting of a precise relationship between the received power and the radar cross section (RCS) in each resolution cell. The classical way to calibrate a radar sensor by measuring the various terms that made up the radar equation [
4] (i.e., antenna gain, transmitted power, losses, etc.) may be lengthy and cumbersome. In the literature [
5,
6], a more cost/effective approach is investigated based on reconstructing the relationship between the RCS of known, standard objects (metallic spheres and corner reflectors) and the received power. This relationship is a linear law above the noise plafond and below the receiver saturation zone.
After calibration, realistic clutter maps can be generated [
7,
8,
9,
10]. They may be used, for example, to design the siting of X-band surveillance radars [
11] for the surveillance of airports (surface movements radar and apron control) and for the detection of unauthorized unmanned aerial vehicles (UAVs also known as drones) in urban and suburban areas. The possibility for the UAVs to carry on payloads makes their use both very interesting and potentially dangerous. The Federal Aviation Administration (FAA) regulations in the USA regarding the use of drones in commercial applications are clearly defined, while hobby or recreation usages are encouraged to follow the safety guidelines [
12]. The problem of drones’ control over a limited area is strongly rising up. The first idea is to, in some way, reproduce the surveillance system in Air Traffic Control (ATC) which uses dedicated radar sensors. However, they are very expensive and not easily deployable: typically, operating at S-band (2–3 GHz), these radars are more suitable for detecting big objects, like an airplane, rather than small drones. On the other hand, drones might be detected by cheaper commercial X-band (9.2–9.5 GHz) marine radars, either with a “solid state” transmitter or with a magnetron transmitter, which are very portable and easily deployable as compared with the large ATC surveillance radars.
The RCS of drones has recently been investigated. In [
13], the numerical and experimental results are shown for a quadcopter drone (DJI Phantom2) at 10 GHz. The measurements of drones, carried out in an anechoic chamber, have shown RCS variations from −30 dBm
2 to −10 dBm
2. In [
14], experimental results, at a range of up to 2 km, with an average altitude (above ground level) of 30 m, are reported using a commercial drone (DJI-Phantom4) and a ubiquitous FMCW radar system working at 8.75 GHz. The authors concluded that the drone echoes show the statistical behavior of Swerling I targets with an average RCS close to −20 dBm
2. In [
15], the radar signals in both the range and Doppler domains are analyzed for a quadcopter drone (DJI Phantom2), as well as three different birds (Hooded Vulture, Eurasian Eagle Owl, and Barn Owl). The authors, using an S-band coherent multistatic radar, found that the micro-drone signatures were comparable (in RCS) to those of birds. In [
16], X-band measurements of flying sea birds have been carried out in order to investigate the potential of extracting the features for target recognition from micro-Doppler signatures.
In spring 2017, at Tor Vergata University in Rome, the preliminary experimental results were obtained using a commercial X-band marine radar to detect drones with the aim to compare the drone’s echo to those coming from known reflectors, in order to derive the radar visibility of the drone [
17]. However, in this experimental research, it was soon clear that back-scattering and multipath, because of both continuous and discrete sources (often found in a suburban area), constitute a serious problem for small aircraft detection. The preliminary results, regarding the multipath analysis as a result of the land surface using an X-band radar, are shown in [
18]. The relevant evaluations are based on the measured RCS variations of a calibrated target taken during the NATO SET-225 (Sensors and Electronic Technology Group-225) field trials in June 2018 at the Fraunhofer Institute FHR in Wachtberg, Germany.
During 2018 and in spring/summer 2019, more trials were carried out at Tor Vergata University. Some of the acquired data will be described in the following so as to highlight the problems related to the back-scattering and to the multipath in the frame of the detection of small targets.
This work is an extended version of the paper “Environmental Effects on Ground-based Radar Measurements” in Proceeding of the Sixth IEEE Workshop Metro Aerospace, Torino, Italy, June 2019 [
19]. The paper is organized as follows.
Section 2 recalls the system calibration procedure, including the estimation of the total radar losses from the antenna to the receiver.
Section 3, considering the effects of terrain, describes the reflection coefficient and the multipath with its possible compensation.
Section 4 and
Section 5 show the evaluation of some aspects of multipath in the Tor Vergata area and in the Wachtberg area using calibrated targets. Using opportunity targets, in
Section 6, the non-stationary nature (in time and in space) of the multipath is highlighted in a suburban environment. Final considerations and perspectives are reported in
Section 7.
Appendix A recalls the multipath.
2. Radar System Calibration and Estimation Losses
For research purposes, the Radar Business Unit (RBU) of the Navico Company has supplied Tor Vergata University of Rome, on loan to use, an International Maritime Organization (IMO) compliant marine radar (model Simrad SRT). The main parameters of this set are shown in
Table 1.
Figure 1 shows the antenna and the Tx/Rx unit on the top of the Information Engineering building at Tor Vergata University in Rome. The radar is placed in the west corner of the upper terrace and is installed on a lifter capable of increasing the height of the radar antenna by up to 1 m. Arbitrary step positions are made possible for data acquisition, so as to evaluate the multipath effects, as explained in the following.
All of the trials shown in this paper have been performed in the short range mode (see
Table 1) with a pulse length of 50 ns (bandwidth B = 20 MHz) and PRF = 3000 Hz. In this set-up, the power of the receiver noise
(
is the noise figure,
is the Boltzmann’s constant, and
is the reference temperature) is equal to
. The range resolution is 7.5 m, with a sampling interval of 3 m. Because of the horizontal beam-width of 1.35°, the nominal azimuth resolution is 12 m and 24 m at distances of 500 and 1000 m, respectively.
Downstream of the analog to digital conversion, for each antenna revolution, the acquired digital signals, coded in 256 levels, hereafter called (0–255), are recorded in a (Range, Azimuth) matrix with 2048 rows (range-bins), corresponding to a maximum range of 6144 m (), and a number of columns (azimuth-bins) nominally equal to 7500 at the rotating speed of 24 RPM and a PRF equal to 3000 Hz. The (Range, Azimuth) matrix allows us to display the plan position indicator (radar PPI).
The received power, expressed in dB above one milliwatt (dBm) is estimated by a linear relationship, as follows:
where the constant
in (dB/level) is associated to each measured ADC level, while the constant
relates to the noise floor level (
). These constants are obtained by radar calibration, without the need for dedicated measurements on the radar set.
The radar calibration is fully described in the literature [
5]. The main steps are summarized as follows:
Three octahedral corner reflectors (CRs) having scaled dimensions of the side (large , medium , and small ), and a coated sphere with a radius of , are put in a suitable position (characterized with terrain masking by suited natural fences) with both land clutter and multipath effects practically absent.
Four independent measurements of the , obtained by averaging eight measurements for each target, are recorded.
The theoretical (in ) of CR’s are normalized to the RCS of the sphere.
From step (ii), on the plane (), four points are plotted.
The slope of the regression line among these points estimates the parameter α, while the parameter is estimated imposing on the straight line to pass on the mean noise point.
The measurements did confirm that in the chosen configuration, the multipath factor [
4] is very close to the unit, thus allowing us to estimate, with an accuracy of the order of 1%, the parameters
and
of Equation (1). It resulted
(dBm/level) and
(dBm). More details on this procedure are presented in [
5].
The received power,
, is theoretically evaluated by the Radar Equation (2), including the overall system losses
up to the analog-to-digital-converter (ADC), as follows:
For a given distance,
, of the target, in the linear range of the receiver response,
is proportional to the
of the target (
) and to the ratio
, where
denotes the losses and
is the contribution due to multipath, analyzed in
Appendix A. Inverting Equation (2) and using Equation (1) for the
expressed in dB, the
,
, of a single target is written as follows:
where
includes the various radar parameters and constants. Using
(dBm/level) and
(dBm; under the hypothesis that in operational conditions, such as the protection by a fence,
), the losses
in Equation (3) can be estimated as the value that minimizes the mean square error between the theoretical and the measured radar cross sections of the dedicated targets (four in this study). This resulted in
. This value is compatible with the manufacturer data, kindly supplied by Navico RBU.
In general, it is not easy to evaluate, in one shot, radar losses with a residual error below 1 dB in typical conditions, and, to the best of the authors’ knowledge, these overall black-box measurements of the radar receiver and overall losses, using standard reflectors and a fence, are novel.
3. Evaluation of Multipath and of Reflection Coefficient
The previous discussion is related to those very particular conditions (e.g., rays blocked by a fence) in which the multipath on the vertical plane is negligible. Conversely, in the presence of a multipath [
1], that is a normal situation in surface radars, the calibration by standard reflectors requires evaluating the multipath factor
in the particular geometry and propagations conditions. The multipath factor is the ratio between the amplitude of the electro-magnetic field on the target and the amplitude in the (theoretical) case of the absence of multipath. This factor depends on the complex reflection coefficient of the surface
, and on the phase shift
, related to the path difference
between the direct and the reflected signals.
Following the developments and the notations in
Appendix A,
is expressed as follows:
For a given polarization, the parameters and depend on the electromagnetic properties of the surface, that is, the real and imaginary part of its dielectric constant , the frequency, the polarization, and on the grazing angle . For the horizontal polarization, .
In addition, considering the roughness of the surface, the global reflection coefficient is , where depend on the grazing angle, , and on the r.m.s. roughness in the vertical direction .
For flat soil (clay) with a moisture content of
, we may use (for the horizontal polarization of the SRT) values of
and
(valid at
[
20]) as a good approximation of the
data. Including the contribution of the land roughness (
),
Figure 2 shows the amplitude,
, of the reflection coefficient for
. For low grazing angles (in the order of
or
degrees), this coefficient is within
.
In the horizontal polarization, the phase of the reflection coefficient for low grazing angles can be supposed to be equal to , as follows: for grazing angles between and , its values range from to (H-polarization, , soil moisture , , and , that is, , is not affected by the variation. The roughness contribution to the reflection coefficient is real.
In
Appendix A,
Table A1 reports the maximum variation in the “apparent”
of the target due to multipath. For example, considering
(corresponding to
grazing angle and
m), a radar target distance
and
target height,
Figure 3 shows a sample of the variation of the multipath factor with the radar antenna height,
, varying from
to
. The distance and the height of the target correspond to a certain target (lamppost L4, see
Section 6) and the antenna height to the installation on the roof of the Information Engineering building at Tor Vergata University. The latter height can be controlled by varying the vertical radar position of the antenna by a mechanical elevator, as shown in
Figure 1. With this set-up, we were able to reconstruct curves, such as in
Figure 3, in order to estimate (with a best fit procedure) the modulus
of the reflection coefficient at the given experimental conditions.
From curves such as the one in
Figure 3, the maximum value
and the minimum value
of the multipath factor are easily estimated. As
is equal to
, and
is equal to
, the modulus of the reflection coefficient is readily estimated, as follows:
In practice, the effect of the multipath in the RCS measurements can be corrected, as described in the following procedure.
- (a)
Set the target under test (distance D from the radar and the height of its main scatter ).
- (b)
Evaluate the period of the factor
versus the antenna height
. It is
(
Appendix A).
- (c)
Define the step to obtain a significant number (e.g., ten) of positions inside the period.
- (d)
Record the received power, proportional to , for each position, and draw the best-fit curve.
- (e)
Find the maximum and the minimum values of .
- (f)
Use Equation (5) to obtain .
- (g)
Use Equation (A7) to calibrate the received power.
4. Multipath Evaluation in the Tor Vergata Area
The first experiments were carried out in 2018 with a fixed radar installation for various heights of the same target (a corner reflector on the top of a pole, see left part of
Figure 4).
Subsequently, in 2019, the radar elevator, shown in
Figure 1, allowed us to vary the antenna height (in addition to the variation of the corner reflector height). The aim was to quantify the effect of the multipath on the
measurements, and to evaluate the reflection coefficient of the terrain (land with vegetation). For this goal, it was decided to keep the reflection point fixed on the ground, in order to be free from unwanted variations of the reflecting surface.
On 16 May 2019, a trihedral CR with an
of
has been posed (as a target) on a tripod at a distance of
from the radar (see
Figure 4). The weather conditions were cloudy sky with sunshine,
, and no wind.
Considering the starting position of the measurements (antenna height
and the target height
), the reflection point on the ground is
away from the radar. To keep the distance of the reflection point unchanged, we varied the height of the antenna and the target, and being
, it was necessary to vary the
and
(Equation (A1)) according to the values shown in
Table 2. Considering the maximum ranges of elevation for the target and for the antenna, five positions were deemed sufficient.
Table 2 reports, for each position, the time of the acquisition and the measured
of the corner reflector. Remember that the theoretical value of the
of this CR is
.
For the five positions, the mean grazing angle is
. Using Equation (1) with
and
, the received power has been evaluated. Hence, the multipath factor has been estimated using the average of the five received power values for each position. Using Equation (5), the amplitude of the reflection coefficient results in
. This value is compatible with vegetated land, with the height of the vegetation in the order of 10 to 15 cm, which is consistent with the observed height of the grass on the meadow.
Figure 5 shows the theoretical multipath factor and the measured one.
The differences between the theoretical model and the measurements of this early experiment have been attributed to various unfavorable elements of the experimental area, including the following: (i) the land not being flat, (ii) the vegetation not being uniform, and (iii) the variation of the grazing angle versus .
In a more regular environment, such as the one described in
Section 5, the multipath analysis has shown a better agreement with the theory [
19].
The ensuing experiments were performed, and their results are presented in the following. One experiment, conducted at Fraunhofer FHR, investigates the influence of variations in the vegetation on the multipath effect, as presented in
Section 5. A second set of experiments, performed over a period of more than one month, and presented in
Section 6, was conducted at Tor Vergata, aiming for an investigation of the temporal influences on the multipath effect.
5. Multipath Effects of Different Vegetations Observed in the Area of Wachtberg
Trials were conducted in June 2018, with a continuous emission (CE) noise radar at the Fraunhofer FHR Institute in Wachtberg, Germany. They addressed different aspects of noise radar, and were carried out by members of the NATO research task group (RTG) SET-225. The noise radar (NR) used during these experiments was developed by FHR, and is referred to as FHR-NR; it operates in the X-Band, more precisely in the standard maritime navigation band ranging from to , and thus, the same frequency band as the Simrad SRT radar, which allows for a direct comparison of the multipath effects observed in Tor Vergata with the ones at the FHR outpost in Wachtberg.
During the multipath trials, the operational signal bandwidth was , corresponding to a range resolution not less than . The transmitted signal was a random waveform with a Gaussian distribution and less than of effective radiated power.
Dedicated experiments were conducted to characterize the effect of different kinds of vegetation on the multipath.
Figure 6 shows the setup of the experiments with two positions of the same reference radar target. FHR-NR was installed inside a shelter at a small facility in Wachtberg, surrounded by fields. On one side of the shelter, there was a field of sweetcorn, with an approximate height of
. On the other side, there was grassland of approximately
in height. The first position of the standard target was located in the sweetcorn field on a block of concrete at a
distance to the radar (path A). The second target position had a radar distance of
, and was located close to the end of the grass field (path B).
The standard target, used for both experiments, was a large trihedral corner reflector (CR) with a radar cross section not less than
. This corner reflector was mounted on a tripod, fully extended during all of the experiments at a constant height of
. The height of the radar transmitting and receiving antenna system, made up by two standard horn antennas on a second tripod, as shown in
Figure 7, varied in intervals of
from
(path A), respectively
(path B), up to a maximum height of
.
A detailed analysis of this experiment can be found in [
19]. The main points are summarized here and compared with the previously discussed results from the experiments performed at the Tor Vergata University of Rome.
Figure 8 shows the measured multipath effect on the RCS measurement of the standard target.
In both graphs of
Figure 8 the periodic multipath effect is clearly seen. The measured values (red dots) are in good agreement with the theoretical values drawn as a blue line. The influence on multipath of the different kinds of vegetation can be easily read from the graphs. The amplitude of the oscillation is higher in the case of propagation over sweetcorn while an oscillation with lower amplitude but higher frequency is observed over grassland.
These results represent a single snapshot of a single day without any significant variation in temperature, humidity and more; anyway, this single comparison between two different kinds of vegetation results in significant observations. Long-term measurements are presented in the ensuing section describing the variations of the multipath effect over time.
6. Long-Term Measurements of Radar Cross Section of Man-Made Targets
In the frame of the 2018 experiments at Tor Vergata campus and their 2019 follow-on, fixed identical and very large lampposts (
high, here referred to as L1, L2, L3, and L4) at different distances from the radar (
,
,
, and
, respectively) were considered as opportunity targets.
Figure 9 displays an overview of the area, showing the position of the radar and of the lampposts.
Each lamppost has six pairs of large lamps arranged on a metallic ring, as shown in the inset of
Figure 9. As the top of each lamppost is similar to a single, small, but strong radar scatterer at
above ground level, the radar echo from these targets are candidates to be significantly affected by multipaths. Using Equation (3) and the (Range, Azimuth) matrix of the ADC levels, an
map of the environment has been generated, as shown in
Figure 10. In these measurement results, two significant target echoes are created by the compound of the Literature Faculty and the long commercial building in via Schiavonetti.
Unexpectedly, during the trials on 17 May 2018, the estimated values of the
, in two close times of this day (11:05 and 12:00) have shown a considerable difference for lampposts L2 and L4. For L2, the measured RCS changed from
to
, while for L4, it changed from
to
(see
Figure 10a,b, where the
in
is shown as a false color in an area including the four lampposts).
The rationale for these significant variations can only be attributed to the propagation effects, in particular, to some time-varying multipath. In fact, in some experiments, we found that the reflecting surface on the ground is characterized by many time-varying scatterers (such as vegetation, foliage of trees, tall grass, etc.), and we know that the multipath can significantly affect the measurements of the .
The sensitivity of the measured
w.r.t. the height of the target (see
Appendix A) increases significantly when this height (in wavelengths) increases. Hence, for a tall scatterer such as a building (e.g., in
Figure 10, the building at Via Schiavonetti, and the ones of the Faculty of Literature), the multipath effects average along the elevation, and, thus, the measured
is stable; the opposite is true for a point target, such as a small flying object or, neglecting the contribution by its pole, a lamppost at a fixed height. Therefore, for large and tall targets and for most buildings, the measured
has steady values as the multipath effect vanishes (
).
To better understand this point, more measurements of the four lampposts were carried out in spring 2019, varying the antenna height. Five 25 cm-step positions, from a maximum elevation (position 1) to a minimum elevation (position 5), were used for radar data acquisition.
Figure 11 shows the reflection points (see point M in
Figure A1) on the ground for lampposts L1, L2, and L4.
For lampposts L1 and L4, the reflection points shifted on the road, by 7−8
, due to the variation of the antenna height. Hence, the reflection points’ positions may move from vegetated land into a road, possibly with some traffic. For the lamppost L2, the reflection point is on the parking lot of the Faculty of Literature, and the multipath may change rapidly if cars or buses move or appear/disappear during the acquisition. This explains the noticed, strong variations in
. The pole of the lamppost L3 is partially masked by the roof of the “Didattica” building (see
Figure 12), hence, multipath effects are considered negligible for the
measurements of L3.
We collected the measurements obtained over the following four different days: 16 May 2019, 4 June 2019, 7 June 2019, and 26 June 2019.
Figure 13 shows the measured
(in dB square meter) of these four lampposts, with a varying antenna elevation (five steps, 0.25 m each).
Lamppost L3, which is the reflection-masked target practically unaffected by the multipath effects, shows the reflection behavior different from the remaining three lampposts. During the different days, the
of the lamppost L3 always decreases from a maximum (step 1) to a minimum (step 5), as shown in
Figure 14, due to an increased occlusion by the “Didattica” building. The temporal variation of the multipath effect is described by the parameter
, which has been evaluated as follows:
where
denotes the number of positions for each day, while
and
are the sample means.
The measurement of the RCS of lamppost L3 does vary a little from day to day, making the correlation coefficient
of the measured
on two different days
close to the unit (
Figure 14).
A different situation arises for the lampposts L1, L2, and L4, where the measured (at the same position on different days) show differences that can reach or more. These significant variations must be attributed to the propagation effects, that is, to a time-varying multipath.
For lamppost L1 and L2,
Figure 15 shows a comparison of the measured
, on different days, varying the antenna position.
For L1, the values show a negative correlation (), while for L2, they are quite de-correlated (). These results are strongly affected by the time variation of the reflector coefficient (moving vehicles on the road, day-to-day variation, etc.)
7. Comments and Conclusions
Field experiments have been performed in suburban and country environments to better understand radar propagation and scattering effects, with particular attention to multipath. To this aim, it has been necessary to calibrate the radar. In general, radar calibration is a cumbersome task, requiring standard reflectors (such as corner reflectors and metal-coated spheres) and the compensation of the backscatter from the underlying surface (i.e., land) and of the multipath in the vertical plane (assuming the horizontal plane free from multipath). In Tor Vergata trials we used a non-coherent radar (with Doppler measurements not available) operating in the marine X-band with horizontal polarization. Some preliminary indications from the experimental results are synthesized as follows.
First, the unwanted reflection and scattering may be mitigated using clutter fences, with the calibration reflectors positioned in order to be seen from the radar “just above the fence”.
Second, the multipath (when not avoided by a fence) may strongly affect the calibration measurements, even when due to an apparently a more “benign” surface (slightly wet land, limited vegetation) that the well-known, strongly reflecting sea surface. The multipath compensation calls for estimating the reflection coefficient by a set of measurements (order of five to ten). In a not uniform environment (e.g., with different vegetation, roads, trees, etc.), it is necessary to vary both the height of the reflector and of the radar antenna, in order to keep the reflection point fixed on ground.
Furthermore, different experiments have proven that the multipath effect is independent of the kind of radar emission, that is, pulsed or continuous emission. Instead, it is the range resolution of the signals that is important. Future investigations might determine the influence of the radar bandwidth on the RCS variations created by the multipath effect.
As a side effect of the calibration procedure after multipath compensation, a significant “bonus”, which is obtained, is the “black box” evaluation of the radar equipment total losses.