Tomographic Imaging of Ionospheric Plasma Bubbles Based on GNSS and Radio Occultation Measurements

Total electron content measurements given by the global navigation satellite system (GNSS) have successfully presented results to capture the signatures of equatorial plasma bubbles. In contrast, the correct reproduction of plasma depletions at electron density level is still a relevant challenge for ionospheric tomographic imaging. In this regard, this work shows the first results of a new tomographic reconstruction technique based on GNSS and radio-occultation data to map the vertical and horizontal distributions of ionospheric plasma bubbles in one of the most challenging conditions of the equatorial region. Twenty-three days from 2013 and 2014 with clear evidence of plasma bubble structures propagating through the Brazilian region were analyzed and compared with simultaneous observations of all-sky images in the 630.0 nm emission line of the atomic oxygen. The mean rate of success of the tomographic method was 37.1%, being more efficient near the magnetic equator, where the dimensions of the structures are larger. Despite some shortcomings of the reconstruction technique, mainly associated with ionospheric scintillations and the weak geometry of the ground-based GNSS receivers, both vertical and horizontal distributions were mapped over more than 30◦ in latitude, and have been detected in instances where the meteorological conditions disrupted the possibility of analyzing the OI 630 nm emissions. Therefore, the results revealed the proposed tomographic reconstruction as an efficient tool for mapping characteristics of the plasma bubble structures, which may have a special interest in Space Weather, Spatial Geodesy, and Telecommunications.


Introduction
Ionospheric plasma bubbles are large-scale structures of density depletion in the Earth's ionosphere.It is well known that the plasma bubbles are generated by plasma instability mechanisms in the lower F region, and grow aligned with the magnetic field up to ~1500 km over the equatorial region [1].Typical ionosonde observations have suggested the irregular plasma density of the bubbles as the main factor causing the spread-F traces in equatorial ionograms [2] and forcing satellite transionospheric signals to fluctuate wildly in phase and amplitude [3,4].Such degradations may cause interference, fading in radio frequency communications and disrupting many applications.Therefore, characterizing the ionospheric plasma bubble extensions and temporal evolutions has been the subject of continued study for decades [5][6][7][8][9].
Temporal evolution and spatial development of the plasma bubbles have been studied extensively using different techniques.Several works have conducted an analysis based on ionosondes measurement, such as Rastogi [10] by reporting seasonal variations of Equatorial Spread F, Saito and Maruyama [11] by studying the day-to-day variability of plasma bubbles, and Wang et al. [12] by using concurrent ionosonde observations with other instruments to correlate the occurrence of scintillations with plasma bubbles.Topside sounders and satellite-based measurements, which can be connected to ground-based observations [13], have also been used to map the plasma bubble distributions, such as for mapping the activity [14] and probability [15] of equatorial spread F, to describe the plasma bubbles' morphology [16][17][18] and specific distributions during storm events [19], and for showing their association with ionospheric irregularities [20] and plasma blobs [21].Furthermore, radio scintillation measurements enabled the monitoring of plasma bubbles in the level of signal fluctuations [22], which makes it possible to capture the occurrence of plasma bubbles [23] and their signatures [24], even in the Brazilian area [25,26].Additionally, all-sky optical instruments providing two-dimensional horizontal images have been presented as relevant tools for mapping instantaneous coverage of the plasma bubble structures [27], as well as for capturing climatological occurrences [28], connections with neutral winds [29], and relations with small-and medium-scale structures [30].
Although many techniques have been developed, a major difficulty in characterizing the plasma bubbles remains due to the limitations of the actual instruments for monitoring the plasma depletions in a sufficiently large area.For instance, the frequently used all-sky optical instruments cover approximately 1600 km of airglow measurements, that is, they are restricted to ~15 • of latitude and longitude in the equatorial region.Such dimensions are sufficiently good to capture the longitudinal variability of the plasma bubbles, which is around 400 km.However, the latitudinal extension of the plasma bubbles can reach distances of 5400 km in the equatorial region, since the geomagnetic conjugate points are located at approximately ±20 • from the equator.In order to overcome such limitations, recent works have made observations of plasma bubbles over a large area using Total Electron Content (TEC) measurements derived from ground-based global navigation satellite system (GNSS) stations [31,32].The ground-based GNSS observations using more than 130 receivers over South America are not dependent on weather conditions and were capable of capturing plasma depletions with a sufficient spatial resolution at TEC level.Therefore, simultaneous analysis of airglow and GNSS measurements has been presented as a useful tool to monitor plasma bubble structures in a more detailed form.
Previous works have shown clear evidence of equatorial plasma bubbles at TEC level [31,32] by comparing the airglow images with the Vertical TEC (VTEC) parameter obtained with GNSS observations.However, these VTEC values have intrinsic shortcomings, due to the approximation of the ionosphere as a thin shell with a constant height at about 350-450 km above the Earth's surface.Consequently, the accuracy of the plasma bubble boundary is limited, and the resulting maps are applicable only in analyzing the horizontal distributions.Otherwise, the GNSS signal intersects many altitudes of the ionosphere, by allowing its three-dimensional reconstruction by using tomographic algorithms [33].The use of tomographic reconstruction techniques for recovering horizontal and vertical distributions of the plasma bubbles is a topic of great interest to the community involved in Space Weather, Spatial Geodesy, and Telecommunications, since much about the plasma bubbles' morphology is not fully understood [34].For instance, the three-dimensional representations could be a useful tool to describe the horizontal and vertical distribution of the plasma bubbles, as well as to contribute to the explanation of their dynamics and evolution.
An ionospheric tomography method was developed by Prol and Camargo [35], with the main purpose of imaging the ionosphere in the Brazilian region.The experimental results provided good estimations of the critical frequency and TEC when compared to conventional ionospheric models.
A further analysis remained to be carried out in order to verify the efficiency of tomography for imaging the ionospheric dynamics over the region.Considering that the spatial variability and evolution of ionospheric plasma bubble irregularities is a topic of major interest in the Brazilian region [36], we present here the possibilities of a new tomographic method for imaging plasma depletions at electron density level.As a result, we show, for the first time, the efficiency of tomographic reconstruction techniques in detecting plasma bubbles over a systematic analysis of several days, using, as reference, simultaneous observations of all-sky images from 630.0 nm atomic oxygen emission, i.e., OI 630 nm all-sky images.Additionally, since the tomographic reconstruction provides three-dimensional representations, the performance of the tomographic algorithm is investigated for analyzing the plasma bubble vertical structures.In this direction, Section 2 shows the proposed method.Section 3 presents the performance of the proposed method to map the horizontal distributions, and Section 4 shows the possibilities of mapping the vertical distributions.Then, Section 5 presents the conclusions of the tomographic method application to map the equatorial plasma bubbles in the Brazilian scenario.

