Issue 
J. Space Weather Space Clim.
Volume 10, 2020
Topical Issue  Scientific Advances from the European Commission H2020 projects on Space Weather



Article Number  29  
Number of page(s)  13  
DOI  https://doi.org/10.1051/swsc/2020030  
Published online  16 July 2020 
Research Article
Impact of mediumscale traveling ionospheric disturbances on network realtime kinematic services: CATNET study case
Research Group on Astronomy and Geomatics (gAGE), Universitat Politècnica de Catalunya (UPC), Jordi Girona 1–3, 08034 Barcelona, Spain
^{*} Corresponding author: cristhian.timote@upc.edu
Received:
2
December
2019
Accepted:
9
June
2020
Mediumscale traveling ionospheric disturbances (MSTIDs) are fluctuations in the plasma density that propagate through the upper layer of the atmosphere at velocities of approximately 100 m/s and periods reaching some tens of minutes. Due to their wavelengths, MSTIDs can degrade the performance of differential positioning techniques, such as realtime kinematics (RTK) or networkRTK (NRTK). This paper defines a novel methodology as a tool for relating the errors in NRTK positioning based on an MSTIDs indicator using the second difference in time of the slant total electron content (STEC). The proposed methodology performs integer ambiguity resolution (IAR) on the undifferenced measurements instead of using doubledifferenced carrierphase measurements, as it is usual in RTK and NRTK. Statistical tests are applied to evaluate the degradation in the position errors caused by the impacts of MSTIDs on RTK and NRTK positioning over a data set spanning one year gathered from the CATNET network; a dualfrequency network of fixed permanent GNSS receivers located at the midlatitudes of northeastern Spain. With the development of the proposed methodology for measuring the position degradation, another results of the present research are the establishment of thresholds for the proposed MSTIDs index, which can be used to monitor the positioning solution and to warn users when the measurements are affected by MSTIDs events, relating the position error to MSTIDs that affect not only the user receivers but also of the reference receivers within the network.
Key words: mediumscale TIDs / ionospheric disturbances / NRTK / integer ambiguity resolution
© C.C. Timoté et al., Published by EDP Sciences 2020
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Traveling ionospheric disturbances (TIDs) are plasma density fluctuations that propagate through the ionosphere with a broad range of velocities and frequencies. Depending on the TID characteristics, some authors (Jacobson et al., 1995) distinguish between largescale TIDs (LSTIDs) with periods greater than 1 hour and moving faster than 0.3 km/s, and mediumscale TIDs (MSTIDs) with shorter periods (less than 1 h) and slower velocities (less than 0.3 km/s). Although the sources of MSTIDs are not unique, MSTIDs are linked to meteorological phenomena such as neutral winds, eclipses and the solar terminator (ST), which produces atmospheric gravity waves (AGW) and manifests them as TIDs at ionospheric heights (HernándezPajares et al., 2006).
The ionosphere represents one of the major signal error sources in positioning using Global Navigation Satellite Systems (GNSS). Without any modelling, the contributions of ionospheric errors to GNSS observables can reach up to tens of metres (Sanz et al., 2013). GNSS doublefrequency receivers can eliminate up to 99.9% of the ionospheric refraction effect by implementing the ionospherefree (IF) combination of code or carrierphase measurements at frequencies f _{1} and f _{2}. In contrast, singlefrequency receivers must apply ionospheric models to account for the ionospheric refraction (RoviraGarcia et al., 2019). The effects of the frequencydependent terms (including the ionospheric refraction) can be estimated using the geometryfree (GF) combination of code (P _{GF} = P _{2} − P _{1}) or carrierphase measurements (L _{GF} = L _{1} − L _{2}), which cancels out the geometric part of the measurements.
In the presence of MSTIDs, the total electron content (TEC) in the ionosphere increases by some TEC units (TECU), where 1 TECU represents an ionospheric delay of approximately 16 cm in the L1/E1 GNSS frequency (1575.42 MHz), which is less than the expected accuracy of the corrections provided by ionospheric models (RoviraGarcia et al., 2016); this error has a minor impact on GNSS standard positioning services based on code measurements. By contrast, highaccuracy positioning services (HAS) based on carrierphase measurements, such as the precise point positioning (PPP), realtime kinematics (RTK), or network RTK (NRTK), require a precise modelling of the ionospheric effects. Indeed, fluctuations reaching several tenths of TECUs in the ionospheric corrections can affect the positioning accuracy of HAS when a precise ionospheric modelling is needed (HernándezPajares et al., 2000).
In the particular case of both RTK and NRTK, dualfrequency receivers can implement the L _{IF} combination to account for ionospheric refraction at the cost of enlarging the noise by a factor of three compared with the noise contributed by each individual signal. This is because the difference of squared frequencies dividing the measurements in the L _{IF} and P _{IF} definition (Sanz et al., 2013). For singlefrequency receivers, the ionospheric delay of GNSS signals is corrected assuming that it is identical to the ionospheric delay experienced by the closest reference receiver (RTK) or to a linear combination of the ionospheric effects experienced by the network of reference receivers (NRTK). MSTIDs can degrade the ionospheric mitigation of singlefrequency users because the baselines of RTK and NRTK are tens of kilometres, which are distances comparable to the typical wavelengths of MSTIDs (Alves et al., 2001; Lachapelle & Alves, 2002).
The present paper presents a new methodology for characterizing the impact of MSTIDs on NRTK positioning, and proposes an index as a tool to alert users about the MSTIDs effect on the GNSS user positioning. The manuscript is organized as follows: Section 1 is the current introduction. Section 2 presents the data set used and the network of NRTK stations to which the study is applied. Section 3 develops a novel methodology for creating the reference measurements from undifferenced measurements for which the carrierphase ambiguities have been fixed. Section 4 presents the implementation of an ionospherefree combination to obtain the reference navigation solution. Section 5 assesses the degradation of the positioning using a single frequency. Section 6 analyses the degradation of the positioning caused by MSTIDs, presents the implementation of a new index to warn RTK and NRTK users about the impacts of MSTIDs effects, and provides statistical results in order to assess the MSTIDs index performance. Finally, Section 7 summarizes the main conclusions of this work.
2 Data: the CATNET NRTK service
The methodology implemented in this study is developed using data from a network of sixteen permanent stations equipped with dualfrequency GNSS receivers. These stations are part of the CATNET NRTK service of the Cartographic and Geologic Institute of Catalonia (ICGC) in Spain (Talaya & Bosch, 1999). Nine stations are selected for the study in the following configuration: three of the stations are used as “service user” or “rover” receivers (PLAN, MARE, and SBAR, depicted in blue in Fig. 1), whereas the remaining six are treated as fixed reference receivers (BEUD, BELL, SONA, GARR, LLIV, and CASE, depicted in red in Fig. 1). Details about the implemented network are presented in Table 1 and Figure 1.
Fig. 1 CATNET NRTK network. Receivers in red: reference receivers; receivers in blue: user receivers. 
CATNET receivers’ specifications.
Additionally, in this study, the reference receiver GARR is used as the reference RTK receiver instead of implementing the typical approach of dynamically selecting the closest reference receiver with respect to the user receiver. This approach is applied to represent the effect of the distance on the position solution for each of the three user receivers. The proposed network configuration does not represent the typical NRTK baselines between the user receiver and reference receivers, as we are implementing fixed stations as user receivers. The data set used in the present study comprises a full year of data from day of year (DoY) 200 in 2017 to DoY 200 in 2018.
3 Methodology for obtaining unambiguous measurements in the CATNET service
The accuracy of a HAS depends on factors such as the satellite geometry and/or the quality of the different HAS corrections (such as the ionospheric corrections). This dependency on several factors makes it difficult to separately attribute the degradation in the navigation solution to any of these factors. RTK and NRTK are based on fixing carrier phase ambiguities; therefore, one of the most relevant topics when computing the network correction is the ability to fix carrierphase ambiguities at the network stations. To exclude this dependency, we have fixed all these carrier phase ambiguities offline. In this section, we explained how we process the data to obtain solutions with fixed carrierphase ambiguities.
3.1 Fixing carrierphase ambiguities
GNSS observables (consisting of the pseudorange P and carrier phase L) emitted by a satellite j and collected by a receiver i, at a frequency f _{ m } can be modelled as (Sanz et al., 2013):(1)where stands for the geometric distance between the satellite and receiver, and c is the speed of light. T _{ i } and T ^{ j } are the receiver and satellite clock offsets, respectively. δ_{mi}, , D _{mi}, and are the instrumental delays (floated numbers), at the frequency f _{ m }, of each of the devices. stands for the tropospheric effect on the GNSS signal. is the ionospheric delay (advance for the phase and delay for the code measurement) proportional to . Finally, for the carrierphase measurement, λ _{ m } represents the carrierphase wavelength, and is the integer ambiguity. Although carrierphase measurements are more precise than code measurements , due to the unknown value of , measurements are ambiguous.
In NRTK positioning, the expressions in equation (1) are simplified by forming double differences (∆∇) with respect to a reference receiver and a reference satellite. One of the advantages of this technique is that the parameters depending upon a single device (either a satellite or a receiver) are cancel out, producing the following observable (Seeber, 2008):(2)
According to equation (2), computations in a NRTK network requires carrierphase measurements with fixed ambiguities between the stations of the network. Estimation of the exact integer value of can be achieved by means of any of the wellknown method for fixing such ambiguities, such as the LAMBDA method (Teunissen, 1995). Therefore, the accuracy of the position estimates depends upon the ability to fix carrierphase ambiguities, more precisely, on how many ambiguities are present and how fast the ambiguities are fixed (Sanz et al., 2013).
To reduce the dependency on the ability to successfully resolve ambiguities (ambiguityfixing success ratio), we fix the integer part of the carrierphase ambiguities of all the measurements of the undifferenced mode (i.e. before Eq. (2)). To achieve this, we have taken advantage of knowing the precise coordinates of each receiver to solve the phase biases, the socalled fractional part of the carrierphase ambiguities (Collins et al., 2008; Ge et al., 2008), to estimate the integer part of the carrierphase ambiguities for all the measurements in the network.
Detailed information about the whole ambiguityfixing process can be found in (Juan et al., 2020); here, we simply present some examples of this process. As a first example of the ambiguityfixing capability, Figure 2a depicts the actual floating values of the ambiguities (modulo 10 of the floating values) for the receiver SBAR during the day 310 in 2017, illustrating that the real estimates of the ambiguities are close to the integer values (modulo 1 of the floating values). Figure 2b presents how different the floating values of the carrierphase ambiguities are from their closest integer values. The differences are typically less than 0.2 cycles, enabling the confident resolution of carrierphase ambiguities to their integer values.
Fig. 2 Fractional part of the L1 carrierphase ambiguities for all the satellites in view from the receivers SBAR on day 359 in 2017. The panels show the fixingambiguity process in which it is possible to fix the carrierphase ambiguities with confidence. 
A second example showing the advantage of fixing undifferenced ambiguities can be found in Figure 3, which graphically presents a comparison between the GF combinations of the pseudoranges and carrier for the user receiver MARE. Figure 3a depicts the values without fixing the carrierphase ambiguities, i.e. the raw measurements. Carrierphase measurements are much more precise (less noisy) than pseudorange measurements. However, due to the presence of uncorrected ambiguities, the carrier phases present biases that make them less accurate than the pseudoranges. Figure 3b depicts the values after fixing the carrierphases ambiguities; after this process, the ionospheric delays obtained from the carrier phases are precise and accurate. In this panel, we highlight the ionospheric delays experienced by the GPS satellite G32 in blue to demonstrate the effect produced by a MSTID. Figure 3c presents a magnified view of the Figure 3b, exposing the previously mentioned disturbance. Finally, Figure 3d shows the advantage of our approach (i.e., fixing the carrierphase ambiguities in an undifferenced manner) with respect to the standard RTK method of fixing the carrierphase ambiguities of the doubledifferenced (DD) measurements. Indeed, this panel depicts the ionospheric delays in the measurements from the satellite G32 and in those collected by the receivers MARE, GARR, and PLAN in blue, red and cyan, respectively. In all three arcs, the presence of the MSTID is clear, as is the propagation of the MSTID. This is not the case for the DD measurements from MARE (green) and PLAN (black) with respect to GARR. Indeed, only the differential effects can be seen, with the amplitude of the MSTID differential effect dependent on the baseline between the receivers.
Fig. 3 Top panels: example of the ionospheric delays measured with the carrier phases (green) and pseudoranges (red) corresponding to the receiver MARE during DoY 310 in 2017. The top left panel (a) displays the raw measurements. The top right panel (b) corresponds to the measurements where the carrierphase ambiguities have been fixed. The ionospheric delay for the satellite G32, affected by a MSTID, is highlighted in blue. The bottom left panel (c) presents a more detailed picture of the MSTID event highlighted in panel (b). The bottom right panel (d) depicts the undifferenced STEC for the satellite G32 measured by receivers MARE (blue), GARR (red), and PLAN (cyan). The double differences of the ionospheric delays with respect to GARR are also depicted in green (MARE) and black (PLAN). 
4 Obtaining the reference navigation solution by means of the ionospherefree combination
Following the methodology presented in the previous section, we fix the carrierphase ambiguities for each of the nine receivers (the user and reference receivers) selected from CATNET. We use these unambiguous data to compute all the navigation solutions. Note that these unambiguous measurements can be modelled as:(3)which is mathematically equivalent to the pseudoranges model in equation (1), except for the sign of the ionospheric effect and the carrierphase biases.
Equation (3) is the standard equation for the navigation problem, i.e. for obtaining the user position . After fixing carrierphase ambiguities, there are several ways to proceed:

