Development of high-resolution Thermosphere-Ionosphere

Tong Dang,Binzheng Zhang,Jiuhou Lei,Wenbin Wang,Alan Burns,Han-li Liu,Kevin Pham,Kareem A. Sorathia

2020-01-01

Abstract:When solving hydrodynamic equations in spherical/cylindrical geometry using explicit finite difference schemes, a 14 major difficulty is that the time step is greatly restricted by the clustering of azimuthal cells near the pole due to the Courant15 Friedrichs-Lewy condition. This paper adapts the azimuthal averaging-reconstruction (Ring Average) technique to finite differ16 ence schemes in order to mitigate the time step constraint in spherical/cylindrical coordinates. The finite-difference Ring Average 17 technique averages physical quantities based on an effective grid and then reconstructs the solution back to the original grid in a 18 piece-wise, monotonic way. The algorithm is implemented in a community upper atmospheric model Thermosphere-Ionosphere 19 Electrodynamics General Circulation Model (TIE-GCM), with horizontal resolution up to 0.625◦ × 0.625◦ in the geographic 20 longitude-latitude coordinates, which enables the capability of resolving critical mesoscale structures within the TIE-GCM. Nu21 merical experiments have shown that the Ring Average technique introduces minimal artifacts in the polar region of the GCM 22 solutions, which is a significant improvement compared to the commonly used low-pass filtering techniques such as the fast 23 Fourier transform method. Since the finite-difference adaption of the Ring Average technique is a post-solver type algorithm, 24 which requires no changes to the original computational grid and numerical algorithms, it has also been implemented in much 25 more complicated models with extended physical/chemical modules such as the coupled Magnetosphere Ionosphere Thermo26 sphere (CMIT) model and the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension 27 (WACCM-X). The implementation of the Ring Average techniques in both models enables CMIT and WACCM-X to perform 28 global simulations with a much higher resolution than that used in the community versions. The new technique is not only a 29 significant improvement in space weather modeling capability, but can also be adapted to more general finite difference solvers 30 for hyperbolic equations in spherical/polar geometries. 31 Key Points: 32 • Ring Average technique is adapted to solve the issue of clustered grid cells in polar/spherical coordinate with finite 33 difference method. 34 • The Ring Average technique is applied to develop a 0.625◦ × 0.625◦ high-resolution TIE-GCM and more complicated 35 geoscientific models. 36 • The high-resolution TIE-GCM shows good capability in resolving mesoscale structures in the I-T system. 37