# Parabolic Equation Accumulated Split Error and Its Correction Method for Tropospheric Radio Propagation

Authors: Peipei Wei; Xiaoyan Du; Xinwei Hu; Changyin Jiang
Journal: International Journal of Antennas and Propagation
 (2018)

---

## Abstract

An increasingly popular approach of modeling electromagnetic wave propagation in the troposphere is to use the slip-step Fourier transform (SSFT) numerical algorithm to solve a parabolic equation (PE). However, available SSFT does not provide perfect solutions by introducing accumulated split error (ASE) for its step-by-step iteration process with a fixed propagation range step. The common analysis on one step split error is only a simple discussion on the split error of PE. Therefore, the main motivation of this study is to provide improvement on the accuracy of the PE by proposing an error correction method for the first time in the literature for PE. This method is on the basis of defining and deriving the accumulated split error. Its performance is then demonstrated with measured radar sea clutter data from Huanghai Sea, China.


---

## Body

## 1. Introduction

Most radio systems have to operate in the lower tropospheric atmosphere. Atmosphere absorbing, scattering, reflection, and refraction could affect radio propagation characteristics. In particular, the occurrence of tropospheric duct can make radio wave beyond-line-of-sight propagation [1], extending communication range and detection distance of radar system. This phenomenon also may cause radio holes and introduce new interfaces, which fundamentally can change the performance of a system that is designed to operate under standard atmospheric conditions [1]. This nonstandard atmosphere has a higher demand on effectively describing propagation characteristics.Over time, radio propagation modeling methods in tropospheric environment mainly focus on ray tracing (RT) [2, 3] and parabolic equation (PE) [4, 5] methods. RT generally neglects the influence of frequency and fails to describe the distribution of field strength at arbitrary location because of restricting to ray numbers. PE is a full-forward-wave approximation calculation method which is derived from wave equation using paraxial approximation [5]. It can deal with electromagnetic propagation problems in arbitrary irregular terrain and inhomogeneous atmosphere with the help of computer numeric calculation. Therefore, PE has been the most popular approach for tropospheric electromagnetic propagation [1–10]. According to propagation angle applied to the Helmholtz wave equation to obtain PE, PE can be divided into wide-angle PE (WPE) and narrow-angel PE (or standard PE (SPE)).At present, there are mainly three algorithms to solve PE: split-step Fourier transform (SSFT) [1], finite difference schemes [11], and finite element method [12], where SSFT has the advantage of higher numerical stability and thus is widely used to solve PE for tropospheric propagation problems. However, PE has approximation error due to paraxial approximation and split error as a result of exponent operator split using SSFT, where approximation error is its inherent drawback. For the past work on the error analysis of PE, approximation error is lower for SPE, and it can be ignored [13]. The expression of one-step split error of SPE is described in [14]. The two kinds of errors for WPE and SPE are discussed with simple calculation in [4]. Some advices are proposed to decrease the split error by studying its relationship with range step in [15]. However, an analytical description of the accumulation of split error (called accumulated split error (ASE)) over many range steps has yet to be successfully developed. In summary, these studies are inefficient to degrade split error. For this problem, we conduct the research on ASE and provide its correction method in this paper. In our work, the accumulated split error of SPE is derived for the first time in the literature on account of the definition of split error [4] and the iteration idea of SSFT. The derivation result shows that the ASE is determined by the atmosphere refractivity, the antenna frequency, the propagation elevation angle, the range step, and the field strength with split error. We simulate the sensitivity of ASE resulting from these factors using two statistical methods in a standard and nonstandard atmosphere. With the research conclusions, the ASE can be degraded by taking small range step in a standard atmosphere. However, this solution does not work well for nonstandard atmosphere on some conditions. For solving this problem, an ASE correction method is provided. This method is involved in complex calculation, including higher order accumulation, multiplication, and power operation. In addition, some supplements are provided for this method to analyze the complexity and efficiency. Finally, this method is verified by measured radar sea clutter from Huanghai Sea, China.The main contributions of this paper are fourfold. Firstly, we derive the accumulated split error of SPE for the first time in the literature. Secondly, the sensitivity analysis of ASE in a standard and nonstandard atmosphere is discussed using two statistical methods. Thirdly, we develop an ASE correction method to improve the accuracy of PE. Lastly, radar sea clutter power calculated by our proposed approach and past PE method is compared with the measured data from Huanghai Sea China, respectively.The remainder of the paper is organized as follows. Section2 presents the deviation of the ASE and its iteration model. Section 3 provides simulation results of sensitivity analysis. In Section 4, the ASE correction method is introduced. In addition, experiment research is provided. Lastly, the conclusions are presented in Section 5.


