GPS slant total electron content accuracy using the single layer model under different geomagnetic regions and ionospheric conditions

The use of observations from the Global Positioning System (GPS) has significantly impacted the study of the ionosphere. As it is widely known, dual-frequency GPS observations can provide very precise estimation of the slant Total Electron Content (sTEC—the linear integral of the electron density along a ray-path) and that the precision level is bounded by the carrier-phase noise and multi-path effects on both frequencies. Despite its precision, GPS sTEC estimations can be systematically affected by errors in the estimation of the satellites and receivers by Inter-Frequency Biases (IFB) that are simultaneously determined with the sTEC. Thus, the ultimate accuracy of the GPS sTEC estimation is determined by the errors with which the IFBs are estimated. This contribution attempts to assess the accuracy of IFBs estimation techniques based on the single layer model for different ionospheric regions (low, mid and high magnetic latitude); different seasons (summer and winter solstices and spring and autumn equinoxes); different solar activity levels (high and low); and different geomagnetic conditions (quiet and very disturbed). The followed strategy relies upon the generation of a synthetic data set free of IFB, multi-path, measurement noise and of any other error source. Therefore, when a data set with such properties is used as the input of the IFB estimation algorithms, any deviation from zero on the estimated IFBs should be taken as indications of the errors introduced by the estimation technique. The truthfulness of this assessment work is warranted by the fact that the synthetic data sets resemble, as realistically as possible, the different conditions that may happen in the real ionosphere. The results of this work show that during the high solar activity period the accuracy for the estimated sTEC is approximately of ±10 TECu for the low geomagnetic region and of ±2.2 TECu for the mid-latitude. During low solar activity the accuracy can be assumed to be in the order of ±2 TECu. For the geomagnetic high-disturbed period, the results show that the accuracy is degraded for those stations located over the region where the storm has the strongest impact, but for those stations over regions where the storm has a moderate effect, the accuracy is comparable to that obtained in the quiet period.     

Climatological study of ionospheric irregularities over the European mid-latitude sector with GPS

High-frequency variability of the ionosphere, or irregularities, constitutes the main threat for real-time precise positioning techniques based on Global Navigation Satellite Systems (GNSS) measurements. Indeed, during periods of enhanced ionospheric variability, GNSS users in the field—who cannot verify the integrity of their measurements—will experience positioning errors that can reach several decimeters, while the nominal accuracy of the technique is cm-level. In the frame of this paper, a climatological analysis of irregularities over the European mid-latitude region is presented. Based on a 10 years GPS dataset over Belgium, the work analyzes the occurrence rate (as a function of the solar cycle, season and local time) as well as the amplitude of ionospheric irregularities observed at a single GPS station. The study covers irregularities either due to space weather events (solar origin) or of terrestrial origin. If space weather irregularities are responsible for the largest effects in terms of ionospheric error, their occurrence rate highly depends on solar activity. Indeed, the occurrence rate of ionospheric irregularities is about 9 % during solar maximum, whereas it drops to about 0 % during medium or low solar activity periods. Medium-scale ionospheric disturbances (MSTIDs) occurring during daytime in autumn/winter are the most recurrent pattern of the time series, with yearly proportions slightly varying with the solar cycle and an amplitude of about 10 % of the TEC background. Another recurrent irregularity type, though less frequent than MSTIDs, is the noise-like variability in TEC observed during summer nighttime, under quiet geomagnetic conditions. These summer nighttime irregularities exhibit amplitudes ranging between 8 and 15 % of the TEC background.     

Application of SWACI products as ionospheric correction for single-point positioning: a comparative study

