Research Article  Open Access
Parabolic Equation Modeling of Propagation over Terrain Using Digital Elevation Model
Abstract
The parabolic equation method based on digital elevation model (DEM) is applied on propagation predictions over irregular terrains. Starting from a parabolic approximation to the Helmholtz equation, a wideangle parabolic equation is deduced under the assumption of forward propagation and the splitstep Fourier transform algorithm is used to solve it. The application of DEM is extended to the Cartesian coordinate system and expected to provide a precise representation of a threedimensional surface with high efficiency. In order to validate the accuracy, a perfectly conducting Gaussian terrain profile is simulated and the results are compared with the shift map. As a consequence, a good agreement is observed. Besides, another example is given to provide a theoretical basis and reference for DEM selection. The simulation results demonstrate that the prediction errors will be obvious only when the resolution of the DEM used is much larger than the range step in the PE method.
1. Introduction
For many years, the parabolic equation (PE) method has been widely used to model the propagation of electromagnetic waves in complex atmospheric, hydrological, and geographical environments. As a forward fullwave analysis method, PE modeling has attracted a lot of interest since it was first proposed by Leontovich and Fock [1] in 1946 to analyze the diffraction of radio waves around the earth. After that, Hardin and Tappert [2] imported the splitstep Fourier transform (SSFT) from underwater acoustic research. Kuttler and Dockery [3] introduced the mixed Fourier transform (MFT), and on this basis, they proposed the discrete mixed Fourier transform (DMFT) [4] in 1996. Since then, this method has been effectively applied.
So far, the three commonly used algorithms of PE are SSFT, the finitedifference method (FDM) [5], and the finiteelement method (FEM) [6], among which SSFT is an analytical method, while FDM and FEM are numerical methods. In terms of the problem of propagation in the troposphere under largescene, longdistance, and lowelevation angle conditions, SSFT has many advantages over FDM and FEM, such as good in stability, small in error, and high in efficiency. In addition, the elevation angle is another critical issue that should be considered, because the accuracy of the PE method could be guaranteed only when the elevation angle is within the scope [7]. The approximation accuracy of the standard PE model is limited to be no more than 10°, so it is also called the narrowangle parabolic equation (NAPE) method. Therefore, many wideangle parabolic equation (WAPE) models have been proposed, which include the Claerbout PE method [8], FeitFleck PE modeling [9], and the Padé PE method [10]. Among these models, the FeitFleck PE modeling is the most popular one and its elevation angle can be expanded to approximately 30°.
In the past decades, a number of scholars [11–18] have studied the problem of propagation over irregular terrain by using the PE method (see [7] for a comprehensive reference list). Ayasli [12] established a terrainmasking model that uniformly zeros the field below the altitude of the boundary at each step. In 1979, Beilis and Tappert [11] developed the terrainshift transform method, also referred to as the shift map. This method performs a general coordinate transformation for NAPE that flattens the irregular surface, thereby solving the problem through SSFT. Barrios [13] adopted this method to present the terrain PE model, which was applied to actual navy propagation calculations. In 1998, Janaswamy [14] improved this technique and introduced grid interpolation to determine the field distribution in height grids. Donohue and Kuttler [15] subsequently extended the shiftmap method to the WAPE, and due to the fact that curvature information may not be available when working with actual terrain problems, they proposed the piecewise linear shift map (PLSM), an alternative method verified as the most accurate approach to handle discretelysampled terrain data.
In recent years, the vectorial version of the PE method has been employed to analyze the field propagation processes in road and railway tunnels where the presence of dissipative and field depolarization processes cannot be neglected. Bernardi et al. [19] used a model based on a hybrid implicit/explicit finitedifference numerical scheme to solve the vector parabolic equation governing the electromagnetic field propagation in straight and curved tunnels. Zhang et al. [20] presented a hybrid channel modeling technique, which combines the ray tracing and vector parabolic equation methods, to enable the modeling of realistic railway scenarios including stations and long guideways within a unified simulation framework. Compared to the twodimensional (2D) one, the vectorial version of the PE method can take the presence of lateral propagation environment and field depolarization processes take into account. However, at present, there is no report about threedimensional (3D) terrain processing method based on SSFT, while irregular terrain modeling is very mature in the 2D PE method.
In this paper, we address the terrain problems in PE modeling and a WAPE method based on the digital elevation model (DEM) is developed using the SSFT algorithm. In fact, as a more advanced modeling approach, the DEM has already been applied in some propagation software programs. However, in the previous models, all the DEMs used are in the geographic coordinate system while in our model, the application of DEM is extended to the 3D Cartesian coordinate system and the effect of the resolution of the DEM on propagation predictions is discussed in order to provide a theoretical basis and reference for DEM selection.
2. The Parabolic Equation Method
2.1. WideAngle Parabolic Equation
In the ensuing theory, a time convention is assumed for the fields and is suppressed throughout. As Figure 1 shows, we are concerned with 2D electromagnetic problems where all the fields can be decomposed into horizontally and vertically polarized components propagating independently. The field component is defined by for horizontal polarization and for vertical polarization. It satisfies the 2D scalar wave equation where is the wave number in free space and is the relative index of refraction of the medium.
We introduce the reduced function associated with the paraxial direction and then the scalar wave equation in terms of can be expressed by
If the backward propagation is ignored, the PE corresponding to forward propagating waves will be where the pseudodifferential operator is defined by
According to the operator splitting proposed by Feit and Fleck [9],
The error is of the order of the coupled term . In a vacuum, expression (8) is exact. It provides a good approximation for radio wave propagation application; however, the elevation angle is limited to be approximately 30°.
If we write as where A and B are defined by then can be approximated by
Thus, the FeitFleck WAPE is obtained.
2.2. SplitStep Fourier Transform Solution
The above WAPE has the formal solution and if an appropriate reduced function is introduced then we have
The Fourier transform of is given by and according to (15) where is the Fourier transform of . Therefore,
Substitute (19) into (16), and we get
It means that the forward propagating field can be obtained at a given range from the field at a previous range, which is quite suitable for numerical computation.
2.3. Terrain Modeling Using Digital Elevation Model
As a digital model was created from terrain elevation data, a DEM is able to provide a 3D representation of a terrain’s surface with high efficiency, especially when the environment is complex enough. The storage format of the DEM data in the 3D Cartesian coordinate system is shown in Figure 2. The terrain elevation values of the study area are stored in a matrix whose size is determined by the resolution of the DEM.
Assuming that the transmitting antenna is located at Tx and the receiving antenna Rx, the terrain elevation data between Tx and Rx in a 2D vertical plane can be easily obtained from the DEM, and according to the PLSM [15], the slope over the piecewise linear terrain is given by where is the angle that the local terrain slope makes with the horizontal. Due to the fact that the tilting of the coordinate system by an angle relative to the surface has locally contracted the horizontal range step, the effective wavenumber in this case may be written as and the WAPE in (13) now becomes which can also be solved by using the SSFT algorithm. Besides, the phase discontinuity when propagating across a boundary is accounted for by
3. Results and Discussions
Two examples are given in this section. In the first one, the DEMbased modeling approach is compared with the shift map to validate its accuracy while the second one deals with the effect of DEM’s resolution on prediction errors. In the following simulations, the transmitting antenna is a horizontally polarized Gaussian antenna with beamwidth located 50 m above the origin. Assuming that , , , and . Besides, the terrain is perfectly conducting and a 3D Gaussian terrain profile is assumed, which has the analytical form where , , , , , , , , and . Its 3D digital map is shown in Figure 3, and the DEM version with a resolution of 100 m is illustrated in Figure 4.
In order to simplify the calculation process and compare the DEMbased method with the shift map, Tx is set to (0 m, 25 km, 50 m) and Rx is set to (50 km, 25 km, 50 m); then the 2D terrain profile between Tx and Rx is given by and thus, the shift map can be carried out.
3.1. Validation of the Proposed Method
As Figure 5 shows, the propagation factor is calculated using the shift map and the DEMbased method, respectively. The resolution of the DEM used is 100 m, for the reason that . It is easy to find that only a tiny difference is observed. Besides, both of them are able to yield good results in the far zone and the interference patterns of the direct wave and the groundreflected wave can be clearly seen. However, the results are not that satisfactory when the propagation angle is too large, because according to the PE method, electromagnetic waves propagate in a conical region along the range axis.
(a)
(b)
The propagation factor at range 50 km is also compared, as shown in Figure 6, to make a quantificational evaluation. A good agreement is seen with a mean error of 0.021 dB and a standard deviation of 0.241 dB.
3.2. Effect of DEM’s Resolution on Prediction Errors
Another issue of interest is the dependence of the prediction errors on the resolution of the DEM used. As shown in Table 1, DEMs with different resolutions are simulated and the average error and standard deviation are computed. All other parameters are the same with the previous example. Figure 7 shows the results in an intuitive way, and we can see that with the increase of the resolution, the standard deviation increases gradually while the average error changes very slowly and maintains close to zero. Besides, a slow increase in the standard deviation, as expected, is observed when the resolution is less than 100 m, which may be explained by the fact that the prediction errors will be obvious only when the resolution is much larger than the range step in the PE method.