Ionospheric Tomography Method
The mathematical formulation of the tomographic reconstruction method was developed in two stages.The first stage uses the climatological model proposed by Prol et al. [37] to build the background ionosphere based on global Radio-Occultation (RO) observations.The second refers to a slight modification of the Multiplicative Algebraic Reconstruction Technique (MART) for an iterative update of the background electron density field, using regional observations of TEC.Both methods were applied using data spanning over twenty-three days from 2013 and 2014 with clear evidence of plasma bubble structures propagating through the Brazilian region, and the following subsections show the mathematical aspects in details.

Ionospheric Background
The background ionosphere was developed in terms of a climatological model, using RO ionospheric profiles derived from the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC)/Formosa Satellite Mission 3 (COSMIC/FORMOSAT-3) mission.Six parameters were estimated to describe climatological patterns of the ionospheric profiles: Electron density at the peak (N m ), peak height (h m ), scale height of the topside referred to h m (H 0 ), scale height gradient at the topside (∂H/∂h), thickness of the bottomside (B 0 ), and scale factor of the bottomside (B 1 ).The peak parameters h m and N m were directly retrieved from ionospheric profiles processed by the University Corporation for Atmospheric Research (UCAR).On the other hand, a least-square fitting was performed by the authors for the topside and bottomside estimation.The following equations were used for the topside: with: and the following equations for the bottomside: being: where n e represents the observed electron densities at the height h for the topside (n t e ) and bottomside (n b e ).
A thirty-day sliding window using RO measurements was defined for the global representation.A spatial bilinear interpolation was carried out for each parameter (h m , N m , H 0 , ∂H/∂h, B 0 , and B 1 ) in terms of local time (LT) and magnetic latitude.The resulting global maps were converted to the Earth-fixed reference frame by projecting the LT maps along with the geomagnetic lines.As a result, two-dimensional grids (2D) were generated with a spatial resolution of 2 • in latitude by 4 • in longitude, and a temporal resolution of 6 min in Universal Time (UT).Therefore, it was possible to interpolate h m , N m , H 0 , ∂H/∂h, B 0 , and B 1 for a specific regular grid in UT, and then to calculate the electron density at any height in the Brazilian grid in order to create the background for regional tomography.

Regional Tomography
Once the ionospheric background is calculated, regional tomography is performed using ground-based dual-frequency receivers from the continuously operating GNSS stations throughout South America.We used data from RBMC (Brazilian Network for Continuous GNSS Monitoring), IGS (International GNSS Service), CALIBRA (Countering GNSS high Accuracy applications Limitations due to Ionospheric disturbances in Brazil), LISN (Low latitude Ionospheric Sensor Network), and RAMSAC (Red Argentina de Monitoreo Satelital Continuo) to derive the tomographic images.Figure 1 shows the location of the GNSS stations operating in the Day of Year (DOY) 3 of 2014.
Remote Sens. 2018, 10, x FOR PEER REVIEW 4 of 16 A thirty-day sliding window using RO measurements was defined for the global representation.A spatial bilinear interpolation was carried out for each parameter (ℎ , , , / ℎ, , and ) in terms of local time (LT) and magnetic latitude.The resulting global maps were converted to the Earth-fixed reference frame by projecting the LT maps along with the geomagnetic lines.As a result, two-dimensional grids (2D) were generated with a spatial resolution of 2° in latitude by 4° in longitude, and a temporal resolution of 6 min in Universal Time (UT).Therefore, it was possible to interpolate ℎ , , , / ℎ, , and for a specific regular grid in UT, and then to calculate the electron density at any height in the Brazilian grid in order to create the background for regional tomography.

Regional Tomography
Once the ionospheric background is calculated, regional tomography is performed using ground-based dual-frequency receivers from the continuously operating GNSS stations throughout South America.We used data from RBMC (Brazilian Network for Continuous GNSS Monitoring), IGS (International GNSS Service), CALIBRA (Countering GNSS high Accuracy applications Limitations due to Ionospheric disturbances in Brazil), LISN (Low latitude Ionospheric Sensor Network), and RAMSAC (Red Argentina de Monitoreo Satelital Continuo) to derive the tomographic images.Figure 1 shows the location of the GNSS stations operating in the Day of Year (DOY) 3 of 2014.Ionospheric tomography was performed using a temporal resolution of 6 minutes and a horizontal resolution of 1° in longitude and 1° in latitude in order to cover the longitudes from 91.5°W to 18°W and the latitudes from 51.5°S to 22°N.The vertical resolution of the 3D grid was constructed with distinct step sizes, varying from 50 km up to 500 km, in order to reduce the number of cells and improve the time process.Each cell of the 3D grid was filled with the initial electron density given by the background.Then, a slightly different method than the conventional MART algorithm was applied in order to iteratively update the initial electron density of each cell using TEC measurements.The general tomographic system can be expressed as below [38]: Ionospheric tomography was performed using a temporal resolution of 6 minutes and a horizontal resolution of 1 • in longitude and 1 • in latitude in order to cover the longitudes from 91.5 • W to 18 • W and the latitudes from 51.5 • S to 22 • N. The vertical resolution of the 3D grid was constructed with distinct step sizes, varying from 50 km up to 500 km, in order to reduce the number of cells and improve the time process.Each cell of the 3D grid was filled with the initial electron density given by the background.Then, a slightly different method than the conventional MART algorithm was applied in order to iteratively update the initial electron density of each cell j using TEC measurements.The general tomographic system can be expressed as below [38]: where y consists of the TEC measurements, x is composed by the unknown vector of electron density, A is the design matrix composed by the path length elements of the GNSS signals inside the boundaries of the voxels, and ε stands for the measurement noises and discretization error.In general, this problem involves many parameters due to the dimensions of the ionosphere.Additionally, the limited number of GNSS viewing angles makes the tomographic inversion an ill-conditioned problem.Therefore, summing up the large dimensions of the matrices and the measurement noise affecting the ill-conditioned system, an iterative process is typically used to perform the estimation procedure.In this regard, the proposed MART algorithm is used with the following equations: (6) and: where x k+1 j represents the j member of the electron density vector x at the iteration k + 1, A ij is the path length of the GNSS signal i inside the boundaries that intersect cell j, A max is the largest path length of the respective signal i, the inner product A i , x k produces the background TEC, w is a weighting parameter that controls the convergence of the algorithm, and y i and y v i are the observed GNSS-based TEC and VTEC, respectively, in which k varies with a stepsize of 2 due to VTEC updates.The iterative process cycles through all ray paths 70 times with a weighting parameter of w =0.2, both parameters being empirically obtained.The slight difference from the conventional MART was defined in this procedure by the incorporation of vertical TEC observations in addition to the slant TEC.In this direction, updates are performed using A v ij and A v max as the vertical distance of the path lengths of A ij and A max , respectively.
TEC observations of Equation ( 5) were obtained with the high precise calibration process developed by Prol et al. [39], and the VTEC observations of Equation (6) were derived by VTEC i = TEC i /M i , where M i stands for the standard mapping function.While the TEC updates correspond to a set of cells in the slant direction of the GNSS signal, the VTEC estimations correspond to the vertical column of voxels above and below the ionospheric pierce point corresponding to the 450 km shell height.This strategy improves the description of the ionospheric horizontal features in regions with lack of slant TEC, without significantly affecting the TEC estimation in regions with a high amount of slant observations.Therefore, the use of both measurements (slant TEC and VTEC) generates a tomographic reconstruction with an increased number of illuminated cells, which is an important aspect for mapping the plasma bubble structures.

