Settimi, Alessandro

This paper describes the main updates characterizing the new version of IONORT (IONOsperic Ray Tracing), a software tool developed at Istituto Nazionale di Geofisica e Vulcanologia to determine both the path of a high frequency (HF) radio wave propagating in the ionospheric medium, and the group time delay of the wave itself along the path. One of the main changes concerns the replacement of a regional three-dimensional electron density matrix, which was previously taken as input to represent the ionosphere, with a global one. Therefore, it is now possible to carry out different ray tracings from whatever point of the Earth’s surface, simply by selecting suitable loop cycles thanks to the new ray tracing graphical user interface (GUI). At the same time, thanks to a homing GUI, it is also possible to generate synthetic oblique ionograms for whatever radio link chosen by the user. Both ray tracing and homing GUIs will be described in detail providing at the same time some practical examples of their use for different regions. IONORT software finds practical application in the planning of HF radio links, exploiting the sky wave, through an accurate and thorough knowledge of the ionospheric medium. HF radio waves users, including broadcasting and civil aviation, would benefit from the use of the IONORT software (version 2023.10).

For the phase refraction index of high frequency (HF) waves in the ionospheric medium exists a well-established theory. However, under the Quasi-Longitudinal (QL) conditions, scientific literature presents various formulas that are not equivalent and that, in some cases, give rise to wrong results. In the present study, further consequences of Booker's rule are discussed, illustrating the validity ranges of the above-mentioned approximate formulas; and the different regimes for applying such QL formulas are described, along with the consequences in simulating the ionospheric HF ray-tracing, oblique and vertical sounding, and absorption.

This paper wants to highlight how the availability of measurements autoscaled at some reference ionospheric stations, and their assimilation by ionospheric models, was of crucial importance in determining, during the solar eclipse conditions occurred on 20 March 2015, a reliable representation of the ionosphere. Even though the solar eclipse falls in the recovery phase of the St. Patrick geomagnetic storm started on 17 March 2015, its influence on the ionospheric plasma seems undeniable. The reference ionospheric stations considered here are those of Rome (41°.8’ N, 12°.5’ E), and Gibilmanna (37°.9’ N,14°.0’ E), Italy. Specifically, in a time interval including that of the eclipse, the electron density profiles autoscaled by the Automatic Real-Time Ionogram Scaler with True-height (ARTIST) system at San Vito (40°.6′ N, 17°.8′ E), Italy, which are here considered as the truth profiles, were compared with both the electron density profiles calculated by the IRI-SIRMUP-Profiles (ISP) model, after assimilating data recorded at Rome and Gibilmanna, and the electron density profiles provided by the IRI-CCIR model. The ISP and IRI-CCIR performances were then evaluated in terms of the root mean square errors made on the whole electron density profiles. The three-dimensional (3-D) electron density mappings of the ionosphere provided by ISP and IRI-CCIR models were also considered as the ionospheric environment by the ray tracing software tool IONORT to calculate quasi-vertical synthesized ionograms over the short radio link San Vito – Brindisi (40°.4′ N, 17°.6′ E), Italy. The corresponding synthesized values of foF2 and fxF2, obtained by IONORT-ISP and IONORT-(IRI-CCIR) procedures, were compared with those autoscaled by ARTIST from the vertical ionograms recorded at the truth site of San Vito. Some examples of IONORT-ISP and IONORT-(IRI-CCIR)synthesized ionograms are shown and discussed. Finally, comparisons in terms of foF2 deduced by long-term prediction and nowcasting maps are also shown. The results achieved in this work demonstrate how the assimilation of autoscaled data into the ionospheric models turned out to be valuable in providing a better representation of the ionospheric electron density under very unusual conditions.

DOME C is located on the East Antarctic Plateau at an altitude of 3233 m above sea level and is the site of the Italian–French base, i.e., the Concordia Station ( 123∘20′ E, 75∘06′ S). It has become an important location for several scientific research studies, including astrophysics, geophysics, glaciology, and climatology, due to the perceived long-term stability and thickness of ice at this location [1], [5], [18], [24]. During the site selection and follow up of the EPICA Dome C coring project, which ultimately provided more than 800 000 years of palaeoclimatic series distributed along 3270 m [6], numerous radio-echo sounding (RES) surveys were undertaken to improve core positioning and subsequently to better leverage data and logistical infrastructures from the coring effort. These studies were conducted at very different scales over the Dome C region and revealed important information about the bedrock physiography and its physical conditions [2]– [4], [8], [9], [19], [20], [30], [31]. This study builds upon these efforts using data from two recently acquired ground-based surveys collected during 2009 and 2011 in a very small area (2.5 × 2.0 km) in the immediate proximity of the EPICA drilling site, to test advances in the acquisition and digitization technology and resolve the basal environment to unprecedented detail.

The ionosphere affects the electromagnetic wave propagation and then its study is important for Earth-Earth, satellite-Earth, and satellite-satellite communication purposes. Diffractive and refractive processes due to irregular electron density structures cause signal fluctuations that can disrupt satellite-ground communications and represent a hazard for navigation systems. The study and the real-time monitoring of the ionosphere are important for Space Weather purposes. The ionospheric vertical sounding is described, together with the automatic scaling of the ionograms.