To work with undifferenced measurements (as in Eq. (3)): in this case, one needs to know the satellite clock offsets (T ^{ j }) and the satellite phase biases (), which are usually provided by several centres, such as the International GNSS Service.

For a common satellite, to make differences between the user measurements and the measurements from a nearreference receiver: in this case, one does not need the satellite parameters.

To make double differences using a common reference satellite (as in Eq. (3)): in this case, one does not need to estimate the receiver clock (T _{ i }).
Once ambiguities are fixed, these three techniques should provide similar results. For simplicity, let us assume we are using double differences: starting from an approximate position () for the user, one can compute a value for the receiversatellite vector (). Therefore, equation (2) can be written as (Sanz et al., 2013):(4)where obtaining the correction vector is the basic goal of the navigation problem. This is accomplished by solving equation (4) using techniques such as least squares or Kalman filter.
Differential tropospheric corrections, , are usually estimated alongside the correction vector using a simple tropospheric model. Therefore, the main unmodelled term in equation (4) is the differential ionospheric correction , which cannot be solved in equation (4). In fact, the way that is treated determines the different types of navigation solutions that we present in this work:

Single frequency RTK solution: For each satellite, the ionospheric delays for both user and reference receiver measurements are the same ( = 0).

Single frequency NRTK solution: The differential ionospheric delays, at the user position can be estimated by linearly interpolating the delays in a set of permanent receivers surrounding the user receiver.

