IONORING: Real-Time Monitoring of the Total Electron Content over Italy

: IONORING (IONOspheric RING) is a tool capable to provide the real-time monitoring and modeling of the ionospheric Total Electron Content (TEC) over Italy, in the latitudinal and longitudinal ranges of 35 ◦ N–48 ◦ N and 5 ◦ E–20 ◦ E, respectively. IONORING exploits the Global Navigation Satellite System (GNSS) data acquired by the RING (Rete Integrata Nazionale GNSS) network, managed by the Istituto Nazionale di Geoﬁsica e Vulcanologia (INGV). The system provides TEC real-time maps with a very ﬁne spatial resolution (0.1 ◦ latitude x 0.1 ◦ longitude), with a refresh time of 10 min and a typical latency below the minute. The TEC estimated at the ionospheric piercing points from about 40 RING stations, equally distributed over the Italian territory, are interpolated using locally (weighted) regression scatter plot smoothing (LOWESS). The validation is performed by comparing the IONORING TEC maps (in real-time) with independent products: (i) the Global Ionospheric Maps (GIM) - ﬁnal product- provided by the International GNSS Service (IGS), and (ii) the European TEC maps from the Royal Observatory of Belgium. The validation results are satisfactory in terms of Root Mean Square Error (RMSE) between 2 and 3 TECu for both comparisons. The potential of IONORING in depicting the TEC daily and seasonal variations is analyzed over 3 years, from May 2017 to April 2020, as well as its capability to account for the effect of the disturbed geospace on the ionosphere at mid-latitudes. The IONORING response to the X9.3 ﬂare event of September 2017 highlights a sudden TEC increase over Italy of about 20%, with a small, expected dependence on the latitude, i.e., on the distance from the subsolar point. Subsequent large regional TEC various were observed in response to related follow-on geomagnetic storms. This storm is also used as a case event to demonstrate the potential of IONORING in improving the accuracy of the GNSS Single Point Positioning. By processing data in kinematic mode and by using the Klobuchar as the model to provide the ionospheric correction, the resulting Horizontal Positioning Error is 4.3 m, lowering to, 3.84 m when GIM maps are used. If IONORING maps are used as the reference ionosphere, the error is as low as 2.5 m. Real-times application and services in which IONORING is currently integrated are also described in the conclusive remarks.