In Global Navigation Satellite Systems (GNSS) using L-band frequencies, the ionosphere causes signal delays that correspond with link related range errors of up to 100 m. In a first order approximation the range error is proportional to the total electron content (TEC) of the ionosphere. Whereas this first order range error can be corrected in dual-frequency measurements by a linear combination of carrier phase- or code-ranges of both frequencies, single-frequency users need additional information to mitigate the ionospheric error. This information can be provided by TEC maps deduced from corresponding GNSS measurements or by ionospheric models. In this paper we discuss and compare different ionospheric correction methods for single-frequency users. The focus is on the comparison of the positioning quality using dual-frequency measurements, the Klobuchar model, the NeQuick model, the IGS TEC maps, the Neustrelitz TEC Model (NTCM-GL) and the reconstructed NTCM-GL TEC maps both provided via the ionosphere data service SWACI (http://swaciweb.dlr.de) in near real-time. For that purpose, data from different locations covering several days in 2011 and 2012 are investigated, including periods of quiet and disturbed ionospheric conditions. In applying the NTCM-GL based corrections instead of the Klobuchar model, positioning accuracy improvements up to several meters have been found for the European region in dependence on the ionospheric conditions. Further in mid- and low-latitudes the NTCM-GL model provides results comparable to NeQuick during the considered time periods. Moreover, in regions with a dense GNSS ground station network the reconstructed NTCM-GL TEC maps are partly at the same level as the final IGS TEC maps.     

Analysis of inversion errors of ionospheric radio occultation

The retrieved electron density profile of ionospheric radio occultation (RO) simulation data can be compared with the background model value during the simulation and the inversion error can be obtained exactly. This paper studies the inversion error of ionospheric RO through simulation. The sources of the inversion errors are analyzed. The impacts of measurement errors, such as the errors in phase measurements and satellite orbits, are very small and can be neglected. The approximation of straight-line propagation introduces errors at the height of the F1 layer under solar maximum condition. The spherical symmetry approximation of the electron density distribution is found to be the main source of the inversion error. The statistical results reveal some characteristics of the inversion errors. (1) The relative error increases with enhanced solar activity. (2) It is larger in winter than in equinox season, and it is smallest in summer. (3) For all seasons, it is smaller at middle latitude than at other latitudes. (4) For all seasons and geomagnetic latitudes, it is smaller at daytime than at other times. The NmF2 of the ROs from COSMIC are compared with the measurements of ionosondes, and the relative differences show the same dependencies on season, geomagnetic latitude and local time, as the relative errors of the simulated ionospheric ROs.

Potentialities of multifrequency ionospheric correction in Global Navigation Satellite Systems

The first-order ionospheric error is reduced in the dual-frequency Global Navigation Satellite Systems (GNSS). In this paper, the possibility of eliminating ionospheric higher-order errors in the multifrequency GNSS is explored. Since the second-order error associated with the geomagnetic field effect on the refractive index can be eliminated in dual-frequency measurements, we explore the possibility of eliminating third-order errors in triple-frequency GNSS in view of phase scintillations. A connection between the possibility of improving the multifrequency GNSS accuracy and diffraction effects in radio signal propagation through the randomly inhomogeneous ionosphere is shown. The numerical simulation has revealed that the systematic, residual ionospheric error is considerably reduced when we pass on from dual-frequency to triple-frequency measurements. The change in the residual error variance during such a transition depends however on the relationship between the inner scale of the turbulent spectrum of ionospheric irregularities and the Fresnel radius. Given the inner scale larger than the Fresnel radius, not only the systematic error, but also the standard deviation reduces when we pass on from dual-frequency to triple-frequency measurements. Otherwise, when the Fresnel radius exceeds the inner scale, the variance increases with increasing number of frequencies in use.     

多GNSS监测下中国区域电离层格网模型可用性分析

多系统全球卫星导航系统(GNSS)的出现为天基增强系统(SBAS)电离层格网模型的性能提升提供了可能,但多系统GNSS测量对电离层格网模型性能提升是有条件的. 为此,利用中国区域GPS观测模拟分析了多系统GNSS测量对中国区域电离层格网模型可用性的影响. 结果表明:多系统GNSS测量可有效提高电离层格网模型的覆盖范围. 中国南方地区存在低纬赤道电离异常(EIA)现象,严重影响SBAS电离层格网模型实现性能,单纯增加GNSS测量不能有效应对低纬电离异常现象影响. 中国北方地区电离层延迟变化平缓,在多系统GNSS测量情况下可以考虑减少地面监测站数量,仍能保持系统原有性能.      

Improving the GNSS positioning stochastic model in the presence of ionospheric scintillation

Ionospheric scintillations are caused by time- varying electron density irregularities in the ionosphere, occurring more often at equatorial and high latitudes. This paper focuses exclusively on experiments undertaken in Europe, at geographic latitudes between ~50°N and ~80°N, where a network of GPS receivers capable of monitoring Total Electron Content and ionospheric scintillation parameters was deployed. The widely used ionospheric scintillation indices S4 and s_j{\sigma_{\varphi}} represent a practical measure of the intensity of amplitude and phase scintillation affecting GNSS receivers. However, they do not provide sufficient information regarding the actual tracking errors that degrade GNSS receiver performance. Suitable receiver tracking models, sensitive to ionospheric scintillation, allow the computation of the variance of the output error of the receiver PLL (Phase Locked Loop) and DLL (Delay Locked Loop), which expresses the quality of the range measurements used by the receiver to calculate user position. The ability of such models of incorporating phase and amplitude scintillation effects into the variance of these tracking errors underpins our proposed method of applying relative weights to measurements from different satellites. That gives the least squares stochastic model used for position computation a more realistic representation, vis-a-vis the otherwise ‘equal weights’ model. For pseudorange processing, relative weights were com- puted, so that a ‘scintillation-mitigated’ solution could be performed and compared to the (non-mitigated) ‘equal weights’ solution. An improvement between 17 and 38% in height accuracy was achieved when an epoch by epoch differential solution was computed over baselines ranging from 1 to 750 km. The method was then compared with alternative approaches that can be used to improve the least squares stochastic model such as weighting according to satellite elevation angle and by the inverse of the square of the standard deviation of the code/carrier divergence (sigma CCDiv). The influence of multipath effects on the proposed mitigation approach is also discussed. With the use of high rate scintillation data in addition to the scintillation indices a carrier phase based mitigated solution was also implemented and compared with the conventional solution. During a period of occurrence of high phase scintillation it was observed that problems related to ambiguity resolution can be reduced by the use of the proposed mitigated solution.     

A two-step ionospheric modeling algorithm considering the impact of GLONASS pseudo-range inter-channel biases

The Global Navigation Satellite System presents a plausible and cost-effective way of computing the total electron content (TEC). But TEC estimated value could be seriously affected by the differential code biases (DCB) of frequency-dependent satellites and receivers. Unlike GPS and other satellite systems, GLONASS adopts a frequency-division multiplexing access mode to distinguish different satellites. This strategy leads to different wavelengths and inter-frequency biases (IFBs) for both pseudo-range and carrier phase observations, whose impacts are rarely considered in ionospheric modeling. We obtained observations from four groups of co-stations to analyze the characteristics of the GLONASS receiver P1P2 pseudo-range IFB with a double-difference method. The results showed that the GLONASS P1P2 pseudo-range IFB remained stable for a period of time and could catch up to several meters, which cannot be absorbed by the receiver DCB during ionospheric modeling. Given the characteristics of the GLONASS P1P2 pseudo-range IFB, we proposed a two-step ionosphere modeling method with the priori IFB information. The experimental analysis showed that the new algorithm can effectively eliminate the adverse effects on ionospheric model and hardware delay parameters estimation in different space environments. During high solar activity period, compared to the traditional GPS + GLONASS modeling algorithm, the absolute average deviation of TEC decreased from 2.17 to 2.07 TECu (TEC unit); simultaneously, the average RMS of GPS satellite DCB decreased from 0.225 to 0.219 ns, and the average deviation of GLONASS satellite DCB decreased from 0.253 to 0.113 ns with a great improvement in over 55%.     

Ionospheric tomography using GNSS: multiplicative algebraic reconstruction technique applied to the area of Brazil

Experimental analysis was performed using multiplicative algebraic reconstruction technique (MART) to map the ionosphere over Brazil. Code and phase observations from the global navigation satellite system (GNSS) together with the international reference ionosphere (IRI) enabled the estimation of ionospheric profiles and total electron content (TEC) over the entire region. Twenty-four days of data collected from existing ground-based GNSS receivers during the recent solar maximum period were used to analyze the performance of the MART algorithm. The results were compared with four ionosondes. It was demonstrated that MART estimated the electron density peak with the same degree of accuracy as the IRI model in regions with appropriate geometrical coverage by GNSS receivers for tomographic reconstruction. In addition, the slant TEC, as estimated with MART, presented lower root-mean-square error than the TEC calculated by ionospheric maps available from the International GNSS Service (IGS). Furthermore, the daily variations of the ionosphere were better represented with the algebraic techniques, compared to the IRI model and IGS maps, enabling a correlation of the elevation of the ionosphere at higher altitudes with the equatorial ionization anomaly intensification. The tomographic representations also enabled the detection of high vertical gradients at the same instants in which ionospheric irregularities were evident.     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

Calibration errors on experimental slant total electron content (TEC) determined with GPS

The Global Positioning System (GPS) has become a powerful tool for ionospheric studies. In addition, ionospheric corrections are necessary for the augmentation systems required for Global Navigation Satellite Systems (GNSS) use. Dual-frequency carrier-phase and code-delay GPS observations are combined to obtain ionospheric observables related to the slant total electron content (sTEC) along the satellite-receiver line-of-sight (LoS). This observable is affected by inter-frequency biases [IFB; often called differential code biases (DCB)] due to the transmitting and the receiving hardware. These biases must be estimated and eliminated from the data in order to calibrate the experimental sTEC obtained from GPS observations. Based on the analysis of single differences of the ionospheric observations obtained from pairs of co-located dual-frequency GPS receivers, this research addresses two major issues: (1) assessing the errors translated from the code-delay to the carrier-phase ionospheric observable by the so-called levelling process, applied to reduce carrier-phase ambiguities from the data; and (2) assessing the short-term stability of receiver IFB. The conclusions achieved are: (1) the levelled carrier-phase ionospheric observable is affected by a systematic error, produced by code-delay multi-path through the levelling procedure; and (2) receiver IFB may experience significant changes during 1 day. The magnitude of both effects depends on the receiver/antenna configuration. Levelling errors found in this research vary from 1.4 total electron content units (TECU) to 5.3 TECU. In addition, intra-day vaiations of code-delay receiver IFB ranging from 1.4 to 8.8 TECU were detected.     

A preliminary evaluation of the performance of multiple ionospheric models in low- and mid-latitude regions of China in 2010–2011

Ionospheric delay is a dominant error source in Global Navigation Satellite System (GNSS). Single-frequency GNSS applications require ionospheric correction of signal delay caused by the charged particles in the earth’s ionosphere. The Chinese Beidou system is developing its own ionospheric model for single-frequency users. The number of single-frequency GNSS users and applications is expected to grow fast in the next years in China. Thus, developing an appropriate ionospheric model is crucially important for the Chinese Beidou system and worldwide single-frequency Beidou users. We study the performance of five globally accessible ionospheric models Global Ionospheric Map (GIM), International Reference Ionosphere (IRI), Parameterized Ionospheric Model (PIM), Klobuchar and NeQuick in low- and mid-latitude regions of China under mid-solar activity condition. Generally, all ionospheric models can reproduce the trend of diurnal ionosphere variations. It is found that all the models have better performances in mid-latitude than in low-latitude regions. When all the models are compared to the observed total electron content (TEC) data derived from GIM model, the IRI model (2012 version) has the best agreement with GIM model and the NeQuick has the poorest agreement. The RMS errors of the IRI model using the GIM TEC as reference truth are about 3.0–10.0 TECU in low-latitude regions and 3.0–8.0 TECU in mid-latitude regions, as observed during a period of 1 year with medium level of solar activity. When all the ionospheric models are ingested into single-frequency precise point positioning (PPP) to correct the ionospheric delays in GPS observations, the PIM model performs the best in both low and mid-latitudes in China. In mid-latitude, the daily single-frequency PPP accuracy using PIM model is ~10 cm in horizontal and ~20 cm in up direction. At low-latitude regions, the PPP error using PIM model is 10–20 cm in north, 30–40 cm in east and ~60 cm in up component. The single-frequency PPP solutions indicate that NeQuick model has the lowest accuracy among all the models in both low- and mid-latitude regions of China. This study suggests that the PIM model may be considered for single-frequency GNSS users in China to achieve a good positioning accuracy in both low- and mid-latitude regions.     

Determination of the optimized single-layer ionospheric height for electron content measurements over China

The ionosphere effective height (IEH) is a very important parameter in total electron content (TEC) measurements under the widely used single-layer model assumption. To overcome the requirement of a large amount of simultaneous vertical and slant ionospheric observations or dense “coinciding” pierce points data, a new approach comparing the converted vertical TEC (VTEC) value using mapping function based on a given IEH with the “ground truth” VTEC value provided by the combined International GNSS Service Global Ionospheric Maps is proposed for the determination of the optimal IEH. The optimal IEH in the Chinese region is determined using three different methods based on GNSS data. Based on the ionosonde data from three different locations in China, the altitude variation of the peak electron density (hmF2) is found to have clear diurnal, seasonal and latitudinal dependences, and the diurnal variation of hmF2 varies from approximately 210 to 520 km in Hainan. The determination of the optimal IEH employing the inverse method suggested by Birch et al. (Radio Sci 37, 2002. doi: 10.1029/2000rs002601) did not yield a consistent altitude in the Chinese region. Tests of the method minimizing the mapping function errors suggested by Nava et al. (Adv Space Res 39:1292–1297, 2007) indicate that the optimal IEH ranges from 400 to 600 km, and the height of 450 km is the most frequent IEH at both high and low solar activities. It is also confirmed that the IEH of 450–550 km is preferred for the Chinese region instead of the commonly adopted 350–450 km using the determination method of the optimal IEH proposed in this paper.     

Higher order ionospheric effects in precise GNSS positioning

With the increasing number of precise navigation and positioning applications using Global Navigation Satellite Systems (GNSS) such as the Global Positioning System (GPS), higher order ionospheric effects and their correction become more and more important. Whereas the first-order error can be completely eliminated by a linear combination of dual- frequency measurements, the second- and third-order residual effects remain uncorrected in this approach. To quantify the second-order residual effect, a simple formula has been derived for GNSS users in Germany. Our proposed correction algorithm reduces the second-order effects to a residual error of fractions of 1 mm up to 2 mm at a vertical total electron content level of 10¹⁸ electrons/m² (100 TECU), depending on satellite azimuth and elevation angles. The correction formula can be implemented in real-time applications as it does not require the knowledge of the geomagnetic field or the electron density distribution in the ionosphere along the signal path. It is expected that the correction will enable more accurate positioning using the line-of-sight carrier-phase measurements.     

Mitigating the impact of ionospheric cycle slips in GNSS observations

Processing of data from global navigation satellite systems (GNSS), such as GPS, GLONASS and Galileo, can be considerably impeded by disturbances in the ionosphere. Cycle-slip detection and correction thus becomes a crucial component of robust software. Still, dealing with ionospheric cycle slips is not trivial due to scintillation effects in both the phase and the amplitude of the signals. In this contribution, a geometry-based approach with rigorous handling of the ionosphere is presented. A detailed analysis of the cycle-slip correction process is also tackled by examining its dependence on phase and code noise, non-dispersive effects and, of course, the ionosphere. The importance of stochastic modeling in validating the integer cycle-slip candidates is emphasized and illustrated through simulations. By examining the relationship between ionospheric bias and ionospheric constraint, it is shown that there is a limit in the magnitude of ionospheric delay variation that can be handled by the cycle-slip correction process. Those concepts are applied to GNSS data collected by stations in northern Canada, and show that enhanced cycle-slip detection can lead to decimeter-level improvements in the accuracy of kinematic PPP solutions with a 30-s sampling interval. Cycle-slip correction associated with ionospheric delay variations exceeding 50 cm is also demonstrated, although there are risks with such a procedure and these are pointed out.     

Higher-order ionospheric effects in GPS time and frequency transfer

In high-precision geodetic time and frequency transfer, which requires precise modeling of code and carrier phase GPS data, the ionosphere-free combinations P ₃ and L ₃ of the codes and carrier phases, made on the two GPS frequencies, are used to remove the first-order ionospheric effect. We quantify the impact of the residual second- and third-order ionospheric effects on geodetic time and frequency transfer solutions for continental and intercontinental baselines. All time transfer computations are done using the ATOMIUM software, developed at the Royal Observatory of Belgium. In order to avoid contamination by some imperfect modeling of the second- and third-order ionospheric effects in the satellite clock products, only single-difference, common-view processing is used, based on code and carrier phase measurements. The results are shown for weak and strong solar activity, as well as for particular epochs of ionospheric storms. Second-order ionospheric delays can lead to corrections up to about 10 ps in the common-view clock solution of intercontinental baselines with very different longitudes. However, realistic values of the geomagnetic field in the ionosphere are required to assess the amplitude of second-order ionospheric effects in time and frequency transfer during an ionospheric storm.     

基于GNSS的电离层模型研究进展

在利用全球导航卫星系统GNSS进行精密定位和导航时,电离层延迟误差是影响其精度和准确度的主要误差源之一,故对电离层模型研究至关重要。本文将电离层模型分成了经典电离层模型和现代电离层模型,并对经典电离层模型进行了比较,重点介绍了目前全世界电离层模型的研究热点、存在问题及研究方向。     

Ionospheric correction using NTCM driven by GPS Klobuchar coefficients for GNSS applications

Global Navigation Satellite Systems (GNSS) require mitigation of ionospheric propagation errors because the ionospheric range errors might be larger than tens of meters at the zenith direction. Taking advantage of the frequency-dispersive property of ionospheric refractivity, the ionospheric range errors can be mitigated in dual-frequency applications to a great extent by a linear combination of carrier phases or pseudoranges. However, single-frequency GNSS operations require additional ionospheric information to apply signal delay or range error corrections. To aid single-frequency operations, the global positioning system (GPS) broadcasts 8 coefficients as part of the navigation message to drive the ionospheric correction algorithm (ICA) also known as Klobuchar model. We presented here an ionospheric correction algorithm called Neustrelitz TEC model (NTCM) which can be used as complementary to the GPS ICA. Our investigation shows that the NTCM can be driven by Klobuchar model parameters to achieve a significantly better performance than obtained by the mother ICA algorithm. Our research, using post-processed reference total electron content (TEC) data from more than one solar cycle, shows that on average the RMS modeled TEC errors are up to 40% less for the proposed NTCM model compared to the Klobuchar model during high solar activity period, and about 10% less during low solar activity period. Such an approach does not require major technology changes for GPS users rather requires only introducing the NTCM approach a complement to the existing ICA algorithm while maintaining the simplicity of ionospheric range error mitigation with an improved model performance.     

层析模型在GNSS探测电离层中的研究进展

在利用全球导航卫星系统GNSS(GPS、GLONASS、Galileo等卫星定位系统)进行精密定位和导航时,电离层延迟误差是影响其精度和准确度的主要误差源之一,故对电离层模型研究至关重要。本文介绍了传统电离层模型的缺点,重点阐述了不同的电离层层析方法,同时分析了现在电离层层析方法所存在的主要问题,最后介绍了目前全世界电离层层析模型的研究展望。     

Mitigation of higher order ionospheric effects on GNSS users in Europe

Current dual-frequency GPS measurements can only eliminate the first-order ionospheric term and may cause a higher-order range bias of several centimeters. This research investigates the second-order ionospheric effect for GNSS users in Europe. In comparison to previous studies, the electron density profiles of the ionosphere/plasmasphere are modeled as the sum of three Chapman layers describing electron densities of the ionospheric F₂, F₁ and E layers and a superposed exponential decay function describing the plasmasphere. The International Geomagnetic Reference Field model is used to calculate the geomagnetic field vectors at numerous points along the incoming ray paths. Based on extended simulation studies, we derive a correction formula to compute the average value of the longitudinal component of the earth’s magnetic field along the line-of-sight as a function of geographic latitude and longitude, and geometrical parameters such as elevation and azimuth angles. Using our correction formula in conjunction with the total electron content (TEC) along the line-of-sight, the second-order ionospheric term can be corrected to the millimeter level for a vertical TEC level of 10¹⁸ electrons/m².