Seasonal variation of inter-hemispheric field-aligned currents deduced from time-series analysis of the equatorial geomagnetic field data during solar cycle 23–24

The east–west component of magnetic field variation (∆D-component) at Davao station (Philippines, geomagnetic latitude: – 2.22˚N) are used to investigate the characteristics of the long-term Inter-Hemispheric Field-Aligned Currents (IHFACs) based on the time-series analysis from August 1998 to July 2020. Recent in situ satellite and ground-based observations have reported that dusk-side current polarity of IHFAC is often opposite to that of the noon IHFAC, being inconsistent with Fukushima's IHFACs model. We investigated the consistency of the dusk-side IHFAC polarity derived from the observations with the polarity expected from Fukushima’s IHFACs model and examined the solar cycle dependence of IHFACs. It was confirmed that the dusk-side IHFACs during June and December solstices flow in the same direction of the noontime IHFACs, which was consistent with the IHFAC polarities suggested by the Fukushima model. The dusk-side IHFACs around March and September–November months disagreed with the Fukushima model. The ∆D variations clearly showed seasonal asymmetry in the dawn and noon sectors, whereas the ∆D variations in the dusk sector demonstrated seasonal symmetry. Solar cycle dependence of IHFACs was exhibited in the dusk sector. For the dawn and noon sectors, the yearly peak-to-peak ∆D amplitude in the later solar cycle SC24 decreased by about 35% in comparison with the earlier solar cycle SC23. In contrast, the dusk-side yearly peak-to-peak ∆D amplitude increased by about 200%. The dusk-side IHFAC yearly amplitude tended to be in inverse proportion to solar activity.


