“…In this study, two nested domains consisting of the first domain (D01) combined with the second domain (D02) were established in the coupled WRF/CMAQ models; details of the domain technical parameters are shown in Table S3 . Specific information for the model configuration, emission inventories, chemical transport mechanisms, initial/boundary conditions, and modelling performance verification can be found in detail in previous studies [ [65] , [66] , [67] , [68] ]. The simulation efficiency of ground-level O 3 concentration by this system was assessed based on a wide range of statistical indicators, including the Nash efficiency index, mean bias (MB), root mean-squared error (RMSE), normalised mean bias (NMB), normalised mean gross errors (NME), and R correlation coefficient from our previous studies [ 67 , 68 ].…”