# Stationary Time Statistical Property of Ionospheric Clutter in High-Frequency Surface-Wave Radar

Authors: Shang Shang; Kangning He; Zhaobin Wang; Xuguang Yang
Journal: Journal of Electrical and Computer Engineering
 (2019)

---

## Abstract

In HFSWR (high-frequency surface-wave radar) system, the detection performance is impacted seriously by ionospheric clutter. Frequency selection is an effective method to avoid the effect of ionospheric clutter. The key to the method is the stationarity of ionospheric clutter over a period of time. This paper mainly researches the stationary time statistical property of the ionospheric clutter. A large number of real data including ionospheric clutter in HFSWR are processed and analyzed. It shows that ionospheric clutter in HFSWR has the characteristics of approximate stationarity within a period of time.


---

## Body

## 1. Introduction

HFSWR is used to achieve overhorizon detection of surface ships and ultralow altitude aircraft by the electromagnetic wave characteristics in which high-frequency band (3–15 MHz) vertical polarization electromagnetic wave spreads along the sea surface [1–3]. HFSWR is considered to be a major technical means to achieve coastal state control and supervision of the offshore exclusive economic zone. HFSWR is regarded as perfect early warning radar. Nowadays, many military and civil fields adopted HFSWR, such as coastal surveillance and early warning, sea-state remote sensing [4], ocean surface pollution monitoring, marine environment monitoring, and resource development [5].However, HFSWR is inevitably affected by ionospheric clutter, sea clutter, radio interference, and other industrial interference. The ionospheric clutter has many characteristics, such as strong energy, long existence time, covering multiple distance units, and Doppler units. All the targets over 100 km may not be detected when ionospheric clutter is powerful [6, 7]. It has become the bottleneck of HFSWR. Therefore, the study of ionospheric clutter suppression is a focus of research on HFSWR.The study of ionospheric clutter characteristics is the base of the suppression. The researchers of various countries have researched its signatures from various angles. In [8], the range, Doppler, and direction signatures of ionospheric clutter have been analyzed. The document pointed that the detection range polluted by ionospheric clutter depends on the altitude of ionosphere which reflects electromagnetic waves. The Doppler frequencies of the electromagnetic wave reflected from various layer are different. The Doppler frequencies of E layer ionospheric clutter are low, and those of F layer are relatively higher. Most ionospheric clutters do not have directionality, or the directionality is not strong. In [9], the directionality of ionospheric clutter has been analyzed in detail, and it indicated that ionospheric clutters are divided into two categories, which are specular reflection clutter and spread clutter. Specular reflection clutter affects few Doppler and range cells. Spread clutter takes up more cells. Specular reflection clutters on the same range cell and from different Doppler cells have little similarity directivity. Jiang et al. [10] pointed that ionospheric clutters which take up fewer cells have the same directivity, but clutters which occupy many cells do not have. In [11], the effects of E layer clutter on HFSWR have been analyzed. The paper indicated that the Es clutter changes with the season and time in the Arctic. Time-frequency characteristics of ionospheric clutter are analyzed in [12]. In [13], the statistical properties of clutter intensity, correlation function, etc., have been analyzed. Plasma theory indicated that the reflection of electromagnetic waves on the ionosphere is related to the electromagnetic wave frequency. In other words, ionospheric clutter characteristics are in connection with the frequency. On some frequencies, the effects of ionospheric clutter are serious and even make radar fail to work. But on some other frequency, it will not happen. According to the conclusion, Shang et al. [14] proposed a method, which is used to avoid ionospheric clutter by choosing frequency. But the basis is time stationarity of ionospheric clutter.According to the above analysis, this paper takes the time statistical characteristics of ionospheric clutter as the research object. First, the ionospheric clutter is preprocessed, and the ionospheric clutter is expressed as a time series. Then, the trend term and stochastic term are separated by wavelet analysis. Finally, a stationary time series analysis is performed on the ionospheric clutter time series.


## 2. Ionospheric Clutter Time Series Preprocessing

HFSWR adopts an ARD (angle-range-Doppler frequency) spectrum to detect targets. The three-dimensional detection platform is established by the coherent integration and digital beam forming of the time-domain data. As shown in Figure1, the system works on N frequencies during different integration times. N ARD spectrums can be obtained during NTs, and Ts is the integration time for each frequency. ARDi is on frequency fi, and fi is one of the alternative frequencies to avoid ionospheric clutter. Time series ARDi,i=1,2,…,N, is obtained by arranging the ARD spectrums in chronological order.Figure 1
Illustration of ARD spectrum time series.The affected scope in the ARD spectrum and average power of ionospheric clutter are regarded as the index to evaluate the changing situation.Vi represents the number of cells polluted by ionospheric clutter in ith ARD spectrum of the time series ARDi. Then, VIC=V1,V2,…,VN is a time series reflecting the ionospheric clutter pollution scope. Similarly, PIC=P1,P2,…,PN is a time series on behalf of ionospheric clutter average power, and Pi is the average power of ionospheric clutter in ARDi. VIC and PIC are normalized for comparison.The time series of the affected scope and average power of ionospheric clutter are disturbed by complexity and randomness of ionospheric clutter, together with the environmental noise. Time seriesVIC and PIC need preprocessing to eliminate or minimize the effect of disturbance. The series VIC (PIC) can be regarded as nonstationary time series. And VIC (PIC) includes the trend term Vd (Pd) and the stochastic term Vr (Pr). They are shown as follows:(1)VIC=Vd+Vr,PIC=Pd+Pr.The trend terms reflect the overall trend of pollution scope and average power of ionospheric clutter. They are the judgment basis of ionospheric clutter to be stationary. The stochastic term reflects the randomness ofVIC and PIC, and they are the disturbance factor. Stationary random series analysis and wavelet transformation are adopted to preprocess the time series of the pollution scope and average power of ionospheric clutter, and the trend terms are extracted. Then, the statistical computations of ionospheric clutter approximate stationary time are performed.