Introduction
The ability to monitor and model the environmental conditions in near-Earth space is of paramount importance for modern technological infrastructures vulnerable to a variety of electromagnetic phenomena of differing spatial and temporal scales, usually originating on the Sun and disrupting the near-Earth space environment via magnetosphere-ionosphere coupling [1]. These phenomena are included in the so-called "Space Weather". Particularly sensitive to the geospace conditions is the ionosphere at high and low latitudes with the formation of electron density irregularities of different scale sizes (from hundreds of km down to cm) affecting radio waves crossing them (see, e.g., [2][3][4]).
The mid-latitude ionosphere may also be a theatre of serious disturbing phenomena propagating from the higher latitudes and combining with the equatorial ionosphere software system comprises a method that is able to de-correlate the estimation of the Differential Code Biases (DCBs) and slant TEC (sTEC) in order to follow the ionospheric disturbances in near real-time. The ROB system has been tested with good results against IGS-GIM products and during several Space Weather events, including the Halloween storm which occurred in 2003 [23].
TEC maps over China and adjacent areas have been generated by using GNSS data from IGS and CMONOC (Crustal Movement Observation Network of China) networks and adopting a data assimilation method based on the Kalman filter [24]. The study highlighted the improvements that can be achieved by data assimilation in providing TEC maps over China with a spatial and temporal resolution of 1 • latitude x 1 • longitude x 5 min, respectively. Mendoza et al., 2019 [25] developed a near real-time TEC mapping over Central and South America, the Caribbean, and the Antarctic Peninsula, with a latency of about 10 min and a spatial resolution of 0.5 • latitude x 0.5 • longitude. The system uses GNSS data from BKG (Germany), IBGE (Brazil), IGN (Argentina), IGS, NASA (USA), SGM (Uruguay), UNAVCO (USA) showing a great collaborative effort for the common goal of the near real-time TEC mapping. Opperman et al., 2007 [26] tested the possibility to adopt the adjusted spherical cap harmonic analysis for TEC mapping over South Africa. The method was developed in the 1990s for mapping the critical frequency of the ionospheric F2 layer from ionosonde data over Europe [27,28]. Opperman and co-authors used TEC data from the Chief Directorate Surveys and Mapping (CDSM), which operated a network of 40 geodetic dual-frequency Global Positioning System (GPS) receivers, and ionosonde data for comparison. They concluded that the method, suitably optimized, could be applied for near real-time TEC mapping over South Africa. Recently, the spherical cap harmonic analysis has been applied over Australia for the near real-time mapping and modeling of TEC, provided by the Australian Regional GPS Network (ARGN). The model makes use of an artificial neural network (ANN) to generate TEC values to fill the grid where experimental data are not available [29]. The model is able to map the TEC over Australia with a Root Mean Square Error (RMSE) of 2-3 TECU and is capable of reproducing its seasonal and hourly features. This paper deals with the real-time monitoring and modeling of TEC over Italy that extends over a geographic area between mid-high latitude (48 • N) to mid-low latitude (35 • N) in a quite confined range of longitude centered at 12 • E. The very fine spatial resolution (0.1 • x 0.1 • latitude, longitude) gives the opportunity to reproduce the ionosphere dynamics under quiet and disturbed conditions as well as its diurnal, seasonal, and solar flux-related variability. In Section 2, the data provided by the dense RING network and the approach adopted to calibrate TEC are detailed. Section 3 deals with the validation of the IONORING real-time product against TEC from GIMs (global and final products) and ROB maps (Europe, real-time product). The performance of IONORING is then evaluated in Section 4, by analyzing the TEC diurnal, seasonal, solar cycle, and latitudinal variations between 2017 and 2020, including also the response to the space weather disturbances that occurred in early September 2017. As a further demonstration of the potential of IONOR-ING, Section 4 also contains an example of how IONORING can support ionospheric error mitigation in GNSS single-frequency positioning. Conclusions are given in Section 5, where the potential of IONORING to support science and end-user services is also highlighted.

