Earth's ionosphere and thermosphere (I-T) system is closely coupled with the magnetosphere, and the electromagnetic energy from the magnetosphere is transferred into the I-T system through field-aligned currents (FACs). The major part of the electromagnetic energy is irreversibly converted into heat through ohmic currents, and such heat is called Joule heating (Cole, 1962;Richmond, 2021;Thayer, 2000). Joule heating can significantly affect the I-T system both locally and globally, especially during geomagnetic storms. For example, the neutral temperature and density increase due to the enhanced Joule heating during geomagnetic storms (e.g., Fuller-Rowell et al., 1994). In addition, Joule heating can effectively change the global circulation within several hours, which markedly alters the thermospheric compositions at different latitudes and can further change the ionospheric electron density (e.g., Buonsanto, 1999;Prölss, 2011). Moreover, gravity waves can be launched due to rapid variations of Joule heating and they can propagate globally, causing large-scale traveling atmospheric disturbances and traveling ionospheric disturbances (e.g., Lu et al., 2016Lu et al., , 2020. A comprehensive review of Joule heating and the I-T response to Joule heating during geomagnetic storms can be found in Richmond (2021).General circulation models (GCMs) of the I-T system are widely used to study variations of the I-T system particularly during geomagnetic storms, and accurate estimations of Joule heating are critical for reproducing observed features. Joule heating in GCMs is calculated from the electric field, conductivities associated with the solar ionization and electron precipitation together with the neutral winds (e.g., Lu et al., 1995). However, accurate estimations of Joule heating is still challenging to date since it is difficult to capture the dynamic variations of the electric field, ionospheric conductivity (mostly associated with the electron precipitation), and neutral winds (e.g.,