# Ionospheric Clutter Model for HF Sky-Wave Path Propagation with an FMCW Source

Authors: Xuguang Yang; Aijun Liu; Changjun Yu; Linwei Wang
Journal: International Journal of Antennas and Propagation
 (2019)

---

## Abstract

A theoretical model of the sky-wave path propagation with frequency modulated continuous wave (FMCW) source for high frequency (HF) radar is proposed in this paper. Based on the modeling of pulsed source, the expression of the received electric field with an FMCW source is derived for the reflection case from the ionospheric irregularities. Subsequently, the ionospheric reflection coefficient with different phase power spectrums for vertical and oblique backscattering propagation paths is incorporated into the ionospheric clutter model. Simulation results show that the peak power of FMCW in average is lower than that of pulsed waveform. Furthermore, different incident angles and magnetic field in mid-latitude can also influence the power density of the backscattering ionospheric clutter. Finally, the data analysis results from the high frequency surface wave radar (HFSWR) and Ionosonde collected in Yellow Sea preliminarily verify the inversion of the variance of the electron density fluctuation and the vertical drift velocity of the irregularities within ionosphere.


---

## Body

## 1. Introduction

In high frequency surface wave radar (HFSWR) system, the electromagnetic wave diffracts along the curvature of the earth to detect ships and low flying aircraft targets over the horizon. Ideally, a perfect conductive plane consisting of sea surface is infinite in coverage area and thus makes transmitted wave entirely traveling along the sea surface. But if considered the affection of poor ground, the overhead nulling of the receiving pattern will shift to left, which causes the whole vertical radiation pattern gain [1]. Simulations by High Frequency Structure Simulator (HFSS) show that transmitting antenna of HFSWR radiates partial energy to elevation especially to the overhead regardless of frequency [2]. The pattern of vertically polarized antenna receiving arrays also has -8dB gain on overhead. Thus partial wave radiated upwards and reflected from the ionosphere is unavoidable. The echo arrived at the receiving antennas from various paths, which is referred to as ionospheric clutter [3, 4].We summarize the path of ionospheric clutter in two kinds: the sky-wave path (include vertical and oblique propagation) and mixed path (via the ocean surface path), as shown in Figure1. Simulations show that the most effective paths are the vertical path and mixed path [1]. However, the oblique sky-wave path should also exist in HFSWR system, since the attenuation of ocean path is much higher than sky path and the detection range of HFSWR can reach 400-600km. Therefore, if the ionosphere-ocean mixed path exists, the oblique propagation path (0.5 jumps, 1 jumps, etc.) should also exist. The backscattering propagation range (0.5 jumps) is far away from the vertical propagation range, but shorter than mixed propagation range.Figure 1
Propagation paths of HFSWR ionospheric clutter.Figure2 shows the comparison between ionogram and RD (range Doppler) spectrum of HFSWR on 10:53 am, June 28, 2018, location: Weihai (37.5°N, 122°E), China. In the ionogram, the green points are locus of received extraordinary frequencies and their corresponding virtual height. The red points are locus of ordinary frequencies and their corresponding virtual height. We can see that the E layer is in 100-150km and the critical frequency is less than 5MHz. But the radar operating frequency is 6MHz and the ionospheric range is about 150-200km. Therefore, it should be the oblique ionospheric clutter in HFSWR.Figure 2
The comparison of ionogram and RD spectrum on 10:53 am, June 28, 2018, location: Weihai (37.5°N, 122°E), China. (a) Ionogram (the horizontal axis is the plasma frequency and the vertical axis is the virtual height). (b) RD spectrum.
(a)
(b)The ionospheric clutter often spread in range, Doppler, and spatial domains, leading to the masking of targets. Most of the clutter is caused by plasma density irregularities in the ionosphere under disturbed ionospheric conditions [5]. Thus it is necessary to develop a mathematic model to quantify the impact of the ionosphere. Many investigations have been made on the characterization of ionospheric clutter and suppression of ionospheric clutter [6–11]. These ionospheric clutter mitigating techniques include many complicated adaptive signal processing algorithms or additional antennas. However, these techniques benefit from a physical model based on scattering mechanism, which demonstrates how the ionospheric clutter manufactured between electromagnetic wave and irregularities caused by plasma instabilities.Walsh et al. [12] laid a solid theoretical foundation to describe the scattering of electromagnetic waves from the ocean surface, which has been extended to model the power spectral density of mixed path propagation for pulsed radar [13] and a frequency modulated continuous wave (FMCW) radar [14]. In those models, the ionospheric reflection coefficient (IRC) is assumed as Gaussian function, which may not accord with the ionosphere physical conditions. In [15], the model of the ionospheric clutter from vertical incidence for HFSWR and Bragg backscatter in the auroral zone for over-the-horizon radar (OTHR) were developed using a geometric optic approach. Chen incorporates the particular ionospheric clutter model into Walsh’s model for mixed path with an FMCW [16, 17] and vertical reflection with pulsed radar [18]. However, the case of vertical propagation for FMCW and oblique sky-wave propagation has not been considered. Since the main waveform in HFSWR system is FMCW or phase coding waveform, the vertical theoretical model with pulsed cannot be validated. In this paper, we investigate the vertical reflection with FMCW source, compared the simulation results with pulsed source and test the experimental results with Ionosonde. Besides, we derive the generalized power spectrum for backscatter sky-wave propagation (0.5 jumps) in mid-latitude.The remaining part of the paper is organized as follows. In Section2, the electric field equation and vertical power spectrum of FMCW source are obtained under the assumption of a continuously excited elementary vertical dipole source. The backscattering power spectrum incorporated with a phase power spectrums of ionospheric irregularities in mid-latitude is derived in Section 3. In Section 4, simulations are conducted for FMCW and compared the experimental results with Ionosonde. Conclusions are summarized in Section 5.


## 2. Vertical Received Power Spectrum for FMCW Source

