IONOLAB-RAY: A wave propagation algorithm for anisotropic and inhomogeneous ionosphere

A high-frequency (HF) radio wave propagation model based on ray tracing for anisotropic, inhomogeneous, and time-varying ionosphere is developed to form the basis for a user-friendly modular algorithm. The algorithm IONOLAB-RAY applies the principles of ray tracing through the spherical 3D grid model of the ionosphere. The physical parameters of the ionosphere in each voxel are obtained from the ionosphere model IRI-Plas-G. The refractive index of each voxel is calculated with the Appleton Hartree formula, which includes electron cyclotron frequency, collision frequency, and anisotropicity due to the geomagnetic field. The developed propagation algorithm allows the input of total electron content (TEC) values to update the ionospheric model to the current conditions. This capability has an important impact on the characterization of the ionosphere under geomagnetic storm conditions. The algorithm can be further developed to represent typical HF transmitter and receiver antenna properties and to calculate the electric field strength, time delay, and polarization rotation at the receiver location.

PDF

___

[1] Budden KG. The Propagation of Radio Waves. Cambridge, UK: Cambridge University Press, 1988.
[2] Chapman S. The atmospheric height distribution of band absorbed solar radiation. P Phys Soc 1939; 51: 93-109.
[3] Shkarofski P. Generalized Appleton-Hartree equation for any degree of ionization and application to the ionosphere. P IRE 1961; 49: 1857-1871.
[4] Bilitza D. International Reference Ionosphere 2000. Radio Sci 2001; 36: 261-273.
[5] Bilitza D, Reinisch BW. International Reference Ionosphere 2007: improvements and new parameters. Adv Space Res 2008; 42: 599-609.
[6] Bilitza D, Altadill D, Zhang Y, Mertens C, Truhlik V, Richards P, McKinnel L, Reinisch B. The International Reference Ionosphere 2012 a model of international collaboration. J Space Weather Spac 2014; 4: 1-12.
[7] Gulyaeva TL, Huang X, Reinisch BW. Plasmaspheric extension of topside electron density profiles. Adv Space Res 2002; 29: 825-831.
[8] Gulyaeva TL, Bilitza D. Towards ISO Standard Earth Ionosphere and Plasmasphere Model. Hauppauge, NY, USA: Nova, 2012.
[9] Arikan F, Sezen U, Cilibas O, Gulyaeva TL. Online, automatic, ionospheric maps: IRI-PLAS-MAP. Adv Space Res 2012; 55: 2106-2113.
[10] Gulyaeva TL, Arikan F, Hernandez-Pajares M, Stanislawska I. GIM-TEC adaptive ionospheric weather assessment and forecast system. J Atmos Sol-Terr Phy 2013; 102: 329-340.
[11] Haselgrove J. Ray theory and a new method of ray tracing. P Phys Soc Lond 1955; 23: 355-364.
[12] Haselgrove CB, Haselgrove J. Twisted ray paths in the ionosphere. P Phys Soc Lond 1960; 75: 357-363.
[13] Jones RM. A Three-Dimensional Ray-Tracing Computer Program. Washington, DC, USA: Government Printing Office, 1966.
[14] Jones RM, Stephenson JJ. A Versatile Three-Dimensional Ray-Tracing Computer Program for Radio Waves in the Ionosphere. Washington, DC, USA: US Department of Communications, 1975.
[15] Varrier NR. Ray tracing analysis for the mid-latitude SuperDARN HF radar at Blackstone incorporating the IRI- 2007 model. MSc, Virginia Tech, VA, USA, 2010.
[16] Huang X. Reinisch BW. Real-time HF ray tracing through a tilted ionosphere. Radio Sci 2006; 41: 1-8.
[17] Li W, Su D, Yan Z. Three-dimensional ray tracing in International Reference Ionosphere. In: IEEE 2010 International Symposium on Antennas and Propagation and EM Theory; 29 November2 December 2010; Guangzhou, China. New York, NY, USA: IEEE, pp. 584-587.
[18] Smith MS. Phase memory in W.K.B. and phase integral solutions of ionospheric propagation problems. P Phys Soc Lond A 1975; 346: 59-79.
[19] Muldrew DB. An ionospheric ray-tracing technique and its application to a problem in a long-distance radio propagation. IRE T Antenn Propag 1959; 7: 393-396.
[20] Norman RJ, Cannon PS. An evaluation of a new two-dimensional analytic ionosphere ray tracing technique: segmented method for analytic ray tracing (SMART). Radio Sci 1999; 34: 489-499.
[21] Norman RJ, Le Marshall J, Carter BA, Wang CS. A new pseudo three-dimensional segment method analytical ray tracing (3-D SMART) technique. IEEE T Antenn Propag 2012; 60: 5818-5824.
[22] Yaxin Y, Niu J, Simpson JJ. A 3D global Earth-ionosphere FDTD model including an anisotropic magnetized plasma ionosphere. IEEE T Antenn Propag 2012; 60: 3246-3256.
[23] Isaakidis SA, Xenos TD. Wave propagation and reflection in the ionosphere: an alternative approach for ray path calculations. Prog Electromagn Res 2004; 45: 201-215.
[24] Azzarone A, Bianchi C, Pezzopane M, Pietrella M, Scotto C, Settimi A. IONORT: A Windows software tool to calculate the HF ray tracing in the ionosphere. Comput Geosci 2012; 42: 57-63.
[25] Settimi A, Pietrella M, Pezzopane M, Bianchi C. The IONORT-ISP-WC system: inclusion of an electron collision frequency model for the D-layer. Adv Space Res 2015; 55: 2114-2123.
[26] Yan Z, Wang G, Tian G, Li W, Su D, Rahman T. The HF channel EM parameters estimation under a complex environment using the modified IRI and IGRF model. IEEE T Antenn Propag 2011; 59: 1178-1783.
[27] Sayin I, Akdogan KE, Arikan F. Investigation of optimum temporal update periods for regional TEC monitoring. In: IEEE 2009 International Conference on Recent Advances in Space Technologies; 1113 June 2009; ˙Istanbul, Turkey. New York, NY, USA: IEEE. pp. 675-679.
[28] Sayin I, Arikan F, Akdogan KE. Optimum temporal update periods for regional ionosphere monitoring. Radio Sci 2010; 45: 1-9.
[29] Qiao J, Rui-Ping M, Ji-Yac X. Variations of the ionospheric conductivity with different solar and geomagnetic conditions. Chinese J Geophys 2006; 49: 1121-1129.
[30] Tanyer SG, Erol CB. Broadcast analysis and prediction in the HF band. IEEE T Broadcast 1998; 44: 226-232.
[31] Finlay C, Maus S, Beggan CD, Bondar TN, Chambodut A, Chernova TA, Chulliat A, Golovkov VP, Hamilton B, Hamoudi M et al. International geomagnetic reference field: the eleventh generation. Geophys J Int 2010; 183: 1216-1230.
[32] Sayin I, Arikan F, Arikan O. Synthetic TEC mapping with ordinary and universal Kriging. In: IEEE 2007 International Conference on Recent Advances in Space Technologies; 1416 June 2007; ˙Istanbul, Turkey. New York, NY, USA: IEEE. pp. 39-43.
[33] Sayin I, Arikan F, Arikan O. Regional TEC mapping with random field priors and kriging. Radio Sci 2008; 43: 1-14.
[34] Erdem E, Arikan F, Deviren MN, Cor I. A model based ray tracing algorithm for anisotropic and inhomogeneous Ionosphere with GIM-TEC assimilation. In: IEEE 2015 International Conference on Recent Advances in Space Technologies; 1619 June 2015; ˙Istanbul, Turkey. New York, NY, USA: IEEE, pp. 477-481.
[35] Wilson D, Mannucci AJ, Edwards CD, Roth T. Global ionospheric maps using global network of GPS receivers. In: International Beacon Satellite Symposium; 612 July 1992; Cambridge, MA, USA. Trieste, Italy: ICTP.
[36] Arikan F, Arikan O, Erol CB. Regularized estimation of TEC from GPS data for certain midlatitude stations and comparison with the IRI model. Adv Space Res 2007; 39: 867-874.