## 2. Parabolic Equation Accumulated Split Error

This section presents the derivation of the ASE on the basis of the one split error of PE.


### 2.1. Parabolic Approximation

In the rectangular coordinate system shown in Figure1, the SPE for the field ux,z is expressed as [5].
(1)∂2ux,z∂z2+2ik∂ux,z∂x+k2n2−1+2zaux,z=0,where k is the wave number and k=2π/λ (λ is the wavelength.), a is the Earth’s radius, and n is the refractive index. The term 2z/a accounts for the spherical shape of the Earth. When it is omitted, (1) describes propagation over a flat Earth.Figure 1
The calculation domain of PE.Only considering the forward wave propagation, (2) is obtained after simple transform of (1)
(2)∂ux,z∂x=A+Bux,z,where the expressions of two operator terms A and B are
(3)A=ik2w,B=i2k∂2∂z2,with w=n2x,z−1+2z/a. The horizontally inhomogeneous refractivity index n is considered here, shown as nx,z.Integrate overx for (2), then the solution is given as
(4)ux+Δx,z=expΔxA+Bux,z,where Δx is the range step.In order to solve this formula, error is incurred by splitting the exponential with the hypothesis thatn is a constant. The simplest split is given by [4].
(5)expΔxA+B≈expΔxAexpΔxB.Then the expression that can be solved is provided by(6)ux+Δx,z=expΔxAexpΔxBux,z.Using (6) and other boundary conditions shown in Figure 1, the field can be described in the propagation domain [5]. It is the commonly used method.The termw is varied with height in the realistic environment, then
(7)∂2wu∂z2≠w∂2u∂z2.Therefore, the two terms splitting in the exponent is not exact in nonhomogeneous medium. The split errorEerror is provided as
(8)Eerror=eB+AΔxux,z−eAΔxeBΔxux,z.This formula is the one-step split error at the rangex assuming that the field ux,z is the same for the right-hand side of (8). However, the field strength values of ux,z are different for the right-hand sides of (4) and (6), except for the initial field u0,z at the range x=0. Therefore, (8) is not applied for correcting the ASE between (4) and (6) over many range steps.


### 2.2. Accumulated Split Error Modeling