Suppose that a dipole transmitting antenna is at coordinate(
0,0
,
0
), the ionospheric reflecting layer is at a height of H
/
2, and the dipole image is at height ofH, as shown in Figure 3. The angle of incidence is θ
i. In the case of vertical reflection, we have θ
i
=
δ
θ, which is the angle error of the antenna pattern. Thus the vertical incidence is a special case of oblique backscattering.Figure 3
Scattering geometry for ionospheric backscattering.Based on image theory [13, 18], the scattered electric field equations for the backscattering ionospheric clutter at (
0,0
,
0
) can be written (when P
(
x
,
y
,
0
) coincides with (
0,0
,
0
)):(1)
E
R
=
j
C
0
R
i
sin
⁡
θ
i
e
-
j
2
k
R
8
π
Rwhere C
0
=
η
0
Δ
l
ω
I
(
ω
)
/
c is the dipole coefficient for antenna of length Δ
l carrying a current I whose radian frequency is ω and wavenumber is k in free space where the intrinsic impedance is η
0, R is the range between the ionospheric irregularities and radar, and R
i
(
x
,
y
) is the IRC at horizontal position (
x
,
y
) of the ionospheric scattering point, which may be considered as a random process that can be written in Fourier form:(2)
R
i
x
,
y
=
∑
K
→
i
P
K
→
i
e
j
K
→
i
·
ρ
→
iwhere P
K
→
i is the Fourier coefficient of the ionospheric irregularities, ρ
i is the diameter of the scattering area varied from 0 to Δ
ρ
i, and K
→
i is the wavenumber of the ionospheric irregularities. The field in (1) can be transformed to the time domain by applying inverse Fourier transform (see (
7
) in [18]):(3)
E
R
t
=
sin
⁡
θ
i
4
π
2
R
Δ
ρ
i
2
F
t
-
1
j
C
0
∗
t
δ
t
-
2
R
C
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
iwhere J
0
(
K
i
ρ
i
) is zero-order Bessel function, F
t
-
1 denotes inverse Fourier transform, ∗
t denotes time convolution, and δ
(
·
) is the Dirac delta function.


### 2.1. The Electric Field for FMCW Waveform

The time domain expression for the FMCW transmitting antenna current waveform within one sweep interval is given as(4)
i
t
=
I
0
e
j
ω
0
t
+
α
π
t
2
h
t
+
T
r
2
-
h
t
-
T
r
2whereI
0 is the peak current, ω
0 is the radian frequency, α
=
B
/
T
r is the frequency sweep rate, B is the sweep bandwidth, T
r is the sweep interval, and h
(
t
) is the Heaviside function. Unlike pulsed radar, the range method for FMCW radar is to use the differential frequency between the transmitting and receiving signal:(5)
R
=
c
Δ
f
2
α
=
c
f
t
-
f
r
2
αwhereR is the target range,c is the light speed, and f
t and f
r are the frequencies of transmitting and receiving signal, respectively.For dipole antenna transmitting FMCW signal, we have the approximate expression in the derivative:(6)
F
t
-
1
j
C
0
=
η
0
Δ
l
c
d
i
t
d
t
≈
j
k
0
η
0
Δ
l
I
0
e
j
ω
0
t
+
α
π
t
2
h
t
+
T
r
2
-
h
t
-
T
r
2where k
0
=
ω
0
/
c. Then (3) becomes(7)
E
R
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
j
ω
0
t
-
2
R
/
c
e
j
α
π
t
-
2
R
/
c
2
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
cAfter the demodulation processing, (7) reduces to(8)
E
R
D
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
-
j
4
π
α
R
2
/
c
2
-
2
k
0
R
e
j
4
π
α
R
/
c
t
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
cwhere the superscriptD represents demodulation.For FMCW, the target range is determined by the difference frequency of the demodulated electric field. Thus the range transform involves Fourier transform the electric field with respect tot for obtaining the difference frequency. Only the last two terms are function of t and their Fourier transforms are(9)
F
e
j
4
π
α
R
/
c
t
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
c
=
∫
2
R
/
c
-
T
r
/
2
2
R
/
c
+
T
r
/
2
e
j
4
π
α
R
/
c
t
e
-
j
ω
t
d
t
=
T
r
e
j
4
π
α
R
/
c
-
ω
2
R
/
c
S
a
T
r
2
ω
-
4
π
α
R
cInserting (9) into (8), the frequency domain expression for FMCW can be obtained:(10)
E
R
D
ω
r
=
F
E
R
D
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
T
r
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
-
j
2
k
0
+
ω
r
/
c
R
e
j
4
π
α
R
2
/
c
2
S
a
T
r
2
ω
r
-
4
π
α
R
cHere,ω has been renamed as ω
r and S
a
(
x
) is the sampling function of sin
⁡
(
x
)
/
x. For a given frequency ω
r, there is a corresponding delay time Δ
t which is the total delay between the ionospheric irregularities and radar. Then the apparent ionosphere reflection range R
0 and range resolution Δ
R can be defined as(11)
R
0
=
c
Δ
t
2
=
c
ω
r
4
π
α
,
Δ
R
=
c
2
BAnd we define(12)
k
B
=
2
π
B
c
,
k
r
=
ω
r
c
,
R
′
=
R
0
-
RThus the variable of the sampling function in (10) becomes(13)
T
r
2
ω
r
-
4
π
α
R
c
=
2
π
B
c
R
0
-
R
=
k
B
R
′Suppose that the sampling function is narrow band and most of the values are from the ranges in the neighborhood of peak, which means those values ofR
′ satisfied k
B
R
′
<
π
/
2. The characteristic of the sampling functionSa implies some range-bin to range-bin interaction, which can be regarded as the varying parameterΔ
r. When Δ
r
=
Δ
R
/
2, it means there is no interaction between the radar range bins. The effect of range bins interaction may be written as(14)
Δ
r
=
Δ
R
2
,
Δ
R
,
⋯
k
Δ
R
2
⋯
,
k
=
1,2
,
⋯Due to random variableR
′
∈
[
-
Δ
r
,
Δ
r
], the interval of k
B
R
′ is [
-
k
π
/
2
,
k
π
/
2
]. Suppose that R
′ is have a uniform distribution, then the expectation ofSa function is given by(15)
E
S
a
k
B
R
′
=
∫
-
Δ
r
Δ
r
S
a
k
B
R
′
d
R
′
2
Δ
r
=
S
i
k
π
/
2
π
≈
1
kwhere S
i
(
x
)
=
∫
0
x
(
sin
⁡
t
/
t
)
d
t. Thus the received electric field equation for FMCW radar can be derived from (10) as(16)
E
R
D
ω
r
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
T
r
4
π
2
R
0
Δ
ρ
i
2
e
-
j
2
k
0
+
k
r
R
e
j
4
π
α
R
2
/
c
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
S
a
k
B
R
′Comparing with the electric field for pulsed radar ((
9
) in [18]), we can find that the main difference is theSa function, which is determined by the range bins interaction parameter k. This term makes the amplitude of ionospheric clutter of FMCW radar lower than pulsed radar.


