We presented the development of the gaseous chemistry adjoint module of the meteorological-chemical model system GRAPES-CUACE (Global/Regional Assimilation and PrEdiction System coupled with CMA Unified Atmospheric Chemistry Environmental Forecasting System) on the basis of the previously constructed aerosol adjoint module. The latest version of the GRAPES-CUACE adjoint model mainly includes the adjoint of the physical and chemical processes, the adjoint of the transport processes, and the adjoint of interface programs, of both gas and aerosol. The adjoint implementation was validated for the full model, and adjoint results showed good agreement with brute force sensitivities. We also applied the newly developed adjoint model to the sensitivity analysis of an ozone episode occurred in Beijing on July 2, 2017, as well as the design of emission-reduction strategies for this episode. The relationships between the ozone concentration and precursor emissions were well captured by the adjoint model. It is indicated that for a case used here, the Beijing peak ozone concentration was influenced mostly by local emissions (6.2–24.3%), as well as by surrounding emissions, including Hebei (4.4–16.8%), Tianjin (1.8–6.6%), Shandong (1.8–2.6%), and Shanxi (<1%). In addition, reduction of NOx, VOCs, and CO emissions in these regions would effectively decrease the Beijing peak ozone concentration. This study highlights the capability of GRAPES-CUACE adjoint model in quantifying “emission-concentration” relationship and in providing guidance for environmental control policy.