<p><strong>Abstract.</strong> Fine particulate matter (PM<sub>2.5</sub>) and surface ozone (O<sub>3</sub>) are major air pollutants in megacities such as Delhi, but the design of suitable mitigation strategies is challenging. Some strategies for reducing PM<sub>2.5</sub> may have the notable side-effect of increasing O<sub>3</sub>. Here, we demonstrate a numerical framework for investigating the impacts of mitigation strategies on both PM<sub>2.5</sub> and O<sub>3</sub> in Delhi. We use Gaussian process emulation to generate a computationally efficient surrogate for a regional air quality model (WRF-Chem). This allows us to perform global sensitivity analysis to identify the major sources of air pollution, and to generate emission-sector based pollutant response surfaces to inform mitigation policy development. Based on more than 100,000 emulation runs during the pre-monsoon period (peak O<sub>3</sub> season), our global sensitivity analysis shows that local traffic emissions from Delhi city region and regional transport of pollutions emitted from the National Capital Region surrounding Delhi (NCR) are dominant factors influencing PM<sub>2.5</sub> and O<sub>3</sub> in Delhi. They together govern the O<sub>3</sub> peak and PM<sub>2.5</sub> concentration during daytime. Regional transport contributes about 80&#8201;% of the PM<sub>2.5</sub> variation during the night. Reducing traffic emissions in Delhi alone (e.g., by 50&#8201;%) would reduce PM<sub>2.5</sub> by 15&#8211;20&#8201;% but lead to a 20&#8211;25&#8201;% increase in O<sub>3</sub>. However, we show that reducing NCR regional emissions by 25&#8211;30&#8201;% at the same time would further reduce PM<sub>2.5</sub> by 5&#8211;10&#8201;% in Delhi and avoid the O<sub>3</sub> increase. This study provides scientific evidence to support the need for joint coordination of controls on local and regional scales to achieve effective reduction on PM<sub>2.5</sub> whilst minimize the risk of O<sub>3</sub> increase in Delhi.</p>