Horizontal Distribution of Plasma Bubbles
The validation of the tomographic method as a tool to map the horizontal distributions of the plasma bubbles was performed comparing the reconstructed VTEC maps with images derived from all-sky optical instruments.The next subsection demonstrates the feasibility of mapping the horizontal distributions of plasma bubbles through a set of comparative examples between the reconstructed VTEC plasma bubbles with the airglow images.Then, Section 3.2 presents the validation of the proposed method by means of a systematic analysis over several days.For both experiments, three all-sky optical instruments, with 180 • wide-angle fisheye lens, were used to detect the red line emission of atomic oxygen at 630 nm in the ionosphere.Two of them are located in the equatorial region, at São João do Cariri (7.4 • S, 36.5 • W, dip latitude 11.11 • S) and Boa Vista (2.8 • N, 60.7 • W, dip latitude 15.14 • N), and the other is located in the southern crest of the Equatorial Ionization Anomaly (EIA), at Cachoeira Paulista (22.7 • S, 45.0 • W, dip latitude 19.7 • S).The OI 630 nm images are freely provided by INPE (Instituto Nacional de Pesquisas Espaciais) and cover a latitudinal extension and a longitudinal extension of 1280 × 1280 km, with a pixel resolution of 2.5 km for the height of 250 km.The airglow images were recorded with a temporal resolution of six minutes, following the same temporal resolution of the VTEC maps.

Feasibility of Mapping Plasma Bubbles by the Tomographic Reconstruction
Figure 2 shows an example of the VTEC map retrieved from the tomographic algorithm (left panel) and the corresponding OI 630 nm images (right panel) in DOY 3 of 2014 at 1.9 h UT (Universal Time).Two well-defined VTEC depletions can be seen aligned with the geomagnetic field lines around 21 and 24 h in Local Time (LT), which are also coincident with signatures of the plasma bubble in the airglow images.At such instants, the plasma that was uplifted to higher altitudes, as a result of the post-sunset enhancement in the equatorial E × B drift, increased the growth rate of the spread-F plasma instabilities.Thus, the well-developed large-scale plasma bubble irregularities expanding to higher latitudes generated depletions observed at VTEC level.Since the emission rate of OI 630.0 nm nightglow is primarily a consequence of a charge exchange involving the ion O + and O 2 and because there is a lower oxygen density (O 2 ) at higher altitudes, the elevation of the F layer generates a large region of ionization depletion seen as darkened flux tubes [40].As expected, we have observed VTEC depletions aligned with a decrease in the intensity of the OI 630 nm images given by the dark brands.Plasma depletions in the VTEC maps were more evident in Cachoeira Paulista (CHPI) and Boa Vista (BOAV) than in São João do Cariri (SJCI).Indeed, SJCI is located at the edge of the reconstructed area, where the geometrical coverage of the GNSS signal is too weak to efficiently perform ionospheric tomography.It is also important to notice that few observations are used in the north-eastern region, above 10 • N, and in the eastern side of 50 • W. In fact, there are some GNSS signals crossing the topside ionosphere, but few observations near the peak height.The VTEC representation of this part is, therefore, mainly associated to the background model and, as expected, some artifacts are observed in the reconstructed images due to the high TEC variation in the transition between the background and the regional reconstruction.
Remote Sens. 2018, 10, x FOR PEER REVIEW 6 of 16 for the height of 250 km.The airglow images were recorded with a temporal resolution of six minutes, following the same temporal resolution of the VTEC maps.