We assume that there is a solution meeting (4), then
(9)uj=expΔxA+Buj−1=F1uj−1,where F1=expΔxA+B, and uj is the field at the jth (0≤j≤Nx−1, Nx is the horizontal grid points) range step, which gives the result using PE without ASE.Letuj′ and uj−1′ obey (6), thus the field can be solved [4, 5], which gives the result using PE with ASE.
(10)uj′=expAΔx⋅expBΔxuj−1′=F2uj−1′,where F2=expAΔx⋅expBΔx.In the following derivations, the Hadamard product is introduced to describe the product and division of elements in the vector, where they are expressed as “∘” and “. /,” respectively.The ASERj at the jth range step is shown as
(11)Rj=uj−uj′=eB+AΔxuj−1−eAΔxeBΔxuj−1′,where uj′ and uj contain Nz vertical grid points, respectively.For the convenience of derivation, we define the correction factor(12)Sj=F1uj′−F2uj′./uj′=eB+AΔxuj′−eAΔxeBΔxuj′./uj′.With the detailed derivation found in the Appendix,Sj is rewritten as
(13)Sj=−Δx24∂nj∂z∘2+nj∘∂2nj∂z2+2iknj∘∂nj∂z+1a11⋮1sinα.Referring to the iteration process of SSFT, we derive the ASE.Given the initial fieldu0, the field of next range step x=Δx is shown as
(14)u1=F1u0,(15)u1′=F2u0.Using (12), the field u1 is given by
(16)u1=F1u0−F2u0./u0∘u0+F2u0=uj′+S0∘u0.For the range stepx=2Δx, u2 satisfies the following relationship
(17)u2−u2′=F1u1′−F2u1′./u1′∘u1′+F1S0∘u0=S1+S0∘u1′+S0∘2∘u0.Then we get(18)u2=u2′+S1+S0∘u1′+S0∘2∘u0.And so on, the field of rangex=jΔx is provided by
(19)uj=∑k=0j∑1Cjk∏p=0kSp∘tp∘uk′,where ∑p=0ktp=j−k,0≤tp≤j−k, tp is the integer, andCjk represents the combination, which can be written as Cjk=k!/j!k−j! using the factorials.We now get(20)Rj=uj−uj′=∑k=0j−1∑1Cjk∏p=0kSp∘tp∘uk′.For horizontally homogeneous atmosphere, the refractivity indexn only varies with z. Then the term Sj in (13) is the same at each range and unified as S0. Considering these relationships given in (14)–(19), Rj is obtained in homogeneous medium
(21)Rj=uj−uj′=∑i=1jCjiS0∘i∘uj−1′.We can conclude from above derivation: with thej increasing, the calculation complexity of uj increases; the ASE Rj is affected by environment parameters (nj, ∂nj/∂z, and ∂2nj/∂z2), antenna parameters (the frequency f), calculation range step Δx, propagation elevation angle α, and uj′; the derivation process is involved in higher order accumulation, multiplication, and power operation. In order to effectively reduce the ASE, we study its sensitivity by varying these parameters.


## 2.1. Parabolic Approximation

In the rectangular coordinate system shown in Figure1, the SPE for the field ux,z is expressed as [5].
(1)∂2ux,z∂z2+2ik∂ux,z∂x+k2n2−1+2zaux,z=0,where k is the wave number and k=2π/λ (λ is the wavelength.), a is the Earth’s radius, and n is the refractive index. The term 2z/a accounts for the spherical shape of the Earth. When it is omitted, (1) describes propagation over a flat Earth.Figure 1
The calculation domain of PE.Only considering the forward wave propagation, (2) is obtained after simple transform of (1)
(2)∂ux,z∂x=A+Bux,z,where the expressions of two operator terms A and B are
(3)A=ik2w,B=i2k∂2∂z2,with w=n2x,z−1+2z/a. The horizontally inhomogeneous refractivity index n is considered here, shown as nx,z.Integrate overx for (2), then the solution is given as
(4)ux+Δx,z=expΔxA+Bux,z,where Δx is the range step.In order to solve this formula, error is incurred by splitting the exponential with the hypothesis thatn is a constant. The simplest split is given by [4].
(5)expΔxA+B≈expΔxAexpΔxB.Then the expression that can be solved is provided by(6)ux+Δx,z=expΔxAexpΔxBux,z.Using (6) and other boundary conditions shown in Figure 1, the field can be described in the propagation domain [5]. It is the commonly used method.The termw is varied with height in the realistic environment, then
(7)∂2wu∂z2≠w∂2u∂z2.Therefore, the two terms splitting in the exponent is not exact in nonhomogeneous medium. The split errorEerror is provided as
(8)Eerror=eB+AΔxux,z−eAΔxeBΔxux,z.This formula is the one-step split error at the rangex assuming that the field ux,z is the same for the right-hand side of (8). However, the field strength values of ux,z are different for the right-hand sides of (4) and (6), except for the initial field u0,z at the range x=0. Therefore, (8) is not applied for correcting the ASE between (4) and (6) over many range steps.


## 2.2. Accumulated Split Error Modeling

