Contributions of Heating and Forcing to the HighLatitude Lower Thermosphere: Dependence on the Interplanetary Magnetic Field
 Author: Kwak YoungSil, Richmond Arthur, Ahn ByungHo, Cho KyungSuk
 Organization: Kwak YoungSil; Richmond Arthur; Ahn ByungHo; Cho KyungSuk
 Publish: Journal of Astronomy and Space Sciences Volume 27, Issue3, p205~212, 15 Sep 2010

ABSTRACT
To better understand the physical processes that maintain the highlatitude lower thermospheric dynamics, we have identified relative contributions of the momentum forcing and the heating to the highlatitude lower thermospheric winds depending on the interplanetary magnetic field (IMF) and altitude. For this study, we performed a term analysis of the potential vorticity equation for the highlatitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions, with the aid of the National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model (NCARTIEGCM). Difference potential vorticity forcing and heating terms, obtained by subtracting values with zero IMF from those with nonzero IMF, are influenced by the IMF conditions. The difference forcing is more significant for strong IMF B_{y} condition than for strong IMF B_{z} condition. For negative or positive B_{y} conditions, the difference forcings in the polar cap are larger by a factor of about 2 than those in the auroral region. The difference heating is the most significant for negative IMF B_{z} condition, and the difference heatings in the auroral region are larger by a factor of about 1.5 than those in the polar cap region. The magnitudes of the difference forcing and heating decrease rapidly with descending altitudes. It is confirmed that the contribution of the forcing to the highlatitude lower thermospheric dynamics is stronger than the contribution of the heating to it. Especially, it is obvious that the contribution of the forcing to the dynamics is much larger in the polar cap region than in the auroral region and at higher altitude than at lower altitude. It is evident that when B_{z} is negative condition the contribution of the forcing is the lowest and the contribution of the heating is the highest among the different IMF conditions.

KEYWORD
highlatitude lower thermosphere , Interplanetary magnetic field , forcing , heating