### 2.2. Vertical Received Power for FMCW

In reality, the ionosphere has a slow time variation. ThusR
i
(
x
,
y
) in (2) should be modified as(17)
R
i
x
,
y
,
t
=
∑
K
→
i
,
ω
i
P
K
→
i
,
ω
i
e
j
K
→
i
·
ρ
→
i
+
ω
i
twhere ω
i is the frequency of the ionospheric irregularities. Supposing ω
i is small enough that time variation during a single sweep can be negligible. Therefore, only sweep to sweep variations are considered. Assuming that R
i
(
x
,
y
,
t
) is a stationary, homogeneous random process, it is useful to investigate the statistical characteristics of the received electric field by its autocorrelation function as(18)
R
τ
=
A
r
2
η
0
1
T
r
2
E
t
+
τ
E
∗
t
=
A
r
η
0
k
0
2
Δ
l
2
I
0
2
s
i
n
2
θ
i
32
π
4
R
0
2
Δ
ρ
i
4
∫
K
→
i
∫
ω
i
S
R
i
K
→
i
,
ω
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
2
S
a
2
k
B
R
′
e
j
ω
i
τ
d
K
→
i
d
ω
iwhere A
r is the receiving antenna effective area. ∗ represents complex conjugation, and S
R
i
(
K
→
i
,
ω
i
) is the spectral density of the ionospheric irregularities which has the wave vector K
→
i and frequency ω
i. The Fourier transform of autocorrelation function R
(
τ
) gives the power spectral density (PSD) as(19)
P
ω
d
=
F
R
τ
=
A
r
η
0
k
0
2
Δ
l
2
I
0
2
s
i
n
2
θ
i
32
π
4
R
0
2
Δ
ρ
i
4
∫
K
→
i
∫
ω
i
S
R
i
K
→
i
,
ω
i
×
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
2
S
a
2
k
B
R
′
δ
ω
d
-
ω
i
d
K
→
i
d
ω
iwhere ω
d is the observed Doppler frequency.Combining the layered ionosphere model proposed in [18] with the vertical ionospheric model of the phase spectrum in [15], S
R
i
(
K
→
i
,
ω
i
) can be expressed as(20)
S
R
i
K
→
i
,
ω
i
=
k
0
1
+
4
z
0
r
e
2
λ
0
2
k
0
log
⁡
z
0
z
0
-
z
0
′
8
π
3
κ
0
N
e
1
2
δ
κ
y
δ
ω
i
-
κ
x
υ
h
-
2
k
0
υ
υ
κ
0
2
+
κ
x
2
3
/
2
+
⋯where z
0 is the ionosphere thickness, z
0
′ is the height of below z
0 on the order of one radar wavelength, r
e is the classical electron radius, λ
0 is the wavelength of radar, κ
x and κ
y are the horizontal components of the density irregularities wave number κ, with x being the east coordinate and y being the north coordinate in geomagnetic coordinate system, κ
0
≈
1
0
-
4
m
-
1 is the “outer” scale length parameter, N
e
1
2 is the variance of the electron density fluctuations at the reflection height, and υ
h and υ
υ are the horizontal and vertical ionospheric plasma drift velocities.Finally, the vertical received PSD for FMCW from (19) is reduced to(21)
P
v
ω
r
,
ω
d
=
k
0
3
η
0
Δ
l
2
I
0
2
A
r
s
i
n
2
δ
θ
32
π
4
R
0
2
Δ
ρ
i
4
S
a
2
k
B
R
′
∫
0
Δ
ρ
i
J
0
ω
d
-
2
k
0
υ
v
υ
h
ρ
i
ρ
i
d
ρ
i
2
×
1
+
4
z
0
r
e
2
λ
0
2
k
0
log
⁡
z
0
z
0
-
z
0
′
8
π
3
κ
0
N
e
1
2
κ
0
2
+
ω
d
-
2
k
0
υ
v
/
υ
h
2
3
/
2
+
·
·
·


## 2.1. The Electric Field for FMCW Waveform

