Model intercomparison studies often report a large spread in simulation results, but quantifying the causes of these differences is hindered by the fact that several processes contribute to the model spread simultaneously. Here we use the Multi-Scale Infrastructure for Chemistry and Aerosols (MUSICA) version 0 to investigate the model resolution dependencies of simulated chemical species, with a focus on the differences between global uniform grid and regional refinement grid simulations with the same modeling framework. We construct two global (ne30 [∼112 km] and ne60 [∼56 km]) and two regional refinement grids over Korea (ne30x8 [∼14 km] and ne30x16 [∼7 km]). The grid resolution can change chemical concentrations by an order of magnitude in the boundary layer, and the importance increases as the species' reactivity increases (e.g., up to 50% and 1,000% changes for ethane and xylenes, respectively). The diurnal cycle of oxidants (OH, O3, and NO3) also varies with the grid resolution, which leads to different oxidation pathways of volatile organic compounds (e.g., the fraction of monoterpenes reacting with NO3 in Seoul around midnight is 90% for ne30, but 65% for ne30x16). The models with high-resolution grids usually do a better job at reproducing aircraft observations during the KORUS-AQ campaign, but not always, implying compensating errors in the coarse grid simulations. For example, ozone is better reproduced by the coarse grid due to the artificial mixing of NOx. When developing new chemical mechanisms and evaluating models over urban areas, the uncertainties associated with model resolution should be considered. © 2023 The Authors. Journal of Advances in Modeling Earth Systems published by Wiley Periodicals LLC on behalf of American Geophysical Union.