If you wish to contribute or participate in the discussions about articles you are invited to contact the Editor

NeQuick Ionospheric Model: Difference between revisions

From Navipedia
Jump to navigation Jump to search
No edit summary
No edit summary
 
(27 intermediate revisions by 5 users not shown)
Line 1: Line 1:
{{Article Infobox2
{{Article Infobox2
|Category=Fundamentals
|Category=Fundamentals
|Authors=J. Sanz Subirana, J.M. Juan Zornoza and M. Hernández-Pajares, Technical University of Catalonia, Spain. [Updated 2017]: R. Prieto-Cerdeira, R. Orús Pérez, European Space Agency, ESA/ESTEC.
|Level=Intermediate
|YearOfPublication=2011, updated 2017
|Title={{PAGENAME}}
|Title={{PAGENAME}}
|Authors=J. Sanz Subirana, JM. Juan Zornoza and M. Hernandez-Pajares, University of Catalunia, Spain.
|Level=Basic
|YearOfPublication=2011
}}
}}
The NeQuick is the ionospheric model proposed to be used by the Galileo single frequency receiver to compute ionospheric corrections. It is based on the original profiler proposed by <ref> [Di Giovanni and Radicella, 1990] Di Giovanni, G. and Radicella, S. M., 1990. An analytical model of the electron density pro_le in the ionosphere. Advances in Space Research. 10(11), pp. 27-30.</ref>.
<b>NeQuick-G</b> is the ionospheric model adopted by the Galileo system in order to help to compute the ionospheric delay corrections for its single-frequency users [Galileo Open Service-Ionospheric Correction Algorithm, 2016] <ref>[http://www.gsc-europa.eu/system/files/galileo_documents/Galileo_Ionospheric_Model.pdf European Union (2016). European GNSS (Galileo) Open Service-Ionospheric Correction Algorithm for Galileo Single Frequency Users. 1.2.] </ref>. The NeQuick-G model is an adaptation for real-time users of the ITU-R NeQuick ionospheric electron density model [ITU-R, 2013] <ref> ITU-R (2013). Ionospheric propagation data and prediction methods required for the design of satellite services and systems. Recommendation ITU-R P. 531-12. </ref> based on the original profiler proposed by [Di Giovanni and Radicella, 1990]<ref> Di Giovanni, G. and Radicella, S. M., 1990. An analytical model of the electron density pro_le in the ionosphere. Advances in Space Research. 10(11), pp. 27-30.</ref> [Hochegger, Nava et al., 2000]<ref>Hochegger, G., B. Nava, S. Radicella and R. Leitinger (2000). “A family of ionospheric models for different uses”. In: Physics and Chemistry of the Earth, Part C: Solar, Terrestrial & Planetary Science 25.4, pp. 307–310. issn: 1464-1917. doi: 10.1016/ S1464-1917(00)00022-2. </ref>.NeQuick a three-dimensional and time dependent ionospheric electron density model based on an empirical climatological representation of the ionosphere, which predicts monthly mean electron density from analytical profiles, depending on solar-activity-related input values: sun spot number or solar flux, month, geographic latitude and longitude, height and Universal Time (UT).  


NeQuick model has been adapted for real-time Galileo single-frequency ionospheric corrections (for convenience referred to as NeQuick-G) in order to derive real-time predictions based on a single input parameter, the Effective ionisation level, Az, which is determined using three coefficients broadcast in the navigation message. As a three-dimensional and time-dependent model, it has the capability to provide vertical or slant ionospheric delay correction by integrating the predicted electron density along the line of sight vector between satellite and receiver.


NeQuick is a tridimensional and time-dependent ionospheric electron density model, which provides electron density in the ionosphere as a function of the position and time. Thence, it allows to compute ionospheric delays (TEC or STEC), as the integrated electron density along any ray path.
In order to take into account both daily variation of the solar activity and the user’s local geomagnetic condition, the NeQuick-G model is driven by and computes its correction by use of the daily solar activity information, so-called effective ionization level (<i>Az</i>) expressed in solar flux unit (<i>10<sup>-22</sup>Wm<sup>-2</sup>Hz<sup>-1</sup></i>)


 
For Galileo single-frequency users, Galileo satellites broadcast three ionospheric coefficients in their navigation message that are used to compute the <i>Az</i> as follows [Galileo - Open Service SIS-ICD, 2016] <ref>[https://www.gsc-europa.eu/system/files/galileo_documents/Galileo-OS-SIS-ICD.pdf Galileo - Open Service - Signal-In-Space Interface Control Document (2016). “European GNSS (Galileo) Open Service, Signal-In-Space Interface Control Document.” Issue 1.3, December 2016. ]</ref>:
NeQuick FORTRAN 77 code was accepted by the \ac{ITU-R} in 2000 and revised in 2002. It is freely available from http://www.itu.int/oth/R0A04000018/en. It is referred to either as version 1 or ITU-R and constitutes the current baseline for GALILEO. This package, includes a comprehensive description of the implementation as well as numerical integration subroutines allowing to compute Vertical and Slant TEC. Online calculations are also available in the Web side http://arpl.ictp.it/nq-online/index.html.
 
 
The input parameters of the model are the position (longitude, latitude and height), the epoch (month and UT) and the solar activity (either expressed <ref  group="footnotes"> The F10.7 index is a measure of the solar activity, i.e., the flux level generated by the sun at at the earth's orbit at a 10.7 cm wavelength. It has been found to correlate well with the sunspot number (Rz). The sunspot number is defined from counts of the number of individual sunspots as well as the number of sunspot groups. The F10.7 index can be measured relatively easily and quickly and has replaced the sunspot number as an index of solar activity for many purposes. The F10.7 and the smoothed sunspot number R12 (12-months moving average) are related by  <math>\mbox{R12}=\sqrt{167273.0+1123.6\,(\mbox{F10.7}-63.7)}-408.99</math>. [ITU-R recommendation].</ref> as  F10.7 or R12). Other internal parameters are the foF2 and M(3000)F2 values, which can be defined according to the ITU-R, amongst other options depending on the purpose.
 
 
The NeQuick model running in the Galileo single frequency receivers will be driven by the ''Effective Ionisation Level'', <math>Az</math> parameter (replacing the solar flux), that is a function of the receiver location.


::<math>
::<math>
Line 30: Line 24:




being <math>I</math> the true magnetic inclination<ref  group="footnotes">{<math>I</math> is <math>0^o</math> at the magnetic equator and <math>90^o</math> at the magnetic poles.</ref>, or ''dip'' in the ionosphere (usually at 300 km), and <math>\phi</math> the geographic latitude of the receiver <ref> [Rawer, 1963] Rawer, K., 1963. Propagation of decameter waves (HF-band), in Meteorological and Astronomical Inuences on Radio Wave Propagation. Ed. Landmark, B. Pergamon Press, New York.</ref>.
being <math>I</math> the true magnetic inclination<ref  group="footnotes"><math>I</math> is <math>0^o</math> at the magnetic equator and <math>90^o</math> at the magnetic poles.</ref>, or ''dip'' in the ionosphere (usually at 300 km), and <math>\Phi</math> the geographic latitude of the receiver [Rawer, 1963]<ref> Rawer, K., 1963. Propagation of decameter waves (HF-band), in Meteorological and Astronomical Inuences on Radio Wave Propagation. Ed. Landmark, B. Pergamon Press, New York.</ref>.
 


The coefficients <math>a_0</math>, <math>a_1</math>,<math>a_2</math> will be broadcasted to the users in the Galileo navigation message and updated at least once a day <ref  group="footnotes"> These parameters will be determined from measured TEC data obtained during the last 24h by the globally distributed Galileo Sensor Stations network.</ref>. Ionospheric Disturbance Flag alerts for five predefined regions (see figure 1) will be also broadcasted to aware the users when the ionospheric correction coming from the Galileo broadcast message might not meet the specified performance. These disturbance flags will be transmitted continuously and updated with the update rate of the Navigation Message (every <math>100</math> minutes).


::[[File: NeQuick_5zones.png |none|thumb|400px|'''''Figure 1:''''' The 5 Regions defined for the Disturbance Flags (in degrees MODIP) -90 to -60, -60 to -30, -30 to 30, 30-60, 60-90. (courtesy of <ref name="Arbesser"> [Arbesser-Rastburg, B., 2006] Arbesser-Rastburg, B., 2006. The Galileo Single Frequency Ionospheric Correction Algorithm. http://sidc.oma.be/esww3/presentations/session4/arbesser.pdf.</ref>).]]
The coefficients <math>a_0</math>, <math>a_1</math>,<math>a_2</math> will be broadcasted to the users in the Galileo navigation message and updated at least once a day <ref group="footnotes"> These parameters are determined from measured TEC data obtained by the globally distributed Galileo Sensor Stations network.</ref>.  


::[[File: NeQuickIonoVTECmap.jpeg |400px|thumb|right|'''''Figure 1:''''' Global vertical Total Electron Content map generated with NeQuick G model at 13h UT for a day in April and Az=193. ]]


The Galileo single frequency receiver algorithm will be based in the following steps (the next paragraph is taken from <ref name="Arbesser"/>):
The Galileo single frequency receiver algorithm will be based in the following steps (the next paragraph is taken from [Arbesser-Rastburg, B., 2006]<ref name="Arbesser">[http://sidc.oma.be/esww3/presentations/Session4/Arbesser.pdf Arbesser-Rastburg, B., 2006. The Galileo Single Frequency Ionospheric Correction Algorithm. European Space Weather Week, Brussels]</ref>):


::'''1.''' <math>Az</math> is evaluated using <math>a_0</math>, <math>a_1</math>,<math>a_2</math> (from navigation message) and MODIP calculated from the data inside the DIPLATS matrix <ref  group="footnotes">In formation on the earth's magnetic field will be stored in the receiver firmware as the DIPLATS matrix. It should be updated every <math>5</math> years, to take account its natural variation. This update shall be considered in the design of the Galileo receivers (MODIP table flash-able in firmware). The ITU-R maps for FoF2 and M3000 (F2) files will be also stored in 12 ASCII files, one for each calendar month.</ref> from the NeQuick model (which depends on estimated receiver position).
::'''1.''' <math>Az</math> is evaluated using <math>a_0</math>, <math>a_1</math>,<math>a_2</math> (from navigation message) and MODIP from the NeQuick model (which depends on estimated receiver position).


::'''2.''' Electron density is calculated for a point along the satellite to receiver path, using the NeQuick model with <math>Az</math> in place of <math>F10.7</math>.
::'''2.''' Electron density is calculated for a point along the satellite to receiver path, using the NeQuick model with <math>Az</math> in place of <math>F10.7</math>.
Line 46: Line 39:
::'''3.''' Steps 1 and 2 are repeated for many discrete points along the  satellite receiver path. The number and spacing of the points will depend on the height and they will be a trade-off between integration error and computational time and power.
::'''3.''' Steps 1 and 2 are repeated for many discrete points along the  satellite receiver path. The number and spacing of the points will depend on the height and they will be a trade-off between integration error and computational time and power.


::'''4.''' All electron density values along the ray are integrated in order to obtain Slant TEC.
::'''4.''' All electron density values along the ray are integrated in order to obtain Slant TEC (STEC).


::'''5.''' Slant TEC is converted to slant delay for correcting pseudo-ranges, by
::'''5.''' Slant TEC, in TECUs, is converted to meters of L1 slant delay for correcting pseudo-ranges, by


::::<math>
::::<math>
I_f = \frac{40.3\cdot 10^{16}}{f^2}TEC \qquad\mbox{(meters)} \quad\mbox{(3)}</math>
I_f = \frac{40.3\cdot 10^{16}}{f^2}TEC \qquad\mbox{(where f is in Hz)} \quad\mbox{(3)}</math>  


An overview of the correction algorithm and initial performance results with broadcast parameters from Galileo satellites have been performed during [[Galileo IOV Satellites|IOV]] in the period April 2013 to March 2014 during solar maximum of Solar Cycle 24 [Prieto-Cerdeira, Orús-Pérez et al, 2014] <ref>[http://gpsworld.com/innovation-the-european-way/ Prieto-Cerdeira, R., R. Orús-Pérez, E. Breeuwer, R. Lucas-Rodríguez and M. Falcone (2014). ‘Innovation: The European Way. Performance of the Galileo Single-Frequency Ionospheric Correction During In-Orbit Validation’. In: GPSworld 25.6, pp. 53–58. ]</ref>.


Notice that, as with the [[Klobuchar Ionospheric Model|Klobuchar]] model, the ionospheric corrections computed by the NeQuick can be used for any GNSS signal (GPS, GLONASS, Galileo…) simply by setting the corresponding frequency in the equation (3).
Notice that, as with the [[Klobuchar Ionospheric Model|Klobuchar]] model, the ionospheric corrections computed by the NeQuick G can be used for any GNSS signal (GPS, GLONASS, Galileo…) simply by setting the corresponding frequency in the equation (3).





Latest revision as of 21:28, 19 September 2017


FundamentalsFundamentals
Title NeQuick Ionospheric Model
Author(s) J. Sanz Subirana, J.M. Juan Zornoza and M. Hernández-Pajares, Technical University of Catalonia, Spain. [Updated 2017]: R. Prieto-Cerdeira, R. Orús Pérez, European Space Agency, ESA/ESTEC.
Level Intermediate
Year of Publication 2011, updated 2017

NeQuick-G is the ionospheric model adopted by the Galileo system in order to help to compute the ionospheric delay corrections for its single-frequency users [Galileo Open Service-Ionospheric Correction Algorithm, 2016] [1]. The NeQuick-G model is an adaptation for real-time users of the ITU-R NeQuick ionospheric electron density model [ITU-R, 2013] [2] based on the original profiler proposed by [Di Giovanni and Radicella, 1990][3] [Hochegger, Nava et al., 2000][4].NeQuick a three-dimensional and time dependent ionospheric electron density model based on an empirical climatological representation of the ionosphere, which predicts monthly mean electron density from analytical profiles, depending on solar-activity-related input values: sun spot number or solar flux, month, geographic latitude and longitude, height and Universal Time (UT).

NeQuick model has been adapted for real-time Galileo single-frequency ionospheric corrections (for convenience referred to as NeQuick-G) in order to derive real-time predictions based on a single input parameter, the Effective ionisation level, Az, which is determined using three coefficients broadcast in the navigation message. As a three-dimensional and time-dependent model, it has the capability to provide vertical or slant ionospheric delay correction by integrating the predicted electron density along the line of sight vector between satellite and receiver.

In order to take into account both daily variation of the solar activity and the user’s local geomagnetic condition, the NeQuick-G model is driven by and computes its correction by use of the daily solar activity information, so-called effective ionization level (Az) expressed in solar flux unit (10-22Wm-2Hz-1)

For Galileo single-frequency users, Galileo satellites broadcast three ionospheric coefficients in their navigation message that are used to compute the Az as follows [Galileo - Open Service SIS-ICD, 2016] [5]:

[math]\displaystyle{ Az=a_0+a_1\mu + a_2 \mu^2 \qquad\mbox{(1)} }[/math]


where [math]\displaystyle{ \mu }[/math] is the modified dip latitude, or MODIP:

[math]\displaystyle{ \tan \mu=\displaystyle \frac{I}{\sqrt{\cos \phi}}\qquad\mbox{(2)} }[/math]


being [math]\displaystyle{ I }[/math] the true magnetic inclination[footnotes 1], or dip in the ionosphere (usually at 300 km), and [math]\displaystyle{ \Phi }[/math] the geographic latitude of the receiver [Rawer, 1963][6].


The coefficients [math]\displaystyle{ a_0 }[/math], [math]\displaystyle{ a_1 }[/math],[math]\displaystyle{ a_2 }[/math] will be broadcasted to the users in the Galileo navigation message and updated at least once a day [footnotes 2].

Figure 1: Global vertical Total Electron Content map generated with NeQuick G model at 13h UT for a day in April and Az=193.

The Galileo single frequency receiver algorithm will be based in the following steps (the next paragraph is taken from [Arbesser-Rastburg, B., 2006][7]):

1. [math]\displaystyle{ Az }[/math] is evaluated using [math]\displaystyle{ a_0 }[/math], [math]\displaystyle{ a_1 }[/math],[math]\displaystyle{ a_2 }[/math] (from navigation message) and MODIP from the NeQuick model (which depends on estimated receiver position).
2. Electron density is calculated for a point along the satellite to receiver path, using the NeQuick model with [math]\displaystyle{ Az }[/math] in place of [math]\displaystyle{ F10.7 }[/math].
3. Steps 1 and 2 are repeated for many discrete points along the satellite receiver path. The number and spacing of the points will depend on the height and they will be a trade-off between integration error and computational time and power.
4. All electron density values along the ray are integrated in order to obtain Slant TEC (STEC).
5. Slant TEC, in TECUs, is converted to meters of L1 slant delay for correcting pseudo-ranges, by
[math]\displaystyle{ I_f = \frac{40.3\cdot 10^{16}}{f^2}TEC \qquad\mbox{(where f is in Hz)} \quad\mbox{(3)} }[/math]

An overview of the correction algorithm and initial performance results with broadcast parameters from Galileo satellites have been performed during IOV in the period April 2013 to March 2014 during solar maximum of Solar Cycle 24 [Prieto-Cerdeira, Orús-Pérez et al, 2014] [8].

Notice that, as with the Klobuchar model, the ionospheric corrections computed by the NeQuick G can be used for any GNSS signal (GPS, GLONASS, Galileo…) simply by setting the corresponding frequency in the equation (3).


Notes

  1. ^ [math]\displaystyle{ I }[/math] is [math]\displaystyle{ 0^o }[/math] at the magnetic equator and [math]\displaystyle{ 90^o }[/math] at the magnetic poles.
  2. ^ These parameters are determined from measured TEC data obtained by the globally distributed Galileo Sensor Stations network.


References

  1. ^ European Union (2016). European GNSS (Galileo) Open Service-Ionospheric Correction Algorithm for Galileo Single Frequency Users. 1.2.
  2. ^ ITU-R (2013). Ionospheric propagation data and prediction methods required for the design of satellite services and systems. Recommendation ITU-R P. 531-12.
  3. ^ Di Giovanni, G. and Radicella, S. M., 1990. An analytical model of the electron density pro_le in the ionosphere. Advances in Space Research. 10(11), pp. 27-30.
  4. ^ Hochegger, G., B. Nava, S. Radicella and R. Leitinger (2000). “A family of ionospheric models for different uses”. In: Physics and Chemistry of the Earth, Part C: Solar, Terrestrial & Planetary Science 25.4, pp. 307–310. issn: 1464-1917. doi: 10.1016/ S1464-1917(00)00022-2.
  5. ^ Galileo - Open Service - Signal-In-Space Interface Control Document (2016). “European GNSS (Galileo) Open Service, Signal-In-Space Interface Control Document.” Issue 1.3, December 2016.
  6. ^ Rawer, K., 1963. Propagation of decameter waves (HF-band), in Meteorological and Astronomical Inuences on Radio Wave Propagation. Ed. Landmark, B. Pergamon Press, New York.
  7. ^ Arbesser-Rastburg, B., 2006. The Galileo Single Frequency Ionospheric Correction Algorithm. European Space Weather Week, Brussels
  8. ^ Prieto-Cerdeira, R., R. Orús-Pérez, E. Breeuwer, R. Lucas-Rodríguez and M. Falcone (2014). ‘Innovation: The European Way. Performance of the Galileo Single-Frequency Ionospheric Correction During In-Orbit Validation’. In: GPSworld 25.6, pp. 53–58.