Penetration electric fields observed at middle and low latitudes during the 22 June 2015 geomagnetic storm

Prompt penetration electric fields were detected during the intense geomagnetic storm on 22 June, 2015 by the HF Doppler sounders at middle and low latitudes; Prague (Czech Republic), Iitate, Oarai, Onna (Japan), and Zhongli (Republic of China). The storm was initiated by the storm sudden commencement (SC) at 1833 UT, immediately followed by the main phase with the minimum SYM-H of − 139 nT at 2017 UT. The convection electric field was found to be up to 5.1 mV/m at low latitudes, which is enough to cause ionospheric disturbances such as the enhancement of the total electron content, equatorial ionization anomaly and so on. The electric fields of the SC and main phase are westward in the night (0330 MLT), while eastward in the evening (1930 MLT). The main phase electric field intensified the eastward equatorial electrojets (EEJ) at Huancayo in South America on the dayside (1330 MLT) and the westward EEJ on the nightside at Guam in the western Pacific (0330 MLT) and Tirunelveli in India (00 MLT). The eastward direction of the evening-time electric field agrees with a feature of the evening anomaly of the penetration electric field, while the intensity is suppressed because of the sunlit condition in the evening in the summer solstice. The sudden northward turning of the interplanetary magnetic field (IMF) caused a transition from the main phase to the recovery phase, when the EEJ turned into the counterelectrojet (CEJ). It was found that the HF Doppler sounders detected an impulsive eastward overshielding electric field at 2005 UT on the nightside during the northward IMF-driven CEJ with the duration of 17 min and intensity of 8.0 mV/m at Iitate (0500 MLT) and 11.8 mV/m at Prague (2130 MLT). The impulsive overshielding electric field was accompanied by the anti-sunward plasma flows at middle latitude (50°–60° GML) equatorward of the dusk sunward convection as observed with SuperDARN, which signifies intensification of the Region-2 field-aligned currents at the flow reversal. Substorm signatures are identified with intensified westward auroral electrojet and magnetic dipolarization as detected by the geosynchronous orbit satellite, ETS-VIII. Consequently, the stormtime substorm may provide strong overshielding electric fields to the middle latitude ionosphere on the nightside.