We assume that there is a solution meeting (4), then
(9)uj=expΔxA+Buj−1=F1uj−1,where F1=expΔxA+B, and uj is the field at the jth (0≤j≤Nx−1, Nx is the horizontal grid points) range step, which gives the result using PE without ASE.Letuj′ and uj−1′ obey (6), thus the field can be solved [4, 5], which gives the result using PE with ASE.
(10)uj′=expAΔx⋅expBΔxuj−1′=F2uj−1′,where F2=expAΔx⋅expBΔx.In the following derivations, the Hadamard product is introduced to describe the product and division of elements in the vector, where they are expressed as “∘” and “. /,” respectively.The ASERj at the jth range step is shown as
(11)Rj=uj−uj′=eB+AΔxuj−1−eAΔxeBΔxuj−1′,where uj′ and uj contain Nz vertical grid points, respectively.For the convenience of derivation, we define the correction factor(12)Sj=F1uj′−F2uj′./uj′=eB+AΔxuj′−eAΔxeBΔxuj′./uj′.With the detailed derivation found in the Appendix,Sj is rewritten as
(13)Sj=−Δx24∂nj∂z∘2+nj∘∂2nj∂z2+2iknj∘∂nj∂z+1a11⋮1sinα.Referring to the iteration process of SSFT, we derive the ASE.Given the initial fieldu0, the field of next range step x=Δx is shown as
(14)u1=F1u0,(15)u1′=F2u0.Using (12), the field u1 is given by
(16)u1=F1u0−F2u0./u0∘u0+F2u0=uj′+S0∘u0.For the range stepx=2Δx, u2 satisfies the following relationship
(17)u2−u2′=F1u1′−F2u1′./u1′∘u1′+F1S0∘u0=S1+S0∘u1′+S0∘2∘u0.Then we get(18)u2=u2′+S1+S0∘u1′+S0∘2∘u0.And so on, the field of rangex=jΔx is provided by
(19)uj=∑k=0j∑1Cjk∏p=0kSp∘tp∘uk′,where ∑p=0ktp=j−k,0≤tp≤j−k, tp is the integer, andCjk represents the combination, which can be written as Cjk=k!/j!k−j! using the factorials.We now get(20)Rj=uj−uj′=∑k=0j−1∑1Cjk∏p=0kSp∘tp∘uk′.For horizontally homogeneous atmosphere, the refractivity indexn only varies with z. Then the term Sj in (13) is the same at each range and unified as S0. Considering these relationships given in (14)–(19), Rj is obtained in homogeneous medium
(21)Rj=uj−uj′=∑i=1jCjiS0∘i∘uj−1′.We can conclude from above derivation: with thej increasing, the calculation complexity of uj increases; the ASE Rj is affected by environment parameters (nj, ∂nj/∂z, and ∂2nj/∂z2), antenna parameters (the frequency f), calculation range step Δx, propagation elevation angle α, and uj′; the derivation process is involved in higher order accumulation, multiplication, and power operation. In order to effectively reduce the ASE, we study its sensitivity by varying these parameters.


## 3. Sensitivity Analysis