Ionospherefree solution: With a receiver having two frequencies, the user can build the L _{IF} combination, as defined in the introduction. This combination cancels out the ionospheric effect on the GNSS signals.
From their definitions, NRTK and RTK navigation solutions will be affected by the errors in the assumptions. These errors will become larger as the baseline between the user and the reference receiver increases. In contrast, ionospherefree measurements are unaffected by the differential ionospheric corrections. Therefore, the corresponding navigation solution will be independent of the baseline between the user and the reference receiver.
Figure 4 depicts an example of the navigation solution obtained with the unambiguous L _{IF} over 24 h for the receiver PLAN, the closest user receiver with respect to the reference receiver (15.2 km from the reference receiver GARR). As it can be seen, the 3D error in the navigation solution is usually below 10 cm, except during some epochs presenting poor satellite geometries, i.e., when the dilution of precision (DOP) is larger than 6 (red points). 3D positioning errors better than 10 cm constitute the nominal performance expected for NRTK positioning (see for instance Cannon et al., 2001; Lachapelle & Alves, 2002; Wielgosz et al., 2005).
Fig. 4 Ionospherefree combination for the navigation solution at the PLAN receiver. From top to bottom: 3D error, up, north, and east components. 
A more general assessment of the 3D positioning error using L _{IF} is performed in Figure 5. This figure depicts the statistic for the three user receivers during the entire data period. This statistic is calculated by means of the complementary of the cumulative distribution function (CCDF, also represented as 1CDF).
Fig. 5 Statistics for the 3D positioning error using the IF combination: the Yaxis depicts the probability of the 3D positioning error being greater than the value on the Xaxis. 
As depicted in Figure 5, 95% of the cases have an error smaller than 8 cm. This value can be established as the nominal reference value at the 95th percentile for the 3D positioning error with L _{IF}.
As expected, the three user receivers report similar statistics regardless of their distance to the reference station GARR. This is because, as mentioned above, L _{IF} measurements are not affected by differential ionospheric delays. Therefore, for each user receiver, this solution can be taken as the reference solution when one wants to assess the degradation in the positioning using the RTK or NRTK techniques, i.e., techniques affected by the correctness of the ionospheric assumptions.
5 Navigation solution using either the RTK or the NRTK techniques with a singlefrequency receiver
When a user navigates solely with L1 measurements (singlefrequency receivers), and therefore it is not able to build the IF combination, the navigation solutions are affected by the way in which the differential ionospheric delays are corrected.
And example of that can be seen in Figure 6, where the navigation solution using L _{IF}, RTK and NRTK are compared. This comparison is done using the CCDFs of the navigation solutions for the three user receivers and using the three HAS techniques defined in the previous section.
Fig. 6 Ionospherefree (L _{IF}) solution in green, NRTK singlefrequency (L1) solution in blue, and RTK singlefrequency (L1) solution in red for three receivers: (a) PLAN, (b) MARE, and (c) SBAR. 
The shortest baseline length of 15.2 km is between PLAN and GARR. In this case, as shown in Figure 6a, up to 15 cm, the 3D position error using L1 (RTK or NRTK) is smaller than that using L _{IF}. This is because the differential ionospheric effects between PLAN and GARR are quite small (see the example in Fig. 3d) and, as mentioned in the introduction, the thermal error of L _{IF} is amplified by a factor of three with respect that of L1 and therefore it dominates the error. On the other hand, the RTK and NRTK singlefrequency solutions are quite equivalent, and the degradation of the solution with respect to that using L _{IF} occurs with relatively low probability (less than 1%).
Figure 6b presents the results for MARE, located at 50.7 km from GARR. The differential ionospheric term in equation (4) becomes greater than in the results for PLAN, and more than 10% of the cases using L1 with RTK positioning have positioning errors larger than 8 cm. However, if the ionospheric delay is corrected using the data from the network (NRTK), the probability of having an error larger than or equal to 8 cm is approximately 8%. This result is not exceedingly different from the percentage found in the case with L _{IF} (approximately 5%).
Finally, panel (c) of Figure 6 depicts the results for the user receiver SBAR (79.3 km from GARR). The degradation of the RTK solution increases considerably. Indeed, the probability of having an error larger than 8 cm in the 3D positioning reaches nearly 40%, while with NRTK positioning, the probability is approximately 10%, i.e., 5% more than the solution using L _{IF}.
Taking into account the previous results, we define the 3D positioning degradation, δ _{3D}, as the modulus of the difference between the singlefrequency navigation (RTK or NRTK) solution with respect to the reference solution :(4)
Figure 7 presents the statistical 3D degradation by means of the CCDF, for the three user receivers with the previously defined data set corresponding to the years 2017 and 2018. These plots confirm the results obtained and presented above. Indeed, for SBAR, as shown in Figure 7c, 10% of the RTK solutions present a degradation larger than 14 cm, while less than 1% of the NRTK cases present a degradation larger than 14 cm. However, if one compares the NRTK results for SBAR with those for PLAN (Fig. 7a), one can see that the 99th percentile is less than 9 cm for PLAN but 12 cm for SBAR. This represents a clear increase in the degradation in the NRTK solution, which is related with to the larger baseline.
Fig. 7 CCDFs of the 3D degradation in the positioning error for the receivers (a) PLAN, (b) MARE, and (c) SBAR during the 1year period of 2017–2018. 
6 Effects of MSTIDs on position degradation: definition of an index linked to MSTIDs
We have shown that the NRTK solutions are closer than RTK solutions to the dualfrequency (L _{IF}) solution. However, as observed in Figure 7, there is some degradation in the positioning for the two receivers with the largest separation distance in the NRTK solution with respect to the IF solution. In this section, we analyse the part of the degradation in the NRTK solutions associated with the presence of MSTIDs. Indeed, NRTK positioning assumes that the differential ionospheric delays between receivers present a linear spatial behaviour. This linearity is broken by the presence of fluctuations with wavelengths (such as MSTIDs wavelengths) comparable to the network baselines. In the first subsection, we will show how MSTIDs can be detected and then propose the definition of an index linked to MSTIDs activity. In the second subsection, we will present some examples of the relationship between the position degradation and the presence of MSTIDs. Finally, in the third subsection, we will perform a statistical analysis of this relationship.
6.1 Detecting MSTIDs and defining a MSTIDs activity index
HernándezPajares et al. (2006) showed that the typical period of a MSTID is on the order of tens of minutes and further defined a method to detect such a fluctuations at this time scale by building the second difference in time (Δ^{2}) of the slant total electron content (STEC) at a given epoch t, defined as:(5)where τ is set in such a way that it is optimal for detecting MSTIDs with a period of 10 min; i.e., when τ is equal to 5 min (HernándezPajares et al., 2006), there is sufficient sensitivity to detect the presence of MSTIDs with periods reaching tens of minutes. Additionally, by setting τ to 5 min, we are able to provide warnings about the presence of MSTIDs activity on a short time scale.
To quantify the activity of a MSTID, for each satellitereceiver pair, we define the MSTIDs index as the integrated amplitude over an interval 2τ (i.e., 10 min, which results in 20 samples when the measurements are sampled at a rate of 30 s):(6)where M (ε) is an obliquity factor for mitigating larger values of Δ^{2} STEC (i) at low elevations.
Figure 8 depicts the Δ^{2} STEC, as defined in equation (5), for the receivers GARR, SONA, and SBAR tracking the same GPS satellite (G01). It is clear that the presence of a MSTID (with an amplitude close to 1 TECU and a period of approximately 1000 s) affects the three receivers. This STEC fluctuation clearly breaks the linearity assumption in NRTK positioning.
Fig. 8 MSTID effects on the second difference of the STEC for three receivers: GARR, SONA, and SBAR. 
Figure 9 presents the results obtained after calculating MSTID_{IDX} (defined in Eq. (6)) for the same satellitereceiver pairs as those in the example depicted in Figure 8. Certainly, it is possible to detect and isolate the MSTID event occurring at approximately 33,000 s of DoY 039, when MSTID_{IDX} reaches a value of 0.5 TECU, which is 0.05 m for L _{GF} (1 TECU = 0.105 L _{GF} m).
Fig. 9 MSTID index definition applied to three receivers on DoY 039 (February 6th) 2014: GARR, SONA, and SBAR. 
6.2 Relationship of MSTID_{IDX} with the error in the ionospheric corrections
Previous sections investigate the relation of the MSTID index with respect to the degradation in the ionospheric corrections and, consequently, in the position accuracy. Regarding the ionospheric corrections, note that this degradation is not influenced only by the MSTIDs activity during the user measurements period. Indeed, because the ionospheric delay is interpolated in NRTK positioning, any activity at any of the reference site contributes to the user positioning degradation. This can be seen in Figure 10, which corresponds to DoY 221 in 2017.
Fig. 10 ∆^{2} STEC for the user receiver SBAR (blue) and reference receiver BEUD (green); the error in the ionospheric correction (red) is clearly correlated with the presence of a MSTID. 
Figure 10 depicts the Δ^{2} STEC experienced by measurements from the GPS satellite G09 for the user receiver SBAR and for the reference receiver BEUD. The resulting error in the ionospheric correction (which is obtained by fitting a linear model from the reference sites) at SBAR for the specific satellite G09 correlates with the presence of a MSTID that propagates towards the southwest. This results in a fluctuation of Δ^{2} STEC observed first at BEUD and then at SBAR. Evidently, a MSTID affecting any of the reference receivers participating in the computation of the ionospheric correction in the network contributes to the error at the user receiver.
In light of previous results, one conservative approach is to define, for each satellite in view, an MSTID index associated with the whole NRTK network. This is accomplished by taking the maximum value of the MSTID_{IDX} experienced by any reference receiver collecting data from that satellite. Taking into account this network of MSTID_{IDX} values, users can be warned about large ionospheric errors in their corrections that could affect their navigation solutions. However, the accuracy of the navigation solution depends on other factors in addition to the quality of the ionospheric corrections, such as the geometry of the observations. For that reason, we can take advantage of having undifferenced and unambiguous carrierphase measurements to assess the errors in the ionospheric corrections at the user positions. Indeed, Figure 11 depicts the ionospheric errors during the studied period for the three user receivers. In this assessment, we consider only observations with elevation angles above 30°; in this way, we exclude large errors at low elevations that are not related to ionospheric activity. Each panel depicts, for the three user receivers, histograms of the 95th and 99th percentiles of the ionospheric errors corresponding to specific values of the MSTID_{IDX}. As it can be seen, as the MSTID_{IDX} increases the probability of having larger ionospheric errors also increases. Therefore, one could select some threshold value for MSTID_{IDX} to exclude observations with large errors. For this purpose, we include the CDF of the overall cases in order to account for the percentage of the cases that will be filtered out with the MSTID_{IDX} threshold. For instance, by setting the MSTID_{IDX} to 0.1 TECU it is possible to guarantee that, for the three user receivers, the 95th percentile of the ionospheric corrections will be approximately less than 0.2 TECU, at the cost of excluding less than a 10% of the total measurements collected.
Fig. 11 Histogram of the error of the ionospheric corrections; 95th percentile (filled red boxes), 99th percentile (empty blue boxes). CDF of the ionospheric corrections (black line). The statistic is done for the three rover receivers: (a) PLAN, (b) MARE, and (c) SBAR. 
6.3 Relationship of MSTID_{IDX} with the degradation in positioning
As mentioned in previous section, for a user that needs to correct for the ionospheric effects on their observations (i.e., users with singlefrequency measurements), the errors in the ionospheric corrections are translated to the position accuracy. Figure 12 presents several examples of the degradation in the 3D position errors for three stations (from top to bottom PLAN, MARE, and SBAR) in contrast to the amplitude of the MSTID provided by the MSTID_{IDX}, on the left and on the right columns respectively. The examples correspond to days 298 and 325 of 2017. A noticeable relationship exists between the degradation in the position error of the receiver and the maximum detected MSTID amplitude in the NRTK network. For a better representation of this relationship, MSTID_{IDX} has been amplified by an arbitrary scale factor from its original amplitude value (in metres). The user receiver PLAN shows less severe degradation of the position error than the two other receivers.
Fig. 12 3D position errors and maximum MSTIDs amplitudes on days 298 and 325 of 2017 for the user receivers (a) PLAN, (b) MARE, and (c) SBAR. 
For results similar to those depicted in Figure 12 encompassing the whole period of time studied in this work, readers are referred to the following public website: http://147.83.47.222/TechTIDE_database/2017/DDD/NRTK_performance_CATNET/Plots/, where “DDD” is the threedigit DoY.
6.4 Statistical analysis
We performed a statistical analysis based on the CCDFs for the position results from DoY 200 in 2017 to DoY 200 in 2018. Figure 13 provides the probability of 3D position degradation for three user receivers (PLAN, MARE, and SBAR). Figure 13 also illustrates a comparison of the degradation using RTK and NRTK positioning. Moreover, to see the relationship between the 3D position degradation and MSTIDs, the NRTK solution is also computed when:

MSTID_{IDX} is lower than 0.15 TECU, for which, in SBAR, 0.1% of the cases present a degradation larger than 10 cm in 3D coordinates.

MSTID_{IDX}, is lower than 0.10 TECU, for which, in SBAR, 0.2% of the cases present a degradation larger than 10 cm in 3D coordinates.
Fig. 13 3D position degradation probabilities for user receivers: (a) PLAN, (b) MARE, and (c) SBAR. The solid lines represent the position solutions for RTK (red) and NRTK (blue) for all MSTIDs amplitudes. The NRTK solutions are also shown with MSTIDs amplitudes of less than 0.15 TECU (black) and with MSTIDs amplitudes of less than 0.10 TECU (green). 
As shown in Figure 13, for the user receivers SBAR and MARE, which are located farther away than PLAN, there is a noticeable reduction in the degradation error when solutions are excluded from the statistics taking into account MSTID_{IDX}. In Tables 2–4, such a reduction in the degradation in the positioning error (by means of the 50th, 95th, and 99th percentiles) is quantified for the same MSTID_{IDX} thresholds presented in Figure 13.
NRTK position degradation (in cm) for the user receiver PLAN with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
NRTK position degradation (in cm) for the user receiver MARE with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
NRTK position degradation (in cm) for the user receiver SBAR with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
From Figure 13 or the aforementioned tables, one can see that the reduction in the degradation is lower for the closest receiver, PLAN (approximately 10%), which is almost unaffected by MSTIDs: the degradation is maintained at approximately 5 cm in the 95th percentile and 8 cm in the 99th percentile.
This is not the case for the other two receivers, where the degradation is approximately 12 cm in the 99th percentile when MSTID_{IDX} is not taken into account. However, if the MSTID_{IDX} is considered, the degradation in the navigation solution is reduced up to 26% for MARE and 23% for SBAR. Furthermore, for the 95th percentile, SBAR and MARE reach levels of reduction similar to PLAN (the closest receiver). Therefore, MSTID_{IDX} can be used as an indicator that warns users of degradation in the navigation solution.
7 Conclusions
The present contribution analyses the impacts of MSTIDs on a network of permanent geodetic receivers located at midlatitudes. The methodology relies on three main steps: First, as detailed in Section 3, reference measurements are obtained by applying a novel method of fixing carrierphase ambiguities in the undifferenced measurements, avoiding the classic use of doubledifferenced carrierphase measurements, where only differential MSTIDs effects can be seen. In this way, we conduct this study independent of the ambiguity fixing success rate. Second, thanks to the previous carrierphase ambiguity fixing, a reference navigation solution based on the ionospherefree navigation combination (L _{IF}) is computed, making it possible to account for only the effects related to the ionosphere. Third, this reference navigation solution is compared with the navigation solution computed using either RTK or NRTK. In this way, we assess the impacts of the ionosphere effects on the 3D positioning.
The relationship between the positioning error and the MSTIDs is also shown, where the presence of an MSTID is a degrading factor for user positioning not only in RTK but also in NRTK. This degradation is related not only to the effect of the TID on the user measurements but also to the measurements of any of the reference receivers.
An MSTID index is defined and implemented as a tool to warn users about possible positioning degradation. The performance is tested with three stations located at 15, 50, and 79 km from the reference receiver. Using this index over a data set covering one year, we show that it is possible to obtain similar accuracies in the three baselines.
As a product related to this work, realtime estimates of MSTID_{IDX}, applied on more than one hundred of worldwide GNSS receivers, are computed and stored in the Warning and Mitigation Technologies for Travelling Ionospheric Disturbances Effects (TechTIDE) open repository, located at the following internet address: http://techtide.space.noa.gr/.
Acknowledgments
This work is sponsored by the European Union’s Horizon 2020 research and innovation programme under grant agreement 776011 (TechTIDE project) and 797461 (NAVSCIN), and the Spanish Ministry of Science, Innovation and Universities RETOS RTI2018094295BI00 (Programa Estatal de I + D + i Orientada a los Retos de la Sociedad). The authors acknowledge the use of data from the CATNET NRTK service from the Cartographic and Geologic Institute of Catalonia, Spain (Institut Cartogràfic i Geològic de Catalunya – ICGC), and data from the International GNSS Service (IGS). The editor thanks Yuichi Otsuka and an anonymous reviewer for their assistance in evaluating this paper.
References
 Alves P, Lachapelle G, Cannon ME, Liu J, Townsend B. 2001. Evaluation of multiplereference DGPS RTK using a large scale network. In: Proceedings of the National Technical Meeting of the Institute of Navigation, ION NTM/2001 (January 2001, Long Beach, USA), pp. 665–671. [Google Scholar]
 Cannon ME, Lachapelle G, Alves P, Fortes LP, Townsend B. 2001. GPS RTK Positioning using a regional reference network: Theory and results. In: Proceedings of the Global Navigation Satellite Systems Conference (May 2001, Seville, Spain), Global Navigation Satellite Systems. [Google Scholar]
 Collins P, Lahaye F, Héroux P, Bisnath S. 2008. Precise point positioning with ambiguity resolution using the decoupled clock model. In: Proceedings of the 21st International Technical Meeting of the Satellite Division of The Institute of Navigation, ION GNSS 2008 (September 2008, Savannah, USA). pp. 1315–1322. [Google Scholar]
 Ge M, Gendt G, Rothacher M, Shi C, Liu J. 2008. Resolution of GPS carrierphase ambiguities in precise point positioning (PPP) with daily observations. J Geodesy 82 (7 ): 389–399. https://doi.org/10.1007/s0019000701874. [CrossRef] [Google Scholar]
 HernándezPajares M, Juan JM, Sanz J, Colombo OL. 2000. Application of ionospheric tomography to realtime GPS carrierphase ambiguities Resolution, at scales of 400–1000 km and with high geomagnetic activity. Geophys Res Lett 27 ( 13 ): 2009–2012. https://doi.org/10.1029/1999GL011239. [CrossRef] [Google Scholar]
 HernándezPajares M, Juan JM, Sanz J. 2006. Mediumscale traveling ionospheric disturbances affecting GPS measurements: Spatial and temporal analysis. J Geophys Res Space Phys 111 ( A07S11 ). https://doi.org/10.1029/2005JA011474. [Google Scholar]
 Jacobson AR, Carlos RC, Massey RS, Wu G. 1995. Observations of traveling ionospheric disturbances with a satellitebeacon radio interferometer: Seasonal and local time behavior. J Geophys Res Space Phys 100 ( A2 ): 1653–1665. https://doi.org/10.1029/94JA02663. [CrossRef] [Google Scholar]
 Juan JM, Sanz J, RoviraGarcia A, GonzálezCasado G, VenturaTraveset J, Cacciapuoti L, Schoenemann E. 2020. A new approach to improve satellite clock estimates, removing the interday jumps. In: Proceedings of the 51st Annual Precise Time and Time Interval Systems and Applications Meeting of the US Institute of Navigation, pp. 279–301. [Google Scholar]
 Lachapelle G, Alves P. 2002. Multiple reference station approach: Overview and current research. J Glob Position Syst 1 ( 2 ): 133–136. [CrossRef] [Google Scholar]
 RoviraGarcia A, Juan JM, Sanz J, GonzálezCasado G, Ibáñez D. 2016. Accuracy of ionospheric models used in GNSS and SBAS: Methodology and analysis. J Geodesy 90 ( 3 ): 229–240. [CrossRef] [Google Scholar]
 RoviraGarcia A, IbáñezSegura D, OrúsPerez R, Juan JM, Sanz J, GonzálezCasado G. 2019. Assessing the quality of ionospheric models through GNSS positioning error: Methodology and results. GPS Solut 24 ( 1 ): 4. https://doi.org/10.1007/s102910190918z. [CrossRef] [Google Scholar]
 Sanz J, Juan Zornoza JM, HernándezPajares M. 2013. GNSS Data Processing. In: Volume I: Fundamentals and algorithms , ESA Communications, ESTEC, Noordwijk, The Netherlands, pp. 145–161. [Google Scholar]
 Seeber G. 2008. Satellite geodesy: Foundations, methods, and applications , Walter de Gruyter, Berlin, Germany. https://doi.org/10.1515/9783110200089. [Google Scholar]
 Talaya J, Bosch E. 1999. CATNET, a permanent GPS network with realtime capabilities. In: Proc. ION GPS99, 12th Int. Tech. Meeting of The Satellite Division of The US Institute of Navigation, pp. 14–17. [Google Scholar]
 Teunissen PJG. 1995. The leastsquares ambiguity decorrelation adjustment: A method for fast GPS integer ambiguity estimation. J Geodesy 70 : 65–82. [CrossRef] [Google Scholar]
 Wielgosz P, Kashani I, GrejnerBrzezinska D. 2005. Analysis of longrange network RTK during a severe ionospheric storm. J Geodesy 79 ( 9 ): 524–531. https://doi.org/10.1007/s001900050003y. [CrossRef] [Google Scholar]