Introduction
The Inter-Hemispheric Field-Aligned Current (IHFAC) is one of the major current systems causing changes in the geomagnetic field at low and middle latitudes. IHFACs flow from the summer hemisphere to the winter hemisphere in the dawn sector and from the winter hemisphere to the summer hemisphere in the noon and dusk sectors (Fukushima, 1994). Van Sabben (1966) first suggested the existence of IHFACs by using the northsouth difference of equivalent Sq current. The IHFACs are caused by the inter-hemispheric imbalance of the ionospheric solar quiet (Sq) current system at the middle-low latitudes (Van Sabben 1966, 1969, 1970, due to asymmetry of the north-south Sq current vortices that establishes a potential difference between the northern and southern hemispheres (Fukushima, 1979;Tomás et al. 2009). Many numerical calculation studies have predicted IHFAC intensity and supported the dependence of IHFAC current polarity changing with local time (Maeda Open Access *Correspondence: fujimoto@ai.kyutech.ac.jp 2 Faculty of Computer Science and Systems Engineering, Kyushu Institute of Technology, 680-4, Kawazu, Iizuka, Fukuoka 820-8502, Japan Full list of author information is available at the end of the article 1974; Schieldge et al. 1973;Stening 1977;Takeda 1982;Van Sabben 1969, 1970).
Fukushima's IHFAC model is characterized by the following features: (1) IHFACs flow from the summer hemisphere to the winter hemisphere in the dawn sector and the opposite current flows in the noon sector; (2) IHFACs polarity between the noon and dusk sectors is in-phase and (3) the absolute intensity of IHFACs is stronger in both dawn and noon sectors than the dusk sector. The current polarity of Fukushima's IHFACs model has been supported by in situ satellite observations and groundbased magnetometer data analysis as well as the abovementioned numerical studies. The first experimental evidence has been provided from Magsat observations (Olsen 1997). The characteristics of IHFACs have been extensively described by many papers, for example, seasonal climatology of IHFACs by the Ørsted satellite (Yamashita and Iyemori 2002), seasonal, longitudinal, and local-time IHFACs variations by the CHAMP satellite (Park et al. 2011) and the Swarm satellite constellation (Lühr et al. 2015(Lühr et al. , 2019Fathy et al. 2019;Park et al. 2020b), ground-based local-time and seasonal dependence (Bolaji et al. 2012;Shinbori et al. 2017;Abidin et al. 2019) and latitudinal dependence (Owolabi et al. 2018).
In recent years, the IHFACs characteristics in the dusk sector, which disagrees with the concept of the Fukushima's model, have been reported by ground-based magnetometer (Shinbori et al. 2017) and satellite observations (Lühr et al. 2015(Lühr et al. , 2019Fathy et al. 2019;Park et al. 2011Park et al. , 2020aPark et al. , 2020b. These papers show the duskside IHFACs flowing southbound irrespective of season, by using the average value of magnetic field obtained from in situ satellite or ground-based magnetometer observations for each month in 1-h local time (LT) bin. Although these analyses revealed the average distribution of IHFAC polarity, little attention has been paid to the relationship of IHFAC polarities between the noon and dusk sectors for individual days. In addition, their analysis periods were almost too limited to demonstrate the solar cycle dependence of IHFACs, namely less than solar cycle period (11 years) (e.g., Park et al. 2011;Owolabi et al. 2018;Lühr et al. 2019). Shinbori et al. (2017) used long-term (59 years) ground-based magnetometer data for the comparison of IHFAC intensity between higher and lower solar activity periods. They examined the monthly average of IHFAC variations, but disregarded to compare the IHFAC polarities in the noon and dusk sectors for individual days. It is essential to investigate the occurrence rate of the dusk-side IHFAC suggested by the Fukushima's model (hereafter, the dusk-side IHFACs of the Fukushima's model type), in order to conclude the consistency between the dusk-side IHFAC polarity derived from the observations and the IHFAC polarity expected from the Fukushima model.
In order to discuss the presence of the dusk-side IHFACs of the Fukushima model type, it is important to conduct a long-term time-series analysis of IHFAC variations, in terms of the comparison of IHFACs polarities between noon and dusk sectors for individual days. In contrast to the in situ satellite observations which provide the global distribution of IHFACs two-dimensional map (longitude-latitude), the ground-based magnetometer data allows us to investigate the daily IHFAC variations. In other words, the in situ observations cannot provide the comparison of IHFAC variations of each LT of a single day, but the ground-based observations enable us to analyze the relationship of IHFAC polarity among the different LT for individual days. The equatorial D-component magnetic fields are used in this paper, since they include more essential effects of IHFAC variation than the equatorial H-component magnetic fields which are dominated by the equatorial electrojet effect (Yamazaki and Maute, 2017).
Field-aligned currents (FACs) observed at the highlatitude regions are excited by the plasma environment and its dynamics in the magnetosphere. Their existence strongly reflects the interaction between the solar wind, magnetosphere, and ionosphere. On the other hand, the IHFAC, which we focus on in this paper, is a current system flowing along the magnetic field line that is excited by reflecting the asymmetry of ionospheric current system between the northern and southern hemispheres. At low and middle latitudes, electromotive forces generating ionospheric currents act as a dynamo caused by the atmospheric wind and the penetrated electric field from the polar to equatorial ionosphere which involves a variety of ranges of spatiotemporal phenomena. Therefore, a close examination of IHFACs is very important for understanding the magnetosphere-ionosphere coupling system excited by the solar wind, the global atmospheric motion driven by sunlight, and the electromagnetic environment of the global near-earth system resulting from their coupling.
Although most of the earlier studies have confirmed the existence of IHFACs, few explanations have fully proposed for the characteristics of dusk-side IHFAC based on the ground-based observational data in terms of the time-series analysis. Since IHFACs are excited along the magnetic field lines at middle and low latitudes, their development is known to be remarkably apparent in the east-west component of the ground magnetic field data. This study aims to understand the long-term variation of the IHFAC using the east-west component of the lowlatitude geomagnetic field data continuously observed for two solar cycles, and to develop the leading edge of systematic IHFAC research.
The main purpose of this paper is to investigate the relationship between noon and dusk of IHFAC polarity for individual days, by using long-term ground-based equatorial magnetometer data from 1998 to 2020. The present study focuses on analyzing the occurrence rate of the dusk-side IHFAC of the Fukushima model type. Namely, we calculate the ratio of the number of days having the same polarity between the noon and dusk sectors to the total number of observational days. In addition, we examine the time-series data of IHFACs for monitoring of the solar cycle dependence of IHFACs. The analysis of the 22-year long-term data enables us to investigate the solar cycle dependence of IHFACs (Fujimoto et al. 2016), which is the secondary purpose of this work. In the following, the section "Observation data and analysis method" briefly describes the observational data, their sources, and data analysis method to determine the variation of IHFACs (∆D magnetic field variations). The section "Results and discussion" demonstrates and discusses characteristics of the long-term ∆D variations, the occurrence rate of the dusk-side IHFACs of the Fukushima model type, and climatology of IHFACs polarity in terms of solar cycle dependence of IHFACs. The section "Conclusion" summarizes the present study.

Observation data and analysis method
We used long-term MAGnetic Data Acquisition System/ Circum-pan Pacific Magnetometer Network Data (MAG-DAS/CPMN) magnetometer 1-h time resolution data at Davao station (Geographical latitude 7˚N, Geographical longitude 124.5˚E, Geomagnetic latitude -2.22˚N, Geomagnetic longitude 197.9˚E, Dip latitude -0.24˚) operated by International Center for Space Weather Science and Education (ICSWSE), Kyushu University (Yumoto et al. 1996(Yumoto et al. , 2001(Yumoto et al. , 2006(Yumoto et al. , 2007, from 1998 August to 2020 July. The D-component (the east-west component) of magnetic field was analyzed for investigating the characteristics of IHFACs, since the northward and southward IHFACs induce the westward and eastward magnetic field variations on the ground, respectively. Since the geomagnetic field at Davao is approximately horizontal on the magnetic equator (Davao station is located near at the magnetic equator), we assume that the D-component is due to the magnetic field variation caused by IHFACs. In this paper, "northward IHFACs" is referred as the current flowing from the southern hemisphere to northern hemisphere. The westward ∆D magnetic field corresponds to the northward IHFACs, whereas the eastward ∆D magnetic field corresponds to the southward IHFACs.
To derive the IHFACs effect from the D-component variation, we first subtracted the base line calculated by using the midnight averaged values from each hour data point. Since the source of the daytime ionospheric wind dynamo is the solar daily radiations, the magnetic effect of this dynamo is generally negligible during the nighttime (Yamazaki and Maute 2017). Fambitakoye and Mayaud (1976) determined the base level by interpolating linearly between two midnights neighboring the day considered. This derivation manner was adapted in this paper.
The daily variation of the D-component, ∆D, was calculated with respect to the midnight-to-midnight baseline (daily baseline). Midnight data points for the daily baseline were selected as hourly values of previous day 23:00 LT (D p23h , the subscript "p" indicates "previous"), previous midnight 00:00 LT (D p00h ), target day 01:00 LT (D 00h ), target day 23:00 LT (D 23h ), next midnight 00:00 LT (D n00h , the subscript "n" indicates "next") and next day 01:00 LT (D n01h ). The baseline was calculated as a linear function which fits for these six data points near midnight. The linear function is represented as the following Eq. (1) and slope a and y-intercept b are constants to be determined.
Domain values for x (D p23h , D p00h , D 00h , D 23h , D n00h , and D n01h ) were determined as (-1, 0, 1, 23, 24, and 25) according to the following function (2) so that the distance between x values of any adjacent hours becomes one and the distance between any x values gives the number of hours between those data points: The linear function constants a and b were calculated using the least squares polynomial fit technique so that the calculated linear function can be as representative as possible of all data points. The projected data points of the daily base line for any given hour were calculated according to the following Eq. (3): ∆D represents the daily variation of the D-component relative to the daily base line, calculating for each day as below: F10.7 cm Radio Flux index (F10.7) and Dst index were used to monitor the solar cycle activity and the magnetospheric disturbances, respectively. Daily F10.7 and Dst (2) index data were obtained from SPDF OMNIWeb (King and Papitashvili 2005; SPDF-OMNIWeb Service 2020). The original source of Dst is WDC geomagnetism Kyoto. The analysis period from 1998 to 2020 includes most of the latest two solar cycles SC23 (1996SC23 ( -2008 and SC24 (2008-2019, NASA officially announced that SC25 has begun in 2019 December [NASA 2020]).

Characteristics of IHFACs in the dusk sector
We examined ∆D variations to demonstrate the transition of IHFAC polarity, especially the dusk-side IHFAC polarity change against the noontime IHFAC polarity. Figure 1 shows examples of ∆D variations as a function of local time (LT) and day of year, in 2010, 2013, 2015 and 2016. These years were selected to show the following typical patterns in different solar activity periods in SC24. The directions of IHFACs among both hemispheres reversed twice around the months of February and October, but not exactly during the March and September equinoxes. During the March-September months, the eastward ∆D (= southward IHFACs) and the westward ∆D (= northward IHFACs) were observed in the dawn (06-08 h, LT) and noon (11-13 h, LT) sectors, respectively. For these two sectors, the directions of IHFACs during October-February months were opposite to those of March-September months. The dusk (16-18 h, LT) IHFACs around June and December solstices had roughly the same polarity as the noon sector. The reversal of IHFAC polarity between the noon and dusk sectors frequently occurred in all months except for the June and December months. We found this tendency in all years analyzed in this paper .
In order to clarify the occurrence rate of the dusk-side IHFAC of the Fukushima model type, we calculated the ratio of the number of days having the same polarity between the noon and dusk sectors to the total number of observational days, for each month in each year bin (Fig. 2a). As shown in Fig. 2a, high occurrences of the dusk-side IHFAC of the Fukushima model type appeared around June and December solstices and the occurrence rate was low during the March-April and September-November months (Fig. 2b). The results indicate that there is a seasonal dependence of the dusk-side IHFAC of the Fukushima model type.
For the dusk-side IHFAC polarity, the earlier studies based on in situ and ground-based observations have shown the consistency of the Fukushima model (Olsen 1997;Yamashita and Iyemori 2002), whereas some recent in situ and ground-based observational studies claimed the opposite direction of dusk-side IHFACs (Lühr et al. 2015(Lühr et al. , 2019Fathy et al. 2019;Park et al. 2011Park et al. , 2020aPark et al. , 2020bShinbori et al. 2017). These previous studies have analyzed the monthly, seasonal (combining 2-3 months) or several-year average values of IHFACs in each LT sector. These approaches, however, have not examined the relationship between the noon and dusk IHFAC polarity for individual days. Consequently, although many studies have been conducted using in situ and ground-based observations, there is little agreement on the dusk-side IHFAC polarity as suggested by Fukushima (1979Fukushima ( , 1994. The occurrence rate of the dusk-side IHFAC of the Fukushima model type is examined by this work for the first time, and our findings reveal the seasonal dependence of occurrence distributions of the dusk-side IHFAC polarity over the long term. Figure 3 shows solar cycle dependences of ∆D variations. Figure 3a is the sample plot of monthly averaged of ∆D variations in the dawn sector in 2019. The dawnside IHFAC polarity flows northbound from November to April and the southward IHFAC occurs from May to October. Figure 3b shows the monthly average value of ∆D for each LT sector and solar cycle activity phase, indicating remarkable asymmetry in seasonal IHFACs variation in both the dawn and noon sectors. Absolute monthly averaged ∆D values in both the dawn and noon sectors were stronger around the September equinox, |∆D|~ 30 nT, and weaker around the December solstice, |∆D|~ 10 nT. The dawn-side and noontime ∆D variations exhibited a remarkable seasonal asymmetry. On the other hand, the dusk-side ∆D variations showed symmetry throughout year, with |∆D|~ 5 nT. Our result of monthly averaged ∆D variations at Davao agrees with the result at Guam presented by Shinbori et al. (2017), but there are weak seasonal asymmetries in dawn and noon sectors at Guam. The asymmetry pattern of ∆D, however, is clearly seen in our result. The magnetic field fluctuations caused by IHFACs could be detected effectively in ∆D variations at Davao, because the Davao station is located nearer to the magnetic equator than the Guam station. The IHFAC model proposed in Fukushima [1979] originates from the summer-winter hemisphere asymmetry of the Sq current system caused by the neutral wind dynamo. Yamazaki et al. (2011) showed that seasonal variations of the equivalent Sq current intensity reproduced from the 210 MM magnetometer chain exhibit north-south asymmetry. The equivalent Sq current in the northern summer has a stronger intensity than that in the southern summer. Since the source of the IHFACs is the north-south asymmetry of Sq current, the seasonal north-south asymmetry of Sq current could cause the asymmetry in seasonal monthly averaged ∆D variations.
Our results show that the dusk-side IHFAC polarity during June and December solstices agreed with the polarity proposed by the Fukushima model, whereas the dusk-side IHFAC polarity disagreed during equinoxes. What is a cause of the dusk-side IHFAC flowing the opposite direction of noontime IHFAC? The period of the dusk-side IHFAC with the opposite polarity of the Fukushima model corresponds to the transition phase when the direction of IHFACs reverses among both hemispheres. Figure 3b presents additional information regarding the month of intersection where IHFACs reversed. For the dusk sector, the intersection months were April and September. The noontime intersections were March and November. It is clear that the noontime and dusk-side IHFACs reversed at the different months. There are known seasonal and solar cycle variabilities of the high-latitude current system, and parts of this system are also highly dependent on external driving such as solar wind conditions (Green et al. 2009;Ridley 2007;Edwards et al. 2017;Christiansen et al. 2002). Considering the above-mentioned different reversal month of the IHFAC polarity, we speculate that the appearance pattern of dusk-side IHFAC results from the present analysis could originate from the superposition of the global current system extending from the polar region. In order to investigate this scenario, it is necessary to proceed with comparative analysis between the long-term development of the global current system from the polar to Integral ∆D [ ∆D [nT] Fig. 3 Solar cycle dependence of ∆D variations. a Illustration of the derivation method of peak-to-peak amplitude and integral intensity of ∆D variations; b monthly average of ∆D variations plotted by the separated groups in terms of solar cycle activity: the ascending-maximum phase (1998)(1999)(2000)(2001)(2002)(2003)(2011)(2012)(2013)(2014)(2015) as shown with solid lines and declining-minimum phase (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2016)(2017)(2018)(2019) as indicated by dashed lines, the error bar shows standard deviations for each data point. c Peak-to-peak amplitude of yearly ∆D variations; d yearly integral intensity of ∆D variations, which is provided by the summation of hatched areas in panel (a). Through (b-d), the blue, red and green colors represent the dawn, noon and dusk sectors, respectively. The data points for 2008, 2009, and 2014 were removed due to the lack of data per year middle and low latitudes and the IHFACs excitation pattern. A further analysis of this point should be conducted.

Solar cycle dependence of IHFACs
In this subsection, we focus on the climatology of IHFACs polarity variation. Thus, the daily IHFAC polarity variations are ignored, and the analysis and discussion are limited to seasonal and solar cycle variation. We applied the 61-day centered moving average to ∆D variations with respect to each local time, in order to monitor the long-term transition of IHFAC polarity. The reason why the smoothing value (= 61 days) was selected is because we wanted to disregard variations with frequencies higher than two months. Figure 4 shows the long-term transition of IHFAC polarity from 1999 to 2019 with F10.7 and Dst index. We found that annual cyclic variations of ∆D in the dawn and noon sectors and multiple periodicities in the dusk sector. Figure 4b-d   Fig. 4 Long-term ∆D variation with F10.7 and Dst index, from 1999 to 2019. a 61-day centered moving-averaged ∆D respect to each local time (LT), gray shade patch indicates the lack of data; b 61-day centered moving-averaged ∆ D in dawn sector (06-08 h, LT); c 61-day centered moving-averaged ∆ D in noon sector (11-13 h, LT); d 61-day centered moving-averaged ∆ D in dusk sector (16-18 h, LT); e 61-day centered moving-averaged daily F10.7 index; f daily Dst index and 61-day centered moving-averaged daily Dst index are plotted with gray and black solid lines, respectively shows the time-series data of ∆D variations in each LT sector (dawn, noon, and dusk). We can see that the relationship of ∆D variations between the dawn and noon sectors is out-of-phase ~ 180˚. As shown in Fig. 4d, the dusk-side IHFAC polarity has an annual variation during the declining phase of the solar cycle and semiannual or terannual variations during the solar cycle's ascendingmaximum phase.
As shown in Fig. 3b, the monthly averaged ∆D variations were classified into two groups in terms of solar cycle activity: the ascending-maximum phase (1998)(1999)(2000)(2001)(2002)(2003)(2011)(2012)(2013)(2014)(2015) and declining-minimum phase (2004)(2005)(2006)(2007)(2008)(2009)(2010)(2016)(2017)(2018)(2019). For the dawn and noon sectors, the absolute monthly averaged ∆D values were larger during the ascending-maximum phase than those during the declining-minimum phase. This result is consistent with the result from the ground-based equatorial magnetometer observations by Shinbori et al. (2017). The features of the dawn and noon sectors in Fig. 3b show an asymmetry that tracks with season, regardless of the solar cycle. The yearly integral intensity of ∆D was calculated to examine the total amount of IHFACs in regards to solar cycle activity (Fig. 3d). The yearly integral intensity of ∆D was derived by the summation of ∆D (hatched areas in Fig. 3a) for each year.
We derived the peak-to-peak amplitude of yearly ∆D in each LT sectors to discuss the solar cycle dependence of yearly IHFACs (Fig. 3c). The peak-to-peak amplitude of yearly ∆D in each sector was defined by the difference between the maximum and minimum of monthly averages of the daily averaged IHFAC values per each year (the derivation method, Fig. 3a). As shown in Fig. 3c, the peak-to-peak amplitude values of yearly ∆D were roughly 2-3 times larger in the dawn and noon sectors than that in the dusk sector. The peak-to-peak amplitude values decreased with increasing year in the dawn and noon sectors, whereas the trend of peak-to-peak amplitude in the dusk sector increased slightly over the long term. This is the opposite trend of the yearly peak-topeak amplitude in dawn-side and noontime. The variation ratio of peak-to-peak amplitude IHFAC for each LT sector was estimated according to the following equation, (∆D 2019 -∆D 1998 )/∆D 1998 . For the dawn and noon sectors, the yearly ∆D amplitude in the later solar cycle SC24 decreased by about 35% in comparison with the earlier solar cycle. In contrast, the yearly dusk-side ∆D amplitude increased by about 200%. Along with the slight overall increase from 1998 to 2019, we see larger yearly amplitudes of dusk IHFACs during the declining phases of both solar cycles (2003-2007, and 2016-2019) than during other phases. Concerning the effect from the magnetosphere, for example, the dusk-side ∆D could be disturbed by the partial ring current caused by geomagnetic storms. Especially, the declining phase has been shown to produce different kinds of solar wind driving conditions than what happens during solar maximum (e.g., Richardson et al. 2001;Echer et al. 2008;Webb 1995), which can result in different behaviors in the magnetospheric storm-time response, such as the duration and intensity of the storm (e.g., Borovsky and Denton 2006).
We found that the distribution of occurrence rate was separated into three groups : 1999-2001, 2002-2018, and 2019 (Fig. 2a, Table 1). The observational year and month of the dusk-side IHFAC types (which are consistent with the Fukushima model or not) are summarized in Table 1. The periods of 1999-2001 and 2012 correspond to the maximum phase of solar cycle in SC23 and SC24, respectively. In these years, the dusk-side IHFAC (bins in blue color, as shown in Fig. 2a), which disagrees with the Fukushima model, had high occurrence rates around not only equinoxes but also around the June solstice. The opposite pattern was found in the year of 2019, which corresponds to the minimum phase of solar cycle SC24. The solar activity in SC24 period was lower than in SC23 (McComas et al. 2013). The F10.7 Flux value during solar maximum phases was larger in SC23 than that in SC24, which may be related to the longer duration of the lower occurrence rate of the dusk-side IHFAC polarity of the Fukushima model type in SC23. The dusk-side IHFAC variations were more related with F10.7 variability, in comparison with the dawn and noon sectors (Fig. 4d-e). Especially, the dusk-side ∆D increased during the declining phase of solar cycle (Figs. 3d and 4b).
As mentioned above in the introduction, other previous studies involving the solar cycle discussion have used the period less than solar cycle period (11 years) (e.g., Park et al. 2011;Owolabi et al. 2018;Lühr et al. 2019). In addition, they have not examined the time-series data of IHFACs in terms of the long-term variation analysis. In the present paper, we used 22-year long-term ∆D data and newly showed the above-mentioned findings as shown in Fig. 4a. In order to better determine the generation mechanism of dusk-side IHFACs from the observational data, we need to investigate the relationship between ∆D and various proxies of magnetospheric disturbances such as Dst and SymH indexes. This point should be done in future work.

Conclusion
We investigated the occurrence rate of the dusk-side IHFAC polarity as suggested by Fukushima's model and examined the solar cycle dependence of IHFACs by using the equatorial ∆D variations at Davao with 1-h time resolution data from August 1998 to July 2020. We found that there were seasonal dependences of the appearance of the dusk-side IHFAC of the Fukushima model type (the same IHFAC polarity between the noon and dusk sectors). The polarities of dawn-side and noontime IHFACs agreed with the Fukushima model irrespective of season: the IHFACs during solstices flow from the summer hemisphere to the winter hemisphere in the dawn sector and from the winter hemisphere to the summer hemisphere around the noon sector ( Figs. 1 and 2). The occurrence rate of the dusk-side IHFAC of the Fukushima model type was high around June and December solstices. In contrast, the dusk-side IHFACs, which are inconsistent with the Fukushima model type, dominated from the September to November months and occurred at slightly rates around the March equinox (Fig. 2, Table 1). The remarkable solar cycle dependence of IHFACs was exhibited in the dusk sector. Especially, large IHFACs flowed in the dusk sector during the declining phase of the solar cycle. It is also clear that the noontime and dusk-side IHFACs polarities reversed at the different intersection months. Seasonal variations of ∆D showed asymmetry in both dawn and noon sectors. In contrast, the seasonal duskside ∆D variations were symmetric. In order to conclude the generation and modulation mechanisms of the duskside IHFAC, we need to conduct further comparative analysis between the long-term development and dynamics of global current systems from polar to middle and low latitudes and the IHFAC excitation pattern.