On account of the derivation results of ASE, its sensitivity involved by various parameters is addressed using computer calculation and statistical analysis methods. Standard and nonstandard atmosphere are examined, where the nonstandard atmosphere takes evaporation duct model as an example.Since the transmission loss is the most commonly used parameter to describe the radio propagation characteristics, we utilize this quantity instead of the field to discuss the sensitivity of ASE. The transmission loss calculated fromuj′ and uj is expressed as Ldepart and Ljoin, respectively. Simulations are given to compare the differences between Ldepart and Ljoin.Antenna parameters are shown in Table1. The distributions of Ldepart and Ljoin in a standard atmosphere are described in Figure 2 with the condition of f = 10 GHz, and Δx=500 m (calculation domain: 60 km × 75 m). Besides, the distributions of Ldepart and Ljoin in an evaporation duct are described in Figure 3 with the same calculation condition with Figure 2, where the simplified Paulus-Jeske (PJ) evaporation duct model [16] is taken with duct height of 20 m.Table 1
Antenna parameters.
Parameters
Values
Frequency
1~12 GHz
Antenna height
10 m
Elevation angle
0°
Beam width
0.5°
Polarization
Horizontal polarizationFigure 2
The transmission loss distribution in a standard atmosphere. (a) The distribution ofLdepart. (b) The distribution of Ljoin.
(a)
(b)Figure 3
The transmission loss distribution for a 20 m evaporation duct. (a) The distribution ofLdepart. (b) The distribution of Ljoin.
(a)
(b)Comparing Figure2(a) to Figure 2(b), the effect of ASE on transmission loss is more prominent at the range greater than 20 km. It is because that the cumulative effect of split error increases with the range increasing in a standard atmosphere. Comparing Figure 3(a) to Figure 3(b), the difference of transmission loss is greater inside the duct, especially nearby the ocean. The reason is that the greater gradient varieties of refractivity index along the height direction in (13) lead to greater ASE.In order to further discuss the effect of environmental parameters,Δx, and f on the ASE, two statistical methods are utilized: one dissimilarity method and one similarity method.The first method defines the mean valueΔL of relative percentage difference which is calculated from the absolute difference between Ljoinx,z and Ldepartx,z, and dividing by Ldepartx,z. Therefore, the mean characteristics of transmission loss difference in the given propagation domain are described.
(22)ΔL=∑x∑zLjoinx,z−Ldepartx,z/Ldepartx,z×100%Nx⋅Nz.The second method applies correlation coefficientR between Ldepart and Ljoin; thus, the bias from the mean value is considered.
(23)R=∑x∑zLjoinx,z−μjoinLdepartx,z−μdepart∑x∑zLjoinx,z−μjoin2∑x∑zLdepartx,z−μdepart2,where μjoin and μdepart are the mean of Ldepart and Ljoin, respectively.With the two methods, the simulation is conducted with the frequencies of L band, S band, C band, and X band, which are commonly used for microwave electromagnetic propagation. In addition, the range steps are taken from 50 m to 500 m.


### 3.1. Sensitivity Analysis in a Standard Atmosphere

We utilize the above-mentioned two methods to analyze the sensitivity in a standard atmosphere.Figure4 presents the distributions of ΔL and R varying with f and Δx in a standard atmosphere. In the four frequency bands, ΔL increases with range step, but R decreases. As noticed, the range step has weak effect on the frequency of L band, making ΔL smaller than 1% and R greater than 0.99. Thus, the ASE can be neglected for PE in this band. In order to control ΔL ≤ 1% and R > 0.99, the range step of smaller than 350 m should be taken for the frequency of S band. We should use the range step smaller than 150 m for the frequency of C band and smaller than 100 m for the frequency of X band.Figure 4
The distribution ofΔL and R in a standard atmosphere. (a) ΔL and (b) R.
(a)
(b)


### 3.2. Sensitivity Analysis in an Evaporation Duct

Two statistical methods are applied for analyzing the sensitivity in an evaporation duct.Figure5 provides the distributions of ΔL and R varying with f and Δx in an evaporation duct, which is obviously different from Figure 4. Comparing Figure 5 to Figure 4, the difference between the maximum and the minimum of ΔL in Figure 5(a) is smaller than that in Figure 4(a), but the difference range of R in Figure 5(b) is greater than that in Figure 4(b). Therefore, the effect of ASE on transmission loss in an evaporation duct is relatively greater to that in a standard atmosphere. In order to control ΔL ≤ 1% and R>0.99, we utilize range step smaller than 300 m for the frequency of L band, smaller than 200 m for the frequency of S band, smaller than 150 m for the frequency of C band, and smaller than 150 m for the frequency of X band.Figure 5
The distribution ofΔL and R for the 20 m evaporation duct. (a) ΔL and (b) R.
(a)
(b)In contrast, examining the distribution of the overall metrics (ΔL and R) for the transmission loss inside the duct is shown in Figure 6. The comparisons of Figure 6, Figure 5, and Figure 4 show that Ldepart in agreement with Ljoin is the least close inside evaporation duct, which leads to the greater difference between the maximum and the minimum of R than that in an entire domain. This is because the first derivative of refractivity index dramatically varieties with heights and the second derivative is nonzero. Therefore, the ASE cannot be decreased by simply taking the smaller range step in the high frequency, such as X band.Figure 6
The distribution ofΔL and R for the region of inside duct (z<20m). (a) ΔL and (b) R.
(a)
(b)


