Developing highly scalable algorithms for global atmospheric modeling is becoming increasingly important as scientists inquire to understand behaviors of the global atmosphere at extreme scales.Nowadays,heterogeneous architecture based on both processors and accelerators is becoming an important solution for large-scale computing.However,largescale simulation of the global atmosphere brings a severe challenge to the development of highly scalable algorithms that fit well into state-of-the-art heterogeneous systems.Although successes have been made on GPUaccelerated computing in some top-level applications,studies on fully exploiting heterogeneous architectures in global atmospheric modeling are still very less to be seen,due in large part to both the computational difficulties of the mathematical models and the requirement of high accuracy for long term simulations.