Session - Modelling the Earth's ionosphere and solutions to counter ionospheric threats to GNSS applications
Marcio Aquino, Cathryn Mitchell, Giorgiana De Franceschi
The Earth's ionosphere is a complex system that is driven by many different factors such as solar radiation, electric and
magnetic fields, and neutral atmosphere dynamics. New models and data to realise the state of the Earth's ionosphere are
of interest to radio system users whose signals are affected by ionospheric propagation, in particular navigation and
communications systems operating below 3 GHz. Global Navigation Satellites Systems (GNSS), which have become essential
in support of a growing number of activities now embedded in modern society, are especially vulnerable to signal propagation
through the ionosphere. In this context this session addresses models and solutions to mitigate ionospheric threats to GNSS
and related applications. The TRANSMIT project, an FP7 funded Marie Curie Initial Training Network,
has focused on the understanding and development of new models that can be tested for their usefulness in
addressing radio system specification with a view to support the development of concepts and operational tools
that could contribute to a service to assist European users in countering GNSS vulnerability to ionospheric phenomena.
Results of this project are expected to be showcased in this session, as well as any other initiatives in this area.
The session welcomes work in the areas of ionospheric tomography and imaging, radio occultation, scintillation and
interference resilient receiver tracking models and implementations, real time positioning algorithms (e.g. for
Precise Point Positioning) to mitigate ionospheric threats affecting legacy and new GNSS signals, scintillation and
TEC prediction models and operational tools, as well as any other related topics. The main aim of the session is to
stimulate discussion and encourage new collaborative work.
Invited Talks and First Class Posters
Tuesday November 18, 09:00 - 10:50
1 |
Oral |
9:00 am |
An Improved Ionospheric Modeling
Technique Using GPS and Empirical-Orthogonal-Function Fits |
|
|
|
Zhang, W1; Langley, R B1; Komjathy, A2; Banville, S3 |
|
|
|
1University
of New Brunswick; 2Jet Propulsion Laboratory & University of New Brunswick; 3Natural Resources Canada
& University of New Brunswick |
|
|
|
The dispersion of the ionosphere
pertaining to Global Positioning System (GPS) signals allows us to gain
information on its morphology in terms of total electron content (TEC) or
electron density along the individual signal ray paths, and the generation of
2D or 3D ionospheric models. Several steps are typically involved in
ionospheric modeling, such as data preprocessing to obtain slant TECs
(STECs), estimating receiver and satellite biases and eliminating them from
the STECs, and fitting the STEC (or converted vertical TEC (VTEC)) values
with a mathematical model. Spherical harmonic functions (SHF), wavelet bases,
and polynomial functions are typically employed in the mathematical model to
depict the horizontal variability of the ionosphere for global and regional
ionospheric modeling. However, while using functions of higher order is a
usual practice for more precise results in terms of lower post-fit residuals,
it may induce model instability, estimated states being unobservable, or a
potential heavy computing burden. This is especially the case for disturbed
ionospheric conditions with irregular “behavior” of the ionosphere. To
overcome any downsides of the function-based models, the Kriging algorithm
was introduced for ionospheric modeling. Its benefit of assigning weights
according to a data-driven weighting function rather than any arbitrary
functions was demonstrated by some related works on the algorithm design for
the Wide Area Augmentation System (WAAS). However, the accuracy of the model
may still be limited by the accuracy of inter-frequency hardware biases
(IHBs), which are typically estimated under some assumptions, such as the
thin shell assumption to convert STECs to VTECs and/or some other
mathematical representation of the ionosphere. In this paper, attempting to
make the estimated parameters more meaningful and effective, we propose an
ionospheric modeling technique that uses data-driven empirical orthogonal
functions (EOFs) to replace arbitrary functions to match the horizontal
variability of the ionosphere and estimate IHBs with the EOF-fit
representation of the ionosphere. We then assess the newly proposed model by
analyzing the post-fit residuals of the estimation process, as well as the
repeatability of estimates of IHBs in both 2D and 3D modeling scenarios. For
the 2D scenario, we use the thin shell mapping function, and for the 3D
scenario, we use EOFs to depict the vertical variability, so that all the
three dimensions in this case are represented by EOFs. |
2 |
Oral |
9:15 am |
Ionospheric Data Assimilation |
|
|
|
Da Dalt, F1; Panicciari, T1; Benton, C1; Smith, N1; Mitchell, C1 |
|
|
|
1University
of Bath |
|
|
|
One of the manifestations of
space weather is the high variability of the ionosphere. Signals propagating
through this part of the atmosphere are adversely affected by its presence
and dynamics, particularly the Global Navigation Satellite Systems (GNSS).
The need for a better understanding of the characteristics and behaviour of
this medium resulted in the realization of MIDAS (Multi-Instrument Data
Analysis System). Developed at the University of Bath, MIDAS is a software
package that provides three dimensional reconstructions of ionospheric plasma
density from GPS (Global Positioning System) measurements collected by
ground-based receivers. These measurements, known as Total Electron content
(TEC) and expressed in TECU (TEC units [1016 particles/m2]), are integrated
electron density values along a specific signal path. MIDAS algorithms gather
and utilise various TEC measurements to compute reconstructions of the
ionosphere through tomographic inversion. Unfortunately, the
satellite-receiver path arrangements do not provide an ideal scan geometry
for the reconstructions. Further, the ground receivers do not cover the Earth
surface homogenously. This translates to a poor data coverage and lack of
information that prevents the inversion from having a unique solution. To
overcome this problem, the inversion in MIDAS is supported by regularization
methods and external information provided by empirical models. The presented
project aims to implement a physics-based ionospheric model into MIDAS by
means of a Data Assimilation (DA) scheme. The objective is twofold: obtain
higher accuracy ionospheric reconstructions using physics-based support and
be used as a forecasting tool. To reach the aforementioned goals, A New
Ionospheric Model (ANIMo) was developed and validated at the University of
Bath. ANIMo is a global model – it
solves the continuity equation for monoatomic oxygen O+ in a three
dimensional grid of altitude, latitude and longitude. It performs reconstructions
in an altitude range from 80 to 600 km and is designed to be used in
mid-latitude regions. ANIMo assumes that the sum of the major ions densities
(O+, NO+ and O2+) is equal to the electron density. It was developed to
conform to specific requirements necessary for its final implementation in
MIDAS. The first-principles nature of ANIMo is very relevant not only for
prediction but also to enable better control over the whole assimilative
procedure. This permits the use of the model to test the importance of
different input parameters or to simulate particular ionospheric conditions.
For similar reasons, the model was intentionally developed to exhibit
simplicity and robustness. The application of ANIMo in MIDAS is mainly
performed through a DA approach known as 4D-Var. Largely used in meteorology,
this type of assimilative scheme typically has an iterative nature. Here,
ANIMo provides the background (or a priori) estimates that, combined with the
GPS observations, give the analysis – an image of the true state of the
ionosphere at a given point in space and time. The analysis can be then used
to initialize ANIMo to compute a new analysis at the next time step or to
perform routine forecasting. The model contribution is weighted by the
background covariance matrix. The latter is constructed in this scheme by
considering time-dependent horizontal correlation distances of electron
density. ANIMo is also used to provide a physical constraint in the
reconstructions of vertical profiles, which represents a critical phase when
imaging the ionosphere through tomography. Promising preliminary results
suggest that the presented assimilative scheme is not only able to provide
accurate ionospheric images but also the opportunity for prediction.
Furthermore, this setup could be used as a benchmark to test the importance
of external forcing factors when modelling the upper atmosphere; therefore
useful for the whole ionospheric forecasting community. |
3 |
Oral |
9:30 am |
Geomagnetic Control of
Equatorial Plasma Bubble Activity Modeled by the TIEGCM with Kp |
|
|
|
Carter, B1; Francis, M2; Retterer, J1; Yizengaw, E1; Groves, K1; Caton, R3; McNamara, L3; Bridgewood, C1; Terkildsen, M2; Norman, R4; Zhang, K4 |
|
|
|
1Institute
for Scientific Research, Boston College; 2IPS Radio and Space Services, Australian Bureau of Meteorology;
3Air Force
Research Laboratory, Kirtland AFB, New Mexico; 4SPACE Research Centre, RMIT University, Melbourne |
|
|
|
Describing the day-to-day
variability of Equatorial Plasma Bubble (EPB) occurrence has remained a
significant challenge over several decades. The importance in pursuing such a
task has increased significantly in recent years due to the increased use and
reliance upon Global Navigation Satellite Systems for navigation, positioning
and timing. In this study we use the Thermosphere Ionosphere Electrodynamics
General Circulation Model (TIEGCM), driven by solar (F10.7) and geomagnetic
(Kp) activity indices, to study daily variations of the linear
Rayleigh-Taylor (R-T) instability growth rate in relation to the measured
scintillation strength at five longitudinally distributed stations. For
locations characterized by generally favorable conditions for EPB growth
(i.e., within the “scintillation season” for that location) we find that the
TIEGCM is capable of identifying days when EPB development, determined from
the calculated R-T growth rate, is suppressed as a result of geomagnetic
activity. Observed and modeled upward plasma drift indicates that the
pre-reversal enhancement best scales linearly with Kp from ~ 3.5 hrs prior,
from which it is concluded that even small Kp changes cause significant
variations in EPB growth, and the likelihood of associated radio
scintillation events on any given day. |
4 |
Oral |
9:45 am |
Zonal Velocity of the Equatorial
Ionospheric Irregularities over São Paulo State, Brazil, during the last
Solar Maximum |
|
|
|
Cesaroni, C1; Spogli, L1; De Franceschi, G1; Alfonsi, L1; Romano, V1; Aquino, M2; Park, J2; Monico, J F G3 |
|
|
|
1INGV;
2University
of Nottingham 3UNESP |
|
|
|
The South American ionosphere is
characterized by the presence of the Equatorial Ionospheric Anomaly (EIA),
which results in two crests of enhanced electron density located at ±15° off
the magnetic equator. Such characterization implies a complex configuration
and dynamics of the local ionospheric plasma, especially during solar maximum
conditions. In the framework of the FP7 project CALIBRA (Countering GNSS high
Accuracy applications LImitations due to ionospheric disturbances in BRAzil)
three dual-frequency multi-constellation GNSS receivers have been installed
in the state of São Paulo, in the region of Presidente Prudente (22°07′32″S,
51°23′20″W), in order to study the dynamics of the irregularities near the
southern crest of the EIA. These receivers (PRU1, PRU2 and PRU3) are located
on the vertices of an equilateral triangle at a distance of about 300 meters
from each other to allow the study of the dynamics of the irregularities of
the dimension of the Fresnel scale for L-band signals. In particular, PRU1
and PRU3 are located along the E-W magnetic direction, while PRU1 and PRU2
along the N-S direction. Following the study carried out by Ledvina (2004),
it is possible to estimate the zonal velocity of the irregularities from the
knowledge of the zonal velocity of the scintillation pattern (vscint) and the
position (xsat) and velocity (vsat) of the satellites. In this paper, we use
a cross-correlation technique, applied to 50 Hz raw-data from the receivers
located in Presidente Prudente, to estimate vscint and the orbital parameters
of the GPS satellites included in the RINEX navigational files to evaluate
xsat and vsat. The present study has been carried out by using data acquired
during a measurement campaign from August 2013 to April 2014, in order to
show the daily variability of the zonal velocity of the irregularities. |
5 |
Oral |
10:00 am |
Current
capability of GNSS ionospheric tomography in Antarctica |
|
|
|
Kinrade, J1; Panicciari, T2; Smith, N2; Mitchell, C2 |
|
|
|
1IEEA;
2University
of Bath; 2University
of Bath |
|
|
|
Global Navigation Satellite
System (GNSS) tomography can provide three-dimensional and time dependent
reconstructions of ionospheric electron density, on global and regional
scales. Tomography is particularly useful at high latitudes, where
ionospheric observations can be sparse and unevenly distributed due to the
respective topography and remote environments of the Arctic and Antarctica.
Following expansion of ground receiver networks over the last decade through
various polar campaigns, tomography is now a proven tool for imaging plasma
morphology at Arctic latitudes. Tomographic reconstructions have been
compared and verified with independent instruments, including Incoherent
Scatter Radar (ISR) and in-situ satellites, revealing features such as the
mid-latitude trough [e.g. Mitchell et al., 1997] and moving patch
enhancements [e.g. Bust & Crowley, 2007; Spencer & Mitchell, 2007].
Tomography played a key role in profiling the ionosphere’s northern
hemisphere response during the Halloween 2003 geomagnetic storm [e.g.
Mitchell et al., 2005]. However, there has been historically poor and uneven
data coverage in Antarctica compared with the Arctic. The first tomographic
images of the Antarctic ionosphere were produced using measurements from two
ground stations [Heaton et al., 1996], but even with a distribution of
approximately 50 stations across the continent today, wide-area tomographic
reconstructions over Antarctica [e.g. Yin et al., 2010; Kinrade et al., 2012]
are limited in area extent and feature resolution compared with those in the
northern hemisphere. This work offers a brief review and current status of
ionospheric tomography in Antarctica, and using case studies highlights the
capability of the University of Bath’s Multi Instrument Data Analysis
Software (MIDAS) [Mitchell & Spencer, 2003]. Current work from the Marie
Curie TRANSMIT Initial Training Network (ITN) aims to adapt the MIDAS
algorithm to achieve multi-scale feature resolution in areas of poor data
coverage, an important step towards better detection and tracking of discrete
ionospheric structures in Antarctica.
References BUST, G. S., and G.
Crowley (2007), Tracking of polar cap ionospheric patches using data
assimilation, J. Geophys. Res., 112, A05307, doi:10.1029/2005JA011597.
HEATON, J.A.T., G.O.L. Jones and L. Kersley (1996), Toward ionospheric
tomography in Antarctica: first steps and comparison with dynasonde
observations, Antarctic Science, 8, pp. 297-302,
doi:10.1017/S0954102096000430. KINRADE, J., C. N. Mitchell, P. Yin, N. Smith,
M. J. Jarvis, D. J. Maxfield, M. C. Rose, G. S. Bust, and A. T. Weatherwax
(2012), Ionospheric scintillation over Antarctica during the storm of 5–6
April 2010, J. Geophys. Res., 117, A05304, doi:10.1029/2011JA017073.
MITCHELL, C. N., I. K. Walker, S. E. Pryse, I. Kersley, I. W. McCrea, and T.
B. Jones (1997), First complementary observations by ionospheric tomography,
the EISCAT Svalbard radar and the CUTLASS HF radar, Ann. Geophys., 16(11),
pp. 1519-1522, doi:10.1007/s00585-998-1519-2. MITCHELL, C. N. and P. S. J.
Spencer (2003), A three-dimensional time-dependent algorithm for ionospheric
imaging using GPS, Ann. Geophys., 46(4), pp. 687-696, doi:10.4401/ag-4373.
MITCHELL, C. N., L. Alfonsi, G. De Franceschi, M. Lester, V. Romano, and A.
W. Wernik (2005), GPS TEC and scintillation measurements from the polar
ionosphere during the October 2003 storm, Geophys. Res. Lett., 32(12),
L12S03. SPENCER, P. S. J. and C. N. Mitchell (2007), Imaging of fast moving
electron-density structures in the polar cap, Annals of Geophysics, 50 (3),
pp. 427-434, doi:10.4401/ag-3074. YIN P., C. N. Mitchell, L. Alfonsi, M.
Pinnock, P. Spencer, G. De Franceschi, V. Romano, P. Newell, P. Sarti, M.
Negusini, A. Capra (2009), Imaging of the Antarctic ionosphere: Experimental
results, J. Atmos. and Sol.-Terr. Phys., 71(17-18), pp. 1757-1765,
doi:10.1016/j.jastp.2009.09.014. |
1 |
First class poster |
10:20 AM |
Development and Investigation of
Applicability of Parameters Describing Performance of GPS Point Positioning
under the Presence of Ionospheric Anomalies |
|
|
|
Wezka, K |
|
|
|
Technische Universität Berlin |
|
|
|
In a real time GNSS single
receiver based processing mode, such as Precise Point Positioning (PPP),
verification and confirmation of a high trust degree towards the estimated
coordinates is a very critical issue. Especially harmful for those
positioning approaches are random and un-modelled errors which can adversely
affect the GPS signal and degrade accuracy and reliability of navigation
solution. Ionospheric perturbations
described as fast and random variability of plasma density in the ionosphere
are difficult to predict, to detect and to model. Occurrence of some strong
ionospheric disturbances can cause, inter alia, degradation and interruption of GNSS
signals or even lead to loss of signal. Under the presence of ionospheric
disturbances a quality monitoring system of PVT (positioning, velocity and
time) solutions can play a crucial role in enhancing of reliability of data
processing results. In the author's previous
reports two sets of system performance parameters have been proposed.
The first set of the parameters has been selected for data quality control
(QC) of the raw GNSS observations. This set is composed of: number of cycle slips detected, number of
not-correctable cycle slips, number of loss of locks of signal, number of
single epoch gaps, length of connected carrier phase arcs. The second one has
been defined for examination of the quality, robustness and performance of
the navigation solutions and of performance of the processing strategies. To
the second set inclusion of precision and accuracy, integrity, availability,
confidence- and significance level, convergence time to reach the request
accuracy have been suggested. In
consideration of perturbations of GNSS signal propagation one of the most
crucial requirement is an effective Receiver Autonomous Integrity Monitoring
(RAIM) approach suitable to track autonomously consistency of the navigation
signals (measurements). Furthermore, important expectation from RAIM systems
is to provide up-to-date and valid warnings information to the user when the
system’s performance exceeds specified tolerance of quality levels. Various
selected RAIM techniques are applied to the data sets recorded in quiet and
perturbed ionospheric conditions. To improve efficiency of RAIM algorithm, an
external information source on state of the ionosphere has been used to
enhance reliability of the controlling process. Based on this analysis
usability of the proposed approaches has been evaluated and verified. The paper presents selected results of the
applicability investigations of the proposed parameters in order to ensure a
high level of reliability and accuracy of GNSS real- and near-real time
positioning under the presence of strong ionospheric anomalies. The data-sets
from continuously operated GNSS stations located at high- and low latitudes,
where ionospheric disturbances occur more frequently, have been used for the
analysis. |
2 |
First class poster |
10:25 AM |
Development Of Ionospheric
Forecast Model Over Malaysia Region |
|
|
|
Bahari, S A1; abdullah, M1; hasbi, A M1 |
|
|
|
1Universiti
Kebangsaan Malaysia |
|
|
|
In the near future, Earth could
be faced with the adverse threat of severe space weather that may
significantly affect satellite operations. The top-down space weather
influence can cause large variations in the ionosphere and penetrate to the
Earth's surface. The study on the ionospheric response to space weather can
serve as an early warning and preventive measure to safeguard and maintain
the accuracy of measurements on board satellites and thereby reduce economic
losses. Since space weather has become more prominent due to the fact that
many countries including Malaysia are moving forward using space-based
instruments, the need to produce high-resolution regional ionosphere models
to support navigation, static and kinematic positioning and space weather
research is essential. This project aims to address the need to better
understand the near equatorial orbit (NEqO) through the flagship program
proposed by the National Space Agency (ANGKASA) under the category of ‘The
Use NEqO Orbit for Scientific and Other Applications’. In this project, an
improved regional model (equatorial model) will be developed based on the
spherical cap harmonic analysis (SCHA) method. This model would be the first
regional ionospheric forecast model in the equatorial region. The SCHA is a
regional modeling technique based on solutions of Laplace’s equation over
constrained cap-like region of the Earth. In this project, the ionospheric
TEC over Malaysia is firstly classified into three ionospheric conditions;
“quiet”, “moderate” and “disturbed” based on existing GPS-TEC data. The
regional model of ionospheric TEC over Malaysia will then be developed. The
regional ionospheric condition over Malaysia based on the developed model
will be produced in near real-time. The output of this project will be highly
beneficial as an early warning system to inform stakeholders and government
agencies such as ANGKASA, ATSB, MINDEF, GNSS users and other NEqO satellite
operators. |
3 |
First class poster |
10:30 AM |
Morphology of Shape and Drift of
Ionospheric Irregularities and GPS Scintillation Statistics |
|
|
|
Stevanović, Đ1; Grzesiak, M2; Materassi, M3; W. Wernik, A2 |
|
|
|
1Space
Research Centre; 2Space Research Centre, Warsaw; 3Istituto dei Sistemi Complessi, CNR, Sesto Fiorentino, Italy |
|
|
|
On its way from the satellite to
the receiver, a trans-ionospheric radio signal is distorted by the
ionospheric irregularity that causes signal delay, changes in signal's phase
and amplitude, hence producing a complicated diffraction pattern on the ground.
The variability of the ionospheric drifting irregularity shape and small
scale statistics, reflects the ionospheric behaviour during quiet and storm
events, influencing the signal propagation.
The purpose of this investigation is to infer the characteristics of
the ionospheric irregularities from those of the scintillation diffraction
pattern. Plasma drift velocity and irregularity geometry are studied through
the eigen-analysis of the cross-correlation matrix of gradients of the
scintillation pattern. As far as the scintillation statistics of the GPS
signal is concerned, the study of probability distributions of the signal
fluctuations at different scales provides important information useful in
applications. This study has been
performed on radio signal phase and power recorded by GPS monitors located at
high latitudes at Hornsund, Svalbard and at near equatorial site in Brazil.
Data covers periods before, during and after some major events in 2010-2011.
We validated the estimated properties of ionospheric irregularities using
additional independent measurements. |
4 |
First class poster |
10:35 AM |
B-spline Model of Ionospheric
Scintillation |
|
|
|
Priyadarshi, S1; Wernik, A W1 |
|
|
|
1Space
Research Centre, Poland |
|
|
|
Using Dynamic Explorer (DE) 2
satellite in-situ data an empirical climatological model for the Northern
Hemisphere high latitude ionosphere is prepared. This model incorporates with
B-spline functions, solar and geomagnetic indices to reproduce amplitude
scintillation index and other ionospheric parameters like turbulence strength
parameter Cs, spectral index. As input to the model Dynamic Explorer 2
satellite retarding potential analyzer plasma density data was utilized with
IRI ionospheric model and phase screen propagation model. Similar model is
prepared for Hornsund (Svalbard) and Warsaw (Poland) using GPS receiver GSV
4004b data. The model for –high latitude scintillation prepared using in situ
data is compared with the ionospheric scintillation model prepared for the
Hornsund (Svalbard). For studying the change in scintillation behavior when
we move from –mid latitude to - high latitude we have compared the B-spline
model for Warsaw (Poland) and Hornsund (Svalbard). The comparison is made on
the basis of seasonal behavior and the behavior of scintillation index for
different geophysical conditions. The GPS data used in the present model have
been corrected for the elevation angle E dependence of scintillation index
using the power law dependence on cosecant (E) derived by Priyadarshi &
Wernik (2013). |
5 |
First class poster |
10:40 AM |
Edge-Preserved Imaging of the
Ionospheric Electron Density |
|
|
|
Panicciari, T1; Smith, N1; Da Dalt, F1; Mitchell, C1 |
|
|
|
1University
of Bath |
|
|
|
The ionosphere is a charged and
dynamic medium. Its interaction with an electromagnetic wave can originate
delay and refraction/diffraction effects. These effects take place at
different locations and are caused by ionospheric structures that occur on multiple
scale sizes. The correct localization and identification of those structures
can significantly improve the study of ionospheric behaviour and its effects
in the electromagnetic wave propagation. The imaging of the ionosphere has
been demonstrated successfully in the last decades by using Computerized
Ionospheric Tomographic (CIT) techniques [1]. CIT allows reconstructing the
ionosphere in terms of electron density, and was formerly proposed by Austen
et al. [2]. Observations are collected by ground-based receivers using Global
Positioning System (GPS) and bring the information of the integrated electron
density along the satellite-receiver path. Geometric limitations, together
with uneven and sparse distribution of observations make CIT problem particularly
difficult to solve. Furthermore, there is a need to accommodate
inconsistencies in observations, e.g. due to representativity errors or
residual dispersive offsets in the observations. Therefore, the problem needs
to be regularized and is more challenging as the structures to be resolved
become smaller, i.e. for high resolution tomography. It is proposed a new
technique for ionospheric tomography based on Total Variation (TV)
regularization. TV techniques can help in preserving important details in the
reconstruction such as edges of ionospheric structures. TV is used in
combination with sparse regularization which uses wavelet decomposition to
promote sparsity in the reconstruction. Sparse regularization is used to
compensate for the non-uniform data coverage and provides the input (initial
solution) for the TV algorithm. TV regularization is expected to improve the
detection and location of sharp gradients, in particular for high resolution
CIT. Sharp gradients may cause disturbances in the signal propagation and
their correct location is an important aspect when a high reliability in the
navigation and communications system service is required. [1] G. S. Bust and
C. N. Mitchell, "History, current state, and future directions of ionospheric
imaging," Reviews of Geophysics, vol. 46, 2008 2008. [2] J. R. Austen,
S. J. Franke, and C. H. Liu, "Ionospheric imaging using computerized
tomography," Radio Science, vol. 23, pp. 299-307, 1988. |
6 |
First class poster |
10:45 AM |
Performance Analysis of
Ionosphere Empirical Models for Different Phases of a Solar Cycle |
|
|
|
Najman, P1; Kos, T1 |
|
|
|
1University
of Zagreb, Faculty of Electrical Engineering and Computing |
|
|
|
Ionospheric empirical models
have many uses such as mitigation of ionospheric error for single frequency
GNSS (Global Navigation Satellite System) receivers, estimation of satellite
and receiver inter-frequency biases or Total Electron Content (TEC) calibration
(Jakowski et al, 2011). In all cases, good modelling performance is
desirable. Even though the ionospheric empirical models are widely used, they
show significant degree of inaccuracy, especially during the periods with
high solar activity and geographic areas with complex ionospheric dynamics
(Radicella, 2013). In this study, we
evaluate TEC modelling performance of four different ionospheric empirical
models: IRI2012 (International Reference Ionosphere), Klobuchar, NeQuick2 and
NTCM (Neustrelitz TEC Model). We assume that since each empirical model was
developed using different data, their TEC modelling performance can
significantly differ. In other words, while one model can be the most
accurate under certain conditions, for other conditions some another model
can make the most accurate TEC estimation.
The evaluation of the ionospheric models was done by comparison of
their modelled data with Global Ionospheric Maps (GIMs) provided by Center
for Orbit Determination in Europe (CODE). Each CODE GIM covers are between 180 western to
180 eastern longitude, and 87.5 northern to 87.5 southern latitude with grid
resolution of 5 degrees in longitude and 2.5 degrees in latitude ( 73x71 grid
points). The results were expressed with respect to universal time, month,
and value of solar radio flux index F10.7. The analysis was done for each of
the GIM grid point separately to show which TEC model is the most accurate in
different areas. The comparison was
done for two periods: 2003-2005 and 2010-2012 to see how the results change
for different phases of a solar cycle. In addition, results were compared
with the TOPEX satellite data to test if using models according to our
comparison leads to improvement in overall TEC modelling. References:
N. Jakowski, C. Mayer, M. M. Hoque and V. Wilken . Total electron
content models and their use in ionosphere monitoring. Radio Science, 2011.
46(6). doi: 10.1029/2010RS004620 S. M.
Radicella. Ionosphere Electron Density Models - Present Trend and Validation
Issues, 7th European Conference on Antennas and Propagation, 2013 pp. 3931 |
|
|