Study and Modelling of the Impact of June 2015 Geomagnetic Storms on the Brazilian Ionosphere

: This study investigated the impact of the June 2015 geomagnetic storms on the Brazilian equatorial and low-latitude ionosphere by analyzing various data sources, including solar wind parameters from the advanced compositional explorer satellite (ACE), global positioning satellite vertical total electron content ( GPS-VTEC ), geomagnetic data, and validation of the SAMI2 model-VTEC with GPS-VTEC . The effect of geomagnetic disturbances on the Brazilian longitudinal sector was examined by applying multiresolution analysis (MRA) of the maximum overlap discrete wavelet transform (MODWT) to isolate the diurnal component of the disturbance dynamo ( Ddyn ), DP2 current fluctuations from the ionospheric electric current disturbance ( Diono ), and semblance cross-correlation wavelet analysis for local phase comparison between the Sq and Diono currents. Our findings revealed that the significant fluctuations in DP2 at the Brazilian equatorial stations (Belem, dip lat: − 0.47 ◦ and Alta Floresta, dip lat: − 3.75 ◦ ) were influenced by IMF Bz oscillations; the equatorial electrojet also fluctuated in tandem with the DP2 currents, and dayside reconnection generated the field-aligned current that drove the DP2 current system. The short-lived positive ionospheric storm during the main phase on 22 June in the Southern Hemisphere in the Brazilian sector was caused by the interplay between the eastward prompt penetration of the magnetospheric convection electric field and the westward disturbance dynamo electric field. The negative iono-spheric storms that occurred during the recovery phase from 23 to 29 June 2015, were attributed to the westward disturbance dynamo electric field, which caused the downward E × B drift of the plasma to a lower height with a high recombination rate. The comparison between the SAMI2 model-VTEC and GPS-VTEC indicates that the SAMI2 model underestimated the VTEC within magnetic latitudes of − 9 ◦ to − 24 ◦ in the Brazilian longitudinal sector from 6 to 17 June 2015. However, it demonstrated satisfactory agreement with the GPS-VTEC within magnetic latitudes of − 9 ◦ to 10 ◦ from 8 to 15 June 2015. Conversely, the SAMI2 model overestimated the VTEC between ± 10 ◦ magnetic latitudes from 16 to 28 June 2015. The most substantial root mean square error (RMSE) values, notably 10.30 and 5.48 TECU , were recorded on 22 and 23 June 2015, coinciding with periods of intense geomagnetic disturbance.