Feasibility of Mapping Plasma Bubbles by the Tomographic Reconstruction
Figure 2 shows an example of the VTEC map retrieved from the tomographic algorithm (left panel) and the corresponding OI 630 nm images (right panel) in DOY 3 of 2014 at 1.9 h UT (Universal Time).Two well-defined VTEC depletions can be seen aligned with the geomagnetic field lines around 21 and 24 h in Local Time (LT), which are also coincident with signatures of the plasma bubble in the airglow images.At such instants, the plasma that was uplifted to higher altitudes, as a result of the post-sunset enhancement in the equatorial E × B drift, increased the growth rate of the spread-F plasma instabilities.Thus, the well-developed large-scale plasma bubble irregularities expanding to higher latitudes generated depletions observed at VTEC level.Since the emission rate of OI 630.0 nm nightglow is primarily a consequence of a charge exchange involving the ion O + and O2 and because there is a lower oxygen density (O2) at higher altitudes, the elevation of the F layer generates a large region of ionization depletion seen as darkened flux tubes [40].As expected, we have observed VTEC depletions aligned with a decrease in the intensity of the OI 630 nm images given by the dark brands.Plasma depletions in the VTEC maps were more evident in Cachoeira Paulista (CHPI) and Boa Vista (BOAV) than in São João do Cariri (SJCI).Indeed, SJCI is located at the edge of the reconstructed area, where the geometrical coverage of the GNSS signal is too weak to efficiently perform ionospheric tomography.It is also important to notice that few observations are used in the north-eastern region, above 10°N, and in the eastern side of 50°W.In fact, there are some GNSS signals crossing the topside ionosphere, but few observations near the peak height.The VTEC representation of this part is, therefore, mainly associated to the background model and, as expected, some artifacts are observed in the reconstructed images due to the high TEC variation in the transition between the background and the regional reconstruction.Although the plasma bubble was not so evident in SJCI at DOY 3/2014, the plasma depletions can be clearly observed in the VTEC maps presented for DOY 359/2013 in Figure 3.It is possible to see the VTEC depletion aligned with the dark brands of the airglow images due to the larger longitudinal dimensions of the plasma bubble of DOY 359/2013 in comparison with the plasma bubble of DOY 3/2014.Thus, even with problems of coverage in the tomographic reconstruction at SJCI, it is possible to detect plasma bubbles in case the dimensions of the ionospheric structure are large enough to be captured.
Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 16 bubble of DOY 3/2014.Thus, even with problems of coverage in the tomographic reconstruction at SJCI, it is possible to detect plasma bubbles in case the dimensions of the ionospheric structure are large enough to be captured.VTEC depletions are not only observed with a similar location and dimensions to the airglow observations, but they also follow the plasma bubble orientation.For instance, Figure 4 shows that the azimuth of the VTEC depletion mapped in CHPI at DOY 3/2014 is changing accordingly with the azimuth of the airglow dark intensities, which also happened for other instruments and days.In all cases with plasma bubble detections, a longitudinal discrepancy can also be observed between the VTEC depletion and the darker brand of the airglow emissions.A similar discrepancy was observed in [31], which might indicate that the bubble is tilted westward with the increasing height [3,41].Additionally, Figure 5 shows a collection of west-east slices of the reconstructed VTEC maps as a form of a keogram.This is presented in order to compare the obtained reconstructions with the keogram presented by Barros et al. [42], which observed equatorial plasma bubbles through a direct interpolation of VTEC values instead of using tomography.The same days and configurations of South America GNSS stations are used here and we are presenting the results for the same 25°S latitudinal section.Considering the zonal velocity as the ratio between Δlongitude/Δtime, we can see three plasma bubble signatures (PB1, PB2, and PB3) in the keogram, with a general motion progressing from west to east throughout the night.Despite the worse resolution of the tomography in comparison to the keogram presented by Barros et al. [42], the interbubble distances around 600 VTEC depletions are not only observed with a similar location and dimensions to the airglow observations, but they also follow the plasma bubble orientation.For instance, Figure 4 shows that the azimuth of the VTEC depletion mapped in CHPI at DOY 3/2014 is changing accordingly with the azimuth of the airglow dark intensities, which also happened for other instruments and days.In all cases with plasma bubble detections, a longitudinal discrepancy can also be observed between the VTEC depletion and the darker brand of the airglow emissions.A similar discrepancy was observed in [31], which might indicate that the bubble is tilted westward with the increasing height [3,41].
Remote Sens. 2018, 10, x FOR PEER REVIEW 7 of 16 bubble of DOY 3/2014.Thus, even with problems of coverage in the tomographic reconstruction at SJCI, it is possible to detect plasma bubbles in case the dimensions of the ionospheric structure are large enough to be captured.VTEC depletions are not only observed with a similar location and dimensions to the airglow observations, but they also follow the plasma bubble orientation.For instance, Figure 4 shows that the azimuth of the VTEC depletion mapped in CHPI at DOY 3/2014 is changing accordingly with the azimuth of the airglow dark intensities, which also happened for other instruments and days.In all cases with plasma bubble detections, a longitudinal discrepancy can also be observed between the VTEC depletion and the darker brand of the airglow emissions.A similar discrepancy was observed in [31], which might indicate that the bubble is tilted westward with the increasing height [3,41].Additionally, Figure 5 shows a collection of west-east slices of the reconstructed VTEC maps as a form of a keogram.This is presented in order to compare the obtained reconstructions with the keogram presented by Barros et al. [42], which observed equatorial plasma bubbles through a direct interpolation of VTEC values instead of using tomography.The same days and configurations of South America GNSS stations are used here and we are presenting the results for the same 25°S latitudinal section.Considering the zonal velocity as the ratio between Δlongitude/Δtime, we can see three plasma bubble signatures (PB1, PB2, and PB3) in the keogram, with a general motion progressing from west to east throughout the night.Despite the worse resolution of the tomography in comparison to the keogram presented by Barros et al. [42], the interbubble distances around 600 Additionally, Figure 5 shows a collection of west-east slices of the reconstructed VTEC maps as a form of a keogram.This is presented in order to compare the obtained reconstructions with the keogram presented by Barros et al. [42], which observed equatorial plasma bubbles through a direct interpolation of VTEC values instead of using tomography.The same days and configurations of South America GNSS stations are used here and we are presenting the results for the same 25 • S latitudinal section.Considering the zonal velocity as the ratio between ∆longitude/∆time, we can see three plasma bubble signatures (PB1, PB2, and PB3) in the keogram, with a general motion progressing from west to east throughout the night.Despite the worse resolution of the tomography in comparison to the keogram presented by Barros et al. [42], the interbubble distances around 600 km and the inclination of the plasma bubble signatures producing a zonal velocity of 160 m/s (PB1), 144 m/s (PB2), and 128 m/s (PB3) are consistent with their results, which confirms the feasibility of the tomographic reconstruction when compared with a robust single-layer interpolation procedure.

Systematic Assessment
In order to give an overview of the efficiency of the tomographic reconstruction for mapping the plasma bubbles, the experimental comparisons of Figures 2-4 were carried out in a systematic analysis over several other events.Twenty-three days were selected during the summer Sun of the solar maximum period of 2013 and 2014, when we found clear evidence of plasma bubbles propagating through the Brazilian region with the optical instruments.The days chosen were DOYs 358-363 of 2013, as well as DOYs 3-6, 8, 9, 11, and 21-30 of 2014.It was possible to analyze the temporal evolution of the VTEC depletions aligned with the signatures of plasma bubbles in the airglow images for many events, and the general rate of success of the tomographic reconstruction was organized in the form of a confusion matrix.The main goal of the confusion matrix was to show the number of detected plasma bubbles with the VTEC and airglow images, as well as the number of non-detected plasma bubbles by each instrument.Since it is sometimes difficult to confirm whether the plasma bubbles were detected, a third piece of information was included in the confusion matrix, named as uncertain bubbles.
The classification of each feature of the confusion matrix was made according to the authors' interpretation, and Figures 6-8  A detailed inspection of the images was conducted by visual analysis of each VTEC and airglow images.Taking Case (a) as an example, we can clearly see the presence of a plasma depletion at VTEC level of Figure 6, and a simultaneous dark brand in Figure 7.When looking at the entire South America scenario of Figure 8, we then confirm that the correspondent VTEC depletion refers to a large-scale irregularity aligned with the geomagnetic field line.Indeed, Figure 8 is also used in the analysis in order to support the decision, mainly because, sometimes, small-scale irregularities, together with few observations in the reconstructed region, can produce substantial VTEC depletions in the observatories and lead to an incorrect classification.For instance, Case (b) shows a small VTEC