Cite this article as: Timoté CC, Juan JM, Sanz J, GonzálezCasado G, RoviraGarcía A, et al. 2020. Impact of mediumscale traveling ionospheric disturbances on network realtime kinematic services: CATNET study case. J. Space Weather Space Clim. 10, 29.
All Tables
NRTK position degradation (in cm) for the user receiver PLAN with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
NRTK position degradation (in cm) for the user receiver MARE with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
NRTK position degradation (in cm) for the user receiver SBAR with respect to MSTID index thresholds (in TECU) by means of the 50th, 90th, and 95th percentiles.
All Figures
Fig. 1 CATNET NRTK network. Receivers in red: reference receivers; receivers in blue: user receivers. 

In the text 
Fig. 2 Fractional part of the L1 carrierphase ambiguities for all the satellites in view from the receivers SBAR on day 359 in 2017. The panels show the fixingambiguity process in which it is possible to fix the carrierphase ambiguities with confidence. 

In the text 
Fig. 3 Top panels: example of the ionospheric delays measured with the carrier phases (green) and pseudoranges (red) corresponding to the receiver MARE during DoY 310 in 2017. The top left panel (a) displays the raw measurements. The top right panel (b) corresponds to the measurements where the carrierphase ambiguities have been fixed. The ionospheric delay for the satellite G32, affected by a MSTID, is highlighted in blue. The bottom left panel (c) presents a more detailed picture of the MSTID event highlighted in panel (b). The bottom right panel (d) depicts the undifferenced STEC for the satellite G32 measured by receivers MARE (blue), GARR (red), and PLAN (cyan). The double differences of the ionospheric delays with respect to GARR are also depicted in green (MARE) and black (PLAN). 