Data and Method
Ionospheric RING (IONORING) leverages the data provided by the Rete Integrata Nazionale GNSS (RING, http://ring.gm.ingv.it/, accessed on 12 July 2021) (INGV RING working group, 2016) owned by the Istituto Nazionale di Geofisica e Vulcanologia. RING represents an important geodetic research infrastructure consisting currently of a distributed network over the Italian territory of more than 200 GNSS sensors (black circles in Figure 1). Almost all the stations are in the free field, thanks to the low-power consumption of the satellite transmission that can be run with solar panels. This significantly reduces the multipath-related errors in the GNSS observables. All the stations are equipped with Remote Sens. 2021, 13, 3290 4 of 21 a professional GPS receiver and a choke ring antenna. These features make the RING receivers particularly suitable for ionospheric studies, as also demonstrated in the recent literature [12,30,31]. Among all the RING receivers, some have been selected to feed IONORING algorithm, by considering the trade-off between a homogenous coverage of the ionospheric sector and a minimum latency of the IONORING products availability. In order to cover the westernmost part of the considered region, the receivers in Ajaccio and Cagliari (red-dotted blue circles in Figure 1) managed by IGS are also included. In Figure 1, the receivers located in Montereale (mter, 42.51 • N 13.21 • E, green dotted) and in INGV headquarter (ingr, 41.83 • N 12.51 • E) are also reported. They are not included in IONORING, but used for validation and positioning application purposes, as described later in the text. IONORING input data are GPS carrier phase observables on L1 and L2 frequencies (1575.42 MHz and 1227.60 MHz, respectively), Clear Acquisition (CA) code observables on L1, and Y-code observables on L2 that are broadcasted at a 1-second rate through Networked Transport of RTCM via Internet Protocol (NTRIP) by the selected receivers (red dots in Figure 1). Additionally, navigational information (i.e., orbital parameters for GPS satellites) are retrieved by NTRIP streaming from the considered receivers or, alternatively, by the IGS navigational information streaming.
From code and phase observables on L1 and L2 frequencies, we derive ionospheric sTEC, that is the integrated electron density along the ray-path connecting satellite and receiver. This is based on the fact that the ionospheric effects on GNSS signal propagation include range delay, being proportional to sTEC according to the following formula: where sTEC is expressed in TEC units (TECu, 1 TECu = 10 16 electrons/m 2 ), I is the ionosphere-induced delay in length unit, and α is a conversion factor to obtain length units from TECu. Thus, the TEC retrieval is obtained by calculating the geometry-free linear combination, that is defined as where subscripts 1 and 2 refer to L1 and L2 GPS frequencies, respectively, while λ refers to the corresponding wavelength, N is the ambiguity on carrier phase measurements, "arc" refers to continuous carrier phase observations (i.e., for which the product λN can be considered constant), τ is the frequency-dependent delay induced by the receiver (R) and satellite (S) hardware, c is the speed of light in vacuum, and ε is the noise in length units on the carrier phase observation combination. From Equations (1) and (2), the following can be derived: α are the so-called "inter-frequency biases" (IFBs) for carrier phase observations related to the receiver and the transmitter, is the ambiguity for the L 1 − L 2 combination and ε L = f 2 α ε is the carrier phase noise in TECu. Analogously, the geometry-free linear combination of the code measurements can be obtained as: where b R and b S are the corresponding IFBs, often called "differential code biases" (DCBs), due to the receiving and the transmitting hardware, and ε P is the noise on code measurements. We let the reader notice that in Equation (4) there is not the ambiguity term present in Equation (2). For every continuous arc of observations, by subtracting the two new observables (one for code measurements and the other one for phase measurements), the average difference between carrier phase and code observables along a single arc can be obtained: (5) in which the noise associated with the phase measurements is neglected [Braasch, 1996]. By subtracting Equation (5) from (3), the following can be obtained: being the carrier-phase ionospheric observable leveled to the code-delay ionospheric observable. The term ε P arc is the effect of the noise and the multipath along the arc on carrier-phase observations, as described in Ciraolo et al., 2007 [32]. One possible approach to estimate sTEC without neglecting ε P arc consists in evaluating a single bias for each arc, instead of different biases for receiver and satellite. Thus, Equation (6) can be written as: where β arc is the sum of receiver and satellite biases (b R + b S ) and multipath and other errors with non-zero mean. This is the equation from which a calibrated TEC can be obtained, by using the method described by Ciraolo et al., 2007 [32]. Hereafter, we refer to this method as the "Gg calibration technique", or simply "Gg technique". In the specific, sTEC values are mapped as a two-dimensional surface by means of the classical thin shell method as where vTEC(φ 1 , φ 2 ) is the unknown describing a surface in the reference frame defined by a couple (φ 1 , φ 2 )· over the thin shell (bi-dimensional) and χ is the angle formed by the ray-path and the perpendicular to the shell at the ionospheric pierce point (IPP) [33,34].
In the case of the Gg technique, the reference frame is local time (LT) and modified dip latitude (Modip) proposed by Rawer et al., 1981 [35]. Then, the vTEC is expanded as a polynomial linear in LT and of the fourth-order in Modip. As a consequence, Equation (8) can be written as a linear relationship: L arc = sec χ ∑ n c n p n (LT, Modip)+β arc (9) where p n is the term of the polynomial and c n the corresponding coefficient. Equation (9) can be solved via the standard least squares method. Although this calibration technique is robust and widely used, especially for midlatitude and low-latitude studies [12,[36][37][38][39][40][41], it has the intrinsic limitation of not being operating in real-time. This is because it assumes the availability of complete arcs of observations. To overcome this limitation and for the purposes of the real-time operation of IONORING, we evaluate β arc for each station-satellite couple and for a given day by computing its mean value on the previous 3 days.
The algorithm for β arc calculation is here summarized: 1. Download of RINEX files for all the GNSS stations for the previous 3 days; 2.
Evaluation of β arc for each satellite-receiver pair for the previous 3 days by applying Gg calibration technique; 3.
Calculation of the mean value over the previous 3 days for each satellite-receiver pair (β arc ); 4.
Creation of a lookup table containing m x nβ arc values, being m the total number of RING stations, n the total number of operating GPS satellites.
The algorithm is run daily so that the lookup table is updated every 24 h. Previous studies have shown that the receiver and satellite DCBs variations between consecutive days are smaller than 0.5 ns (<1 TECu) for GPS satellites and smaller than 1 ns (<2 TECu) for the GPS stations (see, e.g., [42]). This allows considering β arc as stable within 24 h. Assumingβ arc = b R + b S , Equation (4) can be rewritten as: in which the only remaining unknown is ε P . As per Equation (8), the values of vTEC are obtained by projecting sTEC to the vertical at the IPP and by considering the shell height at 350 km. The IPP coordinates are calculated in the WGS84 reference frame by using the positions of the receivers (derived from RINEX files) and of the satellites. The latter are evaluated by leveraging on the navigation information broadcasted by each receiver in its NTRIP streaming or alternately, from the IGS navigational information streaming.
To produce a vTEC map in which the impact of ε P is mitigated, vTEC values estimated at every IPP from all considered satellite-receiver pairs are interpolated through a nonparametric local regression technique named LOWESS [43], which is an interpolation method that includes local data in the computation of the smoothing curve. LOWESS has been proven to be effective for regional TEC mapping purposes (see, e.g., [38]). The interpolation algorithm is first applied on the whole vTEC dataset to identify and reject outliers defined as: where vTEC(φ 1 , φ 2 ) is the value of the smoothed vTEC curve evaluated at the IPP, identified by its coordinates (φ 1 , φ 2 ), while vTEC(φ 1 , φ 2 ) is the actual vTEC at the IPP and RMSE indicates the root mean square error of the dataset. After the outlier rejection, the LOWESS interpolation is reapplied on the clean dataset to provide the final vTEC map that is provided in terms of the smoothed vTEC curve evaluated on a regular grid of 0.1 • latitude

Results
To validate IONORING, we compare the vTEC maps with three independent products by considering 3 years of accumulated data covering the period 1 May 2017 to 30 April 2020. The products are: (i) the vTEC derived from the Montereale (mter, 42.51 • N 13.21 • E, not included in IONORING) receiver (green dot in Figure 1) by applying the Gg calibration technique (post-processing) and described in Section 2; (ii) the vTEC maps over Europe provided by ROB and available at ftp://gnss.oma.be (accessed on 17 August 2021) [23] and (iii) the GIMs provided as final combined solution product ("igsg") by the IGS and available at cddis.gsfc.nasa.gov (accessed on 17 August 2021) [17]. To provide a qualitative comparison with other similar regional TEC products, far from being exhaustive, Figure 5 reports the RMSE resulting from the validation of ION-ORING against ROB and IGS maps (red) and that referred to different regional or local TEC products reported in selected recent literature. Even if a direct comparison among different techniques is out of the scope of this paper, the aim of this figure is to illustrate how the difference of IONORING with the validation datasets is in agreement with what was reported by other authors, validating independent regional TEC mapping techniques. It is worth mentioning that the RMSE values reported in the literature are obtained by considering different inhomogeneous datasets, covering different ionospheric sectors, seasons, and helio-geophysical conditions.   Probability distribution of the difference between IONORING and ROB vTEC and the corresponding Gaussian fit (red), whose expression and parameters are also reported in the textbox (panel c). Boxplots of the hourly (panel d) and monthly (panel e) differences between IONORING and ROB vTEC. In the boxplots, the red line indicates the median, and the bottom and top edges of the box indicate the 25th and 75th percentiles, respectively, while the range marked by the dashed lines indicates the extrema.
By following the same approach adopted to compare with ROB maps, Figure 4 reports the results of the comparison between vTEC from IGS maps and vTEC from IONORING. The IGS GIM maps (final combined solution, "igsg") are provided at a global level with a spatial resolution of 2.5 • latitude x 5 • longitude and every 2 h with a latency of about 11 days [17]. For this comparison, IONORING maps have been down-sampled to 2 h.
The summary of the validation results is reported in Table 1, in which the following parameters are provided: the angular coefficient (a), intercept (b) and coefficient of determination (R 2 ) of the linear regression (panels a of Figures 2-4), the RMSE of the vTEC differences, defined as: where vTEC i IONORI NG is the IONORING vTEC at the time t i and vTEC i REF is the vTEC at the time t i of reference products (REF= Gg, ROB, IGS), as described above. In addition, remind here again that data from mter are not included in IONORING maps and that vTEC over such receiver is retrieved by applying the Gg calibration technique as it is (i.e., in post-processing). Thus, the fit of the comparison with IONORING vTEC over the mter location, featured by a = 1, b~0, reveals: (i) the goodness of the LOWESS interpolation method and (ii) the validity of the assumption that DCBs variability over 24 h is negligible. It is also worth noticing that the differences have no dependence on the year, time of the day, and season (panels b, d, and e of Figure 2, respectively).    Concerning the comparison with ROB and IGS datasets, the coefficient of determination suggests an almost ideal correlation (R 2 = 0.95 and R 2 = 0.96, respectively). However, a mean difference of −1.9 ± 0.8 TECu with the ROB dataset and −2.8 ± 0.8 TECu with the IGS dataset is present.
To provide a qualitative comparison with other similar regional TEC products, far from being exhaustive, Figure 5 reports the RMSE resulting from the validation of IONORING against ROB and IGS maps (red) and that referred to different regional or local TEC products reported in selected recent literature. Even if a direct comparison among different techniques is out of the scope of this paper, the aim of this figure is to illustrate how the difference of IONORING with the validation datasets is in agreement with what was reported by other authors, validating independent regional TEC mapping techniques. It is worth mentioning that the RMSE values reported in the literature are obtained by considering different inhomogeneous datasets, covering different ionospheric sectors, seasons, and helio-geophysical conditions. median, and the bottom and top edges of the box indicate the 25 th and 75 th percentiles, respectively, while the range marked by the dashed lines indicates the extrema.

Discussion
The performance of IONORING has been evaluated by characterizing the daily, seasonal, and latitudinal variations [44] in the same period considered for the validation (1 May 2017 to 30 April 2020). Additionally, the storm-time behavior [45] is also investigated, by considering the geomagnetic events that occurred in September 2017. Figure 6 reports the daily variation of the mean vTEC for Italian low (37°N, blue), mid (42°N, orange) and high (47°N, yellow) latitudes (longitude is 12°E) for the period 1 May 2017 to 30 April 2020. The error bars represent the ±1σ spread around the mean. As expected, the higher the latitude, the lower the TEC, as the latitudinal dependence of the ionization is dependent on the solar zenith angle. It is interesting to note how the peak TEC for low-latitude is at 12:00-13:00 UT (13:00-14:00 LT), while for mid-and high-latitudes it is at 11:00 UT (12:00 LT). This may be due to the fact that the Italian low-latitudes are sensitive to the effect of the increase of ionization in correspondence with the expected position of the northern crest of the EIA. In fact, while at mid-latitudes and at F2-layer peak heights the ionization is strongly driven by the zenith angle that peaks at noon [46], at low-latitudes the EIA is formed mainly from the removal of plasma from around the In red, the results reported in this paper (Figures 3 and 4).

Discussion
The performance of IONORING has been evaluated by characterizing the daily, seasonal, and latitudinal variations [44] in the same period considered for the validation (1 May 2017 to 30 April 2020). Additionally, the storm-time behavior [45] is also investigated, by considering the geomagnetic events that occurred in September 2017. Figure 6 reports the daily variation of the mean vTEC for Italian low (37 • N, blue), mid (42 • N, orange) and high (47 • N, yellow) latitudes (longitude is 12 • E) for the period 1 May 2017 to 30 April 2020. The error bars represent the ±1σ spread around the mean. As expected, the higher the latitude, the lower the TEC, as the latitudinal dependence of the ionization is dependent on the solar zenith angle. It is interesting to note how the peak TEC for low-latitude is at 12:00-13:00 UT (13:00-14:00 LT), while for mid-and high-latitudes it is at 11:00 UT (12:00 LT). This may be due to the fact that the Italian low-latitudes are sensitive to the effect of the increase of ionization in correspondence with the expected position of the northern crest of the EIA. In fact, while at mid-latitudes and at F2-layer peak heights the ionization is strongly driven by the zenith angle that peaks at noon [46], at low-latitudes the EIA is formed mainly from the removal of plasma from around the equator by the upward E × B drift [8], whose peak is expected at 14:00 LT [47], i.e., at a later time with respect to the local noon. To better characterize the average daily vTEC variation provided by IONORING, Figure 7 shows

September 2017 Geomagnetic Storm
The geomagnetic storm that occurred in September 2017 is one of the most investigated of the 24th solar cycle (see e.g., [48][49][50][51][52][53]) and was featured by the strongest flare event in a decade. It started on 6 September 2017 at 11:53 UT, peaked at 12:02 UT, and was classified as an X9.3. Solar irradiance increases rapidly during solar flare events, peaking withiñ 10-30 min, and creating abrupt increases of the ionization in the upper atmosphere [Qian et al., 2019]. The occurrence around UT noon enables to study the TEC response over Italy to such a flare, as, at that UT time, it is located close to the subsolar point (6.48 • N, 0.41 • E), at which the largest TEC gradients driven by solar EUV flux enhancements are expected to take place (see, e.g., [49,53,54]). Within this scope, the top panel of Figure 9 reports the mean value of TEC over Italy for the period 6 September 2017 to 9 September 2017 (red line) and the 27-days (before the X9.3 flare occurrence) running median of mean TEC over Italy (gray dashed line). The mean TEC is calculated by averaging the values over the whole map. The bottom panel of the same figure reports the difference between actual and median TEC (the curves in the top panel, black line) and the disturbance storm time (Dst) index (green line). Blue dashed lines in both panels indicate the starting time of the X9.3 flare. We remind the reader that the TEC-related values in Figure 9 are provided with a 10-minutes cadence. Right after the arrival of the X9.3 flare, a sudden average TEC increase is recorded, passing from 13.6 TECu at 12:00 UT to 15.7 TECu at 12:10 UT, and then it reaches 16.2 TECu at 12:20 UT. Starting from 7 September, two geomagnetic storms occurred (Qian et al., 2019) related to the arrival of Coronal Mass Ejections associated with the X9.3 flare. The Dst index characterizes the geomagnetic storms, with the first and large commencing~20UT on 7th September and the second storm commencing~10UT on 8th September. The ionospheric response over Italy is featured by a typical response of the mid-latitude ionosphere [45]. Specifically, there are positive ionospheric conditions after the flare EUV enhancement on the 6th, with a large daytime enhancement on the 7th before the geomagnetic storms. This may be caused by the transport of ionization at F2 heights into the region as a result of the initial ionospheric disturbance from the flare. The initial geomagnetic storm occurs during the local night and does not cause a TEC enhancement, but the second smaller storm in the following morning of the 8th causes a daytime enhancement positive phase ionospheric storm, though interestingly smaller than that of the 7th before the onset of the geomagnetic storms. Enhanced TEC lasts until 18:00 UT on the 8th and are followed by a negative phase ionospheric disturbance that is maintained until at least 00UT on the 10th, at the end of the period in Figure 9. To further characterize the TEC features depicted by IONORING under effective flaring activity, Figure 10 shows the maps for the times between 11:50 UT and 12:30 UT on 6 September 2017. Panels b and c report the maps right before and right after the arrival of the X9.3 flare. The effect of the flare is an increase in TEC all over the considered area. Highlighting the latitudinal dependence of the flare-induced TEC increase, Figure 11 reports the maps (∆TEC median ) of the percentage difference between median vTEC considering 27 days before the X9.3 flare on 6 September 2017 and vTEC maps at 12:00 UT (panel a), i.e., before the flare, and 12:10 UT (panel b), i.e., after the flare, and the difference ∆TEC f lare = ∆TEC median | 12:10 − ∆TEC median | 12:00 (panel c). By subtracting ∆TEC median | 12:00 , we are removing the deviation from the median behavior due to effects not related to the flare. Such effects are assumed to be constant between 12:00 UT and 12:10 UT, in agreement with what was reported by Cesaroni et al., 2015 [38] for low-latitude ionosphere, which is featured by larger temporal TEC gradients with respect to mid-latitude ones. Hence, ∆TEC flare represents the flare net effect measured as a percentage of extra ionization. In Figure 11c, ∆TEC flare values are of the order of 20% in the mid-and low latitudes, while of about 10% in the northernmost part of the map. This is reasonably due to a twofold effect. First, the dependence of the flare effects on the distance from the subsolar point (see, e.g., [54]). In addition, from the bottom panel of Figure 9, the flare-induced extra ionization lasts for about 4 h, in agreement with what was reported for the Halloween and Bastille Day storms by Tsurutani et al., 2005 [54]. Such duration is compatible with the hypothesis that the extra ionization is mainly due to the EUV effect on the atomic oxygen in the ionospheric F2 layer, which recombines with a time scale of hours [54]. The atomic oxygen concentration decreases with the latitude, likely resulting in the meridional TEC enhancement distribution shown in Figure 11c. Deviations from this behavior are found in the easternmost part of the map in Figure 11c, where a quiet uniform increase of 20% is present. This is likely due to border effects introduced by the interpolation technique and by the different local times between maps in Figure 11a,b, from which ∆TEC flare is calculated, which affect the photoionization in the easternmost part.

Application on Positioning
As an example of application, we test the use of IONORING maps in single-frequency positioning. To the scope, we perform Single Point Positioning (SPP) [55], exploiting the data acquired by the ingr receiver (41.83 • N, 12.51 • E, yellow dot in Figure 1) and processing them in kinematic mode. We remind the reader again that the ingr receiver is not included in IONORING. IONORING-based positioning is compared with the ones obtained by using, as ionospheric information, the Klobuchar model [53] and the IGS GIM final combined solution product [17]. It is worth mentioning that this positioning exercise is performed in order to show a possible application of IONORING products and not for validation purposes in the positioning domain that will be the focus of a companion paper. In fact, validation in the positioning domain requires particular attention in not mixing ionospheric induced errors with the ones introduced by other factors, such as clock and orbit errors, measurement noise, pseudo-range multipath, evaluation metric, and outlier's contributions [56]. For this reason, validation on the positioning domain is usually performed by applying the Precise Point Positioning technique that is able to mitigate some of the errors mentioned above by considering precise orbit and clock information. In this paper, positioning results are obtained by running gLAB software [57], a tool suite developed under a European Space Agency (ESA) Contract by the research group of Astronomy and Geomatics (gAGE) from the Universitat Politecnica de Catalunya (UPC) to process and analyze GNSS data.
To the scope, carrier-smoothing of code pseudo-ranges are used, locations of the satellites are calculated by finding the optimal broadcasting ephemeris according to the observing epoch, satellite clock error is calculated by the broadcasting clock error parameters, tropospheric delay correction is calculated using the University of New Brunswick (UNB) model [58], and relativity correction and earth rotation correction are also calculated with models. We perform the positioning under quiet and disturbed ionospheric conditions. As a disturbed day, we select 6 September 2017, being featured by the already discussed X9.3 flare. As a quiet day, we select 25 September 2017, which is suggested in Alfonsi et al., 2021 [48] as suitable to characterize the quiet ionospheric behavior, being featured by quiet geomagnetic conditions, low auroral activity, and seasonal and solar flux conditions similar to those under which the X9.3 flare occurred. Figure 12 Table 2 summarizes the North, East, and Horizontal Position Error (HPE) and the improvement ratio (IR) with respect to the use of the Klobuchar model, which we define as: The HPE is of the order of meters in all the cases, but the improvement ratio reveals an increase of the accuracy on the horizontal position of 25% (91%) and 67% (25%) when using IONORING (GIM) during quiet and disturbed days, respectively. We remind the reader that GIM are available with a latency of about 11 days, while IONORING is a real-time product, so, despite the fact that GIMs perform better than IONORING in quiet conditions, IONORING is a valuable added value in real-time positioning applications. During disturbed conditions, IONORING seems to perform better than Klobuchar and GIM.

Conclusions
In this paper, we present IONORING: a real-time tool for ionospheric monitoring and modeling over Italy. IONORING exploits the data provided by the RING network managed by INGV, which has been conceived for geodetic purposes. IONORING provides maps of vertical TEC with a cadence of 10 min, a spatial resolution of 0.1 • latitude x 0.1 • longitude, and a latency of less than a minute. The reliability of IONORING has been proved by comparing the vTEC maps produced over a 3-year period (May 2017-April 2020) with independent vTEC measurements (Gg calibration technique, ROB vTEC maps, and GIM). Such comparison has shown a very good accuracy of IONORING products being the RMSEs equal to 1.2 TECu, 2.0 TECu, and 2.9 TECu for Gg vTEC, ROB, and IGS final combined maps, respectively.
In addition, IONORING capability to depict the seasonal, daily, and latitudinal variations as well as the dependence on the solar activity and to depict the ionospheric response to the September 2017 geomagnetic storm has been shown.
The daily variation of mean vTEC reveals a peak of 9.6 ± 0.8 TECu at 11:00 UT (12:00 LT) for higher latitudes and a peak of 14.6 ± 1.6 TECu for lower latitudes at later times, i.e., 13:00 UT (14:00 LT). This highlights the influence of the northern crest of the EIA on vTEC at Italian lower latitudes. The monthly running average of vTEC between May 2017 and April 2020 reproduces the expected seasonal variation (maximizing vTEC in May/June) and presents a decay of about 1 TECu per year, matching the behavior of the solar activity.
The storm-time behavior is also well depicted by IONORING. Specifically, the positive and negative phases of the geomagnetic storm that occurred in September 2017 and the effect of the X9.3 flare that started at 11:53 UT on 6 September are well characterized. In particular, the mean vTEC enhancement over Italy induced by the flare is 2.6 TECu at 12:20 UT, resulting in a net TEC increase ranging between 10% (higher latitudes) and 20% (lower latitudes).
As an example of application, IONORING has been tested as external ionospheric information single frequency standard point positioning showing an improvement in the performance in terms of horizontal error of 25% (91%) and 67% (25%) when using IONORING (GIM) during a quiet and disturbed day, respectively. Despite having been tested on a single day, there is a promising indication that IONORING performs better than the GIM final product under disturbed conditions. Under quiet conditions, IONORING presents a smaller improvement ratio than the GIM final product. However, the latter cannot be used in real-time, likely indicating that IONORING provides a valuable added value in real-time positioning applications. Such features are currently under investigation for further assessments on larger statistics and against also IGS rapid products, for real-time positioning purposes.
Proven such capabilities, IONORING is a precious tool to support either research and services. In the recent past, it has been used to characterize the ionospheric effect of a moderate storm in the mid-low latitude boundary region [12]. Moreover, it is one of the algorithms currently running in the Ionospheric Prediction Service (IPS, https://ionospheric-prediction.jrc.ec.europa.eu, accessed on 13 July 2021) [59] as a nowcasting product and as input to a tool dedicated to the detection and estimation of the main characteristics of Medium-Scale Traveling Ionospheric Disturbances (MSTIDs). Additionally, the maps from IONORING are part of the product portfolio of the PECASUS consortium (http://pecasus.eu/, accessed on 13 July 2021), that is one of the three global centers providing space weather advisories according to International Civil Aviation Organization (ICAO) regulations.
Author Contributions: Conceptualization, all authors; methodology, C.C. and L.S.; software, C.C. and L.S.; validation, C.C. and L.S.; data curation, C.C.; writing-original draft preparation, L.S. and G.D.F.; writing-review and editing, all authors. All authors have read and agreed to the published version of the manuscript.

Data Availability Statement:
The data presented in this study are openly available in eSWua webportal at https://doi.org/10.13127/eswua/tec, accessed on 13 July 2021.