Introduction
Storm time penetration electric fields at low latitude ionosphere have been examined by means of ground-based magnetometers, radars and low Earth orbiting spacecrafts (e.g., Rastogi 1977;Fejer et al. 2007;Kil et al. 2007;Kikuchi et al. 2008;Huang 2015;Astafyeva et al. 2016). When the interplanetary magnetic field (IMF) turns southward, the dynamo that generates the convection electric field and Region-1 field-aligned currents (FACs) is activated around the cusp region in the Open Access *Correspondence: hashi@geosci.jp 1 School of Agriculture, Kibi International University, 370-1, Shichi-Sareo, Minamiawaji, Hyogo, Japan Full list of author information is available at the end of the article outer magnetosphere (Tanaka 1995(Tanaka , 2007. The electric potential and currents are transmitted to the polar ionosphere and then toward the equator by means of the zeroth-order transverse magnetic (TM 0 ) mode waves propagating at the speed of light in the Earth-ionosphere waveguide (Kikuchi et al. 1978;Kikuchi 2014), disseminating the positive and negative potentials to the global ionosphere. The potential difference provides ionospheric electric fields responsible for the motion of plasma in the F-region and DP2-type Hall currents in the E-region.
The equatorial electrojets (EEJs) on the dayside have been examined for many years in studies of the prompt penetration electric fields from high latitudes to low latitudes (Rastogi 1977;Kikuchi et al. 1996;Hashimoto et al. 2011). The penetration electric field suffers severe geometrical attenuation to be less than one tenth of the source intensity at low latitudes. However, the Pedersen currents are intensified at the dayside equator (Kikuchi et al. 1978), because of the anomalously high Cowling conductivity of the equatorial E-region (Hirono 1952). Thus, the magnetometer deployed at the magnetic equator on the dayside well detects the penetration electric field. In other words, the daytime equatorial ionosphere works as a sensitive antenna to the penetration electric field.
During the geomagnetic storms, the H-component of the magnetic field decreases considerably at middleequatorial latitudes, which is caused by the ring current driven by the dawn-to-dusk convection electric field (Ebihara and Ejiri 2000). However, the stormtime magnetic fields are superposed by disturbances due to the ionospheric DP2 currents driven by the dawn-to-dusk convection electric field even at middle latitude (Wilson et al. 2001) and at the equator (Kikuchi et al. 2008;Tsuji et al. 2012). During the storm main phase, the auroral electrojet (AEJ) expands to the middle latitude as low as 55º GML (geomagnetic latitude) (Feldstein et al. 1997;Kikuchi et al. 2008). At the equator, the convection electric field intensifies the eastward EEJ on the dayside during the main phase, while the EEJ decreases to be westward EEJ, i.e., counterelectrojet (CEJ) due to the overshielding during the recovery phase (Kikuchi et al. 2008).
In contrast to the convection electric field associated with the Region-1 FACs, the dusk-to-dawn overshielding electric field originates from the Region-2 FACs that is generated by the dynamo in the inner magnetosphere as a result of development of the magnetospheric plasma convection (Vasyliunas 1972). Convection reduction due to a northward turning of the IMF is one of the possible causes for the overshielding (Kelley et al. 1979). It has also been reported that the overshielding electric field becomes dominant at low latitudes, when the AEJ moves poleward rapidly at the beginning of the storm recovery phase (Kikuchi et al. 2008) and the Region-2 FACs intensify at the substorm expansion (Wei et al. 2009;Hashimoto et al. 2011Hashimoto et al. , 2017. Direct measurements of the stormtime electric fields have been carried out with the incoherent scatter radar (ISR) in Jicamarca (Fejer et al. 2007;Huang et al. 2010) and low Earth orbiting spacecrafts; DMSP (Burke et al. 2000;Huang et al. 2010), C/NOFS and ROCSAT-1 (Kil et al. 2007;Huang 2015;Zhang et al. 2017). In major storms, the ISR detected the electric field with magnitude greater than 3 mV/m at the equator (Fejer et al. 2007). However, only a few observations of stormtime electric fields have been reported from the middle and low latitudes, because of sparse facilities at these latitudes (Blanc et al. 1977). No direct measurements of the penetration electric fields have been reported even for the recent major storm on 22 June, 2015, although quite a few papers have clarified effects of the penetration electric fields on the equatorial plasma bubbles in the postsunset (Cherniak et al. 2019), disturbances in the EEJ (Singh and Sripathi 2017;Astafyeva et al. 2018), and ionospheric and thermospheric density variations (Astafyeva et al. 2016;Ngwira et al. 2019). The present paper reports direct measurements of the penetration electric fields at middle and low latitudes, using the HF Doppler sounders in Europe and western Pacific zones.
The HF Doppler sounders measure Doppler frequency deviations of the HF radio signal reflected from the ionospheric F-region caused by the vertical motion due to the ExB drift of plasma (Davies et al. 1962). The HF Doppler sounder is sensitive to the electric field of short period (1-10 min) disturbances such as the geomagnetic sudden commencement (SC) and geomagnetic pulsations (Davies et al. 1962;Kikuchi et al. 1985Kikuchi et al. , 2016Alperovich et al. 1991;Motoba et al. 2004;Baddeley et al. 2005) and is capable of detecting electric fields of substorms (20 min) (Tsutsui et al. 1988;Hashimoto et al. 2017) and quasi-periodic DP2 fluctuations (30 min) (Reddy et al. 1994;Abdu et al. 1998). Hashimoto et al. (2017) detected eastward overshielding electric field on the nightside that occurred simultaneously with the dayside CEJ during substorms. The results showed that the overshielding electric field is the dusk-to-dawn electric field associated with the Region-2 FACs.
From these observations, we would expect to detect the stormtime electric fields with longer time scales by means of the HF Doppler sounders at middle and low latitudes. In this paper, we report our finding that strong stormtime electric fields were detected simultaneously in Japan and Czech Republic in the post-midnight and evening hours, respectively, during the major storm on 22 June, 2015. During the main phase, the electric field was westward in the post-midnight with the eastward EEJ on the dayside, meeting the dawn-to-dusk convection electric field. In the evening sector, on the other hand, the electric field manifested the evening anomaly with the same direction as on the dayside. In the beginning of the recovery phase, the EEJ turned into the CEJ because of the northward turning of the IMF. We found that an impulsive eastward overshielding electric field with duration of 17 min was observed during the northward IMF-driven CEJ simultaneously in Japan and Czech Republic on the nightside. We will show substorm signatures during the impulsive overshielding and suggest that the stormtime substorm may provide significant overshielding electric fields to the middle latitude ionosphere.
The observed Doppler shifts result from vertical velocities caused by the ExB drift motion of ionospheric plasma. Since the penetration electric field propagates almost instantaneously to the global ionosphere (Kikuchi et al. 1996), we identify the penetration electric field from ionospheric motions with no time difference in two local time zones. In the present paper, the electric fields of geomagnetic storms are observed at multi-points in Japan (MLT = UT + 9), Prague, Czech Republic (MLT = UT + 1), and Zhongli, Republic of China (MLT = UT + 8) as indicated in red in Fig. 1.
The Japanese sounders consist of one transmitter at the Center for Space Science and Radio Engineering, University of Electro-Communications, Chofu, Tokyo, with frequencies of 5.006 and 8.006 MHz, and three receivers at Oarai, Ibaraki; Iitate, Fukushima; and Onna, Okinawa; 120 km, 250 km, and 1500 km away from the transmitter, respectively ( Table 1). The HF Doppler sounders in Prague and Zhongli consist of three transmitters (Tx1, Tx2, Tx3) and one receiver with the sounding frequencies of 3.59 MHz (Prague) and 6.57 MHz (Zhongli) (see details in Laštovička and Chum 2017;Chum et al. 2018). The Doppler frequency data from one of the transmitters, Tx2, are analyzed in this study. The data in Japan are sampled every 10 s and are made smooth through the running average over 60 s to suppress the short-period fluctuations and impulsive noises. The data in Zhongli and Prague are digitized every 30 s.
If the radio wave of HF Doppler reflects in the F-region, where the plasma is magnetized, then the eastward electric field E can be estimated from the observed Doppler shift ∆f caused by vertical plasma motion owing to E × B drift (E: electric field intensity and B: geomagnetic flux density) as:  Table 1 where f, c, I and θ are the radio wave frequency, light speed, inclination angle of the geomagnetic field, and elevation angle of the radio path, respectively (Davies et al. 1962). By applying parameters in Table 1, the electric field is derived from the HF Doppler frequency as below: In order to confirm that the HF Doppler frequency variations are caused by the penetration electric field, we examined the EEJ which is intensified during geomagnetic storms by the magnetospheric electric field penetrating to the equator through the polar ionosphere (Kikuchi et al. 2008;Tsuji et al. 2012 To detect the ionospheric currents flowing in the north-south direction from high latitude to the equator, we used the Y-component of the magnetic field from the INTERMAGNET stations at middle latitudes, Hartland (HAD), UK and Honolulu (HON), Hawaii (Table 2). Furthermore, we used the magnetometer data from Dixon (DIK) and Pebek (PBK) of the THEMIS magnetometer network located at auroral latitudes in Russia (Table 2) to show substorm activities on the nightside. The THEMIS ground magnetometer data are obtained from the coordinated data analysis web (https ://cdawe b.sci.gsfc.nasa.gov).

Solar wind data and geomagnetic indices
In this paper, we focus on the geomagnetic storm that started at 1833 UT on 22 June, 2015 and turned into the recovery at 2007 UT on the same day. The storm (hereafter referred to as 22 June storm) is followed by another development of the ring current on 23 June such that the 22-23 June storm goes into the recovery phase in the beginning of 23 June (Piersanti et al. 2017). However, we focus on the 22 June storm, since the HF Doppler sounders detected the stormtime convection and overshielding electric fields on 22 June. Thus, the storm of our concern goes into the recovery phase at 2007 UT on 22 June. Figure 2a shows three components of the IMF (Bx, By, and Bz), velocity (Vsw), and ion density (Ion N) of the solar wind observed by ACE at the Lagrangian Point L1 during the time interval of the storm. The time of the ACE data is shifted by 35 min so that we can see the correspondence between the sudden increases in the solar wind density and speed and the geomagnetic sudden commencement (SC) on the ground. Figure 2b shows the SYM-H and ASY-H geomagnetic indices from the WDC for Geomagnetism, Kyoto. The sudden increase of the solar wind dynamic pressure caused the SC at 1833 UT on 22 June, 2015. The IMF Bz was slightly negative with − 2.1 nT as averaged over 1 h before the arrival of the ICME. The IMF Bz then dramatically decreased to − 23.6 nT at 1833 UT and reached the minimum value of − 39.5 nT at 1925 UT concurrently with sudden increases in the density and velocity. The average value of the southward IMF was − 29.0 nT over the time interval of 72 min until the northward turning of the IMF Bz at 1945 UT. The solar wind velocity abruptly increased from 475 to 704 km/s and the density jumped from 25 to 44/cc. Subsequently, the density gradually decreased to 25/cc at 1920 UT. As can be seen in Fig. 2b, the SYM-H abruptly increased at 1833 UT by an amount of 110 nT, followed by the steep decrease by an amount of 230 nT until 2015 UT, indicating that the storm ring current developed immediately after the SC, because of the strong southward IMF accompanying the increases in the solar wind speed and density. Figure 3a shows the X-component of the magnetic field at the low latitude, San Juan (SJG) and equator, Huancayo (HUA) in the American zone ( Table 2). The storm occurred when the stations were located in the afternoon  Both the X HUA and X SJG started to increase at 1833 UT, reaching the peak of the amplitude 80 nT at 1836 UT at SJG, whereas reaching the peak of 350 nT at 1843 UT at HUA. The enhanced X HUA was followed by the steep decrease over the period of 1943-2010 UT, when the storm turned into the recovery phase because of the northward turning of the IMF. The EEJ turned into the counterelectrojet (CEJ), indicating the overshielding associated with the convection reduction at the equator (Kelley et al. 1979). It should be noted that the rapid decrease in the ASY-H at 1950 UT (Fig. 2b) is also caused by the northward turning of the IMF. The reduction of X HUA was 580 nT from the peak of the main phase, while the X SJG decreased by 200 nT. The enhancement of the main phase variations at HUA should be due to the EEJ and CEJ, caused by the dawn-to-dusk convection electric field and overshielding electric field associated with the northward turning of the IMF. The EEJ is obtained as the difference between X HUA and X SJG with magnitude of 380 nT under an assumption that the storm ring currents contribute equally to both stations. The main phase EEJ has been detected by the Swarm satellite as the negative magnetic deflections above the EEJ (Astafyeva et al. 2016) and the CEJ at X HUA was shown by Ngwira et al. (2019). It is remarkable that the geomagnetic storm is enhanced at the dayside equator with magnitude three times as large as the storm at low latitude in a similar way to the major storm on 6 November 2001 (Kikuchi et al. 2008). The equatorial enhancement of the geomagnetic storm is caused by both effects of the EEJ and CEJ superposed on the magnetic fields due to the ring current. Therefore, the stormtime EEJ can produce even greater magnetic disturbances than those by the stormtime ring current. Figure 3b shows the X-components from low and equatorial latitudes, Alibag (ABG) and Tirunelveli (TIR) in the Indian zone and Kanoya (KNY) and Guam (GUA) in the western Pacific zone. The storm main phase occurred, when ABG and TIR were located in the midnight (00-02 MLT) and KNY and GUA were located before dawn (03-05 MLT). Unlike the gradual decrease in X SJG at the dayside low latitude (Fig. 3a), the X-component at all of the four stations on the nightside (Fig. 3b) decreased steeply immediately after the SC and reached the minima at 1940 The EEJ HUA denotes the equatorial electrojet defined as the difference between X HUA and X SJG . b Same parameters for the low latitude (Alibag, ABG) and equator (Tirunelveli, TIR) in the Indian zone located in the midnight (00-02 MLT) and for the low latitude (Kanoya, KNY) and equator (Guam, GUA) in the western Pacific zone, located in the pre-dawn (03-05 MLT). c X and Y components of the magnetic field at middle latitudes, Honolulu (HON) in the morning (09 MLT) and at Hartland (HAD) in the evening (19MLT) UT. The asymmetric development is due to the partial ring currents on the nightside that is commonly observed during the storm main phase (Tsuji et al. 2012).

Equatorial electrojets during the storm main phase
The EEJs derived from X TIR and X ABG and from X GUA and X KNY are negative during the main phase, consistent with the EEJ at Tirunelveli analyzed by Singh and Sripathi (2017) and the EEJ in Thailand (Astafyeva et al. 2018;Ngwira et al. 2019). Those results indicate that the westward EEJ is driven even in the nighttime ionosphere at the equator. The nightside EEJs of 30-40 nT are weaker than the dayside EEJ by one order of magnitude, being consistent with the result in Tsuji et al. (2012). It should be noted that the nightside EEJ is not very sensitive to the overshielding electric field associated with the northward turning of the IMF. The similarity in the temporal variations of the day-and night-side EEJs with opposite polarity implies that the EEJs are caused by the dawn-to-dusk penetration electric fields associated with the Region-1 FACs.

X-and Y-component of the magnetic fields at middle latitudes
In order to confirm that the EEJs are driven by the Region-1 FACs, we examine the correlation between the EEJs and the Y-component magnetic field at middle latitudes, which is caused by the Pedersen currents connecting the EEJ with the Region-1 FACs (Motoba et al. 2002Hashimoto et al. 2011. The EEJ on the dayside would be positively correlated with the middle latitude Y-component in the morning sector, while the correlation is negative in the afternoon sector. Figure 3c shows X and Y at Honolulu (X HON , Y HON ) in the morning sector (09 MLT at 19 UT) and X and Y at Hartland (X HAD , Y HAD ) in the dusk sector (19 MLT). The X HON jumped by 50 nT at the SC and monotonically decreased by 200 nT after the SC until 1950 UT, whereas the X HAD remained almost constant with magnitude of about 200 nT. X HON is in the same shape as those at low latitudes, KNY and SJG, which indicates that the X-component is overwhelmed by stormtime ring currents. X HAD , on the other hand, is in a trapezoidal form similar to the EEJ at HUA, which indicates that the X is affected by afternoon/eveningtime DP2 currents flowing at middle latitudes during the storm main phase (Wilson et al. 2001). The strong local time asymmetry in the X may be associated with the high ASY-H index in Fig. 2b.
On the other hand, Y HON increased by 60 nT from the pre-event level in the morning sector, while Y HAD decreased by 80 nT in the evening sector during the storm main phase. It should be noted that the temporal variations of the Y-component are in a trapezoidal form similar to the EEJ and the correlation of Y HON with the dayside EEJ is positive, while negative for Y HAD as predicted from the Pedersen currents flowing from the Region-1 FACs to the EEJs.

HF Doppler observations of the electric field during the main phase
The stormtime electric fields were observed by the HF Doppler sounders at middle and low latitudes in the western Pacific zone, Iitate, Oarai and Onna in Japan and Zhongli in Republic of China, and at middle latitude in the European zone, Prague in Czech Republic. Figure 4 shows electric fields converted from the HF Doppler frequency deviations using Eq. (1). The positive values represent the eastward electric field. When the SC occurred at 1833 UT, the westward electric field sharply increased at Iitate, Oarai, Onna and Zhongli in the post-midnight (0230-0330 MLT), reaching the maximum values within 3 min. The HF Doppler sounder at Prague, on the other hand, observed the eastward electric field in the evening (1930 MLT) that is in the same direction as the dayside EEJ at HUA. The evening-time electric fields have been reported to be in the same direction as in the day for the SC over Japan (Kikuchi et al. 1985(Kikuchi et al. , 2016 and DP2 fluctuation events at the equator (Abdu et al. 1998). The eastward electric field at Prague is indicative of the evening anomaly for the geomagnetic storm main phase. If the evening anomaly is a common property of the penetration electric field (Fejer et al. 2008), it is reasonable to expect eastward electric fields at Prague.
With regard to the intensity of the electric field during the storm main phase, the westward electric field decreases with decreasing latitude to be 3.7 mV/m at Iitate, 5.1 mV/m at Oarai, 3.1 mV/m at Onna, 2.4 mV/m at Zhongli from the pre-event level. Figure 5 (red circles) shows latitudinal profiles of the intensity of the convection electric field at 1836 UT during the main phase, where the positive value refers to the eastward electric field. The electric field of 3.0 mV/m at Prague was less than that at Iitate despite the higher magnetic latitude of 50.2º at Prague. The electric field may have been suppressed by the higher conductivity in the ionosphere over Prague that is sunlit at 1930 MLT in the summer solstice. The simultaneity with the EEJs (Fig. 3) and the latitudinal dependence of the westward electric field (red circles in Fig. 5) indicate that the storm time convection electric fields are transmitted from high latitude to the equator via the middle and low latitude ionosphere on both the day-and night-sides.
It should be stressed that the eastward electric field at Prague is consistent with the positive EEJ at HUA (Fig. 3a) and that the westward electric field at Iitate is consistent with the negative EEJ at GUA and TIR (Fig. 3b), although the duration times of the electric field are shorter than those of the EEJs. The westward electric field in Japan (e.g., Iitate) decreases gradually over 30 min, while the westward EEJ at GUA remains constant in a trapezoidal form during the whole period of the main phase (Fig. 3b). This fact suggests that the reflection height of the HF radio signal tends to cease increasing even under the steady electric field, probably because of the reconfiguration of the F-region ionizations through the production and loss processes (Tsutsui et al. 1988). The time constant of the apparent decay of the electric field remains an issue to be solved in using the HF Doppler sounders as a detector of the stormtime electric field.

Overshielding electric field during the storm recovery phase
In the end of the main phase, the day-and night-side EEJs decreased over the period 1950-2010 UT (Fig. 3). The dayside EEJ turned into the CEJ at HUA because of the northward turning of the IMF (Fig. 3a), while no significant changes are observed on the nightside during the dayside CEJ (Fig. 3b). It is remarkable that the HF Doppler sounders detected impulsive eastward electric fields at Iitate, Oarai (0500 MLT) and Prague (2300 MLT) for about 17 min after 2005 UT (Fig. 4) during the northward IMF-driven CEJ. Figure 5 (blue) shows the intensity of the impulsive electric field at 2014 UT; 8.0, 4.8 and 11.8 mV/m at Iitate, Oarai and Prague, respectively, as measured from the values at the start of the impulsive electric field (2005UT). The impulsive electric field is eastward at all the stations located on the nightside (23-05 MLT), indicating that the electric field is associated with the overshielding during the period of the northward IMF-driven CEJ. It should be noted that the impulsive features are not observed in the EEJs on the nightside as well as on the dayside (Fig. 3a, b). To confirm the impulsive electric field be the overshielding electric field, we examined SuperDARN data in the following section. Figure 6 shows the line of sight (LOS) velocity map constructed with two SuperDARN radars in north America and two radars in Japan looking at the auroral and subauroral latitudes on the dayside (12 MLT to the top) at 1940 UT during the main phase (upper panel) and 2014 UT during the impulsive electric field event (lower panel). The LOS velocity maps are drawn by using on-line plotting tools from the Virginia Tech website (http://vt.super darn.org/). The LOS velocity toward the radar (anti-sunward) is shown in blue and green and away from the radar (sunward) in red color. The upper panel shows the large-scale two-cell convection with sunward flows in the 12-16 MLT sector during the main phase, which is typical for the southward IMF condition (Ruohoniemi and Greenwald 1996). The lower panel shows that the anti-sunward flow (green) appeared at 50°-60° GML equatorward of the sunward convection in the 12-13 MLT that indicates sudden intensification of the southward electric field, which must be provided by the Region-2 FACs (Hashimoto et al. 2011). To show temporal variations of the antisunward flow, we plotted the LOS velocity as a function of time as well as the radar echo intensity detected by the beam 9 of the Blackstone radar in the bottom and upper panels of Fig. 7, respectively. The antisunward plasma flow was observed during the period

Evening anomaly of the main phase electric field
We have shown the stormtime electric fields at middle and low latitudes in the evening and night (Fig. 4), assuming that the electric field is in linear correlation with the HF Doppler frequency through Eq. (1). Although the linearity breaks down for long lasting electric fields, the Fig. 6 Line of sight velocity map constructed with two SuperDARN radars in the northern America and two radars in Japan. The blue and green colors of the velocity refer to toward the radar and red refers to away from the radar. The radars detected the sunward (12 MLT to the top) flow (red) at 1940 UT during the storm main phase and anti-sunward flow (green) equatorward of the sunward flow at 2014 UT during the impulsive electric field event main phase electric field at Prague in the evening (1930 has the eastward polarity same as the eastward EEJ at HUA on the dayside (Fig. 3a). On the other hand, the electric fields over Japan and EEJs at TIR and GUA were westward at the post-midnight (Fig. 3b). Cherniak and Zakharenkova (2016) reported intense plasma bubbles over Europe in the evening that should have been caused by the eastward electric field during the main phase of the 22 June storm that we report in the present paper. The eastward electric field has been known to cause ionospheric disturbances in the evening sector (e.g., Basu et al. 2001). The asymmetric local time features with the evening anomaly has been observed by the HF Doppler sounders during the SC (Kikuchi et al. 1985(Kikuchi et al. , 2016 and DP2 (Abdu et al. 1998) and by spacecrafts during quiet and disturbed periods (Fejer et al. 2008;Huang et al. 2010). The evening anomaly results from the distorted potential distribution caused by the day-night inhomogeneity and Hall effects of the ionospheric conductance (Senior and Blanc 1984;Tsunomura 1999).
At latitudes lower than 30º, the intensity of the convection electric field decreases with decreasing latitude such that 5.1 mV/m at middle latitude (Oarai) is twice stronger than 2.4 mV/m at low latitude (Zhongli) (red circles in Fig. 5). On the other hand, the evening-time electric field over Prague is 3.0 mV/m, which is not subject to the evening enhancement, but is comparable to the intensity at Onna in the night. This feature can be explained with high ionospheric conductivity under sunlit conditions of the ionosphere at Prague in the summer solstice, which may have suppressed the convection electric field. These results lead us to conclude that the evening anomaly of the stormtime electric field expanded to 20 MLT, consistent with that of the SC (Kikuchi et al. 1985).
We have shown that the main phase electric field converted from the HF Doppler frequency deviations decreased monotonically over 30-40 min (Fig. 4), although the EEJ remained almost constant over 90 min during the whole period of the main phase on both the day-and night-sides (Fig. 3a, b). These facts suggest that the production-loss processes in the ionospheric F-region tend to reduce the upward velocity of the reflection height of the HF radio waves, so that the Doppler frequency decreases with time even under a constant electric field driving the motion of plasma (Tsutsui et al. 1988). If we had an incoherent scatter radar, we would have a linear correlation between the electric field and EEJ like those obtained by Anderson et al. (2002) with the incoherent scatter radar and magnetometer at Jicamarca. The quantitative relationship between the electric field derived from the HF Doppler sounder and the EEJ would help establish the HF Doppler sounder system as a tool of the electric field measurement, if we cover the effects of the disequilibrium of the ionization due to the production and loss processes.

Impulsive overshielding electric field during the recovery phase
As shown in Fig. 4, the HF Doppler sounders observed the prominent overshielding electric field for 17 min during the northward IMF-driven CEJ in the beginning of the recovery phase (2005 UT). The impulsive electric field is eastward both at Prague (2100 MLT) and Iitate (0500 MLT) that is opposite to the convection electric field. The blue circles in Fig. 5 represent the intensity of the electric field derived from four HF Doppler sounders except for Zhongli. The intensities of 11.8 mV/m at Prague and 8.0 mV/m at Iitate on the nightside were anomalously strong as compared with a few mV/m of isolated substorms (Hashimoto et al. 2017). Furthermore, the overshielding electric field is more intense than the main phase convection electric field (red circles in Fig. 5) and decreases with decreasing latitude steeply to be almost zero at low latitude (Onna).
The latitude dependence of the overshielding electric field might be due to coupled effects of the decaying Region-1 FACs and growing Region-2 FACs in the transition period from the main phase to the recovery phase. The disturbance dynamo produces significant reversed electric fields at low latitude (Blanc and Richmond 1980). However, the duration (17 min) of the impulsive electric field is too short and the shape is too impulsive for the disturbance dynamo electric field that usually continues for several hours and sometimes may exceed 10 h (e.g., Zhang et al. 2017) Another possible cause of the short-lasting overshielding electric field is an intensification of the Region-2 FACs during the substorm expansion (Tanaka et al. 2010;Hashimoto et al. 2011Hashimoto et al. , 2017. Hashimoto et al. (2017) showed that when the substorm expansion breaks out, the equatorial counterelectrojet often develops on the dayside with sudden increases in the eastward electric field on the nightside and recovers to the pre-event level in 20-30 min. In addition, westward ionospheric currents and anti-sunward plasma flows are often intensified at subauroral latitudes together with the intensification of the eastward currents at the auroral latitudes in the afternoon sector (Hashimoto et al. 2011). As shown in Fig. 6, both the anti-sunward plasma flow at the subauroral latitude and the sunward flow at the auroral latitude as detected by SuperDARN strongly support that the downward Region-2 FACs had developed at the poleward edge of the anti-sunward flow region concurrently with the impulsive overshielding electric field at middle latitudes.
To verify that the impulsive overshielding is associated with the substorm, we examined the auroral electrojet indices (AU, AL) and magnetic H-components at auroral latitudes, Dixon (DIK, 69.06° GML, 02 MLT) and Pebek (PBK, 65.58° GML, 07 MLT) in Fig. 8a, b. The westward AEJ at PBK started to intensify at 1833 UT and both the AU and AL indices enhanced, representing development of the convection electric field during the storm main phase. No appearance of the AEJ at DIK may be due to equatorward shift of the AEJ during the storm main phase (Feldstein et al. 1997;Kikuchi et al. 2008). When the IMF turned northward, the AEJ retreated to the normal auroral latitudes and caused the rapid increase in AU at 1955 UT. In the beginning of the recovery phase, the westward AEJ at DIK and PBK started to increase at 2005 UT concurrently with the decrease in the AL, reaching the minimum of − 1300 nT at 2010 UT. The intensification of the westward AEJs and AL index signifies the onset of substorm (Wiens and Rostoker 1975), starting exactly at the same time as the impulsive overshielding electric field. To ensure the occurrence of the substorm expansion, we further examined a signature of dipolarization of the near-Earth magnetotail magnetic field detected by the Engineering Test Satellite (ETS)-VIII geosynchronous satellite. It is well established that the magnetotail magnetic field in the geosynchronous orbit significantly changes from a taillike configuration to a dipolar configuration as the substorm progresses from the growth phase to the expansion phase (Cumming et al. 1968;McPherron et al. 1973;Kokubun and McPherron 1981). The configuration change of the magnetic field at the substorm expansion onset is called the dipolarization of the near-Earth magnetotail, which has been explained as resulting from diversion of the well-developed cross-tail current into the polar ionosphere (e.g., McPherron et al. 1973;Lui 1978). Consequently, the intensification of the AEJ at a substorm onset is closely related with the dipolarization of the magnetic field in the geosynchronous orbit (Lopez and Rosenvinge 1993).
ETS-VIII is an advanced technology demonstration satellite that JAXA launched on 18 December, 2006 and placed in a geosynchronous orbit at 146º E geographic longitude over the western Pacific. The onboard triaxial fluxgate magnetometer (MAM) observations provide us with signatures of the substorm, i.e., dipolarization of the tail magnetic field (Nose et al. 2014). The Hp (northward) component of MAM is plotted for the 22 June storm in Fig. 8c, which has been corrected in accordance with the Appendix A of Nose et al. (2014).
The Hp gradually decreased with short-term fluctuations as the convection electric field developed during the storm main phase and then substantially jumped up by more than 150 nT at 2001 UT, a few minutes earlier than the onset of the overshielding in the ionosphere. Figure 9 shows zoomed in plots of Hp and He (Earthward) at around the tail magnetic dipolarization. The He component decreases sharply at 20:01:35 UT (vertical dotted line) concurrently with the increase in Hp. The temporal variations of Hp and He match the configuration change from stretched to dipolarized forms. The drastic changes over a few minutes in the near-Earth magnetotail may be associated with the impulsive feature of the overshielding electric field. In spite of the predawn sector of 05 MLT, the magnetic field reconfiguration observed by the ETS-VIII is consistent with the dipolarization signatures in the premidnight hours (Lopez and Rosenvinge 1993;Ohtani et al. 1992;Liou et al. 2002).
The time difference of several minutes between the dipolarization and overshielding may be the time needed for activation of the Region-2 FACs dynamo in the near-Earth magnetotail and propagation from the dynamo to the polar ionosphere. Note that the propagation from the polar ionosphere to low latitude is near-instantaneous by the TM 0 mode wave in the Earth-ionosphere waveguide (Kikuchi and Araki 1979). On the basis of these results, we conclude that the impulsive eastward electric field on the nightside is most likely to be the substorm overshielding electric field found by Hashimoto et al. (2017).
The overshielding electric field detected simultaneously at Prague, Iitate and Oarai shows global penetration of the electric field, which leads us to expect an impulsive EEJ/ CEJ, but we did not see clear signatures as in Fig. 3. As shown in Fig. 5, the overshielding electric field decreases with decreasing latitude steeply to be almost zero at low latitude. The overshielding electric field associated with the Region-2 FACs competes against the convection electric field associated with the Region-1 FACs at the low and equatorial latitudes (Kikuchi 2000;Peymirat et al. 2000;Hashimoto et al. 2011Hashimoto et al. , 2017. Hashimoto et al. (2011) showed that both the Region-1 and Region-2 FACs are intensified during the substorm expansion and that the overshielding electric field often overcomes the convection electric field at low latitudes during isolated substorms. In the storm event, the Region-1 FACs have a large spatial scale extending to middle latitudes and make significant effects at the equator (Kikuchi et al. 2008). If the impulsive overshielding electric field is associated with the substorm, the Region-2 FACs must have a smaller scale extending further equatorward, closer to the HF Doppler sounders at middle latitudes. In this situation, it is likely that the overshielding electric field suffers much severer geometrical attenuation than the convection electric field during its propagation to the equator (Kikuchi and Araki 1979). The severe geometrical attenuation would cause the steep latitudinal profile of the overshielding electric field in Fig. 5 and no significant effects at the

Conclusion
1. Using the HF Doppler sounders at middle and low latitudes and magnetometers at the low and equatorial latitudes, we found that stormtime electric fields penetrated promptly on both the day-and nightsides during the geomagnetic storm on 22 June, 2015. 2. The HF Doppler sounders at Prague, Czech Republic (1930 MLT), Iitate, Oarai, Onna, Japan (0330 MLT) and Zhongli, Republic of China (0230 MLT) successfully detected the eastward electric field in the evening and westward electric field in the night during the storm main phase. The eastward electric field at Prague indicates the evening anomaly of the penetration electric field with the direction same as on the day, which must have caused severe ionospheric disturbances over Europe reported by several other papers. 3. The main phase electric field caused DP2 currents at middle-equatorial latitudes, driving eastward and westward EEJs on the day-and night-sides, respectively. The Y-component magnetic field at middle latitudes detected north-south ionospheric currents connecting the EEJs with the Region-1 FACs developed by the southward IMF. 4. The HF Doppler sounders detected the anomalously strong eastward overshielding electric fields (8.0 mV/m at Iitate and 11.8 mV/m at Prague) on the nightside for 17 min in the beginning of the recovery phase of the storm. The overshielding electric field decreases with decreasing latitude steeply to be zero at low latitude, consistent with no corresponding changes in the EEJ at the equator. 5. The Blackstone SuperDARN radar in the American sector observed anti-sunward plasma flow equatorward of the normal sunward plasma flow during the period of the overshielding electric field, which may indicate development of the Region-2 FACs equatorward of the Region-1 FACs. 6. The impulsive overshielding electric field was found to be associated with substorm signatures; intensification of the westward auroral electrojets on the nightside and dipolarization of the tail magnetic field as detected by the geosynchronous satellite, ETS-VIII. 7. In conclusion, the stormtime convection electric field penetrated to the middle and low latitudes during the main phase, and the stormtime substorm may play a crucial role in providing strong overshielding electric field to the middle latitude ionosphere on the nightside.

Acknowledgements
The results presented in this paper rely on the data collected at Guam, Hartland, Honolulu, Huancayo, Kanoya and San Juan. We thank United States Geological Survey, British Geological survey, Instituto Geofisico del Peru, and Japan Meteorological Agency for supporting its operation and INTERMAGNET for promoting high standards of magnetic observatory practice (http://www.inter magne t.org

Authors' contributions
KKH contributed to the design and implementation of the study, to the analysis of the results and to the writing of the initial draft of the manuscript. TK contributed to interpretation of data and assisted in the preparation of the manuscript. IT and KH performed observations with HF Doppler sounders in Japan. JC and DB performed observations with HF Doppler sounders at Prague and Zhongli. MN contributed to calibration of the magnetometer data onboard the ETS-VIII satellite. KK performed collection of the magnetometer data onboard the ETS-VIII satellite. All authors read and approved the final manuscript.