<p><strong>Abstract.</strong> Here we present results from an evaluation of model simulations from the International Hemispheric Transport of Air Pollution Phase II (HTAPII) and Chemistry Climate Model Initiative (CCMI) model inter-comparison projects against a comprehensive series of ground based, aircraft and satellite observations of ozone mixing ratios made at various locations across India. The study focuses on the recent past (observations from 2008&#8211;2013, models from 2008&#8211;2010) as this is most pertinent to understanding the health impacts of ozone. To our understanding this is the most comprehensive evaluation of these models' simulations of ozone across the Indian sub-continent to date. This study highlights some significant successes and challenges that the models face in representing the oxidative chemistry of the region.</p> <p>The multi-model range in area weighted surface ozone over the Indian subcontinent is 37.26&#8211;56.11&#8201;ppb, whilst the population weighted range is 41.38&#8211;57.5&#8201;ppb. When compared against surface observations from the Modelling Atmospheric Pollution and Networking (MAPAN) network of eight semi-urban monitoring sites spread across India, we find that the models tend to simulate higher ozone than that which is observed. However, observations of NO<sub><i>x</i></sub> and CO tend to be much higher than modelled mixing-ratios, suggesting that the underlying emissions used in the models do not characterise these regions accurately and/or that the resolution of the models is not adequate to simulate the photo-chemical environment about these surface observations. Empirical Orthogonal Function (EOF) analysis is used in order to identify the extent to which the models agree with regards to the spatio-temporal distribution of the tropospheric ozone column, derived using OMI-MLS observations. We show that whilst the models agree with the spatial pattern of the first EOF of observed tropospheric ozone column, most of the models simulate a peak in the first EOF seasonal cycle represented by principle component 1, which is later than the observed peak . This suggest a widespread systematic bias in the timing of emissions or some other unknown seasonal process.</p> <p>In addition to evaluating modelled ozone mixing ratios, we explore modelled emissions of NO<sub><i>x</i></sub>, CO, VOCs, and the ozone response to the emissions. We find a high degree of variation in emissions from non-anthropogenic sources (e.g. lightning NO<sub><i>x</i></sub> and biomass burning CO) between models. Total emissions of NO<sub><i>x</i></sub> and CO over India vary more between different models in the same MIP than the same model used in different MIPs, making it impossible to diagnose whether differences in modelled ozone are due to emissions or model processes. We therefore recommend targeted experiments to pinpoint the exact causes of discrepancies between modelled and observed ozone and ozone precursors for this region. To this end, a higher density of long term monitoring sites measuring not only ozone but also ozone precursors including speciated VOCs, located in more rural regions of the Indian sub-continent, would enable improvements in assessing the biases in models run at the resolution found in HTAPII and CCMI.</p>