Simulation of O3 and NOx in São Paulo street urban canyons with VEIN (v0.2.2) and MUNICH (v1.0)

We evaluate the performance of the Model of Urban Network of Intersecting Canyons and Highways (MUNICH) in simulating ozone (O3) and nitrogen oxides (NOx) concentrations within the urban street canyons in the São Paulo metropolitan area (SPMA). The MUNICH simulations are performed inside the Pinheiros neighborhood (a residential area) and Paulista Avenue (an economic hub), which are representative urban canyons in the SPMA. Both zones have air quality stations maintained by the São Paulo Environmental Agency (CETESB), providing data (both pollutant concentrations and meteorological) for model evaluation. Meteorological inputs for MUNICH are produced by a simulation with the Weather Research and Forecasting model (WRF) over triple-nested domains with the innermost domain centered over the SPMA at a spatial grid resolution of 1 km. Street coordinates and emission flux rates are retrieved from the Vehicular Emission Inventory (VEIN) emission model, representing the real fleet of the region. The VEIN model has an advantage to spatially represent emissions and present compatibility with MUNICH. Building height is estimated from the World Urban Database and Access Portal Tools (WUDAPT) local climate zone map for SPMA. Background concentrations are obtained from the Ibirapuera air quality station located in an urban park. Finally, volatile organic compound (VOC) speciation is approximated using information from the São Paulo air quality forecast emission file and non-methane hydrocarbon concentration measurements. Results show an overprediction of O3 concentrations in both study cases. NOx concentrations are underpredicted in Pinheiros but are better simulated in Paulista Avenue. Compared to O3, NO2 is better simulated in both urban zones. The O3 prediction is highly dependent on the background concentration, which is the main cause for the model O3 overprediction. The MUNICH simulations satisfy the performance criteria when emissions are calibrated. The results show the great potential of MUNICH to represent the concentrations of pollutants emitted by the fleet close to the streets. The street-scale air pollutant predictions make it possible in the future to evaluate the impacts on public health due to human exposure to primary exhaust gas pollutants emitted by the vehicles.

Abstract. We evaluate the performance of the Model of Urban Network of Intersecting Canyons and Highways (MU-NICH) in simulating ozone (O 3 ) and nitrogen oxides (NO x ) concentrations within the urban street canyons in the São Paulo metropolitan area (SPMA). The MUNICH simulations are performed inside the Pinheiros neighborhood (a residential area) and Paulista Avenue (an economic hub), which are representative urban canyons in the SPMA. Both zones have air quality stations maintained by the São Paulo Environmental Agency (CETESB), providing data (both pollutant concentrations and meteorological) for model evaluation. Meteorological inputs for MUNICH are produced by a simulation with the Weather Research and Forecasting model (WRF) over triple-nested domains with the innermost domain centered over the SPMA at a spatial grid resolution of 1 km. Street coordinates and emission flux rates are retrieved from the Vehicular Emission Inventory (VEIN) emission model, representing the real fleet of the region. The VEIN model has an advantage to spatially represent emissions and present compatibility with MUNICH. Building height is estimated from the World Urban Database and Access Portal Tools (WUDAPT) local climate zone map for SPMA. Background concentrations are obtained from the Ibirapuera air quality station located in an urban park. Finally, volatile organic compound (VOC) speciation is approximated using information from the São Paulo air quality forecast emission file and non-methane hydrocarbon concentration measurements. Results show an overprediction of O 3 concentrations in both study cases. NO x concentrations are underpredicted in Pin-heiros but are better simulated in Paulista Avenue. Compared to O 3 , NO 2 is better simulated in both urban zones. The O 3 prediction is highly dependent on the background concentration, which is the main cause for the model O 3 overprediction. The MUNICH simulations satisfy the performance criteria when emissions are calibrated. The results show the great potential of MUNICH to represent the concentrations of pollutants emitted by the fleet close to the streets. The street-scale air pollutant predictions make it possible in the future to evaluate the impacts on public health due to human exposure to primary exhaust gas pollutants emitted by the vehicles.