1. INTRODUCTION
At high latitudes the thermospheric dynamics are governed by various heat and momentum sources: nonuniform solar heating, atmospheric tides and other waves coming from below, and energy and momentum sources associated with magnetosphereionosphere coupling, such as the magnetospheric convective electric field and the precipitation of energetic electrons and ions into the auroral regions.
Recently several modeling studies have been attempted to understand the physical process that control the highlatitude lower thermospheric dynamics. Kwak & Richmond (2007) studied the momentum forcing balance that are mainly responsible for maintaining the highlatitude lower thermospheric wind system by using the National Center for Atmospheric Research Thermosphere Ionosphere Electrodynamics General Circulation Model (NCARTIEGCM). Kwak et al. (2008b) analyzed the divergence and vorticity of the highlatitude neutral wind field in the lower thermosphere during the southern summertime. They showed that the mean neutral wind pattern in the highlatitude lower thermosphere is dominated by rotational flow rather than by divergence flow. Kwak et al. (2008a) performed a term analysis of vorticity equation that describes the driving forces for the rotational component of the horizontal wind in order to determine key processes that causes strong rotational flow in the highlatitude lower thermospheric winds.
We have extended previous works by Kwak & Richmond (2007) and Kwak et al. (2007, 2008a,b), which have been attempted to better understand the physical process maintaining thermospheric dynamics at high latitudes, and here first perform a “term analysis of the potential vorticity equation” for the highlatitude neutral wind field in the lower thermosphere. These analyses can provide separate insight into the relative strength of the momentum forcing and the heating responsible for driving winds at the highlatitude lower thermosphere. The momentum forcing is associated with the viscous force and drag force from convecting ions. The heating is the net heat including the heat transfer by downward molecular and eddy heat conduction, the absorption of solar ultraviolet (UV) and extreme ultraviolet (EUV) radiation, auroral heating by particles, Joule dissipation of ionospheric currents, release of chemical energy by the atomic oxygen recombination, and radiative CO_{2},NO and O infrared emissions.
It is known that highlatitude ionospheric plasma convection and current depend on the direction of the interplanetary magnetic field (IMF) (e.g., Heppner 1972, Foster et al. 1986, Heppner & Maynard 1987, Ruohoniemi & Greenwald 1996, Weimer 1995, 2001). Thus the highlatitude thermospheric wind which are influenced by interaction with the ionospheric plasma, forcing on the wind and Joule/particle heating, are also strongly modulated by the IMF variation (e.g., McCormac & Smith 1984, McCormac et al. 1985, 1991, Killeen et al. 1985, 1995, Meriwether & Shih 1987, Thayer et al. 1987, Rees & FullerRowell 1989, 1990, Sica et al. 1989, Hernandez et al. 1991, Niciejewski et al. 1992, 1994, Won 1994, Richmond et al. 2003, McHarg et al. 2005, Zhang et al. 2005, Kwak & Richmond 2007, Kwak et al. 2007). That is, one can expect there is an intimate relationship between the IMF variation and the relative strength of momentum forcing / heating responsible for maintaining thermospheric dynamics at high latitudes. In this study, we investigate the relative contributions of the momentum forcing and the heating to the highlatitude lower thermospheric dynamics depending on the IMF during the southern summertime by using NCARTIEGCM.
2. MODEL
The NCARTIEGCM used in this study computes selfconsistently the coupled thermospheric/ ionospheric dynamics, the associated dynamo electric fields and currents, and the electrodynamic feedback on the neutral and plasma motions and thermodynamics (Richmond et al. 1992). The nonlinear primitive equations for momentum, energy, continuity, hydrostatics, current density, and the equations of state are solved for the neutrals and the ions. The model calculates global distributions of neutral gas temperature, wind, mass mixing ratios of the major constituents O_{2}, N_{2}, and O, and of the minor constituents N(^{2}D), N(^{4}S) and NO. The Eulerian model of the ionosphere solves for global distributions of electron and ion temperatures and number densities of O^{+}, O_{2}^{+}, NO^{+}, N_{2}^{+}, and N^{+}. The TIEGCM has a 5˚ latitudebylongitude grid. The vertical coordinate of the model consists of 29 constant pressure levels ranging from approximately 97 km to 500 km in altitude, with a vertical resolution of 2 grid points per scale height. In this study we focus on altitudes up to about 200 km. The external inputs required by the model are the solar EUV and UV fluxes, the auroral particle precipitation, the ionospheric convection pattern at high latitudes, and the upward propagating tides from the middle atmosphere.
We analyze the relative contributions of the momentum forcing and the heating to the highlatitude lower thermospheric dynamics for different IMF directions, for seasonal and solar conditions representative of 1993 January 23. This period is a characteristic of conditions for wind observations analyzed by Richmond et al. (2003) and previous analyses of forcing (Kwak & Richmond 2007, Kwak et al. 2007) and vorticity (Kwak et al. 2008a,b). On this date, the daily F10.7 index was recorded as 102.7 × 10^{22} W/m^{2}/Hz. The total hemispheric power (HP) was used to specify the highlatitude auroral particle precipitation pattern and is related with the Kp index by HP (in gigawatts) = 2.78+9.33×
Kp (Maeda et al. 1989). The meanKp index on January 23, 1993 was 1, although the actual timevarying values for this date were used. The empirical electric potential model by Weimer (2001) was adopted to specify the pattern of the ionospheric convection. The prescriptions of the upward propagating diurnal and semi diurnal tides are taken from the Global Scale Wave Model (GSWM) (Hagan & Forbes 2002).To investigate the response of the thermosphere to varying IMF, five TIEGCM simulations for IMF (B_{y},B_{z}) values of (3.2, 0.0), (+3.2, 0.0), (0.0, 2.0), (0.0, +2.0), and (0.0, 0.0) nT are made. The magnitudes of the nonzero reference values of B_{y} and B_{z} , that is, 3.2 nT and 2.0 nT, are their respective rootmeansquare values for the temporally smoothed data set analyzed by Richmond et al. (2003). A 2min time step is used for the entire simulation, and linear interpolations of input values are made at a given time step. The TIEGCM history is recorded hourly for the period that we studied.
3. POTENTIAL VORTICITY EQUATION
A largescale motion in lower atmosphere below the mesosphere has a synoptic time scale and thus is nearly in an adiabatic and nonfriction processes. Therefore, the contour of the potential vorticity on the isentrope surface is nearly material contour and the potential vorticity is conserved along with the isentrope surface. Thus, an analysis of the potential vorticity on the isentrope surface is a useful tool for measuring the flow of air.
Unlike below the mesosphere, since the largescale motion in the thermosphere is in a diabatic process having sources and sinks, the contour of the potential vorticity on the isentrope surface in the thermosphere is no longer material contour. Moreover, the potential vorticity is not conserved along with the isentrope because the molecular diffusion in the thermosphere is so strong and the time scale of the potential vorticity is so short. In the thermosphere, thus, the potential vorticity itself has no meaning for study of the thermospheric dynamics.
A term analysis of the potential vorticity equation, however, provides a quantitative insight into the contribution of the momentum forcing and the heating to the highlatitude lower thermospheric dynamics, because this equation describes the terms that derive the change of the potential vorticity. The potential vorticity equation can be derived from vorticity equation and is given as (Pedlosky 1979)
where Π is the Ertel’s potential vorticity on the isobaric surface. ω_{a} is the vertical component of the absolute vorticity, ρ is atmospheric mass density. λ is some scale fluid property and could be the potential temperature (θ) or ms/R , where
m is the atmospheric molecular mass, s is the entropy, and R is the universal gas constant. ψ is defined as Dλ/D_{t} and is associated the net heat. In Eq. (1), P is pressure, and F_{V} and F_{I} are the viscous and ion drag forces, respectively. The time rate of change of the potential vorticity is determined by the following terms on the right hand side of the Eq. (1): (a) heating term, (b) momentum forcing term associated with the viscous and ion drag forces, and (c) term associated with the pressure gradient.The individual terms in units of m^{2}kg^{1}Ks^{2} at all model grid points for each hour of universal time are obtained by using the potential vorticity Eq. (1) with outputs of the model, where vertical and horizontal finite differencings are performed to calculate terms at a given altitude for each universal time. We then carry out an averaging for 24 hours of each potential vorticity forcing terms. The resultant forcing terms with respect to different altitudes can be then mapped out over magnetic latitude and magnetic local time.
4. RESULTS
In order to examine closely how the potential vorticity forcing and heating terms are influenced by the IMF, we consider the difference potential vorticity forcing and heating terms, obtained by subtracting values with zero IMF from those with nonzero IMF conditions.
Fig. 1 shows the calculated distribution of difference forcing term, which is associated with iondrag and viscosity, contributing to the change of the potential vorticity at 195, 177, 160, 146, 134, 125, 117, 111, and 106 km altitudes over the southern hemisphere in the poleward latitude 50^{˚} for IMF (B_{y}, B_{z}) values of (3.2, 0.0), (3.2, 0.0), (0.0, 2.0), and (0.0, 2.0) nT. These projections are as if one were looking up on the thermosphere from below. Solid and dotted lines represent positive and negative values, respectively. The minimum and maximum values are indicated at the lower corner of each plot. The distributions of the forcing terms at 195, 177, and 160 km altitudes are shown with contour interval of 50 m^{2}kg^{1}Ks^{2}. They are shown with contour interval of 10 m^{2}kg^{1}Ks^{2} at 146, 134, 125 km altitudes, and with contour interval of 0.2 m^{2}kg^{1}Ks^{2} at 117, 111, and 106 km altitudes.
For negative IMF B_{y}, at 195 km strong positive difference forcing occurs in the poleward latitude 75^{˚} (near polar region) with the maximum value of 1038 m^{2}kg^{1}Ks^{2}. There is also negative difference forcing region at 65^{˚}∼75^{˚} latitude (near auroral region) with the minimum value
of 438 m^{2}kg^{1}Ks^{2}. For positive IMF By, the distribution is similar but the sign is opposite with that for negative B_{y} at roughly the same place. For negative IMF B_{z} condition, at 195 km there are positive difference forcing regions on the dusk with the maximum of 434 m^{2}kg^{1}Ks^{2} and negative difference forcing regions on the dawn sides with the minimum of 295 m^{2}kg^{1}Ks^{2}, respectively. When positive IMF B_{z} condition, the sign is opposite with that for negative B_{z} at roughly the same place with the minimum of 885 m^{2}kg^{1}Ks^{2} and the maximum of 311 m^{2}kg^{1}Ks^{2}. It is evident that the difference forcing is more significant for strong IMF B_{y} condition, and the difference forcings in the polar cap are larger by a factor of about 2 than those in the auroral region. It is obvious that the magnitude of the difference forcing decreases rapidly with descending altitudes.
Fig. 2 shows the calculated distribution of difference heating term contributing to the change of the potential vorticity at 195, 177, 160, 146, 134, 125, 117, 111, and 106 km altitudes over the southern hemisphere for IMF (B_{y}, B_{z}) values of (3.2, 0.0), (3.2, 0.0), (0.0, 2.0), and (0.0, 2.0) nT. The distributions of the forcing terms at 195, 177, and 160 km altitudes are shown with contour interval of 5 m^{2}kg^{1}Ks^{2}. They are shown with contour interval of 1 m^{2}kg^{1}Ks^{2} at 146, 134, 125 km altitudes, and with contour interval of 0.1 m^{2}kg^{1}Ks^{2} at 117, 111, and 106 km altitudes.
For negative IMF B_{y}, at 195 km positive difference heating occurs in the poleward latitude 75^{˚} with the maximum value of 19.95 m^{2}kg^{1}Ks^{2}. There is negative difference heating region in the morning auroral region around 65^{˚}∼75^{˚} latitude with the minimum value of 26.32 m^{2}kg^{1}Ks^{2}. For positive IMF B_{y}, the distribution is similar but the sign is opposite with that for negative B_{y} at roughly the same place. For negative IMF B_{z} condition, at 195 km there are positive difference heating regions on the dusk including the polar cap with the maximum of 24.37 m^{2}kg^{1}Ks^{2} and negative difference heating regions in the auroral region on the dawn sides with the minimum of 52.63 m^{2}kg^{1}Ks^{2}, respectively. When positive IMF B_{z} condition, the sign is opposite with that for negative B_{z} at roughly the same place with the minimum of 14.99 m^{2}kg^{1}Ks^{2} and the maximum of 22.91 m^{2}kg^{1}Ks^{2}. It is evident that the difference heating is the most significant for negative IMF B_{z} condition, and the difference heatings in the auroral region are larger by a factor of about 1.5 than those in the polar cap region. It is obvious that the magnitude of the difference heating decreases rapidly with descending altitudes. Although not shown here, term (c) in the potential vorticity Eq. (1) is much smaller than the forcing term (a) or the heating term (c). Indeed, the absolute values of difference term (c) associated with the pressure gradient are less than 5 m^{2}kg^{1}Ks^{2} at 195 km for all IMF conditions.
To see relative contributions of the forcing and heating to the highlatitude lower thermospheric dynamics for different IMF conditions and with altitudes, we examine a forcingtoheating ratio. Fig. 3 shows this ratio with altitudes for different IMF conditions. Solidblack and dottedgrey bars represent forcingtoheating ratio for positive (or maximum) and negative (or minimum) values from Figs. 1 and 2, respectively. For negative B_{y} condition shown in Fig. 3a, the solidblack and dottedgrey bars represent forcingtoheating ratio for in the polar cap and in the auroral region, respectively. For positive B_{y} condition (Fig. 3b), the solidblack and dottedgrey bars represent forcingtoheating ratio for in the auroral region and in the polar cap region, respectively. One can see that the ratio is more than 2 with maximum 65.6 at 195 km for positive B_{y} condition, indicating that the contribution of the forcing to the highlatitude lower thermospheric dynamics is stronger than the contribution of the heating to it. The forcingtoheating ratio is much larger by a factor of 1.2∼5.5 in the polar cap region than in the auroral region above 125 km altitude. It is obvious that the magnitude of the ratio decreases rapidly with descending altitudes. Especially, the ratio in the polar cap region decreases more rapidly with descending altitudes than in the auroral region. For negative B_{z} condition shown in Fig. 3c, the solidblack and dottedgrey bars represent forcingtoheating ratios in the evening auroral region including polar cap and in the morning auroral region, respectively. One can see that the ratio is more than 1.5 with maximum 21.7 at 177 km except at 106 km in the morning auroral region, indicating that the contribution of the forcing to the highlatitude lower thermospheric dynamics is stronger than the contribution of the heating to it. It is shown that the ratio in the evening auroral region including the polar cap is larger by a factor of 1.2∼10 than in the morning auroral region and that the ratio at all altitudes except for 111 km is smaller by a factor of about 0.3 than for other IMF conditions. For positive B_{z} condition (Fig. 3d), the solidblack and dottedgrey bars represent forcingtoheating ratios in the morning auroral region and in the evening auroral region including polar cap, respectively. It is shown that the ratio in the evening auroral region including the polar cap is larger by a factor of 1.5∼5 than in the morning auroral region. One can see again the forcingtoheating ratio decreases rapidly with descending altitudes.
From Fig. 3, it is evident that the contribution of the
forcing to the highlatitude lower thermospheric dynamics is stronger than the contribution of the heating to it. Especially, it is obvious that the contribution of the forcing is much larger at higher altitude than at lower altitude and in the polar cap region than in the auroral region for negative and positive IMF B_{y} conditions. The relative contribution of the forcing in the auroral region for negative B_{z} condition is the lowest among the different IMF conditions.
5. CONCLUSIONS
In this paper, we first present the relative contributions of the momentum forcing and the heating to the highlatitude lower thermospheric dynamics depending on the IMF and altitude. For this study we attempted a “term analysis of the potential vorticity equation” for the highlatitude neutral wind field in the lower thermosphere during the southern summertime for different IMF conditions by using NCARTIEGCM.
By considering the difference potential vorticity forcing and heating terms, obtained by subtracting values with zero IMF from those with nonzero IMF, we can examine closely how the potential vorticity forcing and heating terms are influenced by the IMF. The difference forcing is more significant for strong IMF B_{y} condition than for strong IMF B_{z} condition. For negative or positive B_{y} conditions, the difference forcings in the polar cap are larger by a factor of about 2 than those in the auroral region. The magnitude of the difference forcing decreases rapidly with descending altitudes. The difference heating is the most significant for negative IMF B_{z} condition, and the difference heatings in the auroral region are larger by a factor of about 1.5 than those in the polar cap region. The magnitude of the difference heating decreases rapidly with descending altitudes.
It is confirmed that forcingtoheating ratio is large, indicating that the contribution of the forcing to the highlatitude lower thermospheric dynamics is stronger than the contribution of the heating to it and eventually drag from convecting ions is generally more important than Joule heating for generating dynamical energy in the highlatitude lower thermosphere. Especially, it is obvious that the forcingtoheating ratio is much larger in the polar cap region than in the auroral region and at higher altitude than at lower altitude, indicating that the contribution of the forcing to the dynamics is much larger in the polar cap region than in the auroral region and at higher altitude than at lower altitude. It is evident that when B_{z} is negative condition the contribution of the forcing is the lowest and the contribution of the heating is the highest among the different IMF conditions.
We have presented here only an analysis of the relative contributions of the momentum forcing and the heating to the highlatitude lower thermospheric dynamics. A future study will be devoted to investigating the available energy budget in the highlatitude lower thermosphere.

