We present an application of GPU accelerators in Earth system modeling. We developed software that generates CUDA kernels for numerical integration in the global climate model EMAC, used to study climate change and air quality. We focus on atmospheric chemical kinetics, the most computationally intensive task in climate–chemistry simulations. This approach can serve as the basis for hardware acceleration of numerous geoscientific models that rely on KPP for chemical kinetics applications.