Introduction
Street urban canyons are structures formed by a street and its flanked buildings (Oke et al., 2017). Due to their proximity to emissions from vehicles and their side function as a compartment that limits pollutant dispersion, the street and the associated urban canyons are considered pollutant hotspots (Zhong et al., 2016). As more people start to live in urban areas (United Nations, 2018) and the ubiquity of urban canyons in cities grows, pedestrians, commuters, bikers, and drivers are being exposed to high pollutant concentrations every day (Vardoulakis et al., 2003). Consequently, the study of air pollution inside urban canyons is an important matter when dealing with studies of human health exposure related to traffic emissions.
To estimate the real impact of the pollutants on human health, it is necessary to obtain accurate pollutant concentrations and the lengths of exposure. Most cities are not covered by a high-density network of air quality stations. Even though the measurements provide precise information, it is expensive and also very difficult to cover all of the impacted areas of a city (Zhong et al., 2016). One alternative, that is starting to be contemplated, is the use of numerical modeling to represent the pollutant behavior in urban canyons, which has the advantage of producing pollutant concentration information at high temporal and spatial resolutions.
Computational fluid dynamics (CFD) models are considered to be the best modeling approach to understand air pollutant dispersion inside the urban areas. Due to the limitations of high computational resources, these models cannot be applied for long-time simulation periods nor for a large area (Fellini et al., 2019;Thouron et al., 2019).
A new type of model, the urban-/local-scale operational model, overcomes these limitations by applying simplifications on urban geometry and parameterizations of the mass transfer processes of air pollutants inside the urban canyons. The Operational Street Pollution Model (OSPM) and the Atmospheric Dispersion Model System (ADMS-urban) are two of the most popular operational models which have already been tested for different cities around the world (Berkowicz et al., 1997;McHugh et al., 1997). Their main advantage is that they calculate pollutant concentrations when sources and receptors are in the same street urban canyon, but they present a limited treatment for the pollutant transfer between streets and intersections (Carpentieri et al., 2012).
Street-network models are also operational, having the advantage of dealing with the transport of pollutants in city street intersections. The SIRANE model uses parametric relations to solve advection on the streets, the dispersion in the street intersections, and interchange between the streets and the over-roof atmosphere (Soulhac et al., 2011. Background concentrations at the over-roof atmosphere are estimated using a Gaussian plume model. This estimation method inhibits a comprehensive atmospheric chemistry treatment. Recently, the Model of Urban Network of intersecting Canyons and Highways (MUNICH) was developed by Kim et al. (2018a) using a similar parameterization to that of SIRANE. MUNICH includes improvements in the treatment of the mean wind profile inside the urban canyon and the turbulent vertical mass transfer at the top of the street. It solves pollutant reactions using a chemical mechanism, so it can also simulate the production of ozone inside the urban canyons. MUNICH has been used to simulate ozone (O 3 ) and nitrogen oxides (NO x ) by Wu et al. (2020) in the Tianhe district of Guangzhou city and NO x as part of the Street in Grid (SinG) model in Kim et al. (2018a), Thouron et al. (2019), and Lugon et al. (2020) in the Paris region.
Significant information is required to run this kind of model. It is explained by Vardoulakis et al. (2003) that, in general, these models need at least information from traffic data, emissions, meteorological data, street geometry, and background concentrations. Recently, the Vehicular Emission Inventory (VEIN) model was developed by Ibarra-Espinosa et al. (2018) using information for São Paulo. VEIN is suitable for use in street-network models because it uses the traffic flow, emission factors, and street morphology (i.e., intersection coordinates) to calculate the vehicular emissions. As a matter of fact, due to its architecture, it can be used together with MUNICH.
In Brazil, previous studies of air quality in urban canyons dealt with measurements of black carbon and O 3 inside a street canyon in Londrina's city center (Krecl et al., 2016), and dispersion of NO x was simulated in Curitiba with the ENVI-met model (Krüger et al., 2011). To our knowledge, this is the first study of modeling O 3 and NO x inside street urban canyons in the São Paulo metropolitan area (SPMA), the biggest megacity in South America, where it is very often the exceedance of O 3 state air quality standard (Andrade et al., 2017).
As the management of secondary pollutants remains a challenge in SPMA, we aim to evaluate MUNICH operational street-network model to simulate O 3 and NO x concentration inside urban canyons, coupled with the VEIN emission model, to build a street-level air quality modeling system. This modeling system can be used in air quality and traffic management of the São Paulo neighborhoods, in studies of health effects from traffic emission exposure, in future urban planning, and post-accident analysis.

Data and methods
The experiment consisted of carrying out simulations of O 3 , NO x , NO, and NO 2 concentrations inside the SPMA urban street canyons with the MUNICH model. To evaluate model performance, the model results are compared against the measurements from the São Paulo Environmental Agency (CETESB) air quality network. We choose the Pinheiros urban area to test the model, where there is an air quality station in a mixed residential-commercial area. Once MUNICH and VEIN are calibrated, a study case is prepared by calculating the pollutant concentration inside Paulista Avenue, the economic central area of the city with high canyons. The selected study period covers the week from 6 to 13 October of 2014. This period is chosen before dry weather conditions in SPMA, a period of high O 3 concentrations , the availability of data, and the availability of the emission inventory developed for a typical week in October 2014 (Ibarra-Espinosa et al., 2020).

MUNICH model
MUNICH is conceptually based on the SIRANE model (Soulhac et al., 2011). It has two main components: the street-canyon component, which deals with and solves pollutant concentrations inside the urban-canopy volume, and the intersection component, which calculates the pollutant concentrations inside the intersection volume. MUNICH differs from SIRANE in the treatment of the vertical flux by turbulent diffusion at the roof level (Schulte parameterization;Schulte et al., 2015) and in the mean wind velocity within the street canyon (Lemonsu parameterization;Lemonsu et al., 2004). Currently, MUNICH solves gas-phase pollutants based on the Carbon Bond mechanism version 5 (CB05). Further information is detailed in Kim et al. (2018a).

VEIN emission model
VEIN is an R package (R Core Team, 2020) to estimate vehicular emissions at the street level. VEIN imports functions from the Spatial Features package (Pebesma, 2018), which represent different types of geometries in space and perform geoprocessing tasks, from the data table package (Dowle and Srinivasan, 2019) to perform fast aggregation of databases, and from the units package (Pebesma et al., 2016) to provide binding to the udunits library (https://www.unidata.ucar.edu/ software/udunits/, last access: 28 May 2021). VEIN includes a function to process vehicular flow at each street to generate activity traffic data, different emissions factors, and different sets of emissions calculation and post-processing tools (Ibarra-Espinosa et al., 2018). Specifically, the emissions factors are based on emissions certification tests with dynamometer measurements in laboratories (CETESB, 2015).

MUNICH input data
Urban canyon models required detailed input information, such as building height and street geometry. Their performance depends on the quality of this information (Vardoulakis et al., 2003). In recent years, new tools have been developed to generate this information. Table 1 summarizes the model input used in this simulation experiment.

Emissions and street links coordinates
The vehicular fleet is the principal source of air pollution in SPMA (Andrade et al., , 2017. The particularity of this fleet is the extensive use of biofuels (i.e., gasohol, ethanol, and biodiesel). During 2014, vehicular emissions were responsible for emitting 97 % of CO, 82 % of volatile organic compounds (VOCs), 78 % of NO x , and 40 % of particulate matter (CETESB, 2015). Vehicular emissions inside SPMA streets were estimated using the VEIN emission model (Ibarra-Espinosa et al., 2018).
Street links are segments of roads split at each vertex. Then, a road can be composed of many links. Emission rates inside these street links in the VEIN model are calculated using 104 million GPS vehicles coordinates in southeast Brazil (Ibarra-Espinosa et al., 2019). The GPS dataset is assigned to the OpenStreetMap (2017) dataset and once traffic flow is obtained, the vehicular compositions are generated and assigned to each emission factor reported by CETESB (2015). Emission factors are transformed into speed functions, and then the average speed calculated at each street is used to obtain more representative emissions at each hour of a week. In addition, the estimation was calibrated with fuel consumption for the year 2014. Ibarra-Espinosa et al. (2020a) described all details regarding the emission estimation with the emissions dataset in g h −1 available at https://github.com/ ibarraespinosa/ae1 (last access: 28 May 2021).
The emissions dataset presents two aspects that need to be discussed. The first one is that there are some differences between the traffic flow from travel demand model (TDM) outputs and GPS (Ibarra-Espinosa et al., 2019, 2020a. The ratio between traffic flows from TDM and GPS for our study is 2.22. Regarding the emissions factors used to estimate the emissions, they are based on average measurement of emissions certification tests (CETESB, 2015); therefore, they may underestimate real-drive emissions (Ropkins et al., 2009). For instance, the real-world emission factors derived from tunnel measurements in São Paulo for NO x were 0.3 g km −1 for light vehicles and 9.2 g km −1 for heavy vehicles (Pérez-Martínez et al., 2014), while the respective fleet-weighted CETESB (2015) emission factors are 0.26 and 6.68 g km −1 , as shown in Fig. S1 in the Supplement, resulting in ratios of 1.11 and 1.38. Then, if we consider the mean emission-factor ratio (1.11+1.38)/2 multiplied by the mentioned traffic flow ratio (2.22) results, the NO x emissions might be approximately 2.73 higher than those estimated using pure CETESB (2015) data. Consequently, we expect that air quality simulations for NO x might be lower than observations. Even when VEIN produces hourly emissions for a standard week (Fig. S2 in the Supplement), MUNICH only considers a standard day for weekdays and weekends. We choose Wednesday emission as a typical weekday and Saturday emission for the weekend. Figure 1 shows the mean diurnal profile of NO x and VOCs emission fluxes from street links in the Pinheiros neighborhood.

WRF simulation
Triple-nested domains are set up centered in SPMA. The mother domain has a spatial resolution of 25 km, the second 5 km, and the finest 1 km. The simulation at 1 km provides MUNICH with meteorological information. Initial and boundary conditions are retrieved from Historical Unidata Internet Data Distribution (IDD) Gridded Model Data (https://rda.ucar.edu/datasets/ds335.0/index.html, last access: 28 May 2020). Table 2 shows WRF configuration and Fig. 2 show the WRF domains.
Before using the WRF simulation outputs for MUNICH modeling, a model verification is performed. Model verification was carried out for the same period as MUNICH runs and for the finest domain output (D03). We used meteorolog-  ical information from 16 air quality stations whose locations are shown in Fig. 4. We also use benchmarks suggested by Emery et al. (2001), which were also used in Reboredo et al. (2015) and Pellegati Franco et al. (2019). However, Monk et al. (2019) explained that these benchmarks are suitable for domains in "simple" terrain; they also presented other sets of benchmarks for "complex" terrain, the latter being more suitable for SPMA. The results are detailed in Table 3. The temperature at 2 m (T 2) and relative humidity at 2 m (RH2) reach the simple terrain benchmarks, while wind speed and direction at 10 m (WS10 and WD10, respectively) are very close to them. When compared against complex terrain benchmarks, only the mean bias of WD10 is beyond the benchmark. Finally, T 2, RH2, and WS10 satisfy the good performance criteria of Keyser and Anthes (1977) and Pielke (2013). More details are shown in Tables S1 and S2 in the Supplement.

Building height and street width
Building height is retrieved from the World Urban Database and Access Portal Tools project (WUDAPT) for SPMA ( Fig. 3). WUDAPT classifies urban areas into 17 local climate zones (LCZs). These LCZs are divided into build types, which are LCZs from 1 to 10, and land cover types, which go from A to G. Each of these LCZs presents different thermal, radiative, surface cover, and geometric properties. The building height is the height of roughness elements, which is the geometric average of building heights (Stewart and Oke, 2012). The WUDAPT file for SPMA is a raster with a spatial resolution of 120 m and was previously used in Pellegati The number of lanes is provided by the OpenStreetMap dataset, so the street width is calculated by using 3 m of lane width and by adding 1.9 m to each side of the street as sidewalk width. Most OpenStreetMap streets do not include the number of lanes for this region; therefore, they are hole filled with the average by type of street.   Vardoulakis et al. (2003) explained that the background concentration in street modeling is necessary to include the proportion of air pollutants that are not emitted inside the street. In the SinG model, background concentrations are the concentrations calculated by Polair3D, a mesoscale air quality model (Kim et al., 2018a). Wu et al. (2020) chose measurements from a station located very close to the study zone as the background concentration. Consequently, we consider the concentration outside the MUNICH domain as background concentration. With that in mind, by using the mean wind field from the WRF simulation for the study period, we select Ibirapuera air quality station (AQS) (83 shown in Fig. 4) measurements as background concentration, which, according to the wind field, advect pollutants to Pinheiros station (99) and Cerqueira César (83) as can be seen in Fig. 4. This assumption is only valid during daylight, when ozone concentrations are higher. As seen in Fig. S3 in the Supplement, during nighttime, wind presents a westerly direction.

Background concentration
Measurements of O 3 , NO 2 , and NO in Ibirapuera AQS were used as background concentrations.

Measurements and statistical analysis
Meteorological and air pollutant measurements are retrieved from the CETESB air quality network. To evaluate WRF simulation in the finest domains, observations from 16 AQSs are used. Background concentration comes from the Ibirapuera AQS. The Pinheiros AQS is used to evaluate MUNICH performance in the Pinheiros neighborhood, while Cerqueira César is used to evaluate Paulista Avenue. To evaluate model performance, we follow the recommendations from Emery et al. (2017). We also use the evaluation statistics from Hanna and Chang (2012): fractional bias (FB), normalized mean square error (NMSE), fraction of predictions within a factor of 2 (FAC2), and normalized absolute difference (NAD). The acceptance criteria for urban zones are |FB| < = 0.67, NMSE < = 6, FAC2 > = 0.3, and NAD < = 0.5. We expand the statistical analysis to the background concentration to see the difference against observations and to assess the influence of background concentration in MUNICH simulations.

Model set up
We use MUNICH to simulate two urban areas inside SPMA: the first domain is the Pinheiros neighborhood and the second one is Paulista Avenue. VEIN produces emissions for all the street links in SPMA. This information can be filtered by the neighborhood name of the street links. We subset that information for the Pinheiros neighborhood (Fig. 5a) and for the neighborhoods that contain the Paulista Avenue urban canyon (Fig. 5b). In MUNICH, NO emissions are estimated from NO x and NO 2 emissions. Figure 5 shows MUNICH domain for the Pinheiros neighborhood and Paulista Avenue. The yellow dot represents the location of the air quality stations. The red lines are the street links used by VEIN to calculate the emissions, and the yellow rectangle is the urban canyon selected for comparison against observations.
There are 677 street links for Pinheiros and 535 for Paulista Avenue. In total, nine points of WRF simulation cover the Pinheiros domains, while 12 WRF points represent Paulista Avenue domains. From WUDAPT, we can see that inside Pinheiros there is a variety of buildings with different heights. The Pinheiros AQS is located in an urban canyon that has a mean building height of 5 m (LCZ 6 -open low rise). On the other hand, the Paulista Avenue domain is more uniform, presenting urban canyons with a mean building height of 45 m (LCZ1 -compact high rise).

Results
Here, we present the O 3 and NO x simulations with MUNICH for a week in October 2014. We first calibrated the input emissions by studying the Pinheiros neighborhood to later simulate NO x inside the Paulista Avenue urban canyon. Figure 6 shows the results of MUNICH simulation using the original emissions calculated by VEIN for SPMA. MU-NICH simulations are very close to background concentrations, which leads to an overprediction of O 3 and underpredicted NO and NO x concentrations. This is produced by a dependence of MUNICH on background concentration and by emission underestimation. The emission underestimation is caused by emission factors calculated based on average measurements of emissions certification tests and because emission factors derived from a dynamometer, and cycle measurements do not represent real-drive emissions (Ropkins et al., 2009). It is also probable that the number of vehicles could have been underestimated inside the urban canyon. The underestimation of NO x is caused by the underestimation of NO concentrations. NO 2 concentration magnitude is well represented by MUNICH.

Control case for the Pinheiros neighborhood
The diurnal variations of MUNICH simulation, observation, and background concentrations are shown in Fig. 7. MUNICH coherently simulated the temporal variation of O 3 and NO 2 concentrations inside the urban canyon. For NO and NO x , the temporal variation during the day and until midnight is well simulated, while the morning peak at 06:00 LT is underestimated. After midnight, a higher concentration of NO x occurs with the increase of heavy-duty vehicles at night that mainly run on diesel. In Pinheiros urban canyons, there is a predominant flow of light-duty vehicles, even though it is registered high NO x concentrations that it is transported from the highway. The mean differences between MUNICH simulation and background concentration for O 3 , NO x , NO,  (91), and the blue triangle shows Ibirapuera AQS (83). Circles represent AQS that only measures pollutant concentrations; stars represent AQS that also measures T 2, RH2, WS10, and WD10; diamonds represents AQS that also measures WS10 and WD10; triangles represent AQS that also measures T 2 and RH2. and NO 2 are −13.10, 28.61, 9.25, and 14.43 µg m −3 , respectively.

Emission adjustment
We ran different scenarios with increased NO x and VOC emission from VEIN. The best results were produced when we doubled the NO x and VOC emissions; this scenario is called MUNICH-Emiss. With this adjustment, we achieved an overall improvement of MUNICH simulations. Figure 8 shows the new comparison between the model, background concentration, and observations. O 3 is still overpredicted, which is caused by the higher value of O 3 background concentration together with a low NO background concentration; nevertheless, the simulated O 3 concentration during nighttime is well represented and daily peaks values are closer to observations. NO x and NO simulations are still underpredicted, but NO 2 is of the same magnitude as observations. NO x underprediction is still mainly attributed to the underprediction of NO, especially during 8, 9, and 10 October when high observational values of NO were recorded. NO underestimation is explained by the lower NO background concentration, the underestimation of emissions, and the use of a single-day emission profile to represent all weekdays. Wind speed overestimation also affects this underestimation as it enhances dispersion. However, MUNICH can better represent the observed high concentration during Saturday, October 11, as MUNICH uses the same emission profile for the weekend and weekdays; this high simulated NO concentration resulted from the influence of meteorology. Figure 9 shows the diurnal profiles for this simulation. The new MUNICH-Emiss profiles are closer to observed concentration profiles, with a better representation of the peak concentration magnitude of NO x , NO, and NO 2 . The mean differences over the simulation period between simulated and the background concentrations for O 3 , NO x , NO, and NO 2 are −17.85, −57.26, 23.60, and 21.07 µg m −3 , respectively,   showing bigger differences than the control case previous scenario and the influence of the reaction with NO emissions. Table 4 summarizes the performance statistics for each scenario and background. The performance statistics from the MUNICH-Emiss case show lower values of MB, NMGE, and RMSE for all pollutants, except NO 2 which presents a slight increase in these indicators. They also show high values of R (≥ 0.7) for each pollutant in every case, which indicates that the temporal variations of emission and background concentration are in the same phase as the observations. In general, in both MUNICH simulations, NO 2 and O 3 are better simulated. The MUNICH-Emiss case performs better and also achieves the recommendations of Hanna and Chang (2012) for O 3 , NO 2 NO, and NO x , whereas the MU-NICH control case did not reach these recommendations for NO. Figure 10 shows the mean hourly concentration of O 3 and NO x in the Pinheiros neighborhood; the red diamond points to the location of Pinheiros air quality station. Because the VEIN model can distribute spatially the emissions, there is a variation of concentrations in different street links. For example, the orange diamond shows the location of a traffic light, where traffic jams occur, causing lower O 3 concentrations from higher NO x emissions.
We also perform an additional sensitivity simulation by running the MUNICH scenario using the background concentrations from the Santos AQS (light blue triangle in Fig. 4). Compared to the Ibirapuera AQS site, measured O 3 and NO 2 concentrations are lower, and those of NO concentrations are higher at the Santos AQS. This results in O 3 and NO 2 underprediction and a better simulation of NO concentration magnitude; however, all evaluated pollutants present lower R values and higher NMGE values than MUNICH-Emiss scenario with Ibirapuera AQS as the background concentration. Simulated NO 2 and O 3 follow background concentrations, which indicates that the MUNICH simulations have a strong dependence on the background concentration (see Figs. S4 and S5 in the Supplement).
Lastly, a sensitivity simulation was performed with an only increase of NO x emissions by four and remaining VOC original emission using Ibirapuera background concentration. This results in a better O 3 representation but unrealistic NO x , NO, and NO 2 concentration (see Figs. S6 and S7 in the Supplement). As SPMA has a VOC-limited regime (Andrade et al., 2017), the increment of NO x emission will lead to a reduction of O 3 concentration. Many studies have shown that São Paulo's atmosphere is VOC limited (Schuch et al., 2020) due to the high NO x emission by the heavy-duty vehicles that are under old emissions regulations. The new regulations for diesel engine emissions were established recently and are being implemented according to the recycling of the fleet, which is 20 years of use for diesel trucks (CETESB, 2019).

Application for Paulista Avenue
The MUNICH simulation is performed with calibrated emissions for a domain that contains a well-defined urban canyon: Paulista Avenue. The simulation shows a better representation of NO x , NO, and NO 2 temporal variations and a good representation of concentration magnitude (Fig. 11). Although the MB indicates an overprediction of NO x , NO, and NO 2 (Table 5), Fig. 12 shows that this is caused by an overprediction of these pollutants during night hours, linked to a mismatch of emissions. As in the Pinheiros domain, MU-NICH did not capture the two peaks of NO and NO x during nighttime. This is caused by WRF limitation in representing planetary boundary layer height during nighttime (Hu et al., 2012;McNider and Pour-Biazar, 2020). Also, as shown in Fig. 1a, the NO x emission profile during weekdays presents two peaks during daylight at 07:00 and 16:00 LT (local time) and a smaller emission peak around 23:00 h; it is probable that this nighttime peak was underestimated.
Statistics in Table 5 show an improvement in representing concentration magnitudes of NO x , NO, and NO 2 with mean simulated concentrations close to observations and very low values of MB, NMB, and RMSE. In this case, R values are lower than those in the Pinheiros case but still higher than 0.4 for NO x and NO 2 , confirming that there is a mismatch of simulated concentrations, which is clearer in the MUNICH NO x and NO peak happening before observation. The MUNICH-Emiss simulations achieve Hanna and Chang (2012) performance criteria for NO x and NO 2 . NO 2 is the best simulated species.

Discussion and conclusions
Simulating air pollutants inside urban street canyons is a challenging task. It is even more difficult in cities as heterogeneous as São Paulo, where its urban structure is not always textbook defined. The limited number of air quality stations located inside or near urban canyons, together with the lack of information from detailed emission inventories and urban morphology data, hinders accurate air quality modeling and consequently air quality management.
In this paper, we attempt to fill in this gap by using the MUNICH street-network model together with the VEIN vehicular emissions model. The latter provides temporal and spatially detailed emission fluxes inside the main streets and coordinates and width of the streets (i.e., the street network). The urban morphology is completed by extracting the building height from the WUDAPT database for the São Paulo metropolitan area. The advantages of using MUNICH are that, besides solving pollutant dispersion, it also solves photochemistry reactions and is an operational model that solves pollutant concentration at neighborhood scale considering street intersections.   Results showed that MUNICH simulations that used adjusted emissions can better represent the temporal variation of O 3 , NO x , NO, and NO 2 concentrations inside urban canyons. Nevertheless, the results are highly dependent on background concentrations and emission fluxes. This background concentration dependence is stronger in secondary pollutants such as O 3 , and primary pollutants are more determined by emission fluxes. The reason for the significant contribution of background concentration is that MUNICH is based in SIRANE, and SIRANE also presents a significant contribution from background concentration .
The main cause of O 3 overprediction in our simulation for both tested urban zones is the high value of background O 3 concentration measured in the Ibirapuera AQS. In the Pinheiros neighborhood, the underprediction of NO x concentration is caused by the underprediction of NO concentration in Pinheiros during the second half of the week. This un-  derestimation is caused by the lower NO background concentration together with an emission underestimation. The concentration magnitudes in Paulista Avenue are well represented but there was a mismatch with observed concentration. MUNICH-Emiss scenario fulfills the performance criteria. O 3 concentration simulated in Pinheiros and Paulista Avenue is lower than background concentrations; these same results are reported by Wu et al. (2019). As noted in Krecl et al. (2016), this behavior is caused by the high NO x emissions inside the street urban canyons, which rapidly deplete the formed O 3 and the one from the rooftop (i.e., background concentration).
As the main source of surface NO and NO 2 emissions in São Paulo are vehicles, it is necessary to go deeper into the reasons why the MUNICH-Emiss scenario performs better. The increase of the emissions is necessary because the emissions factors are the average of emission certification tests (CETESB, 2015). It has been shown that emission factors derived from dynamometer and cycle measurements do not represent real-drive emissions (Ropkins et al., 2009). São Paulo does not have an inspection and maintenance (I&M) pro-   gram; therefore, there may exist a fraction of the fleet which are high emitters and do not meet the emission standards; more details can be found in Ibarra-Espinosa et al. (2020a). Furthermore, the comparison of traffic flow between GPS and TDM data for the Pinheiros area showed that TDM traffic flows are 2.22 times higher than GPS. Hence, more representative traffic flows would also improve the emissions compilation. As a conclusion, it is important to develop new and more representative vehicular traffic flow and emission factors for Brazil.
With calibrated emissions (i.e., MUNICH-Emiss scenario), the good performance of MUNICH in representing NO 2 concentrations in both neighborhoods and NO and NO x in Paulista Avenue urban canyon suggests that the VEIN model distributes emissions spatially and temporally efficiently, which proves its potential to be used in other cities. VEIN is being continuously developed and currently offers some utilities to format emissions to the MUNICH model. On the other hand, now Google Earth allows new features such as 3-D view that together with in situ measurements can improve WUDAPT building height estimates. These new features can be used to improve MUNICH input data and therefore the model simulation results. Further, a better estimation of background concentrations from photochemical grid models can potentially improve the model performance.
The results obtained show the promising capability of MUNICH to represent the concentrations of pollutants emitted by the fleet close to the streets. As MUNICH uses the 3264 M. E.  and NO x in São Paulo street urban canyons CB05 gas-phase mechanism, it can also simulate VOCs inside the urban canyon. Measurements of VOCs inside urban canyons are therefore necessary to validate the model in the future. An accurate prediction of street-scale air pollutant concentrations will enable the future assessment of the impacts on human health due to their exposure to air pollutants emitted by the vehicles.    Kim et al., 2018b). VEIN can be installed from CRAN, and it is also available on Zenodo (https://doi.org/10.5281/zenodo.3714187, Ibarra-Espinosa et al., 2020b). Additional information and help are available by contacting the authors.
Author contributions. MEGC performed the simulations and prepared the manuscript with the support of all co-authors. MEGC, MdFA, and YZ designed the experiment. SIE provided the emissions and street morphology information. YK provided support to set up and run MUNICH. MEGC, YZ, MFA, and SIE discussed the results.
Competing interests. The authors declare that they have no conflict of interest.