In the text 
Fig. 4 Ionospherefree combination for the navigation solution at the PLAN receiver. From top to bottom: 3D error, up, north, and east components. 

In the text 
Fig. 5 Statistics for the 3D positioning error using the IF combination: the Yaxis depicts the probability of the 3D positioning error being greater than the value on the Xaxis. 

In the text 
Fig. 6 Ionospherefree (L _{IF}) solution in green, NRTK singlefrequency (L1) solution in blue, and RTK singlefrequency (L1) solution in red for three receivers: (a) PLAN, (b) MARE, and (c) SBAR. 

In the text 
Fig. 7 CCDFs of the 3D degradation in the positioning error for the receivers (a) PLAN, (b) MARE, and (c) SBAR during the 1year period of 2017–2018. 

In the text 
Fig. 8 MSTID effects on the second difference of the STEC for three receivers: GARR, SONA, and SBAR. 

In the text 
Fig. 9 MSTID index definition applied to three receivers on DoY 039 (February 6th) 2014: GARR, SONA, and SBAR. 

In the text 
Fig. 10 ∆^{2} STEC for the user receiver SBAR (blue) and reference receiver BEUD (green); the error in the ionospheric correction (red) is clearly correlated with the presence of a MSTID. 

In the text 
Fig. 11 Histogram of the error of the ionospheric corrections; 95th percentile (filled red boxes), 99th percentile (empty blue boxes). CDF of the ionospheric corrections (black line). The statistic is done for the three rover receivers: (a) PLAN, (b) MARE, and (c) SBAR. 

In the text 
Fig. 12 3D position errors and maximum MSTIDs amplitudes on days 298 and 325 of 2017 for the user receivers (a) PLAN, (b) MARE, and (c) SBAR. 

In the text 
Fig. 13 3D position degradation probabilities for user receivers: (a) PLAN, (b) MARE, and (c) SBAR. The solid lines represent the position solutions for RTK (red) and NRTK (blue) for all MSTIDs amplitudes. The NRTK solutions are also shown with MSTIDs amplitudes of less than 0.15 TECU (black) and with MSTIDs amplitudes of less than 0.10 TECU (green). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.