1. Introduction
Thermospheric neutral winds are critical in the coupling between the ionosphere and the thermosphere [
1,
2,
3,
4,
5,
6,
7]. For instance, in a number of studies, it has been shown that the plasma can be moved along geomagnetic field lines [
1,
3,
4,
6], and in the literature, it is also shown that the ionospheric electric field and currents were generated by the collision between ions and neutrals [
2]. For instance, the equatorial electrojet could be driven westward/eastward by the eastward winds at Hall/Pedersen altitudes, in association with the collision between ions and neutrals.
During geomagnetic quiet time, the thermospheric neutral winds were shown to have significant local time, longitudinal, hemispheric, seasonal, and solar activity dependences [
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19]. Meridional winds are generally equatorward/poleward at nighttime/daytime, and the zonal winds were generally westward/eastward before/after 14 local time (LT) [
4,
8]. Using the CHAllenging Minisatellite Payload (CHAMP) observations and Thermosphere–Ionosphere–Electrodynamic General Circulation Model (TIEGCM) simulations, Zhang et al. [
4] found that the equatorial zonal wind jet at the dip equator blew eastward at 14–06 magnetic local time (MLT) and westward at 06–14 MLT. This windjet has also been studied by Liu et al. [
9], Knodo et al. [
10], and Miyoshi et al. [
11]. They found that the quiet-time windjet is driven by the ion drag, which is related to the electron density and the relative motion between ions and neutrals. Moreover, the equatorial wind jet at 20 MLT increased with the solar activity by approximately 110 m/s and 130 m/s at low and high solar activity. Over the past few decades, a number of researchers have paid attention to the longitudinal pattern of thermospheric winds, e.g., [
12,
13,
14,
15,
16,
17,
18]. In their works, the significant wave structures of thermospheric winds have been found. For instance, based on the Fabry–Pérot interferometer (FPI) measurements, Wu et al. [
12] revealed that zonal winds may behave differently at different longitudes. They showed that zonal winds observed in Boulder turned westward earlier and had a larger diurnal variation than the zonal winds seen at the Chinese stations during geomagnetic quiet conditions. A strong westward wind at fixed longitudes and 50°–60° GLat was obtained in the CHAMP observations [
15]. As reported by Häusler and Lühr [
14], Häusler et al. [
16], Wang et al. [
17], and Wang and Zhang [
18], thermospheric winds had an obvious wave-4 longitudinal structure, caused by the strong wave number 3 (DE3) nonmigrating tidal component. The seasonal and hemispheric dependences of the longitudinal pattern of zonal winds were explored in this study by Zhang et al. [
13]. They found that the longitudinal distributions of zonal winds in the northern hemisphere were almost the opposite of those in the southern hemisphere. Moreover, the longitudinal patterns during the June solstice were significantly different from those in other seasons. By imposing a poor dipole geomagnetic configuration in the TIEGCM, they found that the geomagnetic field configuration was the main cause of the local time, hemispheric asymmetry, and seasonal changes in zonal winds at middle and low latitudes.
In the literature, the behaviors of neutral winds during the geomagnetically disturbed time have also been investigated, for instance, geomagnetic storms, substorms, and subauroral polarization streams [
19,
20,
21,
22,
23,
24,
25,
26,
27,
28]. Dungey [
24] suggested in their study that the interaction between the interplanetary magnetic field and geomagnetic field could lead to energy transfer from the solar wind to the Earth’s upper atmosphere. During geomagnetic storm time, the interaction between the southward IMF Bz and geomagnetic field, the energy deposition could lead to the enhancement of Joule heating and neutral temperature and cause large-scale and medium-scale traveling atmospheric disturbances (TADs) in meridional winds [
20]. Based on the zonal winds observed by CHAMP, Ritter et al. [
24] found that substorm-related disturbed winds increased in the westward with a speed of roughly 50 m/s at midlatitudes around midnight. The universal time (UT) and local time dependences of the substorm effects on zonal winds at high latitudes were explored in this study by Wang et al. [
23]. The disturbed winds were poleward and westward in the dusk sector and equatorward and westward at nighttime. The daytime/nighttime perturbation was related to the ion drag/both the Bz and hemispheric power input. Owing to the low background plasma, the disturbed winds responded somewhat later at nighttime than during the daytime. When the geomagnetic pole moved toward the dayside/nightside, stronger/weaker disturbed winds could be generated. Subauroral polarization streams (SAPS) comprised the strong geomagnetic westward plasma flow at the subauroral latitudes [
21,
26,
27,
28]. SAPS were located in a latitudinally narrow region from dusk to early morning sectors. SAPS were driven by a strong poleward electric field during geomagnetic disturbed and quiet periods and had a speed greater than 500 m/s. Due to the collision between ions and neutrals, the neutrals were moved westward with ions. At subauroral latitudes, the enhanced frictional heating resulted in the upwelling of molecular-rich air from lower altitudes to higher altitudes. Away from the SAPS region, neutral wind convergent flow produced a downwelling of atomic oxygen-rich air [
29]. SAPS-driven nighttime geographic poleward winds at 30°–50° geomagnetic latitudes (MLat) showed obvious UT variations [
27]. Due to the misalignment between geomagnetic and geographic coordinate systems, the strong geomagnetic westward ion drag could be separated into two components: geographic poleward/equatorward and geographic westward. The strong geographic poleward ion drag could drive a poleward wind at nighttime and mid-latitudes. Therefore, the poleward wind changes due to SAPS were stronger at 06 and 18 UT and weaker at 00 and 12 UT. By imposing an empirical SAPS model into the global ionosphere thermosphere model, Wang et al. [
28] found that the SAPS-driven disturbed winds showed a close correlation with the solar zenith angle χ according to cos0.5 χ; hence, with more sunlight, stronger westward winds were generated. The strongest/weakest disturbed winds occurred on 18/04 and 04/16 UT in the northern and southern hemispheres, respectively.
The authors of previous studies have disclosed that temporal oscillations in the IMF Bz with Alfvén waves are frequent in solar wind [
20,
30,
31]. Liu et al. [
30] and Zhang et al. [
20] found that the coupled Magnetosphere–Ionosphere–Thermosphere (MIT) system had the nature of a low-pass filter. That is, with respect to the high-frequency IMF Bz, the coupled MIT system could fully respond to the low-frequency IMF Bz. During the geomagnetic storm on 23–24 April 2023, the temporal oscillations in the IMF Bz were strong. However, it is still unknown as to how thermospheric winds would respond to it. The aim of the present study is to address the potential physical drivers and determine the response of thermospheric horizontal winds at Mohe station. In the literature, it is well known that the neutral winds are controlled by a balance between ion drag, pressure gradient, Coriolis force, centrifugal, and viscosity, e.g., [
32]. The rapid changes of IMF Bz could lead to perturbations in the ionospheric convection, causing disturbances in the neutral winds via several drivers (e.g., pressure gradient and ion drag). However, when they discuss the roles of IMF Bz on neutral winds, the contributions from IMF Bx and By, solar wind speed, and density cannot be completely excluded. In the present work, we performed two cases: one is input with the observed IMF and solar wind density and speed (real case); the other case is specified by a constant Bz with a value of 0 nT, and the other inputs are the same. In comparison with the previous studies, the effects of solar wind density, wind speed, IMF Bx, and By on neutral winds are removed. Therefore, our findings could improve the understanding of thermospheric winds during the temporal oscillations of IMF Bz. Our study is founded on the previous understanding but does not agree in all parts. Thus, new insights into the neutral wind response to IMF Bz oscillations are revealed.
2. Data and Model Description
Swarm satellites have a near-polar orbit with an inclination of 87°, consisting of three identical satellites: Alpha, Bravo, and Charlie (A, B, and C) [
33]. Swarm A and C fly side-by-side at ~450 km, with a 1.4° separation in longitude at the equator. The orbit of Swarm B is ~530 km. The Swarm satellites have an orbital period of ~96 min. The blue lines in
Figure 1 indicate the orbits of Swarm A on 23 April. In the present work, the electron density and neutral density measurements from Swarm A were used to estimate the reliability of the TIEGCM. The GRACE-FO twin satellite mission was launched on 22 May 2018. It has been deployed directly to an initial altitude of approximately 520 km, with a near-polar inclination of 89°. Each day, it passes through the Earth 15.3 times. In
Figure 1, the orbits of GRACE-FO on 23 April have been indicated by the red lines. In the present work, the cross-track winds from GRACE-FO on 23–24 April 2023 have been used to compare with TIEGCM simulations. We used the Fabry–Pérot interferometer (FPI) operated at Mohe station (122.3° geographic longitude (GLon); 53.5° geographic latitude (GLat)) to aid us in understanding the responses of thermospheric winds to the temporal oscillations in the IMF Bz. The location of Mohe station is indicated by the magenta star shown in
Figure 1. The station provides the nighttime wind velocity at around 250 km using the Doppler shift in the airglow in four directions (north, east, south, and west) with an elevation angle of 45°. The FPI observations have a temporal resolution of ~10 min.
The TIEGCM v2.0 is a first principles model of the coupled thermosphere and ionosphere. The drivers include the high-latitude electric field specified by the empirical Heelis model [
34] or Weimer model [
35], solar EUV, and UV spectral fluxes parameterized by the F10.7 index [
36]. In this work, TIEGCM has a horizontal resolution of 2.5° GLat by 2.5° GLon. The vertical resolution is 1/4 scale height, with a bottom/upper boundary of 97/600 km. The lower boundary forcing is specified by either the Global Scale Wave Model (GSWM) [
37,
38] or the derived tides from the Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) and TIDI observations [
12,
13]. In the present work, the 1 min IMF data from OMNI was imposed into the TIEGCM with the high-latitude electric field specified by the Weimer model. Note here that the solar wind data from OMNI are observed by the Advanced Composition Explorer (ACE), with the Time-Shifted to the Nose of the Earth’s Bow Shock. The information on ACE can be found at the URL of
https://science.nasa.gov/mission/ace/ (accessed on 23 April 2023). The high latitude ion convection and auroral particle precipitation in TIEGCM are defined by the empirical Weimer model [
35]. The input parameters for the empirical Weimer model are solar wind density, wind speed, and IMF. The sign of IMF By is changed for the potential pattern in the southern hemisphere. The lower boundary is specified by the migrating and nonmigrating tides from the GSWM model. To reach a diurnally reproducible steady state, the TIEGCM is run for 20 days including the storm event as the last day of the model simulation.
4. Discussion
During the geomagnetic storm, the interaction between southward IMF and the geomagnetic field leads to the occurrence of open field lines [
24]. The charged particles in the solar winds enter the Earth’s upper thermosphere along the open field lines, causing energy deposition and momentum transfer. High-latitude ionospheric convection and Joule heating constitute significant changes. Joule heating could enhance the neutral temperature, and the disturbances in the thermosphere could travel to the middle and low latitudes, causing traveling atmospheric/ionospheric disturbances (TADs/TIDs). The traditional TADs are equatorward. However, unexcepted poleward TADs/TIDs at middle latitudes have been observed by researchers, including Zhang SR et al. [
41]. As described in the literature, the behavior of neutral winds is controlled by a balance between pressure gradient, ion drag, Coriolis force, centrifugal, and viscosity [
4,
11,
32]. The issue is that the potential drivers of neutral wind changes at Mohe station during the strong temporal oscillations in the IMF Bz are still unknown. As we discussed before, the pressure gradient in association with the neutral temperature would be greatly influenced by the IMF Bz oscillations. Furthermore, during IMF Bz oscillations, the prompt penetration electric field (PPEF) and ion drifts at middle and low latitudes are disturbed. Then, the ionospheric plasma density could be greatly perturbed by the disturbed winds, and the ion drifts. The relative motion between ions and neutrals could be also modulated. Therefore, the ion drag effects on neutral winds could be affected during IMF Bz oscillations. The Coriolis force is related to the Coriolis coefficient and the thermospheric winds. At Mohe station, the Coriolis coefficient is constant. During IMF Bz oscillations, the energy deposition at high latitudes causes disturbances in thermospheric winds, driving disturbances in the Coriolis force. In the following paragraphs, the drivers of zonal and meridional wind disturbances will be explored via term analysis in the TIEGCM. The detailed effects from different forces can be separated by the TIEGCM (please refer to Hsu et al. [
32] and the TIEGCM description on the HAO.
To explore the effects of the temporal variations in the IMF Bz on neutral winds, a control case with a constant Bz value of 0 nT is performed. The differences between the neutral winds in the real and control cases are the wind perturbations associated with the strong temporal oscillations of IMF Bz.
Figure 8a depicts the temporal variations of zonal wind responses (ΔUN) to the temporal oscillations in the IMF Bz on 23 April. The equation about ΔUN is shown as follows:
. It can be found that ΔUN is strong westward at 14–16 UT, with an average speed of 150 m/s. The peak westward ΔUN occurs at 15 UT, with a speed of 200 m/s. At the following UTs, the westward ΔUN is rapidly reduced and reverses eastward. The polarity transition is located at 17 UT. The peak of eastward ΔUN at 17–19 UT is 130 m/s at 19 UT. Maybe over 19 UT, the maximum of zonal wind changes might exist but are still unknown. An interesting phenomenon should be noted herein the enhanced ΔUN eastward gets slowed and even weakened at 18 UT. On 24 April (
Figure 6b), a similar temporal variation in the ΔUN to that on 23 April is found. ΔUN is first enhanced westward from −80 m/s at 14 UT to −150 m/s at 15 UT, then eastward to 10 m/s at 17 UT, thereafter westward to −30 m/s at 18.5 UT, and finally eastward to −20 m/s at 19 UT.
Figure 8c,d provides information on the acceleration perturbations due to the total forcing, pressure gradient, ion drag, and Coriolis force on 23–24 April. The effects of other forces (i.e., centrifugal, horizontal advection, and viscosity) are relatively weaker than the above three forces. Therefore, only the acceleration responses associated with the pressure gradient (ΔZL), ion drag (ΔFD_UN), and Coriolis force (ΔCOR_UN) are displayed herein. All of those three factors are related to the collision between ions and neutrals. The plasma collisional heating is the dominant heating mechanism for the neutrals in the upper thermosphere. Therefore, the pressure gradient might be affected by the collision between ions and neutrals. The ion drag is sourced from the plasma density and the relative motion between ions and neutrals. The Coriolis force on the zonal winds is attributed to the Coriolis coefficient and the meridional winds. The meridional winds could be controlled by the pressure gradient and ion drag. Acceleration represents the ratio of speed changes. The negative (positive) acceleration changes indicate the decrease (increase) in eastward winds or the increase (decrease) in westward winds. Therefore, the temporal structures of acceleration might not be the same as those of the neutral winds.
At 14–15 UT on 23 April, the total acceleration changes are weak eastward at the beginning stage and obviously enhanced westward, with a trough of −1.5 cm/s
2. The effects of the pressure gradient are westward from −0.2 cm/s
2 at 14 UT to −2.0 cm/s
2 at 15 UT. A comparison between the accelerations due to the total forcing and the pressure gradient shows that the pressure gradient is the primary driver. The same result can be obtained at the following UTs on 23 April. During the geomagnetic storm, the pressure gradient is related to the neutral temperature changes [
20,
32]. Thus, as shown in
Figure 8a, the westward ΔUN at 14–15 UT might be dominated by the pressure gradient changes. However, the effects of ion drag and Coriolis force are positive, with an average acceleration of 1 and 0.2 cm/s
2, respectively. The positive acceleration indicates the weakening of the westward winds. Hence, at 14–15 UT, ion drag and Coriolis force prevent the formation of ΔUN. As reported by the authors of previous studies [
15,
32], the ion drag is associated with two factors: one is the plasma density and the other is the relative motion between ions and neutrals. The detailed roles of plasma density and relative motion are not the focus of the present study and have thus not been included herein. Instead, we only show the complete effects of ion drag. The effects of Coriolis force on zonal winds are related to the Coriolis coefficient and meridional winds [
27,
42]. The Coriolis force tends to direct westward with an equatorward wind. At Mohe station, ΔVN is poleward at 14–15 UT (
Figure 8a in the following section); hence, the Coriolis force effects are directly eastward, preventing the formation of ΔUN at 14–15 UT. At the phase of the enhanced ΔUN eastward (15–19 UT,
Figure 6a), the acceleration changes due to the pressure gradient are strongly positive, indicating its key and positive role. The effects of ion drag are weak and perturbed at around zero, indicating its insignificant role. In comparison, the effects of Coriolis force are negative, with an average magnitude of −1.5 cm/s2. Thus, the enhanced ΔUN in eastward is also dominated by the pressure gradient, with negative contributions from ion drag and Coriolis force. On 24 April (
Figure 8b,d), a similar conclusion is obtained. That is, ΔUN is controlled by the pressure gradient, with minor contributions from ion drag and negative effects from Coriolis force.
Figure 9 depicts the temporal variations in meridional wind responses (ΔVN) to the IMF Bz on 23–24 April 2023 and the associated accelerations due to different forces. The equation about ΔVN can be expressed as follows:
. It can be seen that on April 23 (
Figure 9a), ΔVN is generally equatorward and enhanced in equatorward. As shown in
Figure 9c, the acceleration changes due to the total force, pressure gradient, ion drag, and Coriolis force are provided. The acceleration changes due to the total forcing share a large degree of similarity with that due to the pressure gradient. The acceleration changes due to ion drag are generally weak and negative, with an average magnitude close to zero. The acceleration perturbations associated with the pressure gradient are strong. They are negative at 14.5–15.5, 16.5–17.5 and 18.5 UT. The Coriolis force effects are negative at 14–17 UT, and positive at the following UTs. In summary, this equatorward ΔVN on 23 April might be related to the combined roles of pressure gradient and Coriolis force, with minor contributions from ion drag. A similar conclusion is achieved for the ΔVN on 24 April.
As shown in
Figure 2e, the IMF Bz at 14–17 UT has an average magnitude of −7.8 nT. The southward IMF Bz occurs continually and lasts for hours. The interaction between southward Bz and the geomagnetic field can lead to the energy deposition and heat of the high-latitude neutrals. Therefore, the nighttime neutral winds at Mohe station in
Figure 9c could be perturbed greatly by the pressure gradient in associated with the enhanced neutral temperature. In
Figure 2e,f, the great oscillations of IMF Bz occur in the period from 18 UT on 23 April to 06 UT on 24 April. The perturbed neutral winds could have existed not only during the above periods but also at the following time (
Figure 8 and
Figure 9). Because the behaviors of the neutral winds are the accumulation of all forces. Furthermore, the disturbances in the neutral winds have a time delay in comparison with IMF Bz. Because the heated air needs time to travel from the high latitudes to the middle latitudes. Therefore, the nighttime neutral wind responses at Mohe station could be related to the great perturbations of IMF Bz. In
Figure 2e,f, IMF Bz has a period ranging from minutes to hours. However, in
Figure 8a,b and
Figure 9a,b, the responses of zonal and meridional winds to IMF Bz are smooth. This is reasonable and related to large-scale traveling atmospheric disturbances (LSTADs). During storm times, the increased Joule heating at high latitudes could significantly enhance the neutral temperature, causing the upwelling of molecular-rich air. The heated air could extend to the lower latitudes through the dynamic processes, that is, LSTADs. Previous studies have demonstrated that the LSTADs have a period of 0.5–3 h [
20]. Therefore, the residual zonal and meridional winds are smooth.
5. Summary
Using Swarm observed electron density and neutral density, FPI measured thermospheric winds, and TIEGCM simulations, the roles of IMF Bz on the thermospheric winds at Mohe station during the storm on 23–24 April are explored in the present work. During the study, a number of interesting results were derived as follows:
1. The meridional winds are strong/weak equatorward at pre-dawn on 23/24 April. The peak/trough of zonal winds occurs at midnight on 23/24 April.
2. The responses of zonal winds to the IMF Bz are westward and eastward at pre-midnight and pre-dawn, respectively. This finding is primarily attributed to the pressure gradient, with contributions from Coriolis force and ion drag.
3. The meridional wind perturbations are strong equatorward (–200 m/s) on 23 April. However, on 24 April, they are generally poleward and peak at midnight. This finding might be the result of the combined roles of both pressure gradient and Coriolis force.