## 3.1. Sensitivity Analysis in a Standard Atmosphere

We utilize the above-mentioned two methods to analyze the sensitivity in a standard atmosphere.Figure4 presents the distributions of ΔL and R varying with f and Δx in a standard atmosphere. In the four frequency bands, ΔL increases with range step, but R decreases. As noticed, the range step has weak effect on the frequency of L band, making ΔL smaller than 1% and R greater than 0.99. Thus, the ASE can be neglected for PE in this band. In order to control ΔL ≤ 1% and R > 0.99, the range step of smaller than 350 m should be taken for the frequency of S band. We should use the range step smaller than 150 m for the frequency of C band and smaller than 100 m for the frequency of X band.Figure 4
The distribution ofΔL and R in a standard atmosphere. (a) ΔL and (b) R.
(a)
(b)


## 3.2. Sensitivity Analysis in an Evaporation Duct

Two statistical methods are applied for analyzing the sensitivity in an evaporation duct.Figure5 provides the distributions of ΔL and R varying with f and Δx in an evaporation duct, which is obviously different from Figure 4. Comparing Figure 5 to Figure 4, the difference between the maximum and the minimum of ΔL in Figure 5(a) is smaller than that in Figure 4(a), but the difference range of R in Figure 5(b) is greater than that in Figure 4(b). Therefore, the effect of ASE on transmission loss in an evaporation duct is relatively greater to that in a standard atmosphere. In order to control ΔL ≤ 1% and R>0.99, we utilize range step smaller than 300 m for the frequency of L band, smaller than 200 m for the frequency of S band, smaller than 150 m for the frequency of C band, and smaller than 150 m for the frequency of X band.Figure 5
The distribution ofΔL and R for the 20 m evaporation duct. (a) ΔL and (b) R.
(a)
(b)In contrast, examining the distribution of the overall metrics (ΔL and R) for the transmission loss inside the duct is shown in Figure 6. The comparisons of Figure 6, Figure 5, and Figure 4 show that Ldepart in agreement with Ljoin is the least close inside evaporation duct, which leads to the greater difference between the maximum and the minimum of R than that in an entire domain. This is because the first derivative of refractivity index dramatically varieties with heights and the second derivative is nonzero. Therefore, the ASE cannot be decreased by simply taking the smaller range step in the high frequency, such as X band.Figure 6
The distribution ofΔL and R for the region of inside duct (z<20m). (a) ΔL and (b) R.
(a)
(b)


## 4. Accumulated Split Error Correction Method

On account of the above-mentioned simulations and analysis, the PE-accumulated split error correction method is proposed. This method is verified by transhorizon tropospheric duct test data.We propose the correction method of ASE contained inuj′ on the basis of its definition and derivation in Section 3. The method can decrease the distribution domain of greater error values, thus improve the calculation accuracy of PE for electromagnetic propagation in the tropospheric environment.


### 4.1. Accumulated Split Error Correction Method