3. Heppner J. P 1972 [JGR] Vol.77 P.4877

5. Hernandez G, McCormac F. G, Smith R. W 1991 [JGR] Vol.96 P.5777

8. Kwak Y. S, Ahn B. H, Kim K. H 2008 [JASS] Vol.25 P.415

9. Kwak Y. S, Lee J. J, Ahn B. H, Hwang J, Kim K. H, Cho K. S 2008 [JASS] Vol.25 P.405

10. Kwak Y. S, Richmond A. D 2007 [JGR] Vol.112 P.A01306

11. Maeda S, FullerRowell T. J, Evans D. S 2007 [JGR] Vol.112 P.A06316

13. McCormac F. G, Smith R. W 1984 [GeoRL] Vol.11 P.935

14. McCormac F. G, Killeen T. L, Gombosi E, Hays P. B, Spencer N. W 1985 [GeoRL] Vol.12 P.155

15. McCormac F. G, Killeen T. L, Thayer J. P 1991 [JGR] Vol.96 P.115

16. McHarg M, Chun F, Knipp D, Lu G, Emery B. A, Ridley A 2005 [JGR] Vol.110 P.A08309

17. Meriwether J. W, Shih P 1987 [AnGeo] Vol.5A P.329

18. Niciejewski R. J, Killeen T. L, Johnson R. M, Thayer J. P 1992 [AdSpR] Vol.12 P.215