Systematic Assessment
In order to give an overview of the efficiency of the tomographic reconstruction for mapping the plasma bubbles, the experimental comparisons of Figures 2-4 were carried out in a systematic analysis over several other events.Twenty-three days were selected during the summer Sun of the solar maximum period of 2013 and 2014, when we found clear evidence of plasma bubbles propagating through the Brazilian region with the optical instruments.The days chosen were DOYs 358-363 of 2013, as well as DOYs 3-6, 8, 9, 11, and 21-30 of 2014.It was possible to analyze the temporal evolution of the VTEC depletions aligned with the signatures of plasma bubbles in the airglow images for many events, and the general rate of success of the tomographic reconstruction was organized in the form of a confusion matrix.The main goal of the confusion matrix was to show the number of detected plasma bubbles with the VTEC and airglow images, as well as the number of non-detected plasma bubbles by each instrument.Since it is sometimes difficult to confirm whether the plasma bubbles were detected, a third piece of information was included in the confusion matrix, named as uncertain bubbles.
The classification of each feature of the confusion matrix was made according to the authors' interpretation, and Figures 6-8  A detailed inspection of the images was conducted by visual analysis of each VTEC and airglow images.Taking Case (a) as an example, we can clearly see the presence of a plasma depletion at VTEC level of Figure 6, and a simultaneous dark brand in Figure 7.When looking at the entire South America scenario of Figure 8, we then confirm that the correspondent VTEC depletion refers to a large-scale irregularity aligned with the geomagnetic field line.Indeed, Figure 8 is also used in the analysis in order to support the decision, mainly because, sometimes, small-scale irregularities, together with few observations in the reconstructed region, can produce substantial VTEC depletions in the observatories and lead to an incorrect classification.For instance, Case (b) shows a small VTEC depletion in Figure 6, which could be misinterpreted as the southern region of a plasma bubble.Therefore, when looking at the entire South America region, we confirm the plasma bubbles' absence.It is relevant to say that the color bar limits were changed accordingly for the comparison of each VTEC depletion.The same analysis of the four cases (a, b, c, and d) was conducted by looking image by image with a temporal resolution of 6 minutes for each time airglow observations were collected.In the end, a total of 739, 709, and 610 images were analyzed in CHPI, SJCI, and BOAV, respectively, to cover the twenty-three days, and the results are summarized in three confusion matrices (Tables 1-3).
Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 16 depletion in Figure 6, which could be misinterpreted as the southern region of a plasma bubble.Therefore, when looking at the entire South America region, we confirm the plasma bubbles' absence.
It is relevant to say that the color bar limits were changed accordingly for the comparison of each VTEC depletion.The same analysis of the four cases (a, b, c, and d) was conducted by looking image by image with a temporal resolution of 6 minutes for each time airglow observations were collected.
In the end, a total of 739, 709, and 610 images were analyzed in CHPI, SJCI, and BOAV, respectively, to cover the twenty-three days, and the results are summarized in three confusion matrices (Tables 1-3).Remote Sens. 2018, 10, x FOR PEER REVIEW 9 of 16 depletion in Figure 6, which could be misinterpreted as the southern region of a plasma bubble.Therefore, when looking at the entire South America region, we confirm the plasma bubbles' absence.
It is relevant to say that the color bar limits were changed accordingly for the comparison of each VTEC depletion.The same analysis of the four cases (a, b, c, and d) was conducted by looking image by image with a temporal resolution of 6 minutes for each time airglow observations were collected.
In the end, a total of 739, 709, and 610 images were analyzed in CHPI, SJCI, and BOAV, respectively, to cover the twenty-three days, and the results are summarized in three confusion matrices (Tables 1-3).Table 1 presents the obtained numbers in Boa Vista, where the 1st line and 1st column of the matrix (or Table) refers to the number of images where the plasma bubbles were mutually detected by VTEC and airglow observations.The number 181, for instance, reveals that 181 plasma bubbles were mutually seen by the instruments, which represents 29% of images used in the analysis in Boa Vista.The last position of the matrix (3rd line and 3rd column) shows the number of images without plasma bubbles in any method.The matrix also shows the relevant information of the 1st line and 3rd column, which is that 117 airglow images detected plasma bubbles not seen on the VTEC maps.This is one of the most important features, since it directly shows the errors obtained with the tomographic method in comparison to the airglow images.On the other hand, the 3rd line and 1st column shows that 21 plasma bubbles were detected by the tomographic method, but not captured by the optical instruments.This is also an important point, because it says that the tomographic method is more accurate in some instances, which is likely true due to the meteorological limitations of the airglow emissions.Table 1 presents the obtained numbers in Boa Vista, where the 1st line and 1st column of the matrix (or Table) refers to the number of images where the plasma bubbles were mutually detected by VTEC and airglow observations.The number 181, for instance, reveals that 181 plasma bubbles were mutually seen by the instruments, which represents 29% of images used in the analysis in Boa Vista.The last position of the matrix (3rd line and 3rd column) shows the number of images without plasma bubbles in any method.The matrix also shows the relevant information of the 1st line and 3rd column, which is that 117 airglow images detected plasma bubbles not seen on the VTEC maps.This is one of the most important features, since it directly shows the errors obtained with the tomographic method in comparison to the airglow images.On the other hand, the 3rd line and 1st column shows that 21 plasma bubbles were detected by the tomographic method, but not captured by the optical instruments.This is also an important point, because it says that the tomographic method is more accurate in some instances, which is likely true due to the meteorological limitations of the airglow emissions.An overall analysis of the results given by Tables 1-3 shows that the highest amount of plasma bubbles detected by optical instruments was obtained in CHPI, since there are more airglow data for this station.The rate of success of the tomographic method was defined by taking the ratio of the mutually detected plasma bubbles (1st line and 1st column) from the sum of the amount of plasma bubbles detected only by the optical images (all columns of the 1st line).In this regard, the rate of success of the tomographic method to detect plasma bubbles was 54.2%, 31.5%, and 25.7% for BOAV, SJCI, and CHPI, respectively.Actually, the lowest rate of success in CHPI suggests that the tomographic method is better at lower latitudes, which is expected because the dimensions of the plasma bubbles are larger in such regions.In addition, CHPI is located in the region with the highest ionospheric variability for the analyzed season.On the other hand, the CHPI region presented 51 images with plasma bubbles in the VTEC maps that were not detected by the airglow images.This is mainly due to the meteorological conditions at CHPI, where clouds of precipitation disrupted the plasma bubbles imaging, or even due to moonlight brightness contamination.
The ionospheric irregularity activity in the Brazilian sector is well known to be more intense during the December solstice months (summer), when, for example, there is a close alignment between the sunset terminator and the magnetic meridian in the longitude of CHPI.This enhancement in irregularity activity is a consequence of the strong evening equatorial vertical plasma drifts, which transport plasma to higher latitudes and increase its density around the crests of the EIA.Consequently, during December 2013 and January 2014, the electron density ambient and the plasma density gradients are more intense to produce ionospheric irregularities, with much higher variabilities in the latitude of CHPI, as confirmed by the airglow detections shown in Tables 1-3.However, the high overall level of ionospheric irregularities is responsible for the generation of many ionospheric scintillations in the GNSS signals crossing the CHPI region, which produces GNSS loss of lock and reinitializations of the TEC calculation.As a result, few observations of TEC are used in the tomographic reconstruction of the CHPI region, due to the high instability of the ionospheric irregularities.Therefore, it needs to be considered that plasma bubble detection may become problematic at instances with high TEC variability, which seems contradictory to the main goal of using the GNSS-based tomographic algorithm for mapping plasma bubbles.