The steps applied to correct the error are provided as follows, and its diagram is shown in Figure7.
(1)
Input the antenna parameters.
(2)
Build propagation environment using refractivity profile.
(3)
Calculate the propagation fielduj′ using PE [1].Figure 7
The ASE correction diagram.On the conditions of (1) and (2),uj′j∈0,Nx is computed using PE [1].
(4)
Compute the ASE correction factor.Using parameters from (1) and a given range stepΔx, ASE correction factor Sj (j∈1,Nx) is calculated by (13) in the environment of (2).
(5)
Correct the ASE.On the basis of (3) and (4), the ASE is corrected using (14)–(19).
(6)
Output the field without the ASE.Using (5) yields the fieldujj∈1,Nx.As noticed in Figure7, the SSFT to solve PE is involved in complex techniques, including the absorbing boundary, surface boundary modeling, and split-step calculation. Besides, the interaction process from (14) to (19) shows that higher order accumulation, multiplication, and power operation are introduced. Therefore, this correction method increases the calculation complexity to improve the accuracy of PE. Taking the standard atmosphere for an example, we provide the complexity and efficiency analysis. The simulation results are shown in Figure 8. These results are based on the two aspects as follows:
(1)
The statistical relationship between calculation step range and execution time is analyzed for PE with ASE method and PE without ASE method, respectively. These results are simulated using the computer CPU of i7-4270HQ and shown in Figure8(a).
(2)
Given the upper threshold of errorΔL, the statistical relationship between frequency and execution time is analyzed and given in Figure 8(b).Figure 8
The complexity and efficiency analysis of PE. (a) The statistical relationship between calculation step range and execution time. (b) The statistical relationship between frequency and execution time.
(a)
(b)As noticed in Figure8(a), the execution time decreases with the range step increasing for the two methods. However, the execution time of PE without ASE method (red color) is much greater than that of PE with ASE method (blue color) on the condition of the same range step, especially for smaller range step. This result agrees with foregoing analysis that the calculation complexity increases by involving in complex correction factor. On the basis of accepted error threshold and execution time in Figure 8(b), the range step can be determined for a given frequency by referring to Figure 8(a). Therefore, Figure 8(b) gives a reference for realistic application.


### 4.2. Experimental Example

The main purpose of this section is to verify the performance of the ASE correction method of PE.The sea clutter data (Figure9(a)) and corresponding data of a refractivity profile (Figure 9(b)) measured in Huanghai, China, August 1, 2014 are considered here. As noticed in Figure 9(a), the north direction is the azimuth 0o, and the sea clutter at 58.5o azimuth is used as an example to analyze the performance of parabolic equation ASE correction method.Figure 9
Measured data distribution. (a) Radar sea clutter strength. (b) Refractivity profile.
(a)
(b)The radar antenna parameters are given in Table2.Table 2
Antenna parameters.
Parameters
Values
Frequency
9.6 GHz
Antenna height
35 m
Elevation angle
0°
Beam width
1.5°
Polarization
Vertical polarizationSince there is a clear weather for this measured data, smooth sea surface and range-independent refractivity profile are considered. By filtering the weather echo, ground echo, noise, and so on of sea clutter strengthX in Figure 9(a), we can obtain the sea clutter power Pobs at range r using
(24)Pobs=0.5X−2−32−20logr−C−30,where C takes into account radar parameters and so on.Applying the measured modified refractivityM into the radar equation [17], the power can be computed by
(25)PcM,r=−2LM,r+10lgr+σ0+C,where LM,r is the transmission loss at range r calculated by PE with ASE and without ASE method; σ0 is the radar cross section of the sea surface.In order to decrease the effect of range step on calculation results, we takeΔx=150 m. The field uj′ is computed as the first step result. Then the field uj can be obtained. On the basis of the field, the theoretical power with ASE Pc′M,r and power without ASE PcM,r is shown in Figure 10, respectively, where the mean of error between Pc′ and Pobs is 12.3 dB, but that between Pc and Pobs is 9.4 dB. Therefore, the parabolic equation ASE correction method can effectively improve the calculation accuracy.Figure 10
Sea clutter power distribution.


## 4.1. Accumulated Split Error Correction Method