The time domain expression for the FMCW transmitting antenna current waveform within one sweep interval is given as(4)
i
t
=
I
0
e
j
ω
0
t
+
α
π
t
2
h
t
+
T
r
2
-
h
t
-
T
r
2whereI
0 is the peak current, ω
0 is the radian frequency, α
=
B
/
T
r is the frequency sweep rate, B is the sweep bandwidth, T
r is the sweep interval, and h
(
t
) is the Heaviside function. Unlike pulsed radar, the range method for FMCW radar is to use the differential frequency between the transmitting and receiving signal:(5)
R
=
c
Δ
f
2
α
=
c
f
t
-
f
r
2
αwhereR is the target range,c is the light speed, and f
t and f
r are the frequencies of transmitting and receiving signal, respectively.For dipole antenna transmitting FMCW signal, we have the approximate expression in the derivative:(6)
F
t
-
1
j
C
0
=
η
0
Δ
l
c
d
i
t
d
t
≈
j
k
0
η
0
Δ
l
I
0
e
j
ω
0
t
+
α
π
t
2
h
t
+
T
r
2
-
h
t
-
T
r
2where k
0
=
ω
0
/
c. Then (3) becomes(7)
E
R
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
j
ω
0
t
-
2
R
/
c
e
j
α
π
t
-
2
R
/
c
2
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
cAfter the demodulation processing, (7) reduces to(8)
E
R
D
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
-
j
4
π
α
R
2
/
c
2
-
2
k
0
R
e
j
4
π
α
R
/
c
t
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
cwhere the superscriptD represents demodulation.For FMCW, the target range is determined by the difference frequency of the demodulated electric field. Thus the range transform involves Fourier transform the electric field with respect tot for obtaining the difference frequency. Only the last two terms are function of t and their Fourier transforms are(9)
F
e
j
4
π
α
R
/
c
t
h
t
+
T
r
2
-
2
R
c
-
h
t
-
T
r
2
-
2
R
c
=
∫
2
R
/
c
-
T
r
/
2
2
R
/
c
+
T
r
/
2
e
j
4
π
α
R
/
c
t
e
-
j
ω
t
d
t
=
T
r
e
j
4
π
α
R
/
c
-
ω
2
R
/
c
S
a
T
r
2
ω
-
4
π
α
R
cInserting (9) into (8), the frequency domain expression for FMCW can be obtained:(10)
E
R
D
ω
r
=
F
E
R
D
t
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
T
r
4
π
2
R
Δ
ρ
i
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
e
-
j
2
k
0
+
ω
r
/
c
R
e
j
4
π
α
R
2
/
c
2
S
a
T
r
2
ω
r
-
4
π
α
R
cHere,ω has been renamed as ω
r and S
a
(
x
) is the sampling function of sin
⁡
(
x
)
/
x. For a given frequency ω
r, there is a corresponding delay time Δ
t which is the total delay between the ionospheric irregularities and radar. Then the apparent ionosphere reflection range R
0 and range resolution Δ
R can be defined as(11)
R
0
=
c
Δ
t
2
=
c
ω
r
4
π
α
,
Δ
R
=
c
2
BAnd we define(12)
k
B
=
2
π
B
c
,
k
r
=
ω
r
c
,
R
′
=
R
0
-
RThus the variable of the sampling function in (10) becomes(13)
T
r
2
ω
r
-
4
π
α
R
c
=
2
π
B
c
R
0
-
R
=
k
B
R
′Suppose that the sampling function is narrow band and most of the values are from the ranges in the neighborhood of peak, which means those values ofR
′ satisfied k
B
R
′
<
π
/
2. The characteristic of the sampling functionSa implies some range-bin to range-bin interaction, which can be regarded as the varying parameterΔ
r. When Δ
r
=
Δ
R
/
2, it means there is no interaction between the radar range bins. The effect of range bins interaction may be written as(14)
Δ
r
=
Δ
R
2
,
Δ
R
,
⋯
k
Δ
R
2
⋯
,
k
=
1,2
,
⋯Due to random variableR
′
∈
[
-
Δ
r
,
Δ
r
], the interval of k
B
R
′ is [
-
k
π
/
2
,
k
π
/
2
]. Suppose that R
′ is have a uniform distribution, then the expectation ofSa function is given by(15)
E
S
a
k
B
R
′
=
∫
-
Δ
r
Δ
r
S
a
k
B
R
′
d
R
′
2
Δ
r
=
S
i
k
π
/
2
π
≈
1
kwhere S
i
(
x
)
=
∫
0
x
(
sin
⁡
t
/
t
)
d
t. Thus the received electric field equation for FMCW radar can be derived from (10) as(16)
E
R
D
ω
r
=
j
k
0
η
0
Δ
l
I
0
sin
⁡
θ
i
T
r
4
π
2
R
0
Δ
ρ
i
2
e
-
j
2
k
0
+
k
r
R
e
j
4
π
α
R
2
/
c
2
∑
K
→
i
P
K
→
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
×
S
a
k
B
R
′Comparing with the electric field for pulsed radar ((
9
) in [18]), we can find that the main difference is theSa function, which is determined by the range bins interaction parameter k. This term makes the amplitude of ionospheric clutter of FMCW radar lower than pulsed radar.


## 2.2. Vertical Received Power for FMCW

In reality, the ionosphere has a slow time variation. ThusR
i
(
x
,
y
) in (2) should be modified as(17)
R
i
x
,
y
,
t
=
∑
K
→
i
,
ω
i
P
K
→
i
,
ω
i
e
j
K
→
i
·
ρ
→
i
+
ω
i
twhere ω
i is the frequency of the ionospheric irregularities. Supposing ω
i is small enough that time variation during a single sweep can be negligible. Therefore, only sweep to sweep variations are considered. Assuming that R
i
(
x
,
y
,
t
) is a stationary, homogeneous random process, it is useful to investigate the statistical characteristics of the received electric field by its autocorrelation function as(18)
R
τ
=
A
r
2
η
0
1
T
r
2
E
t
+
τ
E
∗
t
=
A
r
η
0
k
0
2
Δ
l
2
I
0
2
s
i
n
2
θ
i
32
π
4
R
0
2
Δ
ρ
i
4
∫
K
→
i
∫
ω
i
S
R
i
K
→
i
,
ω
i
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
2
S
a
2
k
B
R
′
e
j
ω
i
τ
d
K
→
i
d
ω
iwhere A
r is the receiving antenna effective area. ∗ represents complex conjugation, and S
R
i
(
K
→
i
,
ω
i
) is the spectral density of the ionospheric irregularities which has the wave vector K
→
i and frequency ω
i. The Fourier transform of autocorrelation function R
(
τ
) gives the power spectral density (PSD) as(19)
P
ω
d
=
F
R
τ
=
A
r
η
0
k
0
2
Δ
l
2
I
0
2
s
i
n
2
θ
i
32
π
4
R
0
2
Δ
ρ
i
4
∫
K
→
i
∫
ω
i
S
R
i
K
→
i
,
ω
i
×
∫
0
Δ
ρ
i
J
0
K
i
ρ
i
ρ
i
d
ρ
i
2
S
a
2
k
B
R
′
δ
ω
d
-
ω
i
d
K
→
i
d
ω
iwhere ω
d is the observed Doppler frequency.Combining the layered ionosphere model proposed in [18] with the vertical ionospheric model of the phase spectrum in [15], S
R
i
(
K
→
i
,
ω
i
) can be expressed as(20)
S
R
i
K
→
i
,
ω
i
=
k
0
1
+
4
z
0
r
e
2
λ
0
2
k
0
log
⁡
z
0
z
0
-
z
0
′
8
π
3
κ
0
N
e
1
2
δ
κ
y
δ
ω
i
-
κ
x
υ
h
-
2
k
0
υ
υ
κ
0
2
+
κ
x
2
3
/
2
+
⋯where z
0 is the ionosphere thickness, z
0
′ is the height of below z
0 on the order of one radar wavelength, r
e is the classical electron radius, λ
0 is the wavelength of radar, κ
x and κ
y are the horizontal components of the density irregularities wave number κ, with x being the east coordinate and y being the north coordinate in geomagnetic coordinate system, κ
0
≈
1
0
-
4
m
-
1 is the “outer” scale length parameter, N
e
1
2 is the variance of the electron density fluctuations at the reflection height, and υ
h and υ
υ are the horizontal and vertical ionospheric plasma drift velocities.Finally, the vertical received PSD for FMCW from (19) is reduced to(21)
P
v
ω
r
,
ω
d
=
k
0
3
η
0
Δ
l
2
I
0
2
A
r
s
i
n
2
δ
θ
32
π
4
R
0
2
Δ
ρ
i
4
S
a
2
k
B
R
′
∫
0
Δ
ρ
i
J
0
ω
d
-
2
k
0
υ
v
υ
h
ρ
i
ρ
i
d
ρ
i
2
×
1
+
4
z
0
r
e
2
λ
0
2
k
0
log
⁡
z
0
z
0
-
z
0
′
8
π
3
κ
0
N
e
1
2
κ
0
2
+
ω
d
-
2
k
0
υ
v
/
υ
h
2
3
/
2
+
·
·
·