4. Conclusion
In this paper, under the assumption of forward propagation, a DEMbased PE method is applied on propagation predictions over irregular terrains using the SSFT algorithm. The terrain is assumed to be perfectly conducting, and a 3D Gaussian terrain profile is simulated by using the DEMbased method and the shift map, respectively. As a result, a good agreement is observed. Furthermore, the simulation results indicate that the prediction errors will be obvious only when the resolution of the DEM used is much larger than the range step in the PE method.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by National Key Laboratory of Electromagnetic Environment and the Foundation for Innovative Research Groups of the National Natural Science Foundation of China (Grant no. 61621005).
References
 M. A. Leontovich and V. A. Fock, “Solution of propagation of electromagnetic waves along the Earth’s surface by the method of parabolic equations,” Journal of Physics, vol. 10, pp. 13–23, 1946. View at: Google Scholar
 R. Hardin and F. Tappert, “Applications of the splitstep Fourier method to the numerical solution of nonlinear and variable coefficient wave equations,” SIAM Review, vol. 15, p. 423, 1973. View at: Google Scholar
 J. R. Kuttler and G. D. Dockery, “Theoretical description of the parabolic approximation/Fourier splitstep method of representing electromagnetic propagation in the troposphere,” Radio Science, vol. 26, no. 2, pp. 381–393, 1991. View at: Publisher Site  Google Scholar
 G. D. Dockery and J. R. Kuttler, “An improved impedanceboundary algorithm for Fourier splitstep solutions of the parabolic wave equation,” IEEE Transactions on Antennas and Propagation, vol. 44, no. 12, pp. 1592–1599, 1996. View at: Publisher Site  Google Scholar
 J. F. Claerbout, Fundamentals of Geophysical Data Processing with Applications to Petroleum Prospecting, Blackwell Scientific Publications, 1985.
 I. Sirkova and H. E. HernandezFigueroa, “Local transparent boundary condition applied to the modeling of tropospheric ducting propagation,” Microwave and Optical Technology Letters, vol. 21, no. 5, pp. 343–346, 1999. View at: Publisher Site  Google Scholar
 P. Zhang, L. Bai, Z. Wu, and L. Guo, “Applying the parabolic equation to tropospheric groundwave propagation: a review of recent achievements and significant milestones,” IEEE Antennas and Propagation Magazine, vol. 58, no. 3, pp. 31–44, 2016. View at: Publisher Site  Google Scholar
 J. Claerbout, “Fundamentals of geophysical data processing,” Geophysical Journal International, vol. 86, no. 1, pp. 217–219, 1986. View at: Publisher Site  Google Scholar
 M. Feit and J. Fleck Jr., “Light propagation in gradedindex optical fibers,” Applied Optics, vol. 17, no. 24, pp. 3990–3998, 1978. View at: Publisher Site  Google Scholar
 M. D. Collins, “A splitstep Padé solution for the parabolic equation method,” The Journal of the Acoustical Society of America, vol. 93, no. 4, pp. 1736–1742, 1993. View at: Publisher Site  Google Scholar
 A. Beilis and F. Tappert, “Coupled mode analysis of multiple rough surface scattering,” The Journal of the Acoustical Society of America, vol. 66, no. 3, pp. 811–826, 1979. View at: Publisher Site  Google Scholar
 S. Ayasli, “SEKE: a computer model for low altitude radar propagation over irregular terrain,” IEEE Transactions on Antennas and Propagation, vol. 34, no. 8, pp. 1013–1023, 1986. View at: Publisher Site  Google Scholar
 A. E. Barrios, “A terrain parabolic equation model for propagation in the troposphere,” IEEE Transactions on Antennas and Propagation, vol. 42, no. 1, pp. 90–98, 1994. View at: Publisher Site  Google Scholar
 R. Janaswamy, “A curvilinear coordinatebased splitstep parabolic equation method for propagation predictions over terrain,” IEEE Transactions on Antennas and Propagation, vol. 46, no. 7, pp. 1089–1097, 1998. View at: Publisher Site  Google Scholar
 D. J. Donohue and J. R. Kuttler, “Propagation modeling over terrain using the parabolic wave equation,” IEEE Transactions on Antennas and Propagation, vol. 48, no. 2, pp. 260–277, 2000. View at: Publisher Site  Google Scholar
 K. Wang and Y. Long, “Propagation modeling over irregular terrain by the improved twoway parabolic equation method,” IEEE Transactions on Antennas and Propagation, vol. 60, no. 9, pp. 4467–4471, 2012. View at: Publisher Site  Google Scholar
 R. Bai, C. Liao, N. Sheng, and Q. Zhang, “Prediction of wave propagation over digital terrain by parabolic equation model,” in 2013 5th IEEE International Symposium on Microwave, Antenna, Propagation and EMC Technologies for Wireless Communications, pp. 458–461, Chengdu, China, 2013. View at: Publisher Site  Google Scholar
 D. D. Wang, X. L. Xi, Y. R. Pu, J. F. Liu, and L. L. Zhou, “Parabolic equation method for LoranC ASF prediction over irregular terrain,” IEEE Antennas and Wireless Propagation Letters, vol. 15, pp. 734–737, 2016. View at: Publisher Site  Google Scholar
 P. Bernardi, D. Caratelli, R. Cicchetti, V. Schena, and O. Testa, “A numerical scheme for the solution of the vector parabolic equation governing the radio wave propagation in straight and curved rectangular tunnels,” IEEE Transactions on Antennas and Propagation, vol. 57, no. 10, pp. 3249–3257, 2009. View at: Publisher Site  Google Scholar
 X. Zhang, N. Sood, J. K. Siu, and C. D. Sarris, “A hybrid raytracing/vector parabolic equation method for propagation modeling in train communication channels,” IEEE Transactions on Antennas and Propagation, vol. 64, no. 5, pp. 1840–1849, 2016. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 XiaoWei Guan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.