19. Niciejewski R. J, Killeen T. L, Won Y 1994 [JATP] Vol.56 P.285

21. Rees D, FullerRowell T. J 1989 [RSPTA] Vol.328 P.139

22. Rees D, FullerRowell T. J 1990 [AdSpR] Vol.10 P.197

25. Ruohoniemi J. M, Greenwald R. A 1996 [JGRA] Vol.101 P.21743

26. Sica R. J, Hernandez G, Emery B. A, Roble R. G, Smith R. W, Rees M. H 1989 [JGRA] Vol.94 P.11921

27. Thayer J. P, Killeen T. L, McCormac F. G, Tschan C. R, Ponthieu J. J, Spencer N. W 1987 [AnGeo] Vol.5A P.363

28. Weimer D. R 1995 [JGRA] Vol.100 P.19595

30. Won Y 1994

31. Zhang X. X, Wang C, Chen Y, Wang Y. L, Tan A, Wu T. S, Germany G. A, Wang W 2005 [JGR] Vol.110 P.A12208

[Fig. 1] Distribution of difference forcing term contributing to the change of the potential vorticity at 195, 177, 160, 146, 134, 125, 117, 111, and 106 km for interplanetary magnetic field (IMF) (By, Bz) values of (3.2, 0.0), (3.2, 0.0), (0.0, 2.0), and (0.0, 2.0) nT. These are obtained by subtracting values with zero IMF from those with nonzero IMF conditions.

[Fig. 2] Distribution of difference heating term contributing to the change of the potential vorticity at 195, 177, 160, 146, 134, 125, 117, 111, and 106 km for interplanetary magnetic field (IMF) (By, Bz) values of (3.2, 0.0), (3.2, 0.0), (0.0, 2.0), and (0.0, 2.0) nT. These are obtained by subtracting values with zero IMF from those with nonzero IMF conditions.

[Fig. 3] The difference forcingtoheating ratio at 195, 177, 160, 146, 134, 125, 117, 111, and 106 km for interplanetary magnetic field (IMF) (By, Bz) values of (a) (3.2, 0.0), (b) (3.2, 0.0), (c) (0.0, 2.0), and (d) (0.0, 2.0) nT. These ratio are obtained from the forcing term by dividing the heating term.