## 3. Backscattering Received Power Spectrum for FMCW Source

The power spectrum of the HF signal phase backscattering from plasma structures in the auroral has been given in [15]. Assumptions are easily satisfied, such as the auroral magnetic field of the Earth which is approximately vertical to the surface, making the integral of phase autocorrelation function much simple. However, most of HFSWRs are located at mid-latitude. Therefore the magnetic field of the HFSWR location dissatisfies those conditions. Here, we rotate the coordinate axis twice, making the x-axis point to the east of geomagnetic field to eliminate thex component and making the z-axis coinciding with the local magnetic fieldB to eliminatez component. Then the integral of the autocorrelation function only leaves the integral abouty. The following is the mathematical modeling process.It is assumed thatα
D and α
I are the magnetic declination and magnetic inclination, respectively (see Figure 4). The coordinate of ray trajectory is X
1
=
x
1
,
y
1
,
z
1 in the geographic coordinate system, with x
1 east, y
1 north, and z
1 vertical. Now fixing z
1-axis, we rotate x
1 and y
1 axis α
D degrees, such that x
2 points to magnetic east and y
2 points to magnetic north. The magnetic field lies in the y
2-z
2 plane. The magnetic coordinate system is denoted byX
2:(22)
X
2
=
cos
⁡
α
D
sin
⁡
α
D
0
-
sin
⁡
α
D
cos
⁡
α
D
0
0
0
1
X
1Figure 4
Scattering geometry for ionospheric backscattering.Then we fixedx
2-axis, and rotated y
2 and z
2 axis π
/
2
-
α
I degrees, such that z
2 coincides with magnetic fieldB. The new coordinate system is denoted byX
3:(23)
X
3
=
1
0
0
0
sin
⁡
α
I
cos
⁡
α
I
0
-
cos
⁡
α
I
sin
⁡
α
I
X
2Thus the Y component inX
3 coordinate system can be expressed as(24)
y
3
=
-
sin
⁡
α
I
sin
⁡
α
D
x
1
+
sin
⁡
α
I
cos
⁡
α
D
y
1
+
cos
⁡
α
I
z
1The horizontal distanceL traveled in the ionosphere is given by(25)
L
=
y
3
-
y
3
′
=
-
sin
⁡
α
I
sin
⁡
α
D
x
1
-
x
1
′
+
sin
⁡
α
I
cos
⁡
α
D
y
1
-
y
1
′
+
cos
⁡
α
I
z
1
-
z
1
′where (
x
1
,
y
1
,
z
1
), (
x
1
′
,
y
1
′
,
z
1
′
) are the locations where the radar signal enters the ionosphere and backscattering takes place, respectively. In the auroral zone, α
D
=
0, α
I
=
π
/
2, then (25) becomes(26)
L
=
y
1
-
y
1
′Therefore this is a special case of (25). Next, we estimate the scope ofL. The ray trajectory has an inverted parabola form within the ionosphere (see (
31
) in [15]):(27)
z
z
0
=
-
1
4
c
s
c
2
θ
ρ
z
0
2
+
c
o
s
2
θwhere θ is the elevation angle with respect to zenith; ρ is the radial range in thex-yplane. The point (
x
,
y
)
=
(
0,0
) is the peak of the parabolic trajectory. At this point, we have(28)
z
=
z
0
c
o
s
2
θThe radar pulse obliquely incident on the ionosphere atz=0, then (27) gives(29)
ρ
=
z
0
sin
⁡
2
θFrom (24), (28), and (29), the scope ofL can be expressed as(30)
0
≤
L
≤
z
0
sin
⁡
2
θ
sin
⁡
α
I
cos
⁡
α
D
-
sin
⁡
α
D
+
c
o
s
2
θ
cos
⁡
α
IThus the approximate expression of the oblique backscattering phase spectrum in mid-latitude can be written as(31)
S
φ
1
κ
x
,
κ
y
,
ω
=
32
π
3
κ
0
z
0
r
e
λ
csc
⁡
θ
2
N
e
1
2
δ
κ
y
δ
ω
-
κ
x
v
d
x
-
2
k
0
v
d
y
×
sin
⁡
2
θ
sin
⁡
α
I
cos
⁡
α
D
-
sin
⁡
α
D
+
c
o
s
2
θ
cos
⁡
α
I
κ
0
2
+
κ
x
2
3
/
2where υ
d
x is the plasma drift velocity which is perpendicular to the radar beam in azimuthal direction; υ
d
y is the plasma drift velocity along the radar beam in azimuthal direction. Due to the scattering geometry, the plasma drift velocities can be represented as(32)
υ
d
x
=
v
h
cos
⁡
θ
+
v
v
sin
⁡
θ
(33)
υ
d
y
=
v
h
sin
⁡
θ
+
v
v
cos
⁡
θFinally, using this specific ionosphere spectral density model, the expression of oblique PSD for FMCW in mid-latitude reduces to(34)
P
o
ω
r
,
ω
d
=
k
0
3
η
0
Δ
l
2
I
0
2
A
r
s
i
n
2
θ
32
π
4
R
0
2
Δ
ρ
i
4
S
a
2
k
B
R
′
∫
0
Δ
ρ
i
J
0
ω
d
-
2
k
0
v
h
sin
⁡
θ
+
v
v
cos
⁡
θ
v
h
cos
⁡
θ
+
v
v
sin
⁡
θ
ρ
i
ρ
i
d
ρ
i
2
×
1
+
32
π
3
κ
0
r
e
λ
0
csc
⁡
θ
i
2
N
e
1
2
sin
⁡
2
θ
sin
⁡
α
I
cos
⁡
α
D
-
sin
⁡
α
D
+
c
o
s
2
θ
cos
⁡
α
I
k
0
κ
0
2
+
ω
d
-
2
k
0
v
h
sin
⁡
θ
+
v
v
cos
⁡
θ
/
v
h
cos
⁡
θ
+
v
v
sin
⁡
θ
2
3
/
2
+
·
·
·