Vertical Distribution of Plasma Bubbles
This Section shows vertical distributions of the reconstructed ionosphere when applying the proposed tomographic algorithm in instances with evident plasma bubbles.To demonstrate the overall scenario of the reconstructed ionosphere, ionospheric profiles during plasma bubble events were analyzed together with the horizontal distributions.Figure 9 shows VTEC maps obtained at DOY 32 of 2014, and Figure 10 shows the correspondent keogram, where it is possible to see a plasma bubble propagating northeastward with an estimated zonal velocity of 149 m/s.A white line was included in the longitudinal sector of 55 • W in order to analyze the vertical distribution at this specific longitude.As can be seen, the plasma bubble becomes noticeable at 0.6 h UT on the left side of the white line, then it becomes more evident at 1 h UT above the line and overpasses the longitudinal sector at 1.8 h UT.The white line is set at 20.7 LT during 0.4 UT and goes up to 22.1 LT at 1.8 UT.  Figure 11 shows the corresponding vertical images of the longitudinal section given by the white line in order to reproduce the vertical distribution of the scenario presented in Figure 9.The reconstructed tomograms, given in terms of latitude by altitude, are similar to those presented by [43], which analyzed a region of Russia in a bidimensional reconstruction.Now, however, the 3D reconstruction allows us to correlate VTEC images with the tomograms.In general, these Brazilian tomograms present typical representations of the north and south crests of EIA in terms of electron density at 0.4 and 0.6 UT, or 20.7 and 20.9 LT.However, when VTEC depletion reaches the 55°W longitudinal sector at 0.8 hours UT (21.1 LT), the electron density of the south crest of EIA reduces substantially, almost disappearing.In such instances, when the radio-frequency signals overpass the diffraction points of the ionospheric irregularities, destructive and constructive interferences are supposed to cause ionospheric scintillations and spread-F signatures on ionograms.However, after 1.6 h UT (21.9 LT), the south crest of EIA reappears and many ionospheric irregularities dissipate.It is worth mentioning that this event was better characterized in the south crest than the north crest, due to the higher amount of GNSS data in the south region.
It is well known that vertical representations of ionospheric tomography are highly dependent on the background [44].Therefore, the intrinsic limitations of ionospheric tomography for the height estimation suggest that the ionospheric irregularities presented in Figure 11 could be at higher or lower altitudes.Additionally, it is typically observed that the electric field in the region of the structure creates bifurcations in the plasma bubbles as it extends towards to higher latitudes.Figure 11 shows the corresponding vertical images of the longitudinal section given by the white line in order to reproduce the vertical distribution of the scenario presented in Figure 9.The reconstructed tomograms, given in terms of latitude by altitude, are similar to those presented by [43], which analyzed a region of Russia in a bidimensional reconstruction.Now, however, the 3D reconstruction allows us to correlate VTEC images with the tomograms.In general, these Brazilian tomograms present typical representations of the north and south crests of EIA in terms of electron density at 0.4 and 0.6 UT, or 20.7 and 20.9 LT.However, when VTEC depletion reaches the 55°W longitudinal sector at 0.8 hours UT (21.1 LT), the electron density of the south crest of EIA reduces substantially, almost disappearing.In such instances, when the radio-frequency signals overpass the diffraction points of the ionospheric irregularities, destructive and constructive interferences are supposed to cause ionospheric scintillations and spread-F signatures on ionograms.However, after 1.6 h UT (21.9 LT), the south crest of EIA reappears and many ionospheric irregularities dissipate.It is worth mentioning that this event was better characterized in the south crest than the north crest, due to the higher amount of GNSS data in the south region.
It is well known that vertical representations of ionospheric tomography are highly dependent on the background [44].Therefore, the intrinsic limitations of ionospheric tomography for the height estimation suggest that the ionospheric irregularities presented in Figure 11 could be at higher or lower altitudes.Additionally, it is typically observed that the electric field in the region of the structure creates bifurcations in the plasma bubbles as it extends towards to higher latitudes.Figure 11 shows the corresponding vertical images of the longitudinal section given by the white line in order to reproduce the vertical distribution of the scenario presented in Figure 9.The reconstructed tomograms, given in terms of latitude by altitude, are similar to those presented by [43], which analyzed a region of Russia in a bidimensional reconstruction.Now, however, the 3D reconstruction allows us to correlate VTEC images with the tomograms.In general, these Brazilian tomograms present typical representations of the north and south crests of EIA in terms of electron density at 0.4 and 0.6 UT, or 20.7 and 20.9 LT.However, when VTEC depletion reaches the 55 • W longitudinal sector at 0.8 hours UT (21.1 LT), the electron density of the south crest of EIA reduces substantially, almost disappearing.In such instances, when the radio-frequency signals overpass the diffraction points of the ionospheric irregularities, destructive and constructive interferences are supposed to cause ionospheric scintillations and spread-F signatures on ionograms.However, after 1.6 h UT (21.9 LT), the south crest of EIA reappears and many ionospheric irregularities dissipate.It is worth mentioning that this event was better characterized in the south crest than the north crest, due to the higher amount of GNSS data in the south region.
It is well known that vertical representations of ionospheric tomography are highly dependent on the background [44].Therefore, the intrinsic limitations of ionospheric tomography for the height estimation suggest that the ionospheric irregularities presented in Figure 11 could be at higher or lower altitudes.Additionally, it is typically observed that the electric field in the region of the structure creates bifurcations in the plasma bubbles as it extends towards to higher latitudes.However, this effect is not reproduced due to the weak geometry of the GNSS signals for tomographic applications.However, despite the problems in correctly mapping the vertical altitudes of the irregularities and bifurcations, the 3D tomographic reconstruction shows a strong potential to analyze the general vertical distribution of ionospheric irregularities, which indicates future possibilities of building ionospheric models to describe the three-dimensional impact of ionospheric irregularities on the transionospheric signals.
Remote Sens. 2018, 10, x FOR PEER REVIEW 13 of 16 However, this effect is not reproduced due to the weak geometry of the GNSS signals for tomographic applications.However, despite the problems in correctly mapping the vertical altitudes of the irregularities and bifurcations, the 3D tomographic reconstruction shows a strong potential to analyze the general vertical distribution of ionospheric irregularities, which indicates future possibilities of building ionospheric models to describe the three-dimensional impact of ionospheric irregularities on the transionospheric signals.