During magnetic storms, two primary physical processes operating on a planetary scale can be observed: (1) the direct penetration of polar cap electric fields to the Equator, giving rise to the DP2 current system [14][15][16], and (2) the disruption of winds caused by auroral joule heating and ion-drag acceleration, leading to the attenuation of the equatorial electrojet in the aftermath of the storm.Mayaud (1982) [17] initially formulated the theory of the direct penetration of polar cap electric fields to equatorial latitudes and predicted the shielding of the electric fields within 30 min [18,19].Over the last three decades, numerous theoretical and experimental studies have focused on the direct penetration of polar cap electric fields into equatorial latitudes [20][21][22][23] and examined the thermospheric response to magnetic storms; subsequently, ref. [24] introduced the ionospheric disturbance dynamo to elucidate the electric field disturbance observed with the incoherent scatter sounder of Saint-Santin in the aftermath of the storm.This disturbance is attributed to the dynamo action of storm winds generated by auroral Joule heating [24][25][26].
The magnetospheric generator builds up excess positive charges on the dawn side of the magnetosphere and negative charges on the dusk side [27].Charges with access to highly conducting magnetic field lines seek to relieve charge build-up by diverging from the source region.Current continuity requires that a divergence in the horizontal current be accompanied by a change in the field-aligned currents [27,28].Magnetic merging immediately leads to field-aligned currents between the magnetosphere and ionosphere [29].Region-1 (R1) currents are largely driven by magnetopause charge separation.These currents exist regardless of the direction of the interplanetary magnetic field [29,30].When IMF Bz is southward, the R1 currents tend to be stronger, consistent with an enhanced merging electric field in the magnetotail.Another set of field-aligned currents, Region-2 (R2) currents, have their sources in the inner magnetosphere.Although their flow patterns are similar to those of the R1 currents, their polarity is opposite, and they intercept the ionosphere at a lower latitude.The R2 current system arises from a charge distribution that is a consequence of the drift of charged particles in Earth's dipole field [1,28,[30][31][32][33].
Under equilibrium conditions with little or no convection, the magnetospheric particle gradient and curvature drift follow the contours of the magnetospheric electric potential Φ [27,28,30].When a cross-tail electric field E from convection enhancement develops, the plasmasheet surges sunward, and a partial westward ring current develops.The dusk edge of the plasmasheet is positively charged and the dawn edge is negatively charged.The excess charge creates a dusk-to-dawn polarization electric field.The superposition of the convection electric field and the polarization electric field tends to shield the near-Earth region from the enhanced dawn-to-dusk convection electric field.The shielding configuration usually prevents widespread surge of the electric field (known as prompt penetration of magnetospheric convection electric field PPMEF) into the inner magnetosphere and lowand mid-latitude ionosphere [14,16,20,21,31,34].
However, the shielding layer and the field did not develop instantaneously.Thus, for short intervals, the inner magnetosphere and low-and mid-latitude ionosphere may be exposed to an enhanced convective electric field.This condition is known as undershielding [1,[29][30][31][32][33]. Undershielding is the temporary penetration of the dawn-to-dusk electric field into the inner magnetosphere during the time of increasing convection [1].During undershielding, the dusk-to-dawn polarization field (R2 region current) is insufficient to cancel the enhanced convective electric field [31].A stronger convection field penetrates the field lines into the ionosphere (R1 region).On the other hand, when the driving convection field suddenly decreases as a result of a northward turning of the IMF Bz, the dusk-to-dawn polarization electric field dominates across the inner magnetosphere, which is referred to as overshielding.Overshielding is the temporary dominance of the dusk-to-dawn polarization field over the dawn-to-dusk convection field [1,19,31,32,[35][36][37][38][39].
During the magnetically quiet period, the ionospheric regular dynamo gives rise to quiet solar (Sq) and equatorial electrojet (EEJ) current systems at middle and equatorial latitudes, respectively [32,40].In contrast, during magnetically disturbed periods, intensified ionospheric electric currents known as auroral electrojets emerge in the auroral zone.These electrojets dissipate energy via the Joule heating effect, thereby influencing the temperature, pressure, and motion of the thermosphere.The perturbed thermospheric winds propagate toward mid-and low-latitudes, generating disturbed ionospheric electric fields, referred to as the disturbance dynamo electric field (DDEF), through the dynamo effect.Blanc and Richmond (1980) [24] introduced the first numerical simulation of an ionospheric disturbance dynamo.
This study focused on the magnetic manifestations of these physical processes on the equatorial electrojet (EEJ) and global positioning satellite vertical total electron content (GPS-VTEC).Mazaudier (1985) [41] observed disturbance winds, a consequence of Joule heating, using incoherent scatter radar data, and their observations aligned with the predictions of the Blanc and Richmond model (1980) [24].
The geomagnetic storms of June 2015 have gained attention from numerous space scientists, such as [42][43][44][45][46][47].Astafyeva et al. (2017Astafyeva et al. ( , 2018) employed a multi-instrument to study the 22-23 June 2015 geomagnetic storms and observed the effect of eastward PPMEF, which caused short-lived positive ionospheric storms during the day, and the eastward DDEF caused long-lived positive ionospheric storms on the nightside of the main phase and the early stage of the recovery phase.Mansilla et al. (2018) [47] also studied the same storm globally and observed increased total electron content (TEC) at high latitudes before the storm's main phase, significant asymmetry in TEC response between the middle and low latitudes of the Northern and Southern Hemispheres, and TEC decrease at equatorial latitudes.The prolonged TEC enhancements closely correlated with an increase in the O/N 2 ratio; however, TEC decreases are not linked to reductions in the O/N 2 ratio, unlike concurrent decrease in electron density.They concluded that eastward PPMEF significantly influenced the equatorial and low-latitude ionospheres during the main phase of the storm.
Paul (2018) examined the latitudinal ionospheric response during the three most prominent geomagnetic storms of 2015 during the solar cycles 24, 16-17 March 2015, 22-23 June 2015, and 19-22 December 2015 using TEC data derived from a latitudinal chain of GPS receivers that extend from 70 • N to 70 • S. They observed a symmetrical hemispheric response during the strongest 17 March (St.Patrick's Day) storm, and an asymmetric hemispheric response of the ionosphere during the 22 June and 20 December storms over the Asian-Australian sector and African equatorial and low-latitude sectors [46].Amaechi (2018) examined the effects of the intense geomagnetic storms of 2015 on the occurrence of large-scale ionospheric irregularities over Africa and reported that irregularities occurred in the post-sunset to midnight period and are associated with TEC depletions and fluctuations due to equatorial plasma bubbles, which are stronger over the equatorial ionospheric anomaly (EIA) crest [48].
Singh (2017) investigated the impact of the 22-23 June 2015 geomagnetic storm on the Indian ionosphere by using ionosondes.They observed suppressed spread F during the westward penetration of the electric field of the overshielding R2 region in the Indian sector.Fluctuations in foF2 with shorter periods are attributed to eastward prompt penetration electric field fluctuations, whereas larger period fluctuations are mainly caused by disturbance wind, traveling ionospheric disturbances (TIDs), and DDEFs.In contrast to their detection in the European sector, plasma bubbles are suppressed in the Indian sector [29].Macho (2020) studied ionospheric dynamics in the South American sector from 21 to 24 June 2015.They used ground ionosonde stations, GNSS receivers, very low-frequency (VLF), and magnetometer data.This study revealed an expansion of the crest of the equatorial ionospheric anomaly (EIA) at mid and high latitudes, mainly due to the prompt penetration of the electric field during the main phase [45].
In this study, we investigated the development of Ddyn in the Brazilian equatorial ionosphere in relation to auroral electrojet intensification due to the interaction between solar winds and the magnetosphere, the response of the equatorial electrojet to DP2 fluctuations, and the westward Ddyn in the Brazilian sector during the main and recovery phases of geomagnetic storms in June 2015.We examined the impact of the Diono current on dayto-day VTEC redistribution and compared the SAMI2 model-VTEC with the GPS-VTEC for the days investigated in this study.
We introduced a novel approach by separating DP2 fluctuations and Ddyn from the Diono current using two innovative techniques.The noteworthy aspect is the inaugural use of multi-resolution analysis of the maximum overlap discrete wavelet transform (MODWT) and semblance cross-correlation wavelet analysis for a detailed exploration of geomagnetic disturbances in Brazilian sector.This paper presents the first study of ionospheric signatures of space weather events in the Brazilian sector over a period of one month.Previous studies have generally focused on individual events.A similar study was conducted in the Indian sector recently [46].
The first section of this research work introduces and emphasizes the research gap that this article fills in the Brazilian sector.The second section presents the data sources and methodology.In the third section, we present the results of this study, and in the fourth section, we discuss and summarize the final results.

Data Source
High-resolution interplanetary parameters in geocentric solar magnetospheric coordinates (GSM) recorded onboard the advanced composition explorer (ACE) satellite are being used to monitor the solar wind conditions.These are the solar wind speed (Vsw), z-component of the interplanetary magnetic field (IMF Bz), y-component of the interplanetary electric field (IEyF), and auroral electrojet indices (auroral upper boundary (AU) and auroral lower boundary (AL)).The data are made available by OMNIWeb data services via http://omniweb.gsfc.nasa.gov/form/omni_min.html.This website was accessed on 28 August 2023.
Equatorial and low-latitude ionosphere responses to magnetic perturbation were examined by using the global navigation satellite system (GNSS-TEC) obtained by processing GNSS-RINEX files made available by the Brazilian Network for Continuous Monitoring of the GNSS System (RBMC) via http://www.ibge.gov.br(accessed on 28 August 2023).A map of the geographic coordinates of the GNSS receivers' stations are shown in Figure 1.

GPS Data
The GPS Receiver Independent Exchange Format (RINEX) files obtained from the Brazilian Network for Continuous Monitoring of GNSS Systems were processed by using GPS-TEC analysis application software developed by Gopi Seemala and Boston College, United States.Details of the GPS-TEC application software can be found in Seemala (2011) [49].
We computed the TEC-MAP using one-minute VTEC resolution data.The VTEC data were grid 1   in geographical latitude and longitude.The dip latitudes were computed using the British Geological Survey website (https://geomag.bgs.ac.uk/cgi-bin/coord_calc (last accessed on 28 August 2023)).Magnetically quiet days were obtained from the International Service of Geomagnetic Indices (ISGI, http://isgi.unistra.fr(las accessed on 28 August 2023)).We computed the storm-time residual VTEC by subtracting the average VTEC of five magnetically quiet days from that of magnetically disturbed days.
Equations ( 1) and (2) expressed the mathematical computation processes for calculating the change in VTEC (∆VTEC), where N denotes the number of magnetically quiet days, Q 5AVR is the average of five magnetically quiet days, and MDs is the magnetically disturbed days.

Computation of Equatorial Electrojet (EEJ)
We computed the horizontal component of Earth's magnetic field from the X and Y components of Earth's magnetic field, as shown in Equation (3).
To eliminate variations in offset values among different magnetometers, we first calculated the nighttime baseline values for the H component using Equation ( 4) for each day.These baseline values were then subtracted from the corresponding magnetometer dataset to obtain the hourly departure of H, which is denoted as δH and can be expressed using Equation (5).The baseline value is defined as the average of the H component nighttime (23: The time variable t represents the hours ranging from 01:00 to 24:00 LT.The hourly departure, δH, was then adjusted for non-cyclic variation using Equation ( 6).This correction method was previously proposed by Rastogi et al. (2004), who defined noncyclic variation as a phenomenon where the value at 01:00 LT differs from that at local midnight (24:00 LT) [50].
S q (t) = δH(t) + (t − 1).∆c (7) where t = 1 to 1440.The values for the solar quiet variation (Sq) can be obtained by correcting the hourly departure of H(δH) for the noncyclic variation based on the magnetometer dataset.This relationship is expressed by Equation (7).
To estimate the vertical drift velocity from the ground-based magnetometer for daytime plasma drift, as given in Equation ( 9) [34,35], we used multiple regression analysis techniques, as suggested by Anderson (2004).

Computation of Ionospheric Electric Current Disturbance (Diono)
The ionospheric electric current disturbance (Diono) was computed by using Equation (10).Diono is the magnetic disturbance caused by the ionospheric electric current system [53][54][55][56][57][58].H is the Earth's horizontal magnetic field measured component, ∅ is the magnetic latitude of the stations, and S H R is the daily regular fluctuation of H caused by the solar quiet (S q ) system.S H R was calculated by averaging five magnetically quiet days provided by the World Data Center Kyoto for Geomagnetism (Kyoto-u.ac.jp).DP2 is the total magnetic disturbance related to the disturbance polar number 2 (DP2) fluctuations during the prompt penetration of the magnetospheric convection electric field (PPMEF), and Ddyn is the magnetic signature of the ionospheric disturbance dynamo caused by Joule heating in the auroral zone (Blanc and Richmond 1980) [24].Joule heating generates disturbed thermospheric winds, which alter the quiet-time global thermospheric circulation.In this study, we refer to Ddyn as the magnetic signature of the ionospheric disturbance dynamo, first described by Le Huy and Amory- Mazaudier (2005Mazaudier ( , 2008) ) [40,59] as shown in Equation (11), ref. [15,16,55,56,60].The ionospheric electric current disturbance (Diono) is the sum of DP2 and Ddyn.We focused on the impact of DP2 current fluctuations and Ddyn on the Brazilian equatorial and low-latitude ionosphere during geomagnetic storms in June 2015.Two techniques were used to isolate the effects of DP2 and Ddyn on the Brazilian equatorial ionosphere.These techniques are: i.
Multiresolution analysis (MRA) of maximum overlap discrete wavelet transforms (MODWT).We extracted the diurnal component of Ddyn, which had a period of approximately 24 h.For DP2 current fluctuations, we extracted the Diono component, which has a period of less than 3 h [61][62][63].ii.
Semblance cross-correlation wavelet analysis [54,64].Interested readers can refer to the article published by Younas (2021) [54,64,65].We performed a semblance analysis by comparing the local phase between the Sq and the Diono currents as a function of time and wavelength.The details of this new method to isolates anti-Sq from Diono can be found in [54].The positive phase is represented by +1 and the negative phase is represented by −1 (which is anti-Sq current).

Computation of Akasofu's Parameter
When the interaction between the solar wind and magnetosphere intensifies, powerful magnetospheric currents are generated, and the magnetic energy density within the magnetosphere increases.This energy leads to the reconfiguration of the magnetosphere, faster movement of plasma in the plasmasheet, heating and acceleration of particles, and amplification of currents [66,67].The movement of solar wind particles on Earth generates kinetic energy flux.The energy flux in the plasma plays a crucial role in shaping the magnetosphere.Additionally, the interplanetary magnetic field influences the magnetosphere through a flux known as the Poynting flux [1,67].
where P = effective solar wind poynting flux W/m 2 µ 0 = permeability of free space 1.26 The Poynting flux is the power delivered by electromagnetic sources.Akasofu (1981) estimated an electromagnetic power flux that represents the solar wind power delivered by magnetic merging, as shown in Equation (12).He assumed that only the tangential component B t of the IMF delivered energy; therefore, the B x component did not play a role.The elements in the square brackets of Equation ( 12) constitute the Poynting flux.The factor outside the bracket accounts for the efficiency of energy entry by the anti-parallel component of the IMF.
For many storm-time applications, space weather scientists are most interested in the power delivered to the cross-section of the dayside magnetosphere.One version of power estimates the Akasofu epsilon parameter (ε), which is obtained by multiplying the effective solar wind Poynting flux by a factor proportional to the cross-sectional area of the magnetopause that intercepts the solar wind, as shown in Equation (13). where Scientists frequently use ε as a reference point when analyzing the energy distribution in magnetic storms.Nevertheless, in the case of severe storms, this estimation may become inaccurate because of significant fluctuations in the cross-sectional area of the magnetopause or the influence of other factors that are not well understood and can affect the efficiency of the energy transfer [1,66,67].

SAMI2 Model
SAMI2 (another model for ionosphere) calculates the evolution of the ionosphere in the low-to mid-latitudes.SAMI2 considers the dynamic plasma and chemical changes of seven ion species (H + , He + , N + , O + , N + 2 , NO + , and O + 2 ) within an altitude range of 85 km to 20,000 km.This altitude range corresponds to a latitudinal span of ± 62.5 • around the magnetic equator [68].
SAMI2 solves the ion continuity and momentum equations for all seven ion species.In addition, it solves the temperature equation for H + , He + , O + and the electrons.The code models the plasma along Earth's geomagnetic field, spanning from one hemisphere to another.For this purpose, it utilizes an offset-tilted dipole field.
The code incorporates a modeled E × B drift of the plasma, which represents the motion of the plasma owing to the interaction between the electric field (E) and magnetic field (B).It also considers the ion inertia in the ion momentum equation to account for the motion along the dipole field line.
SAMI2 employs a fixed, nonorthogonal grid, where one coordinate axis aligns with the geomagnetic field.The neutral species were characterized using the empirical models NRLMSISE00 and HWM14.These models provide information regarding a neutral atmosphere, which is essential for understanding the ionosphere's behavior [68,69].
The input parameters in the SAMI2 model are F10.7A(3 months average F10.7), the daily value of F10.7, the daily value of Ap, the year, and the day.The SAMI2 code is coupled with the NRLMSISE00, HWM14, and EUVAC models developed by Richards [70].By default, it is assumed that the ionosphere's perpendicular vertical E × B drift, which is influenced by the magnetic field, can be either sinusoidal or defined by the E × B drift model [26].
The primary outputs of the SAMI2 model include ion density, ion temperature, ion velocity along the magnetic field, electron temperature, time step in local time, geographic latitude, geographic longitude, and altitude.The details of the SAMI2 model can be found in an article by Hubal et al. (2000) [68].In Figure 3, the 1-min resolution data for the interplanetary parameters, including IMF Bz and IEyF (refer to Figure 3a-c), along with the solar wind speed (Vsw; refer to Figure 3a), are presented.The aim is to examine the concurrent patterns of auroral activity, achieved by comparing 1-min resolution auroral electrojet indices (refer to Figure 3d).These indices were derived from a network of strategically positioned ground-based magnetometers at latitudes where auroras are commonly observed.The AU index (highlighted in red) signifies the envelope of the maximum northward disturbance recorded at all the stations, indicating an eastward auroral electrojet.Conversely, the AL index (depicted in black) represents the envelope of the minimum disturbance (maximum negative value) typically associated with the westward auroral electrojet.These indices offer a comprehensive overview of the auroral activity and prevailing electrojet conditions [71].As part of our analysis, we incorporated 1-min-resolution geomagnetic data, specifically the SYM-H and ASYM-H indices (Figure 3e).The SYM-H and ASYM-H indices are indicators of the symmetric and asymmetric part of the ring current.The asymmetric or partial ring current occurs when hot ions from the tail region are injected into the evening and afternoon sectors, and subsequently drift and close in the ionosphere through fieldaligned currents at dusk and dawn.By examining the enhancements in the ASYM-H index, we can gain insight into the occurrence of substorm injections and field-aligned current systems, which contribute to the asymmetry in the ring current [71,72].

Impact of
To investigate the impact of geomagnetic storms and auroral substorms on the ionosphere over the Brazilian equatorial region, we analyzed the Diono current.By employing the MRA-MODT technique, we separated the DP2 and Ddyn components from the Diono current.The fluctuations in the DP2 current are shown in Figure 3f. Figure 3g depicts the equatorial electrojet (EEJ) in the Brazilian longitudinal sector.The red legend represents the day-to-day variation of the EEJ observed from 6 to 30 June 2015, while the black legend represents the average EEJ during five magnetically quiet days in June 2015.The DP2 current is linked to the prompt penetration of the magnetospheric convection electric field at the magnetic equator during the main phase of geomagnetic storms [14][15][16]18,73].In contrast, Ddyn represents the magnetic signature of ionospheric disturbance dynamo induced by Joule heating in the auroral zone (Blanc and Richmond 1980) [24].In the low-latitude ionosphere, Ddyn exhibits a westward orientation during the day and an eastward orientation during the night.These phenomena have been investigated extensively by various researchers [17,24,54,58,60,65,74].
On 8 June, we observed a geomagnetic storm driven by a co-rotating interaction region (CIR) followed by high-speed solar wind streams (HSSWs).Figure 3a shows the increase in the solar wind speed from ~371 km/s at 00:17, reaching a maximum speed of ~680 km/s at approximately 11:52 UT on 8 June and remaining at an average of 500 km/s until 16 June 2015.Figure 3b shows that IMF Bz oscillated rapidly during the early hours of 8 June and turned southward at 06:25 UT, with a minimum southward incursion of ~−20 nT for 1 h.During the southward incursion of IMF Bz, the y component of the interplanetary electric field (IEyF) turned eastward (see Figure 3c), and it reached a maximum eastward amplitude of ~9 mV/m at 06:25 UT.As shown in Figure 3d, the auroral westward electrojet (AL black legend) activities intensified on 8 June, starting from 06:00 UT, reaching a maximum westward value of −1068 nT at 07:23 UT, and fluctuated continuously until 18 June 2015, under the influence of HSSWs (see Figure 3a).The southward incursion of the IMF Bz on 8 June triggered a moderate geomagnetic storm, as seen in the SYM-H and ASYM-H indices in Figure 3e.The main phase of the 8 June geomagnetic storm lasted 1.26 h.The recovery phase started at 07:40 UT on 8 June and lasted until 20 June 2015.
Another geomagnetic storm commenced from 18:38 UT on 22 June to 04:25 UT on 23 June 2015.Before the commencement of the 22 to 23 June geomagnetic storms, three interplanetary shocks (IS) hit the Earth's magnetic field [42,43,73].
The first IS arrived on 21 June at 17:00 UT, and the solar wind speed was increased from 300 to 400 km/s.This compressed the magnetopause and caused a sudden increase of approximately 40 nT in the SYM-H index, which lasted for 4.5 h (Figure 3d).However, the IMF Bz component remained positive, so no auroral electrojet activities followed, as seen in AU and AL (see Figure 3e) during this period.The second IS arrived on 22 June at approximately 04:45 UT and was accompanied by little enhancement in the solar wind speed from 400 km/s to ~500 km/s, which provoked a short-duration of 20 to 27 nT increase in the SYM-H index (see Figure 3d).Unlike the two previous events, the arrival of the third IS was followed by giant coronal mass ejection (CME) at 18:38 UT on 22 June.The CME was accompanied by a significant increase in solar wind speed (Vsw) from ~450 to ~712 km/s at 18:55 UT (Figure 3a).This CME caused significant oscillations in the IMF Bz component (see Figure 3b), provoking an increase of approximately 88 nT in the SYM-H index at 18:38 UT (see Figure 3d).At the IS3 arrival, the IMF Bz component was intense and directed southward (see Figure 3b), which led to magnetic reconnection and the development of a strong geomagnetic storm, with a minimum SYM-H index value of −208 nT at 04:25 UT on 23 June 2015, and a maximum of 363 nT of the ASYM-H index at 19:50 UT on 22 June 2015.
One of the main features of this storm was the multiple large-amplitude fluctuations of IMF Bz during the main and recovery phases.The IMF Bz was first largely negative (southward) from 18:38 to 19:24 UT with the minimum of ~−39 nT at ~19:24 UT on 22 June, which was the most negative excursion measured in 2015 [42,43] then turned positive (northward) for ~1.5 h with the maximum of ~30 nT at 20:36 UT.From ~20:55 UT to 21:15 UT, IMF Bz was directed southward, and from ~21:25 UT, IMF Bz turned northward and remained until 00:50 UT on 23 June.Between 00:55 and 01.8 UT, there were rapid smallamplitude fluctuations of IMF Bz, and from ~01:55 UT, the Bz component decreased sharply to −28 nT at 02:05 UT on 23 June and remained intensively negative until 05:30 UT.The IMF Bz oscillations continued until ~12:00 UT during the recovery phase, which started at ~04:25 UT on 23 June, as seen from the SYM-H index behavior in Figure 3d.Following these IMF Bz oscillations, the interplanetary electric field (IEyF) component varied between −18 mV/m at 20:36 UT and +27 mV/m at 19:24 UT on 22 June and −21 mV/m at 00:36 UT and 20 mV/m at 02:06 on 23 June (see Figure 3c and Table 1).
Table 1.The minimum and maximum magnitudes of the parameters plotted in Figure 3 and their respective times of occurrence in universal time.To examine the impacts of geomagnetic disturbances on the Brazilian equatorial ionospheric electric field during the June 2015 geomagnetic storms, Figure 3f,g present the DP2 current fluctuations at magnetic stations in Belem (dip lat: −0.47 • in the red legend) and in Alta Floresta (dip lat: −3.75 • in the black legend) and the equatorial electrojet (EEJ: day-to-day variations of EEJ in the red legend and average of five magnetically quiet days presented in black legend in Figure 3g) in the Brazilian longitudinal sector, respectively.

Day
On 8 June, we observed a remarkable increase in fluctuations of DP2 current amplitude in Belem (red legend) during the daytime in the early stage of the recovery phase, but the amplitude of DP2 at Alta Floresta (black legend) was less conspicuous.This observation provides evidence of the prompt penetration of the magnetospheric convection electric field (EPPEF) into the Brazilian equatorial ionosphere, as the Earth's ionosphere was influenced by HSSWs during the recovery phase on 8 June 2015.During the recovery phase of the 8 June storm, no significant decrease in the amplitude of the eastward EEJ was observed.The notable decrease in eastward EEJ amplitude during the daytime on 15 and 17 June 2015, showed the effect of westward disturbance dynamo on EEJ current.
During the main phase of 22-23 June 2015, particularly on the dayside on 22 June, conspicuous fluctuations in the amplitude of DP2 current between 18:25 and 20:00 UT in Belem (dip lat: −0.47 • presented in the red legend) and Alter Floresta (presented in the black legend) signaled the penetration of the PPMEF into the Brazilian equatorial ionosphere.The eastward maximum amplitude of the EEJ was 88.96 nT at 19:26 UT on 22 June, showed the effect of eastward PPMEF in the Brazilian magnetic equator.During the recovery phase on 23 June, EEJ westward flow was caused by westward disturbance dynamo electric field [22,24,59,74].
Table 1 shows the time intervals in Coordinated Universal Time (UTC) for the interplanetary parameters, including IMF Bz, Vsw, and IEyF, to change from minimum to maximum amplitudes during space weather events in June 2015.Additionally, it displays the time intervals for the AU, AL, SYM-H, ASYM-H, EEJ, and DP2 currents to shift from minimum to maximum amplitudes in response to variations in the interplanetary parameters.Days with significant magnetic disturbances are highlighted in red, while all other days are shown in black.As shown in Table 1, there was no relationship between the time when the IEyF reached its maximum eastward amplitude on 8 June and the time when the EEJ reached its peak eastward amplitude.One possible explanation for this is that the main phase of the 8 June geomagnetic storm occurred at night in the Brazilian longitudinal sector.However, a correlation was found between the time at which Vsw reached maximum speed and the periods of DP2 current fluctuations in Belem and Alta Floresta.On 22 June, the IMF B reached the minimum southward incursion of −38.98 nT at 19:24 UT, while the IEyF reached maximum eastward amplitude of 26.57mV/m at the same time.The impact of the IEyF's maximum eastward amplitude on the EEJ was evident, as it reached a maximum eastward value of 88.96 nT at 19:26 UT, just two minutes later.This serves as clear evidence of the eastward PPMEF in the Brazilian magnetic equatorial ionosphere during the dayside main phase on 22 June 2015, which was a geomagnetic storm.The PPMEF also caused fluctuations in the DP2 currents in Belem and Alta Floresta on 22 June, as seen in columns 11 and 12, respectively.Moreover, the westward auroral electrojet reached a maximum negative amplitude of −1508 nT at 18:57 UT, 19 min after the IMF Bz turned southward, whereas the eastward auroral amplitude reached a maximum positive amplitude of 1473 nT at 20:08 UT, 44 min after the IMF Bz reached its lowest southward incursion.

Phase Comparison between the Sq and Diono Currents in Belem from 6 to 30 June 2015
Figure 4 shows a phase comparison between the Sq (Figure 4a) and Diono currents (

Effect of Cowling Conductivity on the Strength of Ddyn
The equatorial electrojet is an eastward current that flows daily along the dip equator and produces a large variation in the magnetic field [51,78].Close to the magnetic equator, the magnetic inclination angle (I) is approximately zero (I ≈ 0).In the E-region, the winds sweep the ions and leave behind the electrons.Electrons are not affected because the collision frequency with the neutral particle is much smaller than the gyro-frequency.Owing to the Hall effect, we have: This current cannot flow through the layer, and charges accumulate on its borders.The accumulated charges generate an upward electric field (E z ).The vertical current cannot flow because the situation is considered to be in a steady state.Then, a current arises to cancel the Hall current, which is called the Pedersen current, as shown in Equation (15).This implies that: The zonal current in Equation ( 18) is calculated by adding the current produced by (E z × B ), the Hall current, and the current E X (Pedersen current).
Equation ( 19) represents the cowling conductivity, which is responsible for the intensification of the electrojet at the magnetic equator.
Figures 5 and 6 display the development of Ddyn in relation to the auroral electrojet activities and effect of cowling conductivity on the Ddyn amplitude variations in the equatorial latitude.In Figure 5, we present the diurnal component of Ddyn (see Figure 5c) at magnetic stations in Belem (dip lat: −0.47 • presented in red) and Alta Floresta (dip lat: −3.75 • presented in the black legend) alongside the Akasofu index (Figure 5a plotted on the left-hand side in the red legend), the AE index (plotted on the right-hand side and presented in the black legend), and the ASYM-H index in Figure 5b.The Akasofu index shows the coupling between the magnetosphere and solar wind during geomagnetic storms and provides insights into the amount of energy poured into the magnetosphere by the solar wind-magnetosphere interaction and the response of auroral electrojet activities (AE index) to variations in the Akasofu index from June 6 to 30, 2015.During the main phase of the geomagnetic storm on June 8, 2015, the AE index responded promptly to the Akasofu index without a time delay.However, during the main phase of the geomagnetic storm on June 22-23, 2015, the AE index response to the Akasofu index variations showed a time delay at 23:00 UT on June 22 owing to the intensification of the storm.This phenomenon suggests decoupling between the ring current and auroral processes during intense storms, with respect to their dependence on solar wind energization [39].The Figure 5b presents the ASYM-H index.The ASYM-H index is a measure of the strength of the partial ring current.The ring current is asymmetric during the main phase of the 8 and 22-23 June 2015 geomagnetic storms.Obviously, the level of disturbance displayed by ASYM-H index during the main phase of the two storms corresponded with amount of energy derived from the solar wind-magnetosphere interaction which determined the amount of the energy input into the auroral ionosphere and magnetosphere as shown in Panel a by the Akasofu index and AE index.The ring current asymmetric magnitude was at the maximum during the main phase of 22-23 June and during the early stage of the recovery phase on 23 June (see Figure 5b), and the disturbance dynamo (Ddyn) presents in Figure 4c showed a remarkable eastward and westward amplitude on 23 June 2015, which corresponded with the period that intense magnetic perturbation was observed as seen on ASYM-H and AE indices.Our observations established that the strength of Ddyn at the low-latitude stations correlates with the magnitude of the auroral electrojet intensification in high-latitude auroral zone [35,54,74,79] (2017) and Astafyeva (2017Astafyeva ( , 2018)).And, on 23 June, the maximum eastward amplitude in Belem was 39 nT at 05:01 UT and maximum westward amplitude was −51.20 nT at 15:00 UT and in Alta Floresta, the maximum eastward amplitude was 20.75 nT at 07:19 UT and the maximum westward amplitude was −15.75 nT at 17:22 UT.
Figure 6 presents a normalized amplitude of Ddyn at two locations: Belem (refer to Figure 6b) and Alta Floresta (refer to Figure 6c).These amplitude values are juxtaposed with the Akasofu parameter (located on the left-hand side of Figure 6a and indicated by a red legend) and the AE index (found on the right-hand side of Figure 6a, denoted by a black legend).
Figure 6b,c represent the anti-Sq current amplitudes presented in Figure 4f.The X-axis is labeled in days, while the Y-axes of Panels b and c display the periods of oscillations of the Ddyn amplitudes in hours.The color index is shown in the leftmost column, representing the normalized amplitude of Ddyn.
A clear observation from Panels b and c is that the dominant periods of Ddyn, marked by a high-power spectrum, typically span approximately 24 h (diurnal).Notably, the highest amplitude was recorded during the strong storm recovery phase on 23 June 2015, which coincided with remarkable fluctuations in the AE and ASYM-H indices.The variations in the magnitude of the amplitude in Belem and Alta Floresta can be attributed to the effect of cowling conductivity, which is strongest in Belem and tends to enhance the amplitude of Ddyn there.
Table 2 displays the minimum and maximum indices of Akasofu, AE, ASYM-H, and Ddyn, as well as the time at which each index reached its individual minimum and maximum values, in response to the space weather event that occurred in June 2015.Days with magnetic disturbances are highlighted in red.The data in Table 2 reveal that the time lag between the maximum Akasofu index and the maximum eastward amplitude of the AE index was 1.40 h on 8 June and 2.02 h on 9 June.This delay may have hindered the development of Ddyn during the recovery phase of the geomagnetic storm on 8 June.On 14 June, we observed that the AE index reached its maximum eastward value of 1418 nT at 09:42 UT, whereas the Akasofu index reached its maximum value of 2.20 × 10 7 GW at 14:49 UT.This indicates that the AE index reached its maximum amplitude 5.07 h before the Akasofu index, suggesting a long duration of coupling between the solar wind and the magnetosphere.This could have contributed to the development of Ddyn at the Brazilian magnetic equatorial latitude from 15 to 17 June 2015.
During the geomagnetic storm of 22-23 June 2015, we observed a rapid response between the Akasofu and AE indices.Table 2 illustrates that on 22 June, the Akasofu index reached peak value of 2.32 × 10 8 GW at 20:13 UT, while the AE index attained maximum eastward amplitude of 2698 nT at 20:10 UT.This signifies that the AE index reached maximum eastward amplitude three minutes after the Akasofu index reached peak value.The significant energy transfer from the solar wind to the magnetosphere due to solar wind-magnetosphere coupling resulted in Joule heating in the auroral high-latitude ionosphere, which drove the disturbed equatorward wind and led to the development of the westward Ddyn.We observed the maximum westward amplitude of −16.  2 reveals that the time lag between the Akasofu maximum index and the AE maximum index was smaller during the 22-23 June geomagnetic storm than during the 8 June geomagnetic storm.Additionally, the data in Table 2 show that the disturbance wind first reached Belem, which is closer to the magnetic equator before Alta Floresta.Zhang et al. (2017) found that disturbance winds tend to approach the magnetic equator more easily in the Northern Hemisphere's Asian sector, because the magnetic equator is located in the Northern Hemisphere and is associated with greater heating rates in the summer hemisphere and prevailing summer-to-winter winds [79].The geomagnetic storm of the 8 June main phase was on the nightside in the Brazilian longitudinal sector (03:20 LT to 04:46 LT; Table 1); therefore, no remarkable geomagnetic activities were observed during the main phase.The Earth's ionosphere was under HSSWs during the recovery phase from 07:45 UT on 8 June to 00:23 on 17 June 2015.During the recovery phase, from 12:00 UT on 8 June, we observed positive ionospheric storms from the Northern Hemisphere to the Southern Hemisphere (Figure 7c-l

SAMI2 Model Run Results
Figures 9 and 10 present a comparative analysis of the daily average values between the SAMI2 model vertical total electron content (SAMI2 model-VTEC) and global positioning system (GPS-VTEC) within the Brazilian region from 6 to 30 June 2015.The Y-axis denotes the VTEC (in TECU), and the x-axis represents the magnetic latitudes in degrees.The SAMI2 model-VTEC is denoted in the red legend and the GPS-VTEC is identified in the black legend.The root mean square errors (RMSE) are indicated within each figure for each day.We input the daily average of the Ap index into the SAMI2 model to describe the geomagnetic activity levels, daily average of F10.7 cm and the three-month average of F10.7 cm was also input into the SAMI2 model.The HWM14 model provides neutral wind velocity, and the nrlmsise00 model provides neutral densities and temperature.The Fejer-Scherliess empirical model [26] was used as an input to the SAMI2 model to describe the ionospheric electric field during the magnetically quiet and auroral substorm periods.The minimum RMSE was documented on 29 June 2015 (1.62 TECU), as the Brazilian longitudinal sector ionosphere recovered from geomagnetic disturbances.The heightened RMSE on 22 June 2015 is attributed to the impact of geomagnetic disturbances arising from the competition between the eastward prompt penetration of magnetospheric convection electric field (EPPMEF) and the westward disturbance dynamo electric field.Similarly, the elevated RMSE on 23 June is linked to the westward disturbance dynamo electric field that induced westward flow of the equatorial electrojet (EEJ) during the daytime.

Effects of IMF Bz Oscillations on DP2 Current Fluctuations in Brazilian Equatorial Latitude
The interaction between the solar wind and the Earth's magnetosphere induces the formation of region 1 (R1) and region 2 (R2) field-aligned current (FAC) systems, along with the associated magnetospheric plasma convection.Generally, the R1 and R2 FACs create an electric field that can generate a high-latitude current system exhibiting quasiperiodic fluctuations, which are closely correlated with variations in the interplanetary magnetic field (IMF) Bz component.The magnetic disturbance observed at the high-latitude ground level is primarily attributed to Hall currents, as the magnetic effects of Pederson currents cancel those of the FAC, but with opposite signs [29,80].Consequently, DP2 fluctuations at high latitudes result from the combined influence of Hall conductivity and the solar wind-driven convection electric field [81][82][83].These DP2 fluctuations rapidly extend to the equatorial ionosphere (refer to Figure 3f) through the transmission of TM0 mode waves in the Earth-ionosphere waveguide [80].11 and 12, respectively.On the night side of the main phase, DP2 currents did not show any remarkable fluctuations.This phenomenon established the observation of Nishida et al. (1968) that the interaction between the geomagnetic field and IMF Bz produced two reconnection locations, at the dayside (sunward side) and nightside (tailside), and that the field-aligned currents that originate from nightside reconnection may not be the cause of the DP2 current system, because the transit time of the event from the tailside reconnection region (away from the ground by tens to a few hundred RE) to the ground is much greater than the DP2 fluctuation time.This implies that field-aligned currents that originate from dayside reconnection drive the DP2 current system [16,80].
The equatorial electrojet also showed a remarkable response to the penetration of the y-component of the interplanetary electric field on 22 June to the Brazilian equatorial station.The EEJ reached the eastward maximum amplitude of ~88 nT (19:26 UT) within 2 min after the IEyF reached maximum eastward amplitude of ~27 mV/m at 19:24 UT on 22 June 2015.This is clear evidence of eastward prompt penetration of magnetospheric convection electric field into the equatorial ionosphere and enhancing the eastward EEJ component (see Figure 3f) [31,44].

Development of Disturbance Dynamo Electric Field during Recovery Phase of June 2015 Geomagnetic Storm and Associated Positive and Negative Ionospheric Storms
The disturbance dynamo induced by storm-time disturbance winds becomes prominent during the recovery phase and persists for extended durations [84,85].Blanc and Richmond (1980) [24] systematically elucidated the disturbance dynamo concept by simulating the impact of mid-latitude disturbance winds on the ionospheric electrodynamics.Their findings revealed equatorward winds propelled by high-latitude Joule heating deflecting westward at middle latitudes due to the Coriolis force.These deflected winds drive an equatorward current, creating a poleward electric field that counterbalances wind-induced current.This electric field induces an eastward Hall current, resulting in positive (negative) charges accumulating at the dusk (dawn) terminator and establishing a dusk-dawn electric field.Richmond et al. (2003) [86] extended these simulations by combining winds driven by high-latitude plasma convection and Joule heating, thereby highlighting the impact of strong nighttime disturbance winds on the global disturbance dynamo.The nighttime westward winds drive an equatorward current, leading to the accumulation of positive charges at approximately midnight in the low-latitude ionosphere, causing an eastward electric field post-midnight and a westward electric field in the pre-midnight sector.  Figure 4 shows the process of isolating the effect of the disturbance dynamo on the equatorial electric field.The effect of the westward disturbance dynamo can be observed on the EEJ during the day of 23 June 2015 (see Figure 3g).As shown in Panel e of Figure 3, the recovery phase of the 22-23 June 2015 geomagnetic storms commenced during postmidnight (01:25 LT) on 23 June 2015.Figure 5a shows that the auroral electrojet current intensified during the main phase and early stage of the recovery phase on 23 June.The equatorward storm-time disturbance wind due to Joule heating at the auroral zone can corotate with the Earth to later local times [84,87,88] and drive a large westward disturbance current at low latitudes, as observed in the Brazilian equatorial latitude on 23 June 2015 (see Figure 5b,c).The remarkable amplitude of Ddyn that developed on 23 June can be associated with two distinct factors; (1) The time that the recovery phase commenced on 23 June (01:25 LT) gave the equatorward storm-time wind an ample time to develop and corotate with Earth to dayside superimposing on the background wind and driving large westward ionospheric electric current.(2) The intensification of the auroral electrojet activities at the high latitude which is capable of depositing heat energy into the thermospheric global circulation as we have seen on 22-23 June 2015 in Figures 5a and 6a.
The westward disturbance dynamo caused downward drift of the E × B drift velocity on 23 June 2015 (see Figure 7b).The downward drift of the vertical E × B drift velocity lowers the ionosphere to the region of high recombination rate and drives negative ionospheric storms, as observed in Figures 7c-l and 8c on 23 June 2015.And, the positive ionospheric storms that characterized the nightside of the main phase and early recovery phase on 23 June can be attributed to the eastward disturbance dynamo electric field [42][43][44].
In conclusion, we studied the impacts of geomagnetic storms of June 2015 on the Brazilian equatorial and low-latitude ionosphere, and our results showed the following: i.
The field-aligned currents that originate from dayside reconnections drive the DP2 current system as we have seen in Figures 11 and 12, respectively.ii.
During the main phase of the 22-23 June geomagnetic storm and in the early stage of the recovery phase, we observed the intensification of auroral activities due to the interaction between the solar wind and the magnetosphere, which caused considerable energy input into the high-latitude auroral ionosphere (see Figures 5a and 6a).The deposition of energy at high latitudes heats the thermosphere, drives equatorward wind surges, and changes global circulation.The change in the global circulation drives downwelling at low latitudes, which decreases molecular species and causes a slight positive ionospheric phase during the main phase on the nightside and early stages of the recovery phase (see Figure 8c).The dynamo effect of the altered wind circulation opposes the normal diurnal variation on 23 June (see Figure 7b), with downward E × B drift velocity during the day and drives the negative ionospheric storm during the daytime on 23 June 2015 (see Figures 7c-l and 8c).
The westward flow of EEJ on 23 June was caused by westward disturbance dynamo electric field.iii.
The interplay between the eastward prompt penetration of magnetospheric electric field and westward disturbance dynamo electric field caused short-lived positive ionospheric storms during the main phase on 22 June 2015 (see Figure 8c) iv.
The positive ionospheric storm observed during the recovery phase of the 8 June geomagnetic storm was not connected to the eastward prompt penetration of the magnetospheric convection electric field and HSSWs, but the possible driver could be solar flux (F10.7 cm).v.
A comparison between the SAMI2 model-VTEC and GSP-VTEC revealed notable dis parities in the context of this study.Specifically, the SAMI2 model underestimated the VTEC from −9 • to −24 • magnetic latitude in the Brazilian longitudinal sector from 6 to 17 June 2015.However, satisfactory agreement with the GPS-VTEC was observed from −9 • to 10 In summary, we observed that the EUVAC model that provides the EUV flux overestimated the ionization during the daytime, particularly from 10:00 to 20:00 LT.These results are not presented to avoid making them cumbersome.Another remarkable observation is that the summer-to-winter hemispheric circulation of meridional neutral wind raised the ionized plasma between ±10 • magnetic latitude and lowered it in the Southern Hemisphere, which might be the reason for the overestimation of VTEC by the SAMI2 model between ±10 • , as shown in Figures 9 and 10.Finally, we observed that the SAMI2 model failed to reproduce the negative ionospheric storm caused by the westward disturbance dynamo on the VTEC during the daytime on 23 June, as we observed when we compared the SAMI2 model-VTEC on 20 June (quiet day) with that on 23 June (disturbed day).However, when we replaced the Fejer-Scherliess empirical model with the ExB drift velocity estimated from the ground-based magnetometer, the SAMI2 model reproduced the effect of the westward disturbance dynamo on the VTEC on 23 June 2015 (we did not present this result in this study to avoid redundancy) [19].Hence, we recommend incorporating disturbed winds into the SAMI2 code to achieve optimal performance during geomagnetic storms.
June 2015 Geomagnetic Storms on Brazilian Ionosphere 3.1.1.The State of the Interplanetary Medium from 6 to 30 June 2015

Figure 3 .
Figure 3. State of the interplanetary medium from 6 to 30 June 2015.(a) Velocity of the solar wind (Vsw).In addition, (b) shows the IMF Bz (interplanetary magnetic field in the Z direction), (c) presents the y-component of interplanetary electric field, (d) portrays the aurora electrojet (AL: aurora lower boundary, AU: aurora upper boundary), (e) demonstrates the H-component symmetry of Earth's magnetic field observed at various low latitudes, and (f) illustrates DP2 currents in Belem (dip lat: −0.47 • ) and Alta Floresta (−3.75 • ).Further, (g) shows the day-to-day variations of EEJ current (red legend) and magnetically five-quiet day average (black legend).
Figure 4 shows a phase comparison between the Sq (Figure 4a) and Diono currents (Figure 4c) from 6 to 30 June 2015.The wavelet-based method, called semblance analysis, compares the local phase relation between Sq and Diono currents as a function of time and wavelength (Figure 4b,d present the continuous wavelet transforms of Sq current and Diono current, respectively).The normalized amplitude of the local phase relationship between the Sq and Diono currents are presented in Figure 4e.The red legend in the color

Figure 4 .
Figure 4. Phase comparison between Sq and Diono currents at Belem from 6 to 30 June 2015.(a) Sq current in Belem from 6 to 30 June 2015.(b) Continuous wavelet transforms (CWT) of the Sq current at Belem, and (c) Diono current at Belem from 6-30 June 2015.(d) Continuous wavelet transforms (CWT) of the Diono current at Belem.(e) Semblance between Sq and Diono current.The value −1 represents the anti-Sq current, and +1 represents the phase correction of the Sq and Diono currents.(f) Amplitude of semblance analysis between Sq and Diono currents.

Figure 5 .
Figure 5. Variation in diurnal component of Ddyn alongside the AE index, Akasofu index, and ASYM-H index from 6 to 30 June 2015.(a) Day-to-day variation of the AE index (right-hand side of (a), presented in the black legend) and Akasofu index left-hand side of (a) presented in the red legend (b) ASYM-H index presented in (b).(c) Diurnal component of Ddyn at magnetic stations in Belem represented in the red legend and Alta Floresta represented in the black legend.

Figure 6 .
Figure 6.Development of anti-Sq amplitude computed from semblance cross-correlation wavelet analysis between Sq and Diono currents in relation to AE index activity and Akasofu index from 6 to 30 June 2015.(a) AE index plotted on the right-hand side (black legend) and Akasofu index plotted on the left-hand side (red legend) in (a).(b) Anti-Sq of semblance analysis of local phase comparison between Sq and Diono current at Belem (dip lat: −0.47 • ).(c) Anti-Sq of semblance analysis of local phase comparison between Sq and Diono currents at Alta Floresta (dip lat: −3.75 • ).
. The Ddyn was observed six days after the main phase of 8 June (from 15-17 June 2015) geomagnetic storm and the maximum eastward amplitude on 15 June was 12 nT at 04:56 UT and maximum westward amplitude was −20 nT at 16:03 UT in Belem and maximum eastward amplitude was 7.45 nT at 06:58 UT and maximum westward amplitude was −12.30 nT in Alta Floresta.But, we started to observe the Ddyn development since the late stage of the main phase at 19:25 UT on 22 June, and the maximum westward amplitude was −16.33 nT in Belem and −20.13 nT at 21:31 UT in Alta Floresta, which established the interplay between the eastward PPMEF and westward DDEF during the main phase on 22 June as first mentioned by Singh et al.
33 nT at 19:25 UT in Belem and −20.13 nT at 21:31 UT in Alta Floresta on 22 June 2015.The maximum westward amplitude of Ddyn was registered on 23 June at 15:00 UT in Belem, with a value of −51.20 nT and at 17:22 UT in Alta Floresta with a value of −15.75 nT.Table

Figure 7
Figure 7 presents the effects of geomagnetic disturbances on vertical E × B drift velocity and ionospheric VTEC in the Brazilian longitudinal sector.Panel a shows the day-to-day variability of the SYM-H index from 6 to 30 June 2015.Panel b illustrates the day-to-day variability of vertical E × B drift velocity (DVDV red legend), plotted along with average of five magnetically quiet days in the month of June 2015.Panels c-l, present day-to-day variability of VTEC from 6 to 30 June 2015 (DVTEC red legend) and the average of five magnetically quiet days in June 2015 (QVTEC black legend).The X-axes represent the days in universal time (UT).

Figure 7 .
Figure 7. Temporal variation of VTEC alongside the SYM-H index and E × B drift velocity from 6 to 30 June 2015.(a) Day-to-day variability of SYM-H index from 6 to 30 June 2015.(b) Day-to-day variability of vertical drift velocity (DVDV: red legend) and average of five magnetically quiet days in June 2015 (QVDV; black legend).(c-l) Day-to-day variations of VTEC from 6 to 30 June 2015 (red legend) and average of five magnetically quiet days in June 2015 (black legend).
) as the value of the DVTEC (red legend) from 8 to 16 June 2015 was greater than the average of five magnetically quiet days (QVTEC black legend).The positive ionospheric storm was strengthened in the Southern Hemisphere EIA crest (from −12.51 • to −23.77 • dip latitude) on 8 June 2015 owing to enhancement in upward drift velocity of the ionospheric plasma, as seen in Figure 7b.Moving to the 22-23 June 2015, the main phase of the geomagnetic storm spanned from 18:38 UT on 22 June to 04:25 UT on 23 June 2015.On the dayside of the main phase, on 22 June, we observed enhancement in the positive ionospheric storm for a short period in the Southern Hemisphere (from −12.51 • to −23.77 • dip latitude).The upward vertical drift of the plasma on the dayside of the main phase can be attributed to the effect of the eastward PPMEF.The interplay of the eastward prompt penetration of the magnetospheric convection electric field and the westward disturbance dynamo electric field might be the cause of the short duration of the positive ionospheric storm during the main phase on 22 June 2015 in the Brazilian longitudinal sector.During the recovery phase, on 23 June 2015, the westward disturbance dynamo electric field caused downward drift of the vertical E × B drift velocity.The downward drift of the E × B drift velocity lowered the ionosphere to a region of high recombination rate, and this effect was observed on 23 June 2015, as the average of five magnetically quiet days of VTEC in June 2015 was greater than the DVTEC on 23 June 2015.

Figure 8
Figure 8 encompasses three panels, elucidating the following: (a) On the left-hand side, it illustrates the storm-time residual of the equatorial electrojet (∆EEJ) from 8 to 30 June 2015, represented in the black legend.The right-hand side of Panel a portrays the Diono, presented in the red legend.(b) Depicts the day-to-day variability of vertical total electron content (VTEC) in the Brazilian longitudinal sector.Additionally, the SYM-H index is illustrated on the right-hand side and indicated by the red curve.(c) The storm-time residual of VTEC (∆VTEC) from 6 to 30 June 2015, featured on the TEC-MAP (the color bars are represented in the left-hand side of Panels c and b), and the Y-axes are the magnetic dip latitude, while the SYM-H index is plotted on the right-hand side and depicted as a red curve.Figure 8a reveals the influence of the Diono current on ∆EEJ.On the dayside, the southward Diono current prompted ∆EEJ to flow westward, particularly on 15, 17, 21, 23, and 29 June 2015.Conversely, the northward Diono current directed the ∆EEJ eastward on 8 and 22 June 2015.In Figure 8c, the ∆VTEC map reveals VTEC enhancement from 8 to 17 June, except for 15 June, when VTEC depletion was observed within the range of −10 • to −20 • , occurring at approximately 13:00 UT and 20:00 UT.VTEC enhancement was noticeable in both the Northern and Southern Hemispheres, ranging from 10 • in the Northern Hemisphere to −18 • in the Southern Hemisphere.However, on 8 and 18 June, substantial enhancement at the EIA crest, ranging from 10 TECU to 20 TECU, was observed.This enhancement occurred between 12:00 UT and 23:00 UT on the respective days.

Figure 8 .
Figure 8. Spatial-temporal variation of VTEC from 8 to 30 June 2015.(a) Change in EEJ (∆EEJ) from 6 to 30 June 2015, plotted on left-hand-side of (a) (black legend) and Diono current from 6-30 June 2015, plotted on right-hand side of Panel a of Figure 8 (red legend).(b) Day-to-day variability of VTEC from 6 to 30 June 2015 shown on TEC-MAP and day-to-day variability of SYM-H index plotted on right-hand side of (b) (red legend).(c) Change in VTEC (∆VTEC) from 6 to 30 June 2015, shown on TEC-MAP and day-to-day variability of SYM-H index plotted on right-hand side (red legend) from 6 to 30 June 2015.On 22 June, we observed positive ionospheric storms that ranged between 10 TECU and 20 TECU in the Southern Hemisphere from −4 • to −24 • dip latitude between 18:38 UT and 19:26 UT due to the effect of the eastward PPMEF.On 23 June, we observed another less intense positive ionospheric storm that ranged from 2 to 8 TECU from 0:00 UT to 11:00 UT in the Northern Hemisphere and Southern Hemisphere (10 • to −10 • dip latitudes), which can be attributed to the eastward disturbance dynamo electric field.During the recovery phase on 23 June, the VTEC exhibited varying degrees of depletion, notably from 23 to 29 June.On 25 June, a period of VTEC enhancement was detected within the range of −4 • to −20 • , coinciding with the intensification of the ring current, as indicated by the SYM-H index.This caused ∆EEJ to shift eastward during the pre-reversal enhancement period.The VTEC depletion from 23 June to 29 June spanned from −20 TECU to −5 TECU

Figure 10 .
Figure 10.SAMI2 model comparison with GPS-VTEC from 20 to 30 June 2015.Comparing the SAMI2 model-VTEC and GPS-VTEC in the two days preceding 8 June 2015, reveals RMSE values of 2.73 TECU and 4.62 TECU (see 6 and 7 June in Figure9).Currently, the SAMI2 model overestimates the VTEC between −7 • and 10 • magnetic latitude and underestimates it between −8 • and −24 • magnetic latitude.But, from 8-9 June, the RMSE values were 4.99 and 5.03 TECU, respectively, and the increase in RSME on 8 and 9 June can be attributed to the effect of geomagnetic storm of 8 June.Subsequently, from 8 to 15 June 2015, the SAMI2 model-VTEC and GPS-VTEC demonstrated favorable agreement between the −7 • and 10 • magnetic latitudes, but the SAMI2 model underestimated the VTEC between the −8 • and −24 • magnetic latitudes (in the Southern Hemisphere).On 16 and 17 June, a recurrence of the pattern observed on 6 and 7 June occurred.June 18 showed a commendable agreement between the SAMI2 model-VTEC and GPS-VTEC from −24 • to ~−10 • magnetic latitude, but the SAMI2 model overestimated VTEC from −10 • to +10 • magnetic latitude; the RMSE recorded on this day was 4.51 TECU.Proceeding to the geomagnetic storm of 22-23 June 2015, the SAMI2 model-VTEC demonstrated concordance with the GPS-VTEC between −8 • and 10 • magnetic latitudes during the magnetically quiet day on 20 June 2015.However, on 22 and 23 June, we recorded the biggest RMSE of 10.30 and 5.48 TECU (see Figure10) respectively.The SAMI2 model overestimated the VTEC from −10 • to +10 • magnetic latitude in the Brazilian sector.

Figures 11 and 12
display the zoomed plots of the DP2 fluctuations alongside the IMF Bz and its effect on EEJ fluctuations on 8 and 22 June 2015 on the dayside (shown by the black rectangular box).The DP2 currents at Belem and Alta Floresta show a remarkable fluctuation during the daytime in response to northward and southward IMF Bz fluctuations.The effect of DP2 fluctuations in the equatorial zone also caused EEJ current to fluctuate, as presented in Panel a of Figures

Table 2 .
Minimum and maximum Akasofu index, AE index, ASYM-index, Ddyn at Belem, and Alta Floresta along with the time of occurrence.

Ddyn (nT) at Alta Floresta (dip lat: 3.75 • )
• magnetic latitudes from 8 to 15 June 2015.In contrast, the SAMI2 model overestimated the VTEC between ±10 • magnetic latitudes from 16 to 28 June 2015.The highest root mean square error (RMSE) values, notably 10.30 and 5.48 TECU, were recorded on 22 and 23 June 2015, coinciding with periods of intense geomagnetic disturbances.The increase in the RSME value during the main phase on 22 June can be attributed to the interplay between the eastward PPMEF and westward DDEF.We recommend that the SAMI2 model should incorporate disturbance winds to optimally perform during intense geomagnetic storms.vi.