## 4. Numerical Simulations and Experiment Results

### 4.1. Numerical Simulations

Firstly we simulate the expectation of the random functionS
a
2
(
k
B
R
′
), which is the main difference between the PSD of pulsed radar and FMCW radar, as shown in Figure 5. We can see that S
a
2
(
k
B
R
′
) changed with the range bins interaction parameterk, which are roughly coincident with (15). When the range resolution Δ
R is 5km (the sweep bandwidth is chosen as 30kHz),k=5 means the range bins interaction is 25km. The corresponding peak power of ionospheric clutter for FMCW waveform is lower than that of pulsed waveform about -
7
d
B. Ask increases, i.e., the radar signal is observed from increasingly larger ionospheric regions, the peak power of ionospheric clutter for FMWC waveform is reduced.Figure 5
The expectation of functionSa
2 with parameterk.In order to study the intensity of the ionosphere clutter and the sea clutter, we use the normalized ionospheric clutter power function defined as the ratio of the ionospheric clutter power to the average first-order sea clutter power. The main parameters are listed in Table1. The parameters are the same as those in [18], for the sake of comparison with pulsed waveform. Figure 6(a) shows the comparison results when FMCW parameters are with sweep bandwidthB= 30 kHz andk=10. It is observed that the peak power of FMCW is lower than that of pulsed waveform about -
10
d
B. The power of the ionospheric clutter for FMCW exceeds the average sea clutter power by about 35dB. Moreover, the main-lobe spectral width of FMCW is obviously narrower than that of the pulsed waveform. Figure 6(b) illustrates the changes in the power of ionospheric clutter with the range bins interaction. We can see that the peak power increases whenk decreases.Table 1
Main parameters in the simulation.
Symbol
Quantity
Parameters
f
0
Radar frequency
4.1MHz
R
0
Ionosphere height
300km
Δ
ρ
i
Scattering region radius
2.5km
z
0
Ionosphere thickness
50km
r
e
Classical electron radius
2.8×10−15m
υ
h
Horizontal ionospheric plasma drift velocity
100m/s
N
e
1
2
Variance of the electron density fluctuations
1018m−6
κ
0
Out scale length parameter
10−4m−1
k
range bins interaction parameter
10
B
sweep bandwidth
30kHzFigure 6
Simulations of vertical normalized ionospheric clutter power: (a) for the pulsed and FMCW waveforms and (b) with differentk.
(a)
(b)Figure7(a) illustrates the PSD of ionospheric clutter for FMCW from oblique backscattering with different incidence angles from 30° to 60°, when the vertical ionospheric velocity v
v is set as -5m/s, 5m/s, 15m/s, respectively. The horizontal ionospheric velocity v
h is set as 50m/s. Magnetic declination is set as 55°, and magnetic inclination is set as -7°, which are the magnetic field of Weihai (38°N, 122°E), China. We can find that the vertical ionospheric motion not only results in a Doppler shift, but also changes the peak power. This is different from the case of vertical reflection [18], where the vertical ionospheric motion only causes Doppler shift. The peak power of the ionospheric clutter with θ
=
π
/
3 exceeds that of θ
=
π
/
4 nearly 17dB.Figure 7
Simulations of oblique backscattering: (a) power with different incidence angles for FMCW. (b) L with different incidence angles.
(a)
(b)Figure7(b) illustrates the horizontal distance in ionosphere with L denoting the magnetic inclination varied from 0° (the equator zone) to 90° (the auroral zone) on the incidence angles. It is observed that L changed not only with magnetic inclinations but also with incidence angles. The maximum distance between different latitudes can reach 30km.Figure8 shows the simulation of the oblique range-Doppler power spectrum for FMCW (background white noise has been added). The main parameters for the three regions are as follows: the ranges are set as 100-150km, 200-250km, and 250-300km, respectively; the incidence angles are set as 30°, 45°, and 60°, respectively; the patch radii are set as 0.5km, 1.5km, and 2.5km, respectively; the variance of the electron density fluctuations is set as 1016m−6, 1018m−6, and 1020m−6, respectively; the vertical ionospheric velocity υ
υ is set as -5m/s, -10m/s, and 5m/s, respectively; the ionospheric velocity of ionospheric velocity υ
h is set as 50m/s, 75m/s, and 100m/s, respectively. The reason for setting these parameters sequentially is to make the peak power and Doppler bandwidth of ionospheric clutter expanding. The oscillating shapes of the power spectrum are caused by the integral over the zero-order Bessel function.Figure 8
Simulated oblique backscattering PSD of ionospheric clutter: (a) R-D spectrum. (b) 3D spectrum.
(a)
(b)By changing each parameter while keeping other parameters invariant, the simulation results show some characteristics as follows: the Doppler frequency of the power peak is mainly determined by the along radar beam velocities; the Doppler bandwidth spreading is mainly determined by the patch radii and the ionospheric velocities of cross radar beam; the peak of power is mainly determined by the variance of electron density fluctuations, radar frequencies, ionosphere ranges, and incidence angles.