Conclusions
A new tomographic reconstruction technique was successfully implemented for imaging the spatial distribution of ionospheric plasma bubbles in the Brazilian region.The mapped horizontal distributions of VTEC depletions were simultaneously compared to the darker brands observed by all-sky airglow imagers in the OI 630 nm emissions.A systematic analysis using 23 days in the peak of the solar cycle 24 was carried out and revealed the rate of success of the tomographic method of 54.2%, 31.5%, and 25.7% for observatories in Boa Vista, São João do Cariri, and Cachoeira Paulista, respectively.The results suggest that the tomographic method is more efficient in the equatorial region, which is expected due to the larger dimensions of the plasma bubbles near the magnetic equator and the higher ionospheric variability in the crests of EIA.However, despite the problems of imaging plasma bubbles with tomographic reconstructions, the latitudes of the VTEC maps covered more than 30° of the plasma bubbles structures and have detected the depletion structures at instances where the meteorological conditions disrupted the possibility of analyzing the OI 630 nm emissions with optical instruments.Additionally, the tomographic reconstructions have shown how the vertical structure of the plasma bubbles can be estimated, which can be useful for further description of the impact of 3D ionospheric irregularities on the transionospheric signals.In future, two analyses are suggested: (1) To validate the plasma bubble vertical reconstruction using the ionosonde data as reference; and (2) to use an accurate method to perform the GNSS-TEC estimation using only single-frequency L1 data.This way, it is expected that fewer TEC gaps will be obtained for plasma bubble imaging, mainly during scintillation events, since the tomographic reconstruction will likely be less sensitive to loss of locks than using dual-frequency TEC estimations.

Conclusions
A new tomographic reconstruction technique was successfully implemented for imaging the spatial distribution of ionospheric plasma bubbles in the Brazilian region.The mapped horizontal distributions of VTEC depletions were simultaneously compared to the darker brands observed by all-sky airglow imagers in the OI 630 nm emissions.A systematic analysis using 23 days in the peak of the solar cycle 24 was carried out and revealed the rate of success of the tomographic method of 54.2%, 31.5%, and 25.7% for observatories in Boa Vista, São João do Cariri, and Cachoeira Paulista, respectively.The results suggest that the tomographic method is more efficient in the equatorial region, which is expected due to the larger dimensions of the plasma bubbles near the magnetic equator and the higher ionospheric variability in the crests of EIA.However, despite the problems of imaging plasma bubbles with tomographic reconstructions, the latitudes of the VTEC maps covered more than 30 • of the plasma bubbles structures and have detected the depletion structures at instances where the meteorological conditions disrupted the possibility of analyzing the OI 630 nm emissions with optical instruments.Additionally, the tomographic reconstructions have shown how the vertical structure of the plasma bubbles can be estimated, which can be useful for further description of the impact of 3D ionospheric irregularities on the transionospheric signals.In future, two analyses are suggested: (1) To validate the plasma bubble vertical reconstruction using the ionosonde data as reference; and (2) to use an accurate method to perform the GNSS-TEC estimation using only single-frequency L1 data.This way, it is expected that fewer TEC gaps will be obtained for plasma bubble imaging, mainly during scintillation events, since the tomographic reconstruction will likely be less sensitive to loss of locks than using dual-frequency TEC estimations.

Figure 1 .
Figure 1.Global navigation satellite system (GNSS) stations operating throughout South America in Day of Year (DOY) 3 of 2014.The colored marks indicate the location of the stations for the different networks, where the dashed line shows the magnetic equator.

Figure 1 .
Figure 1.Global navigation satellite system (GNSS) stations operating throughout South America in Day of Year (DOY) 3 of 2014.The colored marks indicate the location of the stations for the different networks, where the dashed line shows the magnetic equator.

Figure 2 .
Figure 2. Comparative example between the Vertical Total Electron Content (VTEC) depletions and the dark brands of the airglow observations.The white square shows the exact location of the Georeferenced OI 630 nm images and the grayscale is in units of normalized intensity (NI).

Figure 2 .
Figure 2. Comparative example between the Vertical Total Electron Content (VTEC) depletions and the dark brands of the airglow observations.The white square shows the exact location of the Geo-referenced OI 630 nm images and the grayscale is in units of normalized intensity (NI).

Figure 3 .
Figure 3. Example of a plasma bubble propagating through São João do Cariri (SJCI) in DOY 359 of 2013, detected by the VTEC reconstructions and the optical instruments.

Figure 4 .
Figure 4. Example of a plasma bubble propagating through Cachoeira Paulista (CHPI) in DOY 3 of 2014 detected by the VTEC reconstructions and the optical instruments.