The present book chapter discusses the stimulated emission, in strong coupling regime, of an atom embedded inside a one dimensional (1D) Photonic Band Gap (PBG) cavity which is pumped by two counter-propagating laser beams. Quantum electrodynamics is applied to model the atom-field interaction, by considering the atom as a two level system, the e.m. field as a superposition of normal modes, the coupling in dipole approximation, and the equations of motion in Wigner-Weisskopf and rotating wave approximations. In addition, the Quasi Normal Mode (QNM) approach for an open cavity is adopted, interpreting the local density of states (LDOS) as the local density of probability to excite one QNM of the cavity; and therefore rendering this LDOS dependent on the phase difference of the two laser beams. In this book chapter we demonstrate that the strong coupling regime occurs at high values of the LDOS. In accordance with the results of the literature, the emission probability of the atom decays with an oscillatory behaviour, so that the atomic emission spectrum exhibits two peaks (Rabi splitting). The novelty of this book chapter is that the phase difference of the two laser beams can produce a coherent control of both the oscillations for the atomic emission probability and, as a consequence, of the Rabi splitting in the emission spectrum. Possible criteria to design active delay lines are finally discussed.

The IONORT-ISP system (IONOspheric Ray-Tracing – IRI-SIRMUP-PROFILES) was recently developed and tested by comparing the measured oblique ionograms over the radio link between Rome (41.89ºN, 12.48ºE), Italy, and Chania (35.51ºN, 24.02ºE), Greece, with the IONORT-ISP simulated oblique ionograms (Settimi et al., 2013). The present paper describes an upgrade of the system to include: a) electron-neutral collision have been included by using a collision frequency model that consists of a double exponential profile; b) the ISP three dimensional (3-D) model of electron density profile grid has been extended down to the altitude of the D-layer; c) the resolution in latitude and longitude of the ISP 3-D model of electron density profile grid has been increased from 2°x2° to 1°x1°. Based on these updates, a new software tool called IONORT-ISP-WC (WC means with collisions) was developed, and a database of 33 IONORT-ISP-WC synthesized oblique ionograms calculated for single (1-hop paths) and multiple (3-hop paths) ionospheric reflections. The IONORT-ISP-WC simulated oblique ionograms were compared with the IONORT-IRI-WC synthesized oblique ionograms, generated by applying IONORT in conjunction with the International Reference Ionosphere (IRI) 3-D electron density grid, and the observed oblique ionograms over the aforementioned radio link. The results obtained show that (1) during daytime, for the lower ionospheric layers, the traces of the synthesized ionograms are cut away at low frequencies because of HF absorption; (2) during night-time, for the higher ionospheric layers, the traces of the simulated ionograms at low frequencies are not cut off (very little HF absorption); (3) the IONORT-ISP-WC MUF values are more accurate than the IONORT-IRI-WC MUF values.

When applying the ray tracing in ionospheric propagation, the electron density modelling is the main input of the algorithm, since phase refractive index strongly depends on it. Also the magnetic field and frequency collision modelling have their importance, the former as responsible for the azimuth angle deviation of the vertical plane containing the radio wave, the latter for the evaluation of the absorption of the wave. Anyway, the electron density distribution is strongly dominant when one wants to evaluate the group delay time characterizing the ionospheric propagation. From the group delay time, azimuth and elevation angles it is possible to determine the point of arrival of the radio wave when it reaches the Earth surface. Moreover, the procedure to establish the target (T) position is one of the essential steps in the Over The Horizon Radar (OTHR) techniques which require the correct knowledge of the electron density distribution. The group delay time generally gives rough information of the ground range, which depends on the exact path of the radio wave in the ionosphere. This paper focuses on the lead role that is played by the variation of the electron density grid into the ray tracing algorithm, which is correlated to the change of the electron content along the ionospheric ray path, for obtaining a ray tracing as much reliable as possible. In many cases of practical interest, the group delay time depends on the geometric length and the electron content of the ray path. The issue is faced theoretically, and a simple analytical relation, between the variation of the electron content along the path and the difference in time between the group delays, calculated and measured, both in the ionosphere and in the vacuum, is obtained and discussed. An example of how an oblique radio link can be improved by varying the electron density grid is also shown and discussed.

We present an algorithm for the identification of trace characteristics of oblique ionograms allowing determination of the Maximum Usable Frequency (MUF) for communication between the transmitter and receiver. The algorithm automatically detects and rejects poor quality ionograms. We performed an exploratory test of the algorithm using data from a campaign of oblique soundings between Rome, Italy (41.90 N, 12.48 E) and Chania, Greece (35.51 N, 24.01 E) and also between Kalkarindji, Australia (17.43 S, 130.81 E) and Culgoora, Australia (30.30 S, 149.55 E). The success of these tests demonstrates the applicability of the method to ionograms recorded by different ionosondes in various helio and geophysical conditions.

This work will lead to ray theory and ray tracing formulation. To deal with this problem the theory of classical geometrical optics is presented, and applications to ionospheric propagation will be described. This provides useful theoretical basis for scientists involved in research on radio propagation in inhomogeneous anisotropic media, especially in a magneto-plasma. Application in high frequencies (HF) radio propagation, radio communication, over-the-horizon-radar (OTHR) coordinate registration and related homing techniques for direction finding of HF wave, all rely on ray tracing computational algorithm. In this theory the formulation of the canonical, or Hamiltonian, equations related to the ray, which allow calculating the wave direction of propagation in a continuous, inhomogeneous and anisotropic medium with minor gradient, will be dealt. At least six Hamilton’s equations will be written both in Cartesian and spherical coordinates in the simplest way. These will be achieved by introducing the refractive surface index equations and the ray surface equations in an appropriate free-dimensional space. By the combination of these equations even the Fermat’s principle will be derived to give more generality to the formulation of ray theory. It will be shown that the canonical equations are dependent on a constant quantity H and the Cartesian coordinates and components of wave vector along the ray path. These quantities respectively indicated as ri(τ), pi(τ) are dependent on the parameter τ, that must increase monotonically along the path. Effectively, the procedure described above is the ray tracing formulation. In ray tracing computational techniques, the most convenient Hamiltonian describing the medium can be adopted, and the simplest way to choose properly H will be discussed. Finally, a system of equations, which can be numerically solved, is generated.