### 4.2. Experiment Results

The HFSWR and Ionosonde are located at Weihai (37.5°N, 122°E), where is on the edge of the Yellow Sea. Experiment has been processed on 10:55 am, June 20, 2018 (Beijing Time). The data of HFSWR are collected during the coherent integration time (CIT) from 10:55 to 11:00. The radar operating frequency is 6MHz, and the sweep bandwidth is 30 kHz. The ionogram is collected in one minute: 11:02. The ionogram is shown in Figure9(a). It is observed that the critical frequency of E layer is about 5.5MHz, and the range spreads from 100km to 120km. Corresponding spreading in range can be found in Figure 9(b), where it is about 150-200km. This may be considered from the oblique backscattering for HFSWR. The extraordinary reflection in F layer is about 330-380km, which can also be found in RD spectrum. Because the locations of two figures are at a similar range, this may be considered as from the vertical reflection for HFSWR.Figure 9
The comparison of Ionogram and RD spectrum of HFSWR on 10:55 am, June 28, 2018 (Beijing Time), location: Weihai (37.5°N, 122°E): (a) Ionogram; (b) RD spectrum.
(a)
(b)Due to lack of the actual power of sea clutter, we cannot testify the model precisely. However, we can get the distribution shapes of the ionospheric irregularities and the relative intensity of the received power. Figures10(a) and 10(b) illustrate the inversion result of E layer and F layer, respectively.Figure 10
Experimental inversion results: (a) E layer; (b) F layer.
(a)
(b)The selected range of F layer is from 350km to 390km, and the selected Doppler width is from 0.08Hz to 0.23Hz. The F layer critical frequency is 6.3MHz, the corresponding maximum electron density is 4.9∗1011m−3. Therefore we set the average of the electron density fluctuation as 4.9∗ 109m−3, which is 1% of the maximum electron density. The inversion result shows that the main Doppler peak of irregularities is 0.1657Hz, which means the vertical drift velocities of the most fluctuation irregularities is about 4.14m/s. The main peak of irregularities in range distribution is 369km. Noted that the Doppler distribution is stratified and the oscillating shapes, which are consistent with Figure 8.The selected range of E layer is from 160km to 210km, and the selected Doppler width is from -0.5Hz to 0.5Hz. The E layer critical frequency is 5.5MHz, the corresponding maximum electron density is 3.7∗1011m−3. From the figure, we can see that the main peak of irregularities in range distribution is 198km. The range spreading, which should come from multipath reflection, also accords with RD spectrum in Figure 9(b). The peak Doppler is -0.265Hz, which implies the vertical drift velocities of the irregularities is about -6.625m/s. The average of the electron density fluctuation is about 3.9∗ 109m−3, which is almost 10% of the E layer maximum electron density. From Figure 10, the irregularities within ionosphere shows more Doppler characteristics, which means the drift motion of irregularities has a significant influence on ionospheric clutter.


## 4.1. Numerical Simulations

Firstly we simulate the expectation of the random functionS
a
2
(
k
B
R
′
), which is the main difference between the PSD of pulsed radar and FMCW radar, as shown in Figure 5. We can see that S
a
2
(
k
B
R
′
) changed with the range bins interaction parameterk, which are roughly coincident with (15). When the range resolution Δ
R is 5km (the sweep bandwidth is chosen as 30kHz),k=5 means the range bins interaction is 25km. The corresponding peak power of ionospheric clutter for FMCW waveform is lower than that of pulsed waveform about -
7
d
B. Ask increases, i.e., the radar signal is observed from increasingly larger ionospheric regions, the peak power of ionospheric clutter for FMWC waveform is reduced.Figure 5
The expectation of functionSa
2 with parameterk.In order to study the intensity of the ionosphere clutter and the sea clutter, we use the normalized ionospheric clutter power function defined as the ratio of the ionospheric clutter power to the average first-order sea clutter power. The main parameters are listed in Table1. The parameters are the same as those in [18], for the sake of comparison with pulsed waveform. Figure 6(a) shows the comparison results when FMCW parameters are with sweep bandwidthB= 30 kHz andk=10. It is observed that the peak power of FMCW is lower than that of pulsed waveform about -
10
d
B. The power of the ionospheric clutter for FMCW exceeds the average sea clutter power by about 35dB. Moreover, the main-lobe spectral width of FMCW is obviously narrower than that of the pulsed waveform. Figure 6(b) illustrates the changes in the power of ionospheric clutter with the range bins interaction. We can see that the peak power increases whenk decreases.Table 1
Main parameters in the simulation.
Symbol
Quantity
Parameters
f
0
Radar frequency
4.1MHz
R
0
Ionosphere height
300km
Δ
ρ
i
Scattering region radius
2.5km
z
0
Ionosphere thickness
50km
r
e
Classical electron radius
2.8×10−15m
υ
h
Horizontal ionospheric plasma drift velocity
100m/s
N
e
1
2
Variance of the electron density fluctuations
1018m−6
κ
0
Out scale length parameter
10−4m−1
k
range bins interaction parameter
10
B
sweep bandwidth
30kHzFigure 6
Simulations of vertical normalized ionospheric clutter power: (a) for the pulsed and FMCW waveforms and (b) with differentk.
(a)
(b)Figure7(a) illustrates the PSD of ionospheric clutter for FMCW from oblique backscattering with different incidence angles from 30° to 60°, when the vertical ionospheric velocity v
v is set as -5m/s, 5m/s, 15m/s, respectively. The horizontal ionospheric velocity v
h is set as 50m/s. Magnetic declination is set as 55°, and magnetic inclination is set as -7°, which are the magnetic field of Weihai (38°N, 122°E), China. We can find that the vertical ionospheric motion not only results in a Doppler shift, but also changes the peak power. This is different from the case of vertical reflection [18], where the vertical ionospheric motion only causes Doppler shift. The peak power of the ionospheric clutter with θ
=
π
/
3 exceeds that of θ
=
π
/
4 nearly 17dB.Figure 7
Simulations of oblique backscattering: (a) power with different incidence angles for FMCW. (b) L with different incidence angles.
(a)
(b)Figure7(b) illustrates the horizontal distance in ionosphere with L denoting the magnetic inclination varied from 0° (the equator zone) to 90° (the auroral zone) on the incidence angles. It is observed that L changed not only with magnetic inclinations but also with incidence angles. The maximum distance between different latitudes can reach 30km.Figure8 shows the simulation of the oblique range-Doppler power spectrum for FMCW (background white noise has been added). The main parameters for the three regions are as follows: the ranges are set as 100-150km, 200-250km, and 250-300km, respectively; the incidence angles are set as 30°, 45°, and 60°, respectively; the patch radii are set as 0.5km, 1.5km, and 2.5km, respectively; the variance of the electron density fluctuations is set as 1016m−6, 1018m−6, and 1020m−6, respectively; the vertical ionospheric velocity υ
υ is set as -5m/s, -10m/s, and 5m/s, respectively; the ionospheric velocity of ionospheric velocity υ
h is set as 50m/s, 75m/s, and 100m/s, respectively. The reason for setting these parameters sequentially is to make the peak power and Doppler bandwidth of ionospheric clutter expanding. The oscillating shapes of the power spectrum are caused by the integral over the zero-order Bessel function.Figure 8
Simulated oblique backscattering PSD of ionospheric clutter: (a) R-D spectrum. (b) 3D spectrum.
(a)
(b)By changing each parameter while keeping other parameters invariant, the simulation results show some characteristics as follows: the Doppler frequency of the power peak is mainly determined by the along radar beam velocities; the Doppler bandwidth spreading is mainly determined by the patch radii and the ionospheric velocities of cross radar beam; the peak of power is mainly determined by the variance of electron density fluctuations, radar frequencies, ionosphere ranges, and incidence angles.


