Ionospheric irregularities as plasma bubbles occur in the ionosphere, and generally, they are characterized by the low plasma density regions. The Rate Of TEC Index (ROTI) was defined in terms of the Total Electron Content (TEC) variation and is used to characterize these plasma bubbles. It is essential to evaluate the ROTI behavior since the ionospheric irregularities can cause fluctuations in the radio signal, interfering in the ionospheric data analysis. Therefore, we performed in this work a comparative study of five different methods available to calculate the ROTI to evaluate the most suitable over the Brazilian region. The ROTI was calculated over three GNSS stations at different latitudes: São Luís (SALU, 2°31′S, 44°16′W; dip: −3.8°), Cachoeira Paulista (CHPI, 22°40'S, 44°59'W; dip: −36.4°), and Santa Maria (SMAR, 29°41′S, 53°48′W, dip: −37°). The results show that the most viable method for calculating ROTI in the Brazilian region is based on the Slant TEC equation as defined by Cherniak et al. (2018), https://doi.org/10.1007/s10291‐018‐0730‐1. Our results are supported by the comparison between the ROTI with TEC maps, ionograms collected at Fortaleza (FZA0M, 3°43'S, 38°32'W, dip: −15.8°), SALU, and CHPI and All‐Sky imagers collected at the São João do Cariri (SJCA, 7°23'S, 36°31'W, dip: −11°) and CHPI. Additionally, we observe Equatorial Plasma Bubbles (EPBs) in the Brazilian region using the ROTI index map.

During disturbed periods, E region electric fields can cause anomalous Es layer behavior, which is observed in the digital ionosonde data. To investigate the influence of these electric fields in the Es layer development, we analyzed a set of 20 magnetic storms from 2015 to 2018 over Boa Vista (BV, 2.8°N, 60.7°W, dip ∼18°), São Luís (SLZ, 2.3°S, 44.2°W, dip ∼8°), and Cachoeira Paulista (CXP, 22.41°S, 45°W, dip ∼35°). The electric field zonal components during the main and recovery phases of each magnetic storm are computed to study the corresponding characteristics of these Es seen in ionograms. Additionally, a numerical model (MIRE, Portuguese acronym for E Region Ionospheric Model) is used to analyze the Es layer dynamics modification around disturbed times. Using observation data and simulations, we were able to establish a threshold value for the electric field intensity for each region that can affect the Es layer formation. The results sustain that the strong Es layer in BV can be an indicator of the disturbed dynamo event. At SLZ, on the other hand, the Es layers are affected by the competition mechanisms of their formation, as equatorial electrojet irregularities and winds, during the main phase of the magnetic storm. Over CXP, the Es layer dynamics are dominated by the wind shear mechanism. Finally, this study provides new insights into the real impact of the electric field in the Es layer development over the Brazilian sector. Thus, our results lead to a better understanding of the underlying mechanisms related to the Es layer formation and dynamics.