Figure 3 .
Figure 3. Example of a plasma bubble propagating through São João do Cariri (SJCI) in DOY 359 of 2013, detected by the VTEC reconstructions and the optical instruments.

Figure 3 .
Figure 3. Example of a plasma bubble propagating through São João do Cariri (SJCI) in DOY 359 of 2013, detected by the VTEC reconstructions and the optical instruments.

Figure 4 .
Figure 4. Example of a plasma bubble propagating through Cachoeira Paulista (CHPI) in DOY 3 of 2014 detected by the VTEC reconstructions and the optical instruments.

Figure 4 .
Figure 4. Example of a plasma bubble propagating through Cachoeira Paulista (CHPI) in DOY 3 of 2014 detected by the VTEC reconstructions and the optical instruments.
Remote Sens. 2018, 10, x FOR PEER REVIEW 8 of 16km and the inclination of the plasma bubble signatures producing a zonal velocity of 160 m/s (PB1), 144 m/s (PB2), and 128 m/s (PB3) are consistent with their results, which confirms the feasibility of the tomographic reconstruction when compared with a robust single-layer interpolation procedure.

Figure 5 .
Figure 5. Keogram of the VTEC values at the latitudinal section of 25°S in terms of longitude and hours at the night of DOYs 358 and 359 from 2013.
show a few examples of the classification.They show the main classes of the confusion matrix: (a) Images with plasma bubbles mutually detected by the tomographic method and the optical instruments; (b) plasma bubbles detected by the airglow images, but not detected in the VTEC reconstructions; (c) images with detected plasma bubbles at VTEC level, but not detected in the airglow images; and (d) plasma bubbles not detected by any instrument.All of these examples are presented for the observatory of Cachoeira Paulista, in which Case (a) refers to DOY 22 of 2014 at 1.5 UT (22.5 LT at the central longitude), Case (b) to DOY 24 of 2014 at 0.8 UT (21.8 LT at the central longitude), Case (c) to DOY 358 of 2013 at 0.5 UT (21.5 LT at the central longitude), and Case (d) to DOY 360 of 2013 at 0 UT (21 LT at the central longitude).

Figure 5 .
Figure 5. Keogram of the VTEC values at the latitudinal section of 25 • S in terms of longitude and hours at the night of DOYs 358 and 359 from 2013.
show a few examples of the classification.They show the main classes of the confusion matrix: (a) Images with plasma bubbles mutually detected by the tomographic method and the optical instruments; (b) plasma bubbles detected by the airglow images, but not detected in the VTEC reconstructions; (c) images with detected plasma bubbles at VTEC level, but not detected in the airglow images; and (d) plasma bubbles not detected by any instrument.All of these examples are presented for the observatory of Cachoeira Paulista, in which Case (a) refers to DOY 22 of 2014 at 1.5 UT (22.5 LT at the central longitude), Case (b) to DOY 24 of 2014 at 0.8 UT (21.8 LT at the central longitude), Case (c) to DOY 358 of 2013 at 0.5 UT (21.5 LT at the central longitude), and Case (d) to DOY 360 of 2013 at 0 UT (21 LT at the central longitude).

Figure 6 .
Figure 6.Four cases to build the confusion matrix using the VTEC reconstructions in CHPI.

Figure 7 .
Figure 7. Four cases to build the confusion matrix using the airglow observations in CHPI.

Figure 6 .
Figure 6.Four cases to build the confusion matrix using the VTEC reconstructions in CHPI.

Figure 6 .
Figure 6.Four cases to build the confusion matrix using the VTEC reconstructions in CHPI.

Figure 7 .
Figure 7. Four cases to build the confusion matrix using the airglow observations in CHPI.Figure 7. Four cases to build the confusion matrix using the airglow observations in CHPI.

Figure 7 .
Figure 7. Four cases to build the confusion matrix using the airglow observations in CHPI.Figure 7. Four cases to build the confusion matrix using the airglow observations in CHPI.

Figure 8 .
Figure 8. Four cases to build the confusion matrix using the VTEC reconstructions over all the South America region.The white square shows the coverage of the CHPI images.

Figure 8 .
Figure 8. Four cases to build the confusion matrix using the VTEC reconstructions over all the South America region.The white square shows the coverage of the CHPI images.

Figure 9 .
Figure 9. Reconstructed VTEC maps during the propagation of an ionospheric plasma bubble over the South America region.The white line is designed at the longitudinal sector of 55°W.

Figure 10 .
Figure 10.Keogram of VTEC values in terms of longitude and hour at the latitudinal section of 24°S during DOY 32 from 2014.

Figure 9 . 16 Figure 9 .
Figure 9. Reconstructed VTEC maps during the propagation of an ionospheric plasma bubble over the South America region.The white line is designed at the longitudinal sector of 55 • W.

Figure 10 .
Figure 10.Keogram of VTEC values in terms of longitude and hour at the latitudinal section of 24°S during DOY 32 from 2014.

Figure 10 .
Figure 10.Keogram of VTEC values in terms of longitude and hour at the latitudinal section of 24 • S during DOY 32 from 2014.

Figure 11 .
Figure 11.Reconstructed ionospheric profiles during the propagation of an ionospheric plasma bubble over the longitudinal sector of 55°W.

Author Contributions:
This research was performed by F.S.P.; supervised by P.d.O.C.; M.T.d.A.H.M.; and M.H.-P.These authors contributed extensively and equally to prepare this paper.Funding: This work was jointly funded by Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP Grants: 2015/15027-7 and

Figure 11 .
Figure 11.Reconstructed ionospheric profiles during the propagation of an ionospheric plasma bubble over the longitudinal sector of 55 • W.

Table 1 .
Confusion matrix showing the number and percentage of plasma bubbles detected by the tomographic reconstruction and optical instruments in the observatory of Boa Vista.

Table 2 .
Confusion matrix showing the number and percentage of plasma bubbles detected by the tomographic reconstruction and optical instruments in the observatory of São João do Cariri.

Table 1 .
Confusion matrix showing the number and percentage of plasma bubbles detected by the tomographic reconstruction and optical instruments in the observatory of Boa Vista.

Table 2 .
Confusion matrix showing the number and percentage of plasma bubbles detected by the tomographic reconstruction and optical instruments in the observatory of São João do Cariri.

Table 3 .
Confusion matrix showing the number and percentage of plasma bubbles detected by the tomographic reconstruction and optical instruments in the observatory of Cachoeira Paulista.