## 4.2. Experiment Results

The HFSWR and Ionosonde are located at Weihai (37.5°N, 122°E), where is on the edge of the Yellow Sea. Experiment has been processed on 10:55 am, June 20, 2018 (Beijing Time). The data of HFSWR are collected during the coherent integration time (CIT) from 10:55 to 11:00. The radar operating frequency is 6MHz, and the sweep bandwidth is 30 kHz. The ionogram is collected in one minute: 11:02. The ionogram is shown in Figure9(a). It is observed that the critical frequency of E layer is about 5.5MHz, and the range spreads from 100km to 120km. Corresponding spreading in range can be found in Figure 9(b), where it is about 150-200km. This may be considered from the oblique backscattering for HFSWR. The extraordinary reflection in F layer is about 330-380km, which can also be found in RD spectrum. Because the locations of two figures are at a similar range, this may be considered as from the vertical reflection for HFSWR.Figure 9
The comparison of Ionogram and RD spectrum of HFSWR on 10:55 am, June 28, 2018 (Beijing Time), location: Weihai (37.5°N, 122°E): (a) Ionogram; (b) RD spectrum.
(a)
(b)Due to lack of the actual power of sea clutter, we cannot testify the model precisely. However, we can get the distribution shapes of the ionospheric irregularities and the relative intensity of the received power. Figures10(a) and 10(b) illustrate the inversion result of E layer and F layer, respectively.Figure 10
Experimental inversion results: (a) E layer; (b) F layer.
(a)
(b)The selected range of F layer is from 350km to 390km, and the selected Doppler width is from 0.08Hz to 0.23Hz. The F layer critical frequency is 6.3MHz, the corresponding maximum electron density is 4.9∗1011m−3. Therefore we set the average of the electron density fluctuation as 4.9∗ 109m−3, which is 1% of the maximum electron density. The inversion result shows that the main Doppler peak of irregularities is 0.1657Hz, which means the vertical drift velocities of the most fluctuation irregularities is about 4.14m/s. The main peak of irregularities in range distribution is 369km. Noted that the Doppler distribution is stratified and the oscillating shapes, which are consistent with Figure 8.The selected range of E layer is from 160km to 210km, and the selected Doppler width is from -0.5Hz to 0.5Hz. The E layer critical frequency is 5.5MHz, the corresponding maximum electron density is 3.7∗1011m−3. From the figure, we can see that the main peak of irregularities in range distribution is 198km. The range spreading, which should come from multipath reflection, also accords with RD spectrum in Figure 9(b). The peak Doppler is -0.265Hz, which implies the vertical drift velocities of the irregularities is about -6.625m/s. The average of the electron density fluctuation is about 3.9∗ 109m−3, which is almost 10% of the E layer maximum electron density. From Figure 10, the irregularities within ionosphere shows more Doppler characteristics, which means the drift motion of irregularities has a significant influence on ionospheric clutter.


## 5. Conclusions

This paper has developed an ionospheric clutter model for the vertical and backscatter reflection with FMCW source. The theory contained the development of the electric field equations and the derivation of the power spectrum for ionospheric clutter. The main difference of PSD between FMCW and pulsed waveform is the Sa function, which makes the peak power of FMCW in average lower than pulsed waveform, and the Doppler bandwidth of FMCW is much narrower than that of the pulsed waveform. The generalized power spectrum for backscatter sky-wave propagation (0.5 jumps) in mid-latitude was also derived.Numerical simulations illustrated that the Doppler bandwidth was spreading and the peak was heightening with the elevation angle increases. The 2D and 3D of RD spectrum were also simulated. At last, we processed the measured HFSWR data and inversed the variance of the electron density fluctuation and the vertical velocity of the irregularities within ionosphere by the Ionosonde data. Such researches could be beneficial to the better understanding of the characteristic of ionospheric clutter and suppressing ionospheric clutter techniques for HFSWR.



---

*Source: 1782942-2019-05-08.xml*