The steps applied to correct the error are provided as follows, and its diagram is shown in Figure7.
(1)
Input the antenna parameters.
(2)
Build propagation environment using refractivity profile.
(3)
Calculate the propagation fielduj′ using PE [1].Figure 7
The ASE correction diagram.On the conditions of (1) and (2),uj′j∈0,Nx is computed using PE [1].
(4)
Compute the ASE correction factor.Using parameters from (1) and a given range stepΔx, ASE correction factor Sj (j∈1,Nx) is calculated by (13) in the environment of (2).
(5)
Correct the ASE.On the basis of (3) and (4), the ASE is corrected using (14)–(19).
(6)
Output the field without the ASE.Using (5) yields the fieldujj∈1,Nx.As noticed in Figure7, the SSFT to solve PE is involved in complex techniques, including the absorbing boundary, surface boundary modeling, and split-step calculation. Besides, the interaction process from (14) to (19) shows that higher order accumulation, multiplication, and power operation are introduced. Therefore, this correction method increases the calculation complexity to improve the accuracy of PE. Taking the standard atmosphere for an example, we provide the complexity and efficiency analysis. The simulation results are shown in Figure 8. These results are based on the two aspects as follows:
(1)
The statistical relationship between calculation step range and execution time is analyzed for PE with ASE method and PE without ASE method, respectively. These results are simulated using the computer CPU of i7-4270HQ and shown in Figure8(a).
(2)
Given the upper threshold of errorΔL, the statistical relationship between frequency and execution time is analyzed and given in Figure 8(b).Figure 8
The complexity and efficiency analysis of PE. (a) The statistical relationship between calculation step range and execution time. (b) The statistical relationship between frequency and execution time.
(a)
(b)As noticed in Figure8(a), the execution time decreases with the range step increasing for the two methods. However, the execution time of PE without ASE method (red color) is much greater than that of PE with ASE method (blue color) on the condition of the same range step, especially for smaller range step. This result agrees with foregoing analysis that the calculation complexity increases by involving in complex correction factor. On the basis of accepted error threshold and execution time in Figure 8(b), the range step can be determined for a given frequency by referring to Figure 8(a). Therefore, Figure 8(b) gives a reference for realistic application.


## 4.2. Experimental Example

The main purpose of this section is to verify the performance of the ASE correction method of PE.The sea clutter data (Figure9(a)) and corresponding data of a refractivity profile (Figure 9(b)) measured in Huanghai, China, August 1, 2014 are considered here. As noticed in Figure 9(a), the north direction is the azimuth 0o, and the sea clutter at 58.5o azimuth is used as an example to analyze the performance of parabolic equation ASE correction method.Figure 9
Measured data distribution. (a) Radar sea clutter strength. (b) Refractivity profile.
(a)
(b)The radar antenna parameters are given in Table2.Table 2
Antenna parameters.
Parameters
Values
Frequency
9.6 GHz
Antenna height
35 m
Elevation angle
0°
Beam width
1.5°
Polarization
Vertical polarizationSince there is a clear weather for this measured data, smooth sea surface and range-independent refractivity profile are considered. By filtering the weather echo, ground echo, noise, and so on of sea clutter strengthX in Figure 9(a), we can obtain the sea clutter power Pobs at range r using
(24)Pobs=0.5X−2−32−20logr−C−30,where C takes into account radar parameters and so on.Applying the measured modified refractivityM into the radar equation [17], the power can be computed by
(25)PcM,r=−2LM,r+10lgr+σ0+C,where LM,r is the transmission loss at range r calculated by PE with ASE and without ASE method; σ0 is the radar cross section of the sea surface.In order to decrease the effect of range step on calculation results, we takeΔx=150 m. The field uj′ is computed as the first step result. Then the field uj can be obtained. On the basis of the field, the theoretical power with ASE Pc′M,r and power without ASE PcM,r is shown in Figure 10, respectively, where the mean of error between Pc′ and Pobs is 12.3 dB, but that between Pc and Pobs is 9.4 dB. Therefore, the parabolic equation ASE correction method can effectively improve the calculation accuracy.Figure 10
Sea clutter power distribution.


## 5. Conclusion

In this paper, we derived the ASE due to exponent split using SSFT numerical algorithm to solve PE. In addition, sensitivity analysis of ASE in standard and nonstandard atmosphere is addressed using two statistical methods. On account of above studies, we propose an ASE correction method for PE. We also verify that PE without ASE method has a greater accuracy by using the measured data from Huanghai, China. Although the ASE of SPE is only derived and simulated, the study idea we developed is also applied to WPE. This is also our future work.



---

*Source: 5765913-2018-04-02.xml*