Application of radio tomographic imaging to HF oblique incidence ray tracing

Radio tomography is a technique for generating images of the spatial structure of ionospheric electron density over a wide area. This paper assesses the potential use of radio tomography in HF oblique propagation and ray tracing applications. Synthetic ionograms produced by ray tracing through tomographic images and ionospheric models have been compared with experimental oblique ionograms from six paths lying close to the image plane in the United Kingdom. In particular, study has been made of the effects of various types of input information used to constrain the vertical electron density structure in the tomographic reconstructions. It was found that use of a fine height resolution (5 km) and incorporation of information from one vertical ionosonde in the reconstruction process makes significant improvements to the overall reliability of the tomographic image. As expected, E layer propagation is better defined using a climatological model than by tomography. However, in comparison with three ionospheric models, use of tomographic images can significantly reduce the RMS error in the determination of the F2 layer maximum usable frequency

Tomographic imaging of ionospheric plasma bubbles based on GNSS and radio occultation measurements

Total electron content measurements given by the global navigation satellite system (GNSS) have successfully presented results to capture the signatures of equatorial plasma bubbles. In contrast, the correct reproduction of plasma depletions at electron density level is still a relevant challenge for ionospheric tomographic imaging. In this regard, this work shows the first results of a new tomographic reconstruction technique based on GNSS and radio-occultation data to map the vertical and horizontal distributions of ionospheric plasma bubbles in one of the most challenging conditions of the equatorial region. Twenty-three days from 2013 and 2014 with clear evidence of plasma bubble structures propagating through the Brazilian region were analyzed and compared with simultaneous observations of all-sky images in the 630.0 nm emission line of the atomic oxygen. The mean rate of success of the tomographic method was 37.1%, being more efficient near the magnetic equator, where the dimensions of the structures are larger. Despite some shortcomings of the reconstruction technique, mainly associated with ionospheric scintillations and the weak geometry of the ground-based GNSS receivers, both vertical and horizontal distributions were mapped over more than 30° in latitude, and have been detected in instances where the meteorological conditions disrupted the possibility of analyzing the OI 630 nm emissions. Therefore, the results revealed the proposed tomographic reconstruction as an efficient tool for mapping characteristics of the plasma bubble structures, which may have a special interest in Space Weather, Spatial Geodesy, and Telecommunications.Peer ReviewedPostprint (published version

Radio tomographic imaging of sporadic-E layers during SEEK-2

During the SEEK-2 Rocket Campaign in August 2002, a Dual Band Beacon (DBB) transmitting to Ground Receivers provided unique data on E-Region electron densities. Information from two rocket beacons and four ground receivers yielded multiple samples of E-region horizontal and vertical variations. The radio beacon measurements were made at four sites (Uchinoura, Tarumizu, Tanegashima, Takazaki) in Japan for two rockets (S310-31 and S310-32) launched by the Institute of Space and Aeronautical Science (ISAS). Analysis was completed for four sets of beacon data to provide electron density images of sporadic-E layers. Signals from the two-frequency beacons on the SEEK-2 rockets were processed to yield total electron content (TEC) data that was converted into electron density measurements. Wide variations in layer structures were detected. These included horizontal sporadic-E variations, vertical profiles of double, single, and weak layers. The radio beacon measurements were shown to be in agreement with the in-situ SEEK-2 sensors. The first tomographic image of a sporadic-E layer was produced from the data. The rocket beacon technique was shown to be an excellent tool to study sporadic-E layers because absolute TEC accuracy of 0.01 TEC Units can be easily obtained and, with proper receiver placement, electron density images can be produced using computerized ionospheric tomography with better than 1km horizontal and vertical resolution

Three-dimensional, high resolution, computerized ionospheric tomographic imaging and computational modeling of an artificial ionospheric cavity

Tomographic techniques use line integral measurements to reconstruct local values of the measured parameter. These techniques have been applied to the ionosphere by using radio transmissions to measure the integral of electron density between a satellite and a chain of ground-based receiving stations. The resultant reconstructions form a two-dimensional map of the electron density in the plane of the satellite/receiver chain. Insufficient quantity of receivers and not having a complete range of accessible look angles with the available receivers are responsible for the non-uniqueness in the reconstructions. The limited look angle problem can be alleviated by making use of alternative data sources such as incoherent scatter radars ( ISR) that provide information on the vertical structure of the electron density. The non-optimal receiver placement problem can be alleviated through the use of sophisticated reconstruction algorithms. The computerized ionospheric tomography (CIT) technique has recently been used to image the artificially modified ionosphere above the Arecibo Observatory (AO) HF heating facility. A total of nine radio beacon receivers forming a three by three matrix were deployed across the entire island of Puerto Rico. The arrangement maximizes the likelihood that several of the receiver-to-satellite lines of sight would intersect the heated region of the ionosphere. Several satellite passes associated with such an intersection occurred, and the resultant tomographic inversions indicate the existence of an electron density cavity approximately 45 kilometers in latitudinal extent. The reduction of electron density in this cavity was typically on the order of 20%. The experimental observations were supported by theoretical work using the open-source SAMI2 ionospheric model. Methods were developed to model both the ohmic heating of the electrons and the suprathermal heating caused by nonlinear effects. Modeled ohmic heating values of 941 K/s match the observed heated temperature profiles. Modeled suprathermal electrons effects on the vibrational temperature of N2 indicate a vibrational/translation temperature differential of 308°K

Three-dimensional ionospheric tomography algebraic reconstruction technique

An improved algebraic reconstruction technique (IART) is presented for the tomographic reconstruction of ionospheric electron density (IED). This method applies the total electron content (TEC) measurements to invert the spatial distribution of the IED from a set of apriori IED distributions. In this new method, a data-driven adjustment of the relaxation parameter is performed to improve the computation efficiency and image quality of the classical algebraic reconstruction technique (ART). In addition, the new algorithm is also combined with ionospheric space discretization technique to simplify the inversion of IED, and it applies CHAMP occultation data to improve the vertical resolution. A numerical simulation experiment is carried out to validate the reliability of the new method. It is then applied to the inversion of IED from real GPS data. Inverted results show that the IART algorithm has better accuracy and efficiency than the conventional ART algorithm. The reliability of the IART algorithm is also validated by ionosonde data recorded at Wuhan station

A hybrid reconstruction algorithm for 3-D ionospheric tomography

In this paper, a hybrid reconstruction algorithm (HRA) is presented to solve the ill-posed inverse problem associated with 3-D ionospheric stochastic tomography. In this new method, the ionospheric electron density (IED) can be inverted by using two steps. First, a truncated singular value decomposition (TSVD) method, whose value is independent on any initial estimation, is used to resolve the ill-posed problem of the tomography system. Second, taking into account the "approximation" of its solution, an iterative improvement process of the solution is then implemented by utilizing the conventional algebraic reconstruction algorithm (ART). The HRA, therefore, offers a more reasonable approach to choose an initial approximate for the ART and to improve the quality of the final reconstructed image. A simulated experiment demonstrates that the HRA method is superior to the TSVD or the ART alone for the tomographic inversion of IED. Finally, the HRA is used to perform GPS-based tomographic reconstruction of the IED at mid- and low-latitude regions

Bayesian approach to ionospheric imaging with Gaussian Markov random field priors

Ionosfääri on noin 60–1000 kilometrin korkeudella sijaitseva ilmakehän kerros, jossa kaasuatomien ja -molekyylien elektroneja on päässyt irtoamaan auringon säteilyn ja auringosta peräisin olevien nopeiden hiukkasten vaikutuksesta. Näin syntyneillä ioneilla ja vapailla elektroneilla on sähkö- ja magneettikenttien kanssa vuorovaikuttava sähkövaraus. Ionosfäärillä on siksi merkittävä rooli radioliikenteessä. Se voi mahdollistaa horisontin yli tapahtuvat pitkät radiolähetykset heijastamalla lähetetyn sähkömagneettisen signaalin takaisin maata kohti. Toisaalta ionosfääri vaikuttaa myös sen läpäiseviin korkeampitaajuuksisiin signaaleihin. Esimerkiksi satelliittipaikannuksessa ionosfäärin vaikutus on parhaassakin tapauksessa otettava huomioon, mutta huonoimmassa se voi estää paikannuksen täysin. Näkyvin ja tunnetuin ionosfääriin liittyvä ilmiö lienee revontulet. Yksi keskeisistä suureista ionosfäärin tutkimuksessa on vapaiden elektronien määrä kuutiometrin tilavuudessa. Käytännössä elektronitiheyden mittaaminen on mahdollista mm. tutkilla, kuten Norjan, Suomen ja Ruotsin alueilla sijaitsevalla EISCAT-tutkajärjestelmällä, sekä raketti- tai satelliittimittauksilla. Mittaukset voivat olla hyvinkin tarkkoja, mutta tietoa saadaan ainoastaan tutkakeilan suunnassa tai mittalaitteen läheisyydestä. Näillä menetelmillä ionosfäärin tutkiminen laajemmalla alueella on siten vaikeaa ja kallista. Olemassa olevat paikannussatelliitit ja vastaanotinverkot mahdollistavat ionosfäärin elektronitiheyden mittaamisen alueellisessa, ja jopa globaalissa mittakaavassa, ensisijaisen käyttötarkoituksensa sivutuotteena. Satelliittimittausten ajallinen ja paikallinen kattavuus on hyvä, ja kaiken aikaa kasvava, mutta esimerkiksi tarkkoihin tutkamittauksiin verrattuna yksittäisten mittausten tuottama informaatio on huomattavasti vähäisempää. Tässä väitöstyössä kehitettiin tietokoneohjelmisto ionosfäärin elektronitiheyden kolmiulotteiseen kuvantamiseen. Menetelmä perustuu matemaattisten käänteisongelmien teoriaan ja muistuttaa lääketieteessä käytettyjä viipalekuvausmenetelmiä. Satelliittimittausten puutteellisesta informaatiosta johtuen työssä on keskitytty etenkin siihen, miten ratkaisun löytymistä voidaan auttaa tilastollisesti esitetyllä fysikaalisella ennakkotiedolla. Erityisesti työssä sovellettiin gaussisiin Markovin satunnaiskenttiin perustuvaa uutta korrelaatiopriori-menetelmää. Menetelmä vähentää merkittävästi tietokonelaskennassa käytettävän muistin tarvetta, mikä lyhentää laskenta-aikaa ja mahdollistaa korkeamman kuvantamisresoluution.Ionosphere is the partly ionised layer of Earth's atmosphere caused by solar radiation and particle precipitation. The ionisation can start from 60 km and extend up to 1000 km altitude. Often the interest in ionosphere is in the quantity and distribution of the free electrons. The electron density is related to the ionospheric refractive index and thus sufficiently high densities affect the electromagnetic waves propagating in the ionised medium. This is the reason for HF radio signals being able to reflect from the ionosphere allowing broadcast over the horizon, but also an error source in satellite positioning systems. The ionospheric electron density can be studied e.g. with specific radars and satellite in situ measurements. These instruments can provide very precise observations, however, typically only in the vicinity of the instrument. To make observations in regional and global scales, due to the volume of the domain and price of the aforementioned instruments, indirect satellite measurements and imaging methods are required. Mathematically ionospheric imaging suffers from two main complications. First, due to very sparse and limited measurement geometry between satellites and receivers, it is an ill-posed inverse problem. The measurements do not have enough information to reconstruct the electron density and thus additional information is required in some form. Second, to obtain sufficient resolution, the resulting numerical model can become computationally infeasible. In this thesis, the Bayesian statistical background for the ionospheric imaging is presented. The Bayesian approach provides a natural way to account for different sources of information with corresponding uncertainties and to update the estimated ionospheric state as new information becomes available. Most importantly, the Gaussian Markov Random Field (GMRF) priors are introduced for the application of ionospheric imaging. The GMRF approach makes the Bayesian approach computationally feasible by sparse prior precision matrices. The Bayesian method is indeed practicable and many of the widely used methods in ionospheric imaging revert back to the Bayesian approach. Unfortunately, the approach cannot escape the inherent lack of information provided by the measurement set-up, and similarly to other approaches, it is highly dependent on the additional subjective information required to solve the problem. It is here shown that the use of GMRF provides a genuine improvement for the task as this subjective information can be understood and described probabilistically in a meaningful and physically interpretative way while keeping the computational costs low

Electron Densities in the Lower Thermosphere From GUVI 135.6 nm Tomographic Inversions in Support of SpreadFEx

The SpreadFEx campaign was conducted with the goal of investigating potential neutral atmospheric dynamics influences in seeding plasma instabilities and bubbles extending to higher altitudes from September to November 2005, with primary measurements in Brazil. In this paper, we present the results of space-based UV and ground-based optical observations in support of this campaign. Specifically, we present multi-dimensional electron density images obtained tomographically from the 135.6 nm emissions measured by the GUVI instrument aboard the TIMED satellite that result from radiative recombination of O+ and compare those with the corresponding 630.0 nm OI images recorded in the Brazilian sector. The GUVI results provide altitude vs. longitude information on depleted regions in the ionospheric plasma density that are complementary to the single-height latitude-longitude images obtained with the airglow imager