### 2.1. One-Dimensional Wavelet Transform

Essentially, WT (wavelet transform) analyzes signals by primary function with constantly changing analysis width to realize different resolution ratios in different ranges of frequency [15]. On the basis of multiresolution theory and filter bank theory, the one-dimensional wavelet decomposition is that original signal Y passes through a dual-band bandpass filter group constituted by a high-pass filter and a low-pass filter bank, and high-frequency signal and low-frequency signal are extracted. Then, the low-frequency signal is further decomposed by passing through a dual-band bandpass channel filter, and this process is repeated until the low-frequency signal is decomposed to the mth layer. Signals obtained after the final decomposition are reconstructed by wavelet transform, and the original signal X can be expressed as(2)X=Lm+H1+H2+⋯+Hm,where Lm is the reconstruction result of the mth layer low-frequency signal and Hi is the reconstruction result of the ith layer high-frequency signal, and i=1,2,…,m.During the analysis of ionospheric clutter time series, low-frequency signal reconstruction resultLm reflects the variation trend as the trend term, and high-frequency signal reconstruction result H1+H2+⋯+Hm reflects the interference factors as the stochastic term.


### 2.2. Stationary Time Series Analysis

The only way we can regard the high-frequency signal as interference is to guarantee that it is a stable zero-mean random process. That is to say that the random item does not include trend items. So, during the extraction of trend items, it is necessary to take stationarity and zero-mean test.The common methods adopted for the stationary test include root test, autocorrelation function, partial autocorrelation functions, and parameter test. This paper adopts relatively the objective root test method. The basic idea is to set up time series mode first and then to solve the characteristic root of the characteristic equation established by the model parameter. Finally, the stationarity of time series is judged by determining whether the characteristic roots are in the unit circle. If all the characteristic tootsλ<1, that is, mean in unit circle, then the time series is regarded as stationary.Zero‐mean test is the hypothesis testing carried to judge whether the series mean is zero, that is, to testH0:μ=0. First variable U is defined as(3)U=x¯−μsn,where x¯ is the mean of the series sample, s is the standard deviation of the series sample, and n is the length of the series sample.When the hypothesisH0 is tenable, U approximately subjects to the normal distribution based on the law of large numbers and central-limit theorem. With regard to the confidence level α (here α=0.05), critical value μα/2 can be obtained by table look-up or calculation. If x¯−μ≤μα/2, the hypothesis that the mean of the series is zero is accepted.


## 2.1. One-Dimensional Wavelet Transform

Essentially, WT (wavelet transform) analyzes signals by primary function with constantly changing analysis width to realize different resolution ratios in different ranges of frequency [15]. On the basis of multiresolution theory and filter bank theory, the one-dimensional wavelet decomposition is that original signal Y passes through a dual-band bandpass filter group constituted by a high-pass filter and a low-pass filter bank, and high-frequency signal and low-frequency signal are extracted. Then, the low-frequency signal is further decomposed by passing through a dual-band bandpass channel filter, and this process is repeated until the low-frequency signal is decomposed to the mth layer. Signals obtained after the final decomposition are reconstructed by wavelet transform, and the original signal X can be expressed as(2)X=Lm+H1+H2+⋯+Hm,where Lm is the reconstruction result of the mth layer low-frequency signal and Hi is the reconstruction result of the ith layer high-frequency signal, and i=1,2,…,m.During the analysis of ionospheric clutter time series, low-frequency signal reconstruction resultLm reflects the variation trend as the trend term, and high-frequency signal reconstruction result H1+H2+⋯+Hm reflects the interference factors as the stochastic term.


## 2.2. Stationary Time Series Analysis

The only way we can regard the high-frequency signal as interference is to guarantee that it is a stable zero-mean random process. That is to say that the random item does not include trend items. So, during the extraction of trend items, it is necessary to take stationarity and zero-mean test.The common methods adopted for the stationary test include root test, autocorrelation function, partial autocorrelation functions, and parameter test. This paper adopts relatively the objective root test method. The basic idea is to set up time series mode first and then to solve the characteristic root of the characteristic equation established by the model parameter. Finally, the stationarity of time series is judged by determining whether the characteristic roots are in the unit circle. If all the characteristic tootsλ<1, that is, mean in unit circle, then the time series is regarded as stationary.Zero‐mean test is the hypothesis testing carried to judge whether the series mean is zero, that is, to testH0:μ=0. First variable U is defined as(3)U=x¯−μsn,where x¯ is the mean of the series sample, s is the standard deviation of the series sample, and n is the length of the series sample.When the hypothesisH0 is tenable, U approximately subjects to the normal distribution based on the law of large numbers and central-limit theorem. With regard to the confidence level α (here α=0.05), critical value μα/2 can be obtained by table look-up or calculation. If x¯−μ≤μα/2, the hypothesis that the mean of the series is zero is accepted.


## 3. The Algorithm of Trend Term Extraction of Ionospheric Clutter Time Series

X (here it represents VIC or PIC) is obtained by real data normalization processing. The purpose of trend term extraction is to eliminate interference composition at most. Trend information should retain integrally. Multilayer wavelet transform is adopted to decompose X. The layer number of WT is decided by stationary-zero‐mean testing for detail signals. Figure 2 shows the concrete flow. The high-frequency signal, which is from the signal decomposition of signal X by the first layer WT, is tested by stable zero‐mean. If the high-frequency signal which stands for the detail is a stable zero‐mean series, then the low-frequency signal should be carried on for second time wavelet decomposition. By that analogy, the number of WT is denoted by m when the high-frequency signal dissatisfies stable zero-mean test. The scale signal got from the m−1th WT decomposition is the trend term.Figure 2
Illustration of the trend extraction algorithm.Take the continuous observation data from 12:35 noon to 15:55 afternoon as example. The coherent integration timeTCIT=1.9min, and ionospheric clutter is identified from the ARD spectrum during the observation time. The scope VIC and mean power PIC of normalized ionospheric clutter are shown by the dotted line, and the trend term of them are shown by the solid line in Figures 3(a) and 3(b).Figure 3
Result of trend extraction. (a) Time seriesVIC. (b) Time series PIC.
(a)
(b)


## 4. The Approximate Stability Time

Ionospheric clutter fluctuates with time. When the affected scope and average power both do not exceed the certain threshold, it is considered approximately stable. The time during which ionospheric clutter keeps an approximate stable station is called the ionospheric clutter approximate stationary time.The approximate stability time of ionospheric clutter is the time during which the trend terms of normalized ionospheric clutter pollution scope time seriesVIC and average power time series PIC both do not exceed the threshold value δ (here δ=0.1). In Figure 4(a), STVi denotes the time during which the variability scope of trend term of VIC does not exceed δ begin with ti. In Figure 4(b), STPi denotes the time during which the variability scope of trend term of PIC does not exceed δ begin with ti. Then, the approximate stability time of ionospheric clutter beginning from ti STi=minSTVi,STPi.Figure 4
Illustration of ionospheric clutter nearly static time. (a) Nearly static time ofVIC. (b) Nearly static time of PIC.
(a)
(b)During the observation time, the ionospheric clutter stability time beginning with various time periods is shown as histogram in Figure5(a). The normal fitted curve is drawn on the normal probability paper shown in Figure 5(b), and it is almost a straight line. The approximate stability time of ionospheric clutter conforms to the normal distribution, and the mean value μ=61.02min and standard deviation σ=2.21min.Figure 5
Normality test of ionospheric clutter nearly static time. (a) Density of ionospheric clutter nearly static time. (b) Probability of ionospheric clutter nearly static time.
(a)
(b)For normal distribution, the ratio that the approximate stability time is within two standard deviations (approximate stability time is within 56.60–65.44 min) is 95.4%. Then, the ratio that the approximate stability time is longer than 56.60 min is bigger than 95.4%. Therefore, for the method to avoid ionospheric clutters by choosing frequency, the ionospheric clutters on various frequencies within 25 min are real-time monitored and the effects on HFSWR are analyzed. The frequency on which the influence of ionospheric clutter is slight is selected, and the result of frequency chosen applies to the next 25 min.


## 5. Conclusion

Approximate stability property and approximate stationary time of ionospheric clutter are defined based on the analysis of real data of HFSWR, and statistical analysis is performed.The ionospheric clutter pollution scope and average power time series are preprocessed by stationary time series analysis and WT. The statistical computing of the approximate stability time of ionospheric clutter is taken by the extraction of trend term. The threshold value of ionospheric clutter approximate stability is chosen asδ=0.1. The analysis of a large number of measured data shows that the rate that the approximate stability time is longer than 56.60 min is larger than 95.4%. It means that on some frequency, ionospheric clutter is approximately stable over 50 min.It provides the theoretical basis to avoid ionospheric clutter by choosing frequency. According to the analysis results, it is reasonable that ionospheric clutter is monitored and analyzed within 25 min, and it is the basis for choosing frequency during next 25 min.



---

*Source: 2450191-2019-11-03.xml*