<figure><img src="image_1.png" /></figure>

HAL
open science

# Modelling concentration heterogeneities in streets using the street-network model MUNICH

Thibaud Sarica, Alice Maison, Yelva Roustan, Matthias Ketzel, Steen Solvang

Jensen, Youngseob Kim, Christophe Chaillou, Karine Sartelet

► To cite this version:

Thibaud Sarica, Alice Maison, Yelva Roustan, Matthias Ketzel, Steen Solvang Jensen, et al. Modelling concentration heterogeneities in streets using the street-network model MUNICH. Geoscientific Model Development, 2023, 16 (17), pp.5281-5303. 10.5194/gmd-16-5281-2023. hal-04457354

HAL Id: hal-04457354

https://hal.inrae.fr/hal-04457354

Submitted on 14 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Modelling concentration heterogeneities in streets using the street-network model MUNICH

Thibaud Sarica, Alice Maison, Yelva Roustan, Matthias Ketzel, Steen Solvang

Jensen, Youngseob Kim, Christophe Chaillou, Karine Sartelet

► To cite this version:

Thibaud Sarica, Alice Maison, Yelva Roustan, Matthias Ketzel, Steen Solvang Jensen, et al. Modelling concentration heterogeneities in streets using the street-network model MUNICH. Geoscientific Model Development, 2023, 16 (17), pp.5281-5303. 10.5194/gmd-16-5281-2023. hal-04457354

HAL Id: hal-04457354

https://hal.inrae.fr/hal-04457354v1

Submitted on 14 Feb 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Modelling concentration heterogeneities in streets using the street-network model MUNICH

Thibaud Sarica<fnref n="1" />, Alice Maison<fnref n="1" />,<fnref n="2" />, Yelva Roustan<fnref n="1" />, Matthias Ketzel<fnref n="3" />, Steen Solvang Jensen<fnref n="3" />, Youngseob Kim<fnref n="1" />, Christophe Chaillou<fnref n="4" />, and Karine Sartelet<fnref n="1" />

**Correspondence:** Thibaud Sarica (thibaud.sarica@enpc.fr) and Karine Sartelet (karine.sartelet@enpc.fr)

Received: 29 March 2023 – Discussion started: 17 April 2023

Revised: 18 July 2023 – Accepted: 16 August 2023 – Published: 14 September 2023

**Abstract.** Populations in urban areas are exposed to high local concentrations of pollutants, such as nitrogen dioxide and particulate matter, because of unfavourable dispersion conditions and the proximity to traffic. To simulate these concentrations over cities, models like the street-network model MUNICH (Model of Urban Network of Intersecting Canyons and Highways) rely on parameterizations to represent the air flow and the concentrations of pollutants in streets. In the current version, MUNICH v2.0, concentrations are assumed to be homogeneous in each street segment. A new version of MUNICH, where the street volume is discretized, is developed to represent the street gradients and to better estimate peoples' exposure. Three vertical levels are defined in each street segment. A horizontal discretization is also introduced under specific conditions by considering two zones with a parameterization taken from the Operational Street Pollution Model (OSPM). Simulations are performed over two districts of Copenhagen, Denmark, and one district of greater Paris, France. Results show an improvement in the comparison to observations, with higher concentrations at the bottom of the street, closer to traffic, of pollutants emitted by traffic (NOₓ, black carbon, organic matter). These increases reach up to 60 % for NO₂ and 30 % for PM₁₀ in comparison to MUNICH v2.0. The aspect ratio (ratio between building height and street width) influences the extent of the increase of the first-level concentrations compared to the average of the street. The increase is higher for wide streets (low aspect ratio and often higher traffic) by up to 53 % for NOₓ and 18 % for PM₁₀. Finally, a sensitivity analysis with regard to

the influence of the street network highlights the importance of using the model MUNICH with a network rather than with a single street.

## 1 Introduction

Pollution is estimated to have been responsible for approximately 9 million premature deaths in 2015 (Landrigan et al., 2018). This figure remains valid in 2019 despite an improvement in the types of pollution associated with extreme poverty (e.g. household air pollution and water pollution) (Fuller et al., 2022). This is partly due to an increase in the number of premature deaths attributable to ambient air pollution. The consequences of air pollution are particularly substantial in urban areas, where individuals are exposed to local high concentrations of air pollutants due to unfavourable dispersion conditions and proximity to traffic. As more than half of the world's population already lives in urban areas, with this being expected to rise to 68 % by 2050 (United Nations, 2019), it is crucial to estimate as accurately as possible the exposure of the population to atmospheric pollutant concentrations in urban areas. For many years, various modelling approaches have been developed to contribute to the understanding of the phenomena that drive the concentrations of pollutants in the atmosphere and to provide decision support tools (Collett and Oduyemi, 1997; Vardoulakis et al., 2003; El-Harbawi, 2013; Conti et al., 2017; Khan and Quamrul, 2021).

<footnote n="1">¹CEREA, École des Ponts ParisTech, EDF R&D, IPSL, Marne-la-Vallée, 77455, France</footnote>

<footnote n="2">²Université Paris-Saclay, INRAE, AgroParisTech, UMR EcoSys, Palaiseau, 91120, France</footnote>

<footnote n="3">³Department of Environmental Science, Aarhus University, Roskilde, Denmark</footnote>

<footnote n="4">⁴Aramco Fuel Research Center, Aramco Overseas Company, Rueil-Malmaison, 92500, France</footnote>
# Modelling concentration heterogeneities in streets using the street-network model MUNICH

Thibaud Sarica<fnref n="1" />, Alice Maison<fnref n="1" />,<fnref n="2" />, Yelva Roustan<fnref n="1" />, Matthias Ketzel<fnref n="3" />, Steen Solvang Jensen<fnref n="3" />, Youngseob Kim<fnref n="1" />, Christophe Chaillou<fnref n="4" />, and Karine Sartelet<fnref n="1" />

**Correspondence:** Thibaud Sarica (thibaud.sarica@enpc.fr) and Karine Sartelet (karine.sartelet@enpc.fr)

Received: 29 March 2023 – Discussion started: 17 April 2023

Revised: 18 July 2023 – Accepted: 16 August 2023 – Published: 14 September 2023

**Abstract.** Populations in urban areas are exposed to high local concentrations of pollutants, such as nitrogen dioxide and particulate matter, because of unfavourable dispersion conditions and the proximity to traffic. To simulate these concentrations over cities, models like the street-network model MUNICH (Model of Urban Network of Intersecting Canyons and Highways) rely on parameterizations to represent the air flow and the concentrations of pollutants in streets. In the current version, MUNICH v2.0, concentrations are assumed to be homogeneous in each street segment. A new version of MUNICH, where the street volume is discretized, is developed to represent the street gradients and to better estimate peoples' exposure. Three vertical levels are defined in each street segment. A horizontal discretization is also introduced under specific conditions by considering two zones with a parameterization taken from the Operational Street Pollution Model (OSPM). Simulations are performed over two districts of Copenhagen, Denmark, and one district of greater Paris, France. Results show an improvement in the comparison to observations, with higher concentrations at the bottom of the street, closer to traffic, of pollutants emitted by traffic (NOₓ, black carbon, organic matter). These increases reach up to 60 % for NO₂ and 30 % for PM₁₀ in comparison to MUNICH v2.0. The aspect ratio (ratio between building height and street width) influences the extent of the increase of the first-level concentrations compared to the average of the street. The increase is higher for wide streets (low aspect ratio and often higher traffic) by up to 53 % for NOₓ and 18 % for PM₁₀. Finally, a sensitivity analysis with regard to

the influence of the street network highlights the importance of using the model MUNICH with a network rather than with a single street.

## 1 Introduction

Pollution is estimated to have been responsible for approximately 9 million premature deaths in 2015 (Landrigan et al., 2018). This figure remains valid in 2019 despite an improvement in the types of pollution associated with extreme poverty (e.g. household air pollution and water pollution) (Fuller et al., 2022). This is partly due to an increase in the number of premature deaths attributable to ambient air pollution. The consequences of air pollution are particularly substantial in urban areas, where individuals are exposed to local high concentrations of air pollutants due to unfavourable dispersion conditions and proximity to traffic. As more than half of the world's population already lives in urban areas, with this being expected to rise to 68 % by 2050 (United Nations, 2019), it is crucial to estimate as accurately as possible the exposure of the population to atmospheric pollutant concentrations in urban areas. For many years, various modelling approaches have been developed to contribute to the understanding of the phenomena that drive the concentrations of pollutants in the atmosphere and to provide decision support tools (Collett and Oduyemi, 1997; Vardoulakis et al., 2003; El-Harbawi, 2013; Conti et al., 2017; Khan and Quamrul, 2021).

<footnote n="1">¹CEREA, École des Ponts ParisTech, EDF R&D, IPSL, Marne-la-Vallée, 77455, France</footnote>

<footnote n="2">²Université Paris-Saclay, INRAE, AgroParisTech, UMR EcoSys, Palaiseau, 91120, France</footnote>

<footnote n="3">³Department of Environmental Science, Aarhus University, Roskilde, Denmark</footnote>

<footnote n="4">⁴Aramco Fuel Research Center, Aramco Overseas Company, Rueil-Malmaison, 92500, France</footnote>
Regional-scale chemistry transport models, such as Po-lair3D (Mallet et al., 2007; Sartelet et al., 2018), CHIMERE (Menut et al., 2021; Falasca and Curci, 2018), and CMAQ (Wong et al., 2012; de la Paz et al., 2015), represent the urban background concentrations by solving the chemistry trans-port equation for spatial resolutions down to 1 km². How-ever, they cannot represent street concentrations, which are often higher than background concentrations for pollutants such as NO₂ and particles (Lugon et al., 2020). To represent these concentrations, local-scale models are thus developed with different approaches of variable complexity and compu-tational cost. Models based on computational fluid dynamics (CFD), such as code_saturen (Archambeau et al., 2004; Gao et al., 2018), OpenFoam (Lin et al., 2023), and PALM (Wolf et al., 2020), are able to represent the dispersion of pollutants and the physicochemical processes taking place in urban dis-tricts and streets with a fine spatial resolution by solving the Navier-Stokes equations and mass conservation equations for pollutants. However, they suffer from high computational cost as they use fine meshes to describe the morphology of buildings and streets. Other models use approaches that are less accurate than CFD but that run faster. They are typically based on a Gaussian or an Eulerian approach (Vardoulakis et al., 2003; Liang et al., 2023). Among these, the following can be mentioned: ADMS-Urban (McHugh et al., 1997; Hood et al., 2021), SIRANE (Soulhac et al., 2011, 2023), AERMOD (Cimorelli et al., 2004; Rood, 2014), EPISODE (Karl et al., 2019), and CALIOPE-Urban (Benavides et al., 2019). Either they consider each street independently of the others with exchanges between the street and the background concentra-tions above the street (Berkowicz, 2000b), or they consider a street network with incoming and/or outcoming flows be-tween streets at intersections (Soulhac et al., 2009; Kim et al., 2022). The Operational Street Pollution Model (OSPM) couples a Gaussian plume model for traffic emissions and a box model for the recirculation in the street (Berkowicz, 2000b). It is thus able to represent concentration heterogeneities in the street but cannot include complex chemistry. The Model of Urban Network of Intersection Canyons and Highways (MUNICH) uses solely an Eulerian box model approach (Lugon et al., 2020). As MUNICH is coupled with the SSH-aerosol model (Sartelet et al., 2020), the formation and age-ing of primary and secondary gas and particles in streets are represented. In the current version of MUNICH (v2.0) (Kim et al., 2022), concentrations are considered to be homoge-neous in each street segment. However, as shown by several on-site and modelling studies, the concentrations are very heterogeneous in streets with traffic emissions (Xie et al., 2003; Vardoulakis et al., 2011; Lateb et al., 2016; Sanchez et al., 2016; Amato et al., 2019; Lin et al., 2023). They are higher near the ground than at the top of the street, especially for primary pollutants.

A heterogeneous version of MUNICH is developed in this study, with the aim of representing the concentration hetero-geneities in the street while keeping the Eulerian approach of

MUNICH to retain its ability to accurately model chemistry and aerosol dynamics. The street volume is discretized verti-cally in three subvolumes. Traffic emissions are not instanta-neously diluted in the whole street volume, as in MUNICH v2.0, but only in the first subvolume, i.e. the one closest to the ground. With this discretization we do not aim to repro-duce finely the vertical profile of concentrations. The main objective is to improve the representation of concentrations close to the ground by avoiding the excessive dilution as-sociated with the homogeneity assumption. To represent the streets' horizontal heterogeneities, a recirculation zone in the shape of a trapeze is defined based on a parameterization of OSPM. It depends on the meteorological conditions and the street morphology, and it is applied under specific conditions in MUNICH that are described later.

A description of the differences between the homogeneous version of MUNICH (v2.0) and the new heterogeneous ver-sion is presented in Sect. 2. The applications to two street networks in Copenhagen, Denmark, with comparisons to ob-servations of NO₂ and CO, and to concentrations simulated by OSPM are discussed in Sect. 3. MUNICH is applied in Sect. 4 to the street network near Paris, France, used to val-idate MUNICH v2.0 (Kim et al., 2022), and the impacts of the heterogeneous version on concentrations of NO₂ and par-ticles are studied. To estimate the impact of the modelling hy-pothesis on the transport of pollutants between streets in the new version of MUNICH, a sensitivity analysis with regard to the presence of a street network is performed in Sect. 5.

## 2 Model description

The homogeneous version of MUNICH (v2.0) and the new heterogeneous version (see Fig. 1) are briefly described in this section, with a focus on the differences between the two versions. The heterogeneous version was developed from MUNICH v2.0. In the following, the homogeneous and heterogeneous versions of MUNICH are referred to as MUNICH-homo and MUNICH-hete, respectively. For a complete description of MUNICH, please refer to Kim et al. (2018, 2022) and Lugon et al. (2020, 2021a).

In order to solve the evolution equation of the street con-centrations C_str, a first-order operator splitting between trans-port and chemistry is performed:

$$ \frac{dC_{\text{str}}}{dt} = \left. \frac{dC_{\text{str}}}{dt} \right|_{\text{tr}} + \left. \frac{dC_{\text{str}}}{dt} \right|_{\text{ch}} \quad (1) $$

The transport term includes advection from one street to another and vertical transport between the street and the background above it. The chemistry includes gas-phase chemistry, as well as aerosol dynamics (coagulation, con-densation, and evaporation). As deposition and resuspension processes have minor effects compared to transport and chemistry (Lugon et al., 2021b; Kim et al., 2022), they are omitted in the rest of this study.
Regional-scale chemistry transport models, such as Po-lair3D (Mallet et al., 2007; Sartelet et al., 2018), CHIMERE (Menut et al., 2021; Falasca and Curci, 2018), and CMAQ (Wong et al., 2012; de la Paz et al., 2015), represent the urban background concentrations by solving the chemistry trans-port equation for spatial resolutions down to 1 km². How-ever, they cannot represent street concentrations, which are often higher than background concentrations for pollutants such as NO₂ and particles (Lugon et al., 2020). To represent these concentrations, local-scale models are thus developed with different approaches of variable complexity and compu-tational cost. Models based on computational fluid dynamics (CFD), such as code_saturen (Archambeau et al., 2004; Gao et al., 2018), OpenFoam (Lin et al., 2023), and PALM (Wolf et al., 2020), are able to represent the dispersion of pollutants and the physicochemical processes taking place in urban dis-tricts and streets with a fine spatial resolution by solving the Navier-Stokes equations and mass conservation equations for pollutants. However, they suffer from high computational cost as they use fine meshes to describe the morphology of buildings and streets. Other models use approaches that are less accurate than CFD but that run faster. They are typically based on a Gaussian or an Eulerian approach (Vardoulakis et al., 2003; Liang et al., 2023). Among these, the following can be mentioned: ADMS-Urban (McHugh et al., 1997; Hood et al., 2021), SIRANE (Soulhac et al., 2011, 2023), AERMOD (Cimorelli et al., 2004; Rood, 2014), EPISODE (Karl et al., 2019), and CALIOPE-Urban (Benavides et al., 2019). Either they consider each street independently of the others with exchanges between the street and the background concentra-tions above the street (Berkowicz, 2000b), or they consider a street network with incoming and/or outcoming flows be-tween streets at intersections (Soulhac et al., 2009; Kim et al., 2022). The Operational Street Pollution Model (OSPM) couples a Gaussian plume model for traffic emissions and a box model for the recirculation in the street (Berkowicz, 2000b). It is thus able to represent concentration heterogeneities in the street but cannot include complex chemistry. The Model of Urban Network of Intersection Canyons and Highways (MUNICH) uses solely an Eulerian box model approach (Lugon et al., 2020). As MUNICH is coupled with the SSH-aerosol model (Sartelet et al., 2020), the formation and age-ing of primary and secondary gas and particles in streets are represented. In the current version of MUNICH (v2.0) (Kim et al., 2022), concentrations are considered to be homoge-neous in each street segment. However, as shown by several on-site and modelling studies, the concentrations are very heterogeneous in streets with traffic emissions (Xie et al., 2003; Vardoulakis et al., 2011; Lateb et al., 2016; Sanchez et al., 2016; Amato et al., 2019; Lin et al., 2023). They are higher near the ground than at the top of the street, especially for primary pollutants.

A heterogeneous version of MUNICH is developed in this study, with the aim of representing the concentration hetero-geneities in the street while keeping the Eulerian approach of

MUNICH to retain its ability to accurately model chemistry and aerosol dynamics. The street volume is discretized verti-cally in three subvolumes. Traffic emissions are not instanta-neously diluted in the whole street volume, as in MUNICH v2.0, but only in the first subvolume, i.e. the one closest to the ground. With this discretization we do not aim to repro-duce finely the vertical profile of concentrations. The main objective is to improve the representation of concentrations close to the ground by avoiding the excessive dilution as-sociated with the homogeneity assumption. To represent the streets' horizontal heterogeneities, a recirculation zone in the shape of a trapeze is defined based on a parameterization of OSPM. It depends on the meteorological conditions and the street morphology, and it is applied under specific conditions in MUNICH that are described later.

A description of the differences between the homogeneous version of MUNICH (v2.0) and the new heterogeneous ver-sion is presented in Sect. 2. The applications to two street networks in Copenhagen, Denmark, with comparisons to ob-servations of NO₂ and CO, and to concentrations simulated by OSPM are discussed in Sect. 3. MUNICH is applied in Sect. 4 to the street network near Paris, France, used to val-idate MUNICH v2.0 (Kim et al., 2022), and the impacts of the heterogeneous version on concentrations of NO₂ and par-ticles are studied. To estimate the impact of the modelling hy-pothesis on the transport of pollutants between streets in the new version of MUNICH, a sensitivity analysis with regard to the presence of a street network is performed in Sect. 5.

## 2 Model description

The homogeneous version of MUNICH (v2.0) and the new heterogeneous version (see Fig. 1) are briefly described in this section, with a focus on the differences between the two versions. The heterogeneous version was developed from MUNICH v2.0. In the following, the homogeneous and heterogeneous versions of MUNICH are referred to as MUNICH-homo and MUNICH-hete, respectively. For a complete description of MUNICH, please refer to Kim et al. (2018, 2022) and Lugon et al. (2020, 2021a).

In order to solve the evolution equation of the street con-centrations C_str, a first-order operator splitting between trans-port and chemistry is performed:

$$ \frac{dC_{str}}{dt} = \left. \frac{dC_{str}}{dt} \right|_{tr} + \left. \frac{dC_{str}}{dt} \right|_{ch} \qquad (1) $$

The transport term includes advection from one street to another and vertical transport between the street and the background above it. The chemistry includes gas-phase chemistry, as well as aerosol dynamics (coagulation, con-densation, and evaporation). As deposition and resuspension processes have minor effects compared to transport and chemistry (Lugon et al., 2021b; Kim et al., 2022), they are omitted in the rest of this study.
<figure><img src="image_4.png" /><figcaption>Figure 1. Representation of the processes in the homogeneous version (a) and the heterogeneous version (b) of MUNICH. The dotted red lines represent the street volume. In (a), the street canopy is represented by a single volume, whereas in (b), it is divided into three subvolumes delimited by the dotted blue lines. The rose arrows represent the traffic emissions (including brake, tyre, and road wear), the light-green arrows represent the fluxes entering the street via the upwind intersection, and the dark-green arrows represent the fluxes leaving the street via the downwind intersection. The yellow arrows symbolize the vertical turbulent exchanges with the background, and in (b), the orange arrows symbolize the vertical exchanges among the subvolumes.</figcaption></figure>

## 2.1 Homogeneous approach

Using a box model approach, the concentrations are assumed to be homogeneous in the whole street volume, and the effects of the processes on the concentrations are represented by the following equation:

$$ \frac{dC_{\text{str}}}{dt} \bigg|_{\text{tr}} = \frac{1}{V} (Q_{\text{em}} + Q_{\text{inflow}} + Q_{\text{outflow}} + Q_{\text{vert}}), \quad (2) $$

where V is the volume of the rectangular cuboid street, $Q_{\text{em}}$ is the traffic emission flux, $Q_{\text{inflow}}$ is the flux entering the street via the upwind intersection, $Q_{\text{outflow}}$ is the flux leaving the street via the downwind intersection, and $Q_{\text{vert}}$ is the vertical turbulent flux between the background and the street.

The street volume is defined as $V = HWL$, with H being the average building height over the street segment, W being the mean street width, and L being its length. MUNICH considers it to be that buildings on each side of the street are continuous, thus not representing inflow and/or outflow that could be induced by gaps in the buildings. The inflow term $Q_{\text{inflow}}$ is obtained from the computation of the fluxes at the upwind intersection (Kim et al., 2018; Soulhac et al., 2009) (see Sect. 5 for a brief description). The outflow term $Q_{\text{outflow}}$ is expressed as follows:

$$ Q_{\text{outflow}} = H W u_{\text{str}} C_{\text{str}}, \quad (3) $$

with $u_{\text{str}}$ being the mean horizontal wind speed in the street.

The vertical turbulent flux $Q_{\text{vert}}$ between the street and the overlying atmosphere is

$$ Q_{\text{vert}} = q_{\text{vert}} W L \frac{(C_{\text{str}} - C_{\text{bkgd}})}{H}, \quad (4) $$

with $q_{\text{vert}}$ being the vertical transfer coefficient and $C_{\text{bkgd}}$ being the background concentration.

The vertical transfer coefficient and the horizontal wind speed are the key parameters representing the dispersion of concentrations. As their formulation differs between the homogeneous and the heterogeneous versions of MUNICH, they are now detailed.

### 2.1.1 Vertical transfer coefficient for turbulent flux

Three parameterizations are implemented in MUNICH to determine the vertical transfer coefficient $q_{\text{vert}}$ between the street and the overlying atmosphere (Maison et al., 2022). However, currently, only the parameterization adapted from Wang (2014) is designed to provide vertical profiles for both wind speed and mixing length within the street. We therefore limit our analysis to the latter.

In MUNICH-homo, the vertical transfer coefficient at the roof level is expressed as follows:

$$ q_{\text{vert}} = \sigma_{\text{w}} l_{\text{m}}(z = H) \quad \text{with} \quad l_{\text{m}}(z) = \frac{\kappa z l_{\text{c}}}{l_{\text{c}} + \kappa z}, \quad (5) $$

where $\sigma_{\text{w}}$ is the standard deviation of the vertical wind velocity at roof level, and $l_{\text{m}}$ is the mixing length defined as a harmonic mean between two length scales (Coceal and Belcher, 2004), namely (i) $\kappa z$, with the Von Kármán constant ($\kappa = 0.42$), and (ii) $l_{\text{c}}$, which is a characteristic length of the street chosen to be equal to 0.5W (Maison et al., 2022).

### 2.1.2 Mean horizontal wind speed in the street

As for the calculation of the vertical transfer coefficient, three parameterizations are proposed in MUNICH to determine the mean horizontal wind speed in the street (Maison et al., 2022). In the Wang parameterization, the mean horizontal
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Representation of the processes in the homogeneous version (a) and the heterogeneous version (b) of MUNICH. The dotted red lines represent the street volume. In (a), the street canopy is represented by a single volume, whereas in (b), it is divided into three subvolumes delimited by the dotted blue lines. The rose arrows represent the traffic emissions (including brake, tyre, and road wear), the light-green arrows represent the fluxes entering the street via the upwind intersection, and the dark-green arrows represent the fluxes leaving the street via the downwind intersection. The yellow arrows symbolize the vertical turbulent exchanges with the background, and in (b), the orange arrows symbolize the vertical exchanges among the subvolumes.</figcaption></figure>

## 2.1 Homogeneous approach

Using a box model approach, the concentrations are assumed to be homogeneous in the whole street volume, and the effects of the processes on the concentrations are represented by the following equation:

$$ \frac{dC_{\text{str}}}{dt} \bigg|_{\text{tr}} = \frac{1}{V} (Q_{\text{em}} + Q_{\text{inflow}} + Q_{\text{outflow}} + Q_{\text{vert}}), \quad (2) $$

where V is the volume of the rectangular cuboid street, $Q_{em}$ is the traffic emission flux, $Q_{inflow}$ is the flux entering the street via the upwind intersection, $Q_{outflow}$ is the flux leaving the street via the downwind intersection, and $Q_{vert}$ is the vertical turbulent flux between the background and the street.

The street volume is defined as $V = HWL$, with H being the average building height over the street segment, W being the mean street width, and L being its length. MUNICH considers it to be that buildings on each side of the street are continuous, thus not representing inflow and/or outflow that could be induced by gaps in the buildings. The inflow term $Q_{inflow}$ is obtained from the computation of the fluxes at the upwind intersection (Kim et al., 2018; Soulhac et al., 2009) (see Sect. 5 for a brief description). The outflow term $Q_{outflow}$ is expressed as follows:

$$ Q_{\text{outflow}} = H W u_{\text{str}} C_{\text{str}}, \quad (3) $$

with $u_{str}$ being the mean horizontal wind speed in the street.

The vertical turbulent flux $Q_{vert}$ between the street and the overlying atmosphere is

$$ Q_{\text{vert}} = q_{\text{vert}} W L \frac{(C_{\text{str}} - C_{\text{bkgd}})}{H}, \quad (4) $$

with $q_{vert}$ being the vertical transfer coefficient and $C_{bkgd}$ being the background concentration.

The vertical transfer coefficient and the horizontal wind speed are the key parameters representing the dispersion of concentrations. As their formulation differs between the homogeneous and the heterogeneous versions of MUNICH, they are now detailed.

### 2.1.1 Vertical transfer coefficient for turbulent flux

Three parameterizations are implemented in MUNICH to determine the vertical transfer coefficient $q_{vert}$ between the street and the overlying atmosphere (Maison et al., 2022). However, currently, only the parameterization adapted from Wang (2014) is designed to provide vertical profiles for both wind speed and mixing length within the street. We therefore limit our analysis to the latter.

In MUNICH-homo, the vertical transfer coefficient at the roof level is expressed as follows:

$$ q_{\text{vert}} = \sigma_w l_m(z = H) \quad \text{with} \quad l_m(z) = \frac{\kappa z l_c}{l_c + \kappa z}, \quad (5) $$

where $\sigma_w$ is the standard deviation of the vertical wind velocity at roof level, and $l_m$ is the mixing length defined as a harmonic mean between two length scales (Coceal and Belcher, 2004), namely (i) $\kappa z$, with the Von Kármán constant ($\kappa = 0.42$), and (ii) $l_c$, which is a characteristic length of the street chosen to be equal to 0.5W (Maison et al., 2022).

### 2.1.2 Mean horizontal wind speed in the street

As for the calculation of the vertical transfer coefficient, three parameterizations are proposed in MUNICH to determine the mean horizontal wind speed in the street (Maison et al., 2022). In the Wang parameterization, the mean horizontal
wind speed in the street $u_{str}$ is equal to

$$ 
\begin{aligned}
u_{str} &= \frac{1}{H - z_{0s}} \int_{z_{0s}}^{H} u_{Y}(z) dz \\
&= \frac{u_{H}|\cos(\varphi)|}{(H - z_{0s})} \int_{z_{0s}}^{H} [J_{1}I_{0}(g(z)) + J_{2}K_{0}(g(z))] dz,
\end{aligned} 
\quad (6) $$

with $u_H|\cos(\varphi)|$ being the wind speed at roof level in the direction of the street and $z_{0s}$ being the wall and ground roughness length in the street (fixed to 0.01 m). $I_0$ and $K_0$ are the first- and second-kind modified Bessel functions of order 0. $J_1$ and $J_2$ are integration coefficients equal to

$$ 
\begin{aligned}
J_1 &= \frac{1}{I_0(g(H)) - I_0(g(z_{0s}))K_0(g(H))/K_0(g(z_{0s}))} \quad \text{and} \\
J_2 &= -\frac{J_1 I_0(g(z_{0s}))}{K_0(g(z_{0s}))}.
\end{aligned} 
\quad (7) $$

The function $g(z)$ is calculated as follows:

$$ g(z) = 2\sqrt{C_B a_r \frac{z}{l_m(H)}}, \quad (8) $$

with $a_r = H/W$ being the aspect ratio of the street and $C_B$ being a coefficient that is dependent on the wind angle with the street and the aspect ratio (Maison et al., 2022).

## 2.2 Heterogeneous approach

In MUNICH-hete, the street is divided into three vertical levels to limit the artificial dilution of the traffic emissions and the concentrations in the whole street volume (see Figs. 1b and 2). Levels are ordered from the ground to the top of the street. The first level ($i=1$) contains the traffic emissions. The thickness $h_1$ is taken as 2 m, which corresponds to a zone where the traffic producing turbulence mixes and dilutes traffic emissions (Solazzo et al., 2008). This traffic-induced turbulence is not explicitly considered in the model. The second-level ($i=2$) thickness $h_2$ is also of 2 m. It acts as a buffer zone between the first level, where traffic emissions are, and the third level, where exchanges with the background take place. Starting at 4 m, the third level ($i=3$) goes to the roof level ($h_3 = (H-4)$ m). The minimum street height considered in the model is set at 6 m. The three levels of the heterogeneous version are referred to as MUNICH-hete-11, MUNICH-hete-12, and MUNICH-hete-13. Each level $i$ is thus associated with a specific volume $V_i$, and the evolution equations may be written as follows:

$$ \frac{dC_{str}^i}{dt} \bigg|_{tr} = \frac{1}{V_i} (Q_{em}^i + Q_{inflow}^i + Q_{outflow}^i + Q_{vert}^{i,i+1} + Q_{vert}^{i-1,i}), \quad (9) $$

<figure><img src="image_5.png" /><figcaption>Figure 2. Schematic of the discretization in MUNICH-hete. The red triangle represents the ventilation zone that is present under specific conditions, and the white trapeze represents the recirculation zone. $z_i^m$ is the middle height of level $i$.</figcaption></figure>

where $Q_{em}^i$ is the traffic emission flux (only in the first level ($i=1$)), $Q_{inflow}^i$ is the flux entering the level via the upwind intersection, $Q_{outflow}^i$ is the flux leaving the level via the downwind intersection, $Q_{vert}^{i,i+1}$ is the vertical turbulent flux between the levels $i$ and $i+1$ (for $i=3$, it is exchanged with the background), and $Q_{vert}^{i-1,i}$ is the vertical turbulent flux between the levels $i-1$ and $i$ (equal to zero if $i=1$).

Note that more than three vertical levels could be defined in MUNICH-hete as the vertical variations within the streets of winds and mixing lengths are parameterized when discretizing the streets. However, the first vertical level at the bottom of the street should not be too thin because of mixing due to traffic turbulence. A minimum height for the first layer of 1.5 m seems reasonable.

In OSPM, the flow that develops into a vortex in the street between buildings is represented by a recirculation zone. It occupies the whole street volume for narrow streets, and it has the shape of a trapeze for wider streets (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). When the recirculation zone does not occupy the whole street volume, there is a ventilation zone (see Fig. 2) where concentrations of pollutants emitted by traffic are usually lower than in the recirculation zone. In MUNICH-homo, the recirculation zone is not explicit, whatever the street ratio (H/W) is. In MUNICH-hete, the volume of the recirculation zone is computed as in OSPM, as detailed in Appendix A. For now, concentrations in the ventilation zone are considered to be homogeneous and equal to the background concentrations, i.e. concentrations above the street. The ventilation zone is thus taken into account when it is not affected by traffic emissions. In practice, this means that the width of the trapeze base can only be equal to or larger than the width $W$ of the street (see Appendix A). The width $W_i$ of the level $i$ can thus be inferior to the width $W$ of the street, reducing the level volume (see Fig. 2). Appendix A presents the algorithm implemented in
wind speed in the street $u_{str}$ is equal to

$$ 
\begin{aligned}
u_{str} &= \frac{1}{H - z_{0s}} \int_{z_{0s}}^{H} u_{Y}(z) dz \\
&= \frac{u_{H}|\cos(\varphi)|}{(H - z_{0s})} \int_{z_{0s}}^{H} [J_{1}I_{0}(g(z)) + J_{2}K_{0}(g(z))] dz,
\end{aligned} 
\quad (6) $$

with $u_H|\cos(\varphi)|$ being the wind speed at roof level in the direction of the street and $z_{0s}$ being the wall and ground roughness length in the street (fixed to 0.01 m). $I_0$ and $K_0$ are the first- and second-kind modified Bessel functions of order 0. $J_1$ and $J_2$ are integration coefficients equal to

$$ 
\begin{aligned}
J_1 &= \frac{1}{I_0(g(H)) - I_0(g(z_{0s}))K_0(g(H))/K_0(g(z_{0s}))} \quad \text{and} \\
J_2 &= -\frac{J_1 I_0(g(z_{0s}))}{K_0(g(z_{0s}))}.
\end{aligned} 
\quad (7) $$

The function $g(z)$ is calculated as follows:

$$ g(z) = 2\sqrt{C_B a_r \frac{z}{l_m(H)}}, \quad (8) $$

with $a_r = H/W$ being the aspect ratio of the street and $C_B$ being a coefficient that is dependent on the wind angle with the street and the aspect ratio (Maison et al., 2022).

## 2.2 Heterogeneous approach

In MUNICH-hete, the street is divided into three vertical levels to limit the artificial dilution of the traffic emissions and the concentrations in the whole street volume (see Figs. 1b and 2). Levels are ordered from the ground to the top of the street. The first level ($i = 1$) contains the traffic emissions. The thickness $h_1$ is taken as 2 m, which corresponds to a zone where the traffic producing turbulence mixes and dilutes traffic emissions (Solazzo et al., 2008). This traffic-induced turbulence is not explicitly considered in the model. The second-level ($i = 2$) thickness $h_2$ is also of 2 m. It acts as a buffer zone between the first level, where traffic emissions are, and the third level, where exchanges with the background take place. Starting at 4 m, the third level ($i = 3$) goes to the roof level ($h_3 = (H - 4)$ m). The minimum street height considered in the model is set at 6 m. The three levels of the heterogeneous version are referred to as MUNICH-hete-11, MUNICH-hete-12, and MUNICH-hete-13. Each level $i$ is thus associated with a specific volume $V_i$, and the evolution equations may be written as follows:

$$ \frac{dC_{str}^i}{dt}|_{tr} = \frac{1}{V_i}(Q_{em}^i + Q_{inflow}^i + Q_{outflow}^i + Q_{vert}^{i,i+1} + Q_{vert}^{i-1,i}), \quad (9) $$

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Schematic of the discretization in MUNICH-hete. The red triangle represents the ventilation zone that is present under specific conditions, and the white trapeze represents the recirculation zone. $z_i^m$ is the middle height of level $i$.</figcaption></figure>

where $Q_{em}^i$ is the traffic emission flux (only in the first level ($i = 1$)), $Q_{inflow}^i$ is the flux entering the level via the upwind intersection, $Q_{outflow}^i$ is the flux leaving the level via the downwind intersection, $Q_{vert}^{i,i+1}$ is the vertical turbulent flux between the levels $i$ and $i+1$ (for $i = 3$, it is exchanged with the background), and $Q_{vert}^{i-1,i}$ is the vertical turbulent flux between the levels $i-1$ and $i$ (equal to zero if $i = 1$).

Note that more than three vertical levels could be defined in MUNICH-hete as the vertical variations within the streets of winds and mixing lengths are parameterized when discretizing the streets. However, the first vertical level at the bottom of the street should not be too thin because of mixing due to traffic turbulence. A minimum height for the first layer of 1.5 m seems reasonable.

In OSPM, the flow that develops into a vortex in the street between buildings is represented by a recirculation zone. It occupies the whole street volume for narrow streets, and it has the shape of a trapeze for wider streets (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). When the recirculation zone does not occupy the whole street volume, there is a ventilation zone (see Fig. 2) where concentrations of pollutants emitted by traffic are usually lower than in the recirculation zone. In MUNICH-homo, the recirculation zone is not explicit, whatever the street ratio (H/W) is. In MUNICH-hete, the volume of the recirculation zone is computed as in OSPM, as detailed in Appendix A. For now, concentrations in the ventilation zone are considered to be homogeneous and equal to the background concentrations, i.e. concentrations above the street. The ventilation zone is thus taken into account when it is not affected by traffic emissions. In practice, this means that the width of the trapeze base can only be equal to or larger than the width $W$ of the street (see Appendix A). The width $W_i$ of the level $i$ can thus be inferior to the width $W$ of the street, reducing the level volume (see Fig. 2). Appendix A presents the algorithm implemented in
MUNICH-hete to consider the volume reduction of the ventilation zone. Further work is needed to differentiate the two zones for cases where the ventilation zone develops further into the street.

To quantify mass transfer through the intersections, the fluxes are assumed to be vertically homogeneous and remain determined as proposed by Soulhac et al. (2009) as the shapes of the intersections may differ from one to another, and turbulence is not quantified. The flux $Q_{inflow}$ entering the street is assumed to be the same for each vertical level. The flux $Q_{outflow}$ leaving a street is a surface-weighted average of the fluxes leaving each vertical level of the street:

$$Q_{\text{outflow}} = \sum_{i=1}^{i=3} Q_{\text{outflow}}^{i} = \sum_{i=1}^{i=3} S_{i}^{v} u_{i} C_{\text{str}}^{i}, \quad (10)$$

with $S_i^v$ being the vertical surface of the level i, as presented in Appendix A, and $u_i$ being the mean horizontal wind speed of the level i (see Sect. 2.2.2).

## 2.2.1 Vertical turbulent fluxes

To compute the vertical turbulent transfer $Q_{\text{vert}}^{i,i+1}$ at the interface of the vertical levels i and i + 1, Eq. (4) is modified to represent the vertical exchanges between the levels i and i + 1:

$$Q_{\text{vert}}^{i,i+1} = q_{\text{vert}}^{i} W_{i} L \frac{(C_{i} - C_{i+1})}{\Delta z_{i}^{m}}, \quad (11)$$

where $C_i$ and $C_{i+1}$ are the concentrations in the levels i and i + 1 respectively; $W_i$ is the width of the level i that is inferior or equal to the width W of the street; and $\Delta z_i^m$ is the difference in altitude between the middles of the levels i + 1 and i, which are as noted $z_{i+1}^m$, $z_i^m$, and $z_3^m$ (see Fig. 2). For i = 3, i.e. the highest level, $z_{i+1}^m$ is taken as the symmetrical of $z_3^m$ to the roof level and is noted as $z_{\text{bkgd}}^m$. It gives an approximation of the volume that effectively exchanges with the third level.

Among the three parameterizations implemented in MUNICH to determine the vertical transfer coefficient $q_{\text{vert}}$, only the Wang parameterization is adapted to the discretization. It is thanks to its explicit vertical dependency and validity for a wide range of street canyon and wind characteristics (Maison et al., 2022). To compute the vertical turbulent transfer $Q_{\text{vert}}^{i,i+1}$, the vertical transfer coefficient is taken at the height of the interface between the two vertical levels considered, and Eq. (5) is now written as follows:

$$q_{\text{vert}}^{i} = \sigma_{w} \kappa z_{i} \frac{l_{c}}{l_{c} + \kappa z_{i}}, \quad (12)$$

where $z_i$ is the height of the interface between levels i and i + 1. The influence of the atmospheric stability on vertical mixing is taken into account by modifying the standard deviation of the vertical wind velocity at roof level and thus the vertical transfer rate, depending on the length of Monin-Obukhov, as in MUNICH-homo.

By combining Eqs. (11) and (12), the vertical turbulent transfer $Q_{i,i+1}$ can be written for each level as follows:

$$Q_{\text{vert}}^{1,2} = \sigma_{w} \kappa z_{1} \frac{l_{c}}{l_{c} + \kappa z_{1}} W_{i} L \frac{C_{1} - C_{2}}{z_{2}^{m} - z_{1}^{m}}, \quad (13)$$

$$Q_{\text{vert}}^{2,3} = \sigma_{w} \kappa z_{2} \frac{l_{c}}{l_{c} + \kappa z_{2}} W_{i} L \frac{C_{2} - C_{3}}{z_{3}^{m} - z_{2}^{m}}, \quad (14)$$

$$Q_{\text{vert}}^{3, \text{bkgd}} = \sigma_{w} \kappa H \frac{l_{c}}{l_{c} + \kappa H} W_{i} L \frac{C_{3} - C_{\text{bkgd}}}{z_{\text{bkgd}}^{m} - z_{3}^{m}}, \quad (15)$$

with $C_1$, $C_2$, $C_3$, and $C_{\text{bkgd}}$ being the concentrations of the three levels and the background respectively.

## 2.2.2 Mean horizontal wind speed in the street

As for the vertical turbulent flux, the Wang parameterization is preferred among the three parameterizations available in MUNICH to compute the mean horizontal wind speeds in the street. This is thanks to its explicit vertical dependency and the no-slip condition at the ground that is always satisfied ($u(z=0) = 0$) (Maison et al., 2022). Therefore, the mean horizontal wind speed can be computed at each level in the street by modifying Eq. (6) to integrate vertically between the level heights:

$$u_1 = \frac{1}{z_1 - z_{0s}} \int_{z_{0s}}^{z_1} u_Y(z) dz, \quad (16)$$

$$u_2 = \frac{1}{z_2 - z_1} \int_{z_1}^{z_2} u_Y(z) dz, \quad (17)$$

$$u_3 = \frac{1}{H - z_2} \int_{z_2}^{H} u_Y(z) dz, \quad (18)$$

with $z_1$ and $z_2$ being the limits of the first two levels, as presented in Fig. 2.

# 3 Application to street networks in Copenhagen with comparison to OSPM

This section presents two applications of MUNICH-hete to assess its capabilities compared to MUNICH-homo and OSPM. Simulations are performed over the year 2019 to generate hourly concentrations for two street networks in Copenhagen, Denmark. The first street network is centred around H. C. Andersens Boulevard, and the second is centred around Jagtvej street. They are respectively named HCAB and JGTV in the following. These street networks have been selected as observational data of CO, NO₂, NOₓ, and O₃ are available for HCAB, and observational data of NO₂ and NOₓ are available for JGTV, both over the whole year. OSPM simulations were also performed to compare model performances. As OSPM does not represent air fluxes at intersections, OSPM
MUNICH-hete to consider the volume reduction of the ventilation zone. Further work is needed to differentiate the two zones for cases where the ventilation zone develops further into the street.

To quantify mass transfer through the intersections, the fluxes are assumed to be vertically homogeneous and remain determined as proposed by Soulhac et al. (2009) as the shapes of the intersections may differ from one to another, and turbulence is not quantified. The flux $Q_{inflow}$ entering the street is assumed to be the same for each vertical level. The flux $Q_{outflow}$ leaving a street is a surface-weighted average of the fluxes leaving each vertical level of the street:

$$Q_{\text{outflow}} = \sum_{i=1}^{i=3} Q_{\text{outflow}}^{i} = \sum_{i=1}^{i=3} S_{i}^{v} u_{i} C_{\text{str}}^{i}, \qquad (10)$$

with $S_i^v$ being the vertical surface of the level $i$, as presented in Appendix A, and $u_i$ being the mean horizontal wind speed of the level $i$ (see Sect. 2.2.2).

### 2.2.1 Vertical turbulent fluxes

To compute the vertical turbulent transfer $Q_{\text{vert}}^{i,i+1}$ at the interface of the vertical levels $i$ and $i+1$, Eq. (4) is modified to represent the vertical exchanges between the levels $i$ and $i+1$:

$$Q_{\text{vert}}^{i,i+1} = q_{\text{vert}}^{i} W_{i} L \frac{(C_{i} - C_{i+1})}{\Delta z_{i}^{m}}, \qquad (11)$$

where $C_i$ and $C_{i+1}$ are the concentrations in the levels $i$ and $i+1$ respectively; $W_i$ is the width of the level $i$ that is inferior or equal to the width $W$ of the street; and $\Delta z_i^m$ is the difference in altitude between the middles of the levels $i+1$ and $i$, which are as noted $z_{i+1}^m$, $z_2^m$, and $z_3^m$ (see Fig. 2). For $i=3$, i.e. the highest level, $z_{i+1}^m$ is taken as the symmetrical of $z_3^m$ to the roof level and is noted as $z_{\text{bkgd}}^m$. It gives an approximation of the volume that effectively exchanges with the third level.

Among the three parameterizations implemented in MUNICH to determine the vertical transfer coefficient $q_{\text{vert}}$, only the Wang parameterization is adapted to the discretization. It is thanks to its explicit vertical dependency and validity for a wide range of street canyon and wind characteristics (Maison et al., 2022). To compute the vertical turbulent transfer $Q_{\text{vert}}^{i,i+1}$, the vertical transfer coefficient is taken at the height of the interface between the two vertical levels considered, and Eq. (5) is now written as follows:

$$q_{\text{vert}}^{i} = \sigma_{w} \kappa z_{i} \frac{l_{c}}{l_{c} + \kappa z_{i}}, \qquad (12)$$

where $z_i$ is the height of the interface between levels $i$ and $i+1$. The influence of the atmospheric stability on vertical mixing is taken into account by modifying the standard deviation of the vertical wind velocity at roof level and thus the vertical transfer rate, depending on the length of Monin-Obukhov, as in MUNICH-homo.

By combining Eqs. (11) and (12), the vertical turbulent transfer $Q_{i,i+1}$ can be written for each level as follows:

$$Q_{\text{vert}}^{1,2} = \sigma_{w} \kappa z_{1} \frac{l_{c}}{l_{c} + \kappa z_{1}} W_{i} L \frac{C_{1} - C_{2}}{z_{2}^{m} - z_{1}^{m}}, \qquad (13)$$

$$Q_{\text{vert}}^{2,3} = \sigma_{w} \kappa z_{2} \frac{l_{c}}{l_{c} + \kappa z_{2}} W_{i} L \frac{C_{2} - C_{3}}{z_{3}^{m} - z_{2}^{m}}, \qquad (14)$$

$$Q_{\text{vert}}^{3,\text{bkgd}} = \sigma_{w} \kappa H \frac{l_{c}}{l_{c} + \kappa H} W_{i} L \frac{C_{3} - C_{\text{bkgd}}}{z_{\text{bkgd}}^{m} - z_{3}^{m}}, \qquad (15)$$

with $C_1$, $C_2$, $C_3$, and $C_{\text{bkgd}}$ being the concentrations of the three levels and the background respectively.

### 2.2.2 Mean horizontal wind speed in the street

As for the vertical turbulent flux, the Wang parameterization is preferred among the three parameterizations available in MUNICH to compute the mean horizontal wind speeds in the street. This is thanks to its explicit vertical dependency and the no-slip condition at the ground that is always satisfied ($u(z=0) = 0$) (Maison et al., 2022). Therefore, the mean horizontal wind speed can be computed at each level in the street by modifying Eq. (6) to integrate vertically between the level heights:

$$u_1 = \frac{1}{z_1 - z_{0s}} \int_{z_{0s}}^{z_1} u_Y(z) dz, \qquad (16)$$

$$u_2 = \frac{1}{z_2 - z_1} \int_{z_1}^{z_2} u_Y(z) dz, \qquad (17)$$

$$u_3 = \frac{1}{H - z_2} \int_{z_2}^{H} u_Y(z) dz, \qquad (18)$$

with $z_1$ and $z_2$ being the limits of the first two levels, as presented in Fig. 2.

## 3 Application to street networks in Copenhagen with comparison to OSPM

This section presents two applications of MUNICH-hete to assess its capabilities compared to MUNICH-homo and OSPM. Simulations are performed over the year 2019 to generate hourly concentrations for two street networks in Copenhagen, Denmark. The first street network is centred around H. C. Andersens Boulevard, and the second is centred around Jagtvej street. They are respectively named HCAB and JGTV in the following. These street networks have been selected as observational data of CO, NO₂, NOₓ, and O₃ are available for HCAB, and observational data of NO₂ and NOₓ are available for JGTV, both over the whole year. OSPM simulations were also performed to compare model performances. As OSPM does not represent air fluxes at intersections, OSPM
simulations are performed only for the streets where there are observations. Thanks to its coupled approach of a Gaussian plume model and a box model, OSPM is able to calculate concentrations on the two sides of the street (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). Two receptors are used to compare with the observed concentrations on each side of the street: OSPM-R1 and OSPM-R2. The height of the receptors is taken as 2 m to correspond roughly to the height at which observations were performed. They are compared to the concentrations simulated at the first two levels of MUNICH-hete, which are representative of the concentrations at 1 and 3 m.

MUNICH and OSPM use the same input data to estimate the street concentrations. The meteorological parameters originate from simulations performed with the Weather Research and Forecasting model (WRF, Skamarock et al., 2008). The background concentrations are simulated using the Urban Background Model (UBM, Berkowicz, 2000a), which is “a multiple source model that applies a Gaussian approach for horizontal dispersion and a linear approach for vertical dispersion up to the boundary layer” (Jensen et al., 2016). Traffic emissions are generated using the procedure implemented in the local-scale Gaussian air pollution model OML-Highway (Olesen et al., 2015), allowing for precise information for each street segment. Traffic data (average daily traffic, travel speed, and share of heavy-duty vehicles) are used to generate emissions by use of the European emission model COPERT IV. Traffic emissions include exhaust emissions of gases and particles and non-exhaust emissions of particles. Non-exhaust emissions consist of brake, tyre, and road wear. The street parameters (building height and street width) used in MUNICH originate from the OSPM setups.

OSPM represents NO₂ and O₃ chemical transformations using a system of two reactions (Berkowicz et al., 1997; Berkowicz, 2000b). The first one describes the production of NO₂ due to the reaction of NO with O₃, and the second one describes the photodissociation of NO₂ that leads to the reproduction of NO and O₃. For a fair comparison, MUNICH is configured to run with a simple chemistry scheme, the Leighton photostationary state for O₃ (Leighton, 1961; Kim et al., 2018):

$$NO_2 + h\nu \rightarrow NO + O(^3P), \qquad (R1)$$

$$O(^3P) + O_2 + M \rightarrow O_3 + M, \qquad (R2)$$

$$NO + O_3 \rightarrow NO_2 + O_2. \qquad (R3)$$

## 3.1 H. C. Andersens Boulevard

H. C. Andersens Boulevard is a wide, densely trafficked boulevard, with an aspect ratio $a_r = H/W$ of about 0.2 ($H = 9.9$ m, $W = 50$ m). It is open on one side with trees instead of buildings. This configuration can be represented in OSPM. However, in MUNICH, a mean building height is defined for each street. Here, it is estimated by averaging the building and the tree heights. The simulated street network is com-

posed of 86 street segments centred around the street where the observation station is located (see brown cross in Fig. 4). Figure 3 presents the monthly averaged concentrations of CO and NO₂ from the OSPM and MUNICH simulations compared to observations. Appendix C contains monthly averaged concentrations of NOₓ and O₃ and statistical indicators of the comparison for all four pollutants.

OSPM-R1 is the receptor that is close to the measurement station and is thus better suited to be compared to observations. Differences in the concentrations between OSPM-R1 and OSPM-R2 highlight the importance of taking into account the horizontal heterogeneities in the street. Observations lie between the concentrations simulated at the two receptors, except in the case of CO, for which OSPM slightly underestimates concentrations in the first half of the year. This underestimation could be linked to underestimation of sources other than traffic, e.g. biomass burning, at the regional scale but could also be due to the absence of volatile organic compounds in the simulation. Overall, the concentrations are well estimated with OSPM, with errors between 26% and 33% for CO and between 34% and 46% for NO₂.

For CO, NO₂, and NOₓ, which are emitted by traffic in the bottom of the street, the concentrations are higher in the first level, MUNICH-hete-11, near the bottom and are lower in the third level, MUNICH-hete-13, near the roof level. For O₃, the opposite behaviour is observed as it is mainly imported by the atmosphere above the street, and it is titrated by NO near the ground. The first two levels, MUNICH-hete-11 and MUNICH-hete-12, have higher concentrations of CO, NO₂, and NOₓ than munich-homo, thus improving the comparison to observations and OSPM concentrations. For CO, the error is improved from 36% in MUNICH-homo to 28% in MUNICH-hete-11, and for NO₂, it is improved from 48% in MUNICH-homo to 35% in MUNICH-hete-11. Although the concentrations of NO₂ and CO are underestimated in MUNICH-homo compared to observations, they are well modelled in the first level, MUNICH-hete-11, with low errors and biases (see Appendix C). For CO, the concentrations simulated in the first level, MUNICH-hete-11, lie between the two OSPM receptors for most of the year. This is the case for NO₂ as well, except between April and August, when NO₂ concentrations are underestimated in MUNICH-hete compared to observations and OSPM concentrations. The formation of NO₂ in the lower levels (see Reaction R3) could be limited with the current model setup because volatile organic compounds are not taken into account in the Danish cases. However, they may participate in O₃ and NO₂ formation through HO₂ and RO₂ radicals (Atkinson, 2000; Kwak and Baik, 2012; Zhong et al., 2017; Dai et al., 2021). Furthermore, the vertical discretization is coarse, limiting O₃ transport deep into the street.

Over the whole street network, as presented in Fig. 4, the differences between the concentrations simulated in the first vertical level of MUNICH-hete and those simulated in MUNICH-homo vary. For wide streets and avenues with
simulations are performed only for the streets where there are observations. Thanks to its coupled approach of a Gaussian plume model and a box model, OSPM is able to calculate concentrations on the two sides of the street (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). Two receptors are used to compare with the observed concentrations on each side of the street: OSPM-R1 and OSPM-R2. The height of the receptors is taken as 2 m to correspond roughly to the height at which observations were performed. They are compared to the concentrations simulated at the first two levels of MUNICH-hete, which are representative of the concentrations at 1 and 3 m.

MUNICH and OSPM use the same input data to estimate the street concentrations. The meteorological parameters originate from simulations performed with the Weather Research and Forecasting model (WRF, Skamarock et al., 2008). The background concentrations are simulated using the Urban Background Model (UBM, Berkowicz, 2000a), which is “a multiple source model that applies a Gaussian approach for horizontal dispersion and a linear approach for vertical dispersion up to the boundary layer” (Jensen et al., 2016). Traffic emissions are generated using the procedure implemented in the local-scale Gaussian air pollution model OML-Highway (Olesen et al., 2015), allowing for precise information for each street segment. Traffic data (average daily traffic, travel speed, and share of heavy-duty vehicles) are used to generate emissions by use of the European emission model COPERT IV. Traffic emissions include exhaust emissions of gases and particles and non-exhaust emissions of particles. Non-exhaust emissions consist of brake, tyre, and road wear. The street parameters (building height and street width) used in MUNICH originate from the OSPM setups.

OSPM represents NO₂ and O₃ chemical transformations using a system of two reactions (Berkowicz et al., 1997; Berkowicz, 2000b). The first one describes the production of NO₂ due to the reaction of NO with O₃, and the second one describes the photodissociation of NO₂ that leads to the reproduction of NO and O₃. For a fair comparison, MUNICH is configured to run with a simple chemistry scheme, the Leighton photostationary state for O₃ (Leighton, 1961; Kim et al., 2018):

$$NO_2 + h\nu \rightarrow NO + O(^3P), \qquad (R1)$$

$$O(^3P) + O_2 + M \rightarrow O_3 + M, \qquad (R2)$$

$$NO + O_3 \rightarrow NO_2 + O_2. \qquad (R3)$$

## 3.1 H. C. Andersens Boulevard

H. C. Andersens Boulevard is a wide, densely trafficked boulevard, with an aspect ratio a_r = H/W of about 0.2 (H = 9.9 m, W = 50 m). It is open on one side with trees instead of buildings. This configuration can be represented in OSPM. However, in MUNICH, a mean building height is defined for each street. Here, it is estimated by averaging the building and the tree heights. The simulated street network is com-

posed of 86 street segments centred around the street where the observation station is located (see brown cross in Fig. 4). Figure 3 presents the monthly averaged concentrations of CO and NO₂ from the OSPM and MUNICH simulations compared to observations. Appendix C contains monthly averaged concentrations of NOₓ and O₃ and statistical indicators of the comparison for all four pollutants.

OSPM-R1 is the receptor that is close to the measurement station and is thus better suited to be compared to observations. Differences in the concentrations between OSPM-R1 and OSPM-R2 highlight the importance of taking into account the horizontal heterogeneities in the street. Observations lie between the concentrations simulated at the two receptors, except in the case of CO, for which OSPM slightly underestimates concentrations in the first half of the year. This underestimation could be linked to underestimation of sources other than traffic, e.g. biomass burning, at the regional scale but could also be due to the absence of volatile organic compounds in the simulation. Overall, the concentrations are well estimated with OSPM, with errors between 26% and 33% for CO and between 34% and 46% for NO₂.

For CO, NO₂, and NOₓ, which are emitted by traffic in the bottom of the street, the concentrations are higher in the first level, MUNICH-hete-11, near the bottom and are lower in the third level, MUNICH-hete-13, near the roof level. For O₃, the opposite behaviour is observed as it is mainly imported by the atmosphere above the street, and it is titrated by NO near the ground. The first two levels, MUNICH-hete-11 and MUNICH-hete-12, have higher concentrations of CO, NO₂, and NOₓ than munich-homo, thus improving the comparison to observations and OSPM concentrations. For CO, the error is improved from 36% in MUNICH-homo to 28% in MUNICH-hete-11, and for NO₂, it is improved from 48% in MUNICH-homo to 35% in MUNICH-hete-11. Although the concentrations of NO₂ and CO are underestimated in MUNICH-homo compared to observations, they are well modelled in the first level, MUNICH-hete-11, with low errors and biases (see Appendix C). For CO, the concentrations simulated in the first level, MUNICH-hete-11, lie between the two OSPM receptors for most of the year. This is the case for NO₂ as well, except between April and August, when NO₂ concentrations are underestimated in MUNICH-hete compared to observations and OSPM concentrations. The formation of NO₂ in the lower levels (see Reaction R3) could be limited with the current model setup because volatile organic compounds are not taken into account in the Danish cases. However, they may participate in O₃ and NO₂ formation through HO₂ and RO₂ radicals (Atkinson, 2000; Kwak and Baik, 2012; Zhong et al., 2017; Dai et al., 2021). Furthermore, the vertical discretization is coarse, limiting O₃ transport deep into the street.

Over the whole street network, as presented in Fig. 4, the differences between the concentrations simulated in the first vertical level of MUNICH-hete and those simulated in MUNICH-homo vary. For wide streets and avenues with
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Monthly averaged concentrations (in µg m⁻³) of CO (a) and NO₂ (b) at HCAB monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represents the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The solid red line represents the OSPM receptor that is close to the measurement station, and the solid orange line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** CO and NO₂ time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the HCAB street network. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for CO and NO₂ are shown in the upper- and lower-right panels (b, d) respectively for the HCAB street network. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Monthly averaged concentrations (in µg m⁻³) of CO (a) and NO₂ (b) at HCAB monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represents the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The solid red line represents the OSPM receptor that is close to the measurement station, and the solid orange line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4.** CO and NO₂ time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the HCAB street network. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for CO and NO₂ are shown in the upper- and lower-right panels (b, d) respectively for the HCAB street network. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
dense traffic, the concentrations are higher in MUNICH-hete-11 than in MUNICH-homo, with an increase of up to 23% for CO and 30% for NO₂. This increase is lower in more narrow and less frequented streets. Although the concentrations in MUNICH-hete-11 are always higher than those in MUNICH-homo for CO, for NO₂, in narrow streets, the concentrations are lower in MUNICH-hete-11. These lower NO₂ concentrations are probably related to the limited transport of O₃ from the background to the bottom of the street, limiting the titration of NO.

## 3.2 Jagtvej

Jagtvej is a conventional street canyon with an aspect ratio of about 0.8 ($H = 20.3$ m, $W = 26.2$ m). The simulated street network is composed of 265 street segments centred around the street where the observation station is located (see brown cross in Fig. 6). The monthly averaged concentrations of CO and NO₂ from the OSPM and MUNICH simulations compared to observations for NO₂ only are presented in Fig. 5. Appendix D contains monthly averaged concentrations of NOₓ and O₃ and statistical indicators of the comparison for NO₂ and NOₓ.

OSPM-R2 is the receptor that is close to the measurement station and is thus better suited to be compared to observations. The differences between the concentrations of the two OSPM receptors are lower in JGTV than in HCAB (see Fig. 3). JGTV is narrower, with higher buildings on both sides, thus limiting the ventilation zone and the penetration into the street of background concentrations that would reduce the concentrations at the downwind receptor. OSPM tends to slightly overestimate NO₂ and NOₓ concentrations (with errors between 41% and 53% and between 42% and 63% respectively).

Concentrations from the homogeneous version of MUNICH are close to OSPM concentrations for CO. They are lower than OSPM concentrations for NO₂, but they compare well to observations, with a bias of -3% and an error of 46%. These lower NO₂ concentrations are linked to lower NOₓ concentrations (with a bias of -19% compared to observations and an error of 53%). As for HCAB, the MUNICH-hete concentrations decrease from the bottom to the top of the street. For CO and NOₓ, only the first level has concentrations higher than MUNICH-homo, while for NO₂, all three levels have concentrations lower than MUNICH-homo. O₃ concentrations are also higher for the three levels (see Fig. D1). The NOₓ concentrations of MUNICH-hete-11 compare slightly better to observations than MUNICH-homo, while MUNICH-homo is slightly better for NO₂. However, the statistics of the two models are quite close (Appendix D).

Conclusions are similar for the whole street network (see Fig. 6). The CO concentrations are slightly higher in the first level of MUNICH-hete than in MUNICH-homo. However, NO₂ concentrations at the bottom of the street in the first

level of MUNICH-hete tend to be lower than in MUNICH-homo. In some specific street segments of the network, the differences in the concentrations for both CO and NO₂ are higher in MUNICH-hete-11 than in MUNICH-homo. This is due to a mix of different traffic emissions and street morphologies which favour the transport of pollutants.

## 4 Application to a street network in greater Paris

The impacts of the discretization on gas and particle concentrations are evaluated over the street network near Paris, France, which was used to validate MUNICH v2.0 (Kim et al., 2022) and in several sensitivity studies (Lugon et al., 2021b; Sarica et al., 2022, 2023b). The street network represents a district of Le Perreux-sur-Marne, a suburb 13 km east of Paris, France. It is composed of 577 street segments (see Fig. 8). The street parameters (building height and street width) were obtained from the BD TOPO database (https://geoservices.ign.fr/bdtopo, last access: 12 September 2023). Simulations were performed from 22 March to 15 June 2014 to generate hourly concentrations with input data (emissions, including exhaust emissions and brake, tyre, and road wear; meteorological parameters; and background concentrations) from the reference simulation SCNO of Sarica et al. (2023b). For this case, MUNICH is coupled to SSH-aerosol (Sartelet et al., 2020) to represent gas-phase chemistry and aerosol dynamics. Observational data are available for the whole simulation period for NO₂, NOₓ, PM₂.₅, PM₁₀, and black carbon (BC) for a segment of Boulevard d'Alsace Lorraine (see brown cross in Fig. 8). The simulations were performed at a height of about 2 m, and they are thus to be compared to the concentrations of the first two levels ($i = 1$ and $i = 2$) of MUNICH-hete. The segment has an aspect ratio of about 0.3 ($H = 8.6$ m, $W = 26$ m).

The average daily profiles of NOₓ and PM₁₀ concentrations are shown in Fig. 7, and NO, NO₂, and the statistical indicators are shown in Appendix E1. As in the Danish streets, for pollutants emitted by traffic, the concentrations are higher in MUNICH-hete at the bottom of the street, decreasing to the roof level. The concentrations of PM₁₀, NO₂, NOₓ, and BC are higher on average by 12%, 21%, 40%, and 30% respectively in MUNICH-hete-11 than in MUNICH-homo. The lower concentration difference for PM₁₀ than for the other compounds reflects the fact that non-traffic sources are more important for PM₁₀, and inversely, they are small for BC. Despite the higher concentrations in MUNICH-hete-11, BC concentrations remain strongly underestimated compared to observations, in agreement with the CFD simulations of Lin et al. (2023). For PM₁₀, NO₂, and NOₓ, the concentrations compare well to observations; e.g. the error is 33% for NO₂ and 37% for PM₁₀, with slightly better statistics using MUNICH-hete-11 than MUNICH-homo. For all pollutants, the concentrations of the intermediate level ($i = 2$) are very close to the concentrations of MUNICH.
dense traffic, the concentrations are higher in MUNICH-hete-11 than in MUNICH-homo, with an increase of up to 23% for CO and 30% for NO₂. This increase is lower in more narrow and less frequented streets. Although the concentrations in MUNICH-hete-11 are always higher than those in MUNICH-homo for CO, for NO₂, in narrow streets, the concentrations are lower in MUNICH-hete-11. These lower NO₂ concentrations are probably related to the limited transport of O₃ from the background to the bottom of the street, limiting the titration of NO.

## 3.2 Jagtvej

Jagtvej is a conventional street canyon with an aspect ratio of about 0.8 ($H = 20.3$ m, $W = 26.2$ m). The simulated street network is composed of 265 street segments centred around the street where the observation station is located (see brown cross in Fig. 6). The monthly averaged concentrations of CO and NO₂ from the OSPM and MUNICH simulations compared to observations for NO₂ only are presented in Fig. 5. Appendix D contains monthly averaged concentrations of NOₓ and O₃ and statistical indicators of the comparison for NO₂ and NOₓ.

OSPM-R2 is the receptor that is close to the measurement station and is thus better suited to be compared to observations. The differences between the concentrations of the two OSPM receptors are lower in JGTV than in HCAB (see Fig. 3). JGTV is narrower, with higher buildings on both sides, thus limiting the ventilation zone and the penetration into the street of background concentrations that would reduce the concentrations at the downwind receptor. OSPM tends to slightly overestimate NO₂ and NOₓ concentrations (with errors between 41% and 53% and between 42% and 63% respectively).

Concentrations from the homogeneous version of MUNICH are close to OSPM concentrations for CO. They are lower than OSPM concentrations for NO₂, but they compare well to observations, with a bias of -3% and an error of 46%. These lower NO₂ concentrations are linked to lower NOₓ concentrations (with a bias of -19% compared to observations and an error of 53%). As for HCAB, the MUNICH-hete concentrations decrease from the bottom to the top of the street. For CO and NOₓ, only the first level has concentrations higher than MUNICH-homo, while for NO₂, all three levels have concentrations lower than MUNICH-homo. O₃ concentrations are also higher for the three levels (see Fig. D1). The NOₓ concentrations of MUNICH-hete-11 compare slightly better to observations than MUNICH-homo, while MUNICH-homo is slightly better for NO₂. However, the statistics of the two models are quite close (Appendix D).

Conclusions are similar for the whole street network (see Fig. 6). The CO concentrations are slightly higher in the first level of MUNICH-hete than in MUNICH-homo. However, NO₂ concentrations at the bottom of the street in the first

level of MUNICH-hete tend to be lower than in MUNICH-homo. In some specific street segments of the network, the differences in the concentrations for both CO and NO₂ are higher in MUNICH-hete-11 than in MUNICH-homo. This is due to a mix of different traffic emissions and street morphologies which favour the transport of pollutants.

## 4 Application to a street network in greater Paris

The impacts of the discretization on gas and particle concentrations are evaluated over the street network near Paris, France, which was used to validate MUNICH v2.0 (Kim et al., 2022) and in several sensitivity studies (Lugon et al., 2021b; Sarica et al., 2022, 2023b). The street network represents a district of Le Perreux-sur-Marne, a suburb 13 km east of Paris, France. It is composed of 577 street segments (see Fig. 8). The street parameters (building height and street width) were obtained from the BD TOPO database (https://geoservices.ign.fr/bdtopo, last access: 12 September 2023). Simulations were performed from 22 March to 15 June 2014 to generate hourly concentrations with input data (emissions, including exhaust emissions and brake, tyre, and road wear; meteorological parameters; and background concentrations) from the reference simulation SCNO of Sarica et al. (2023b). For this case, MUNICH is coupled to SSH-aerosol (Sartelet et al., 2020) to represent gas-phase chemistry and aerosol dynamics. Observational data are available for the whole simulation period for NO₂, NOₓ, PM₂.₅, PM₁₀, and black carbon (BC) for a segment of Boulevard d'Alsace Lorraine (see brown cross in Fig. 8). The simulations were performed at a height of about 2 m, and they are thus to be compared to the concentrations of the first two levels ($i = 1$ and $i = 2$) of MUNICH-hete. The segment has an aspect ratio of about 0.3 ($H = 8.6$ m, $W = 26$ m).

The average daily profiles of NOₓ and PM₁₀ concentrations are shown in Fig. 7, and NO, NO₂, and the statistical indicators are shown in Appendix E1. As in the Danish streets, for pollutants emitted by traffic, the concentrations are higher in MUNICH-hete at the bottom of the street, decreasing to the roof level. The concentrations of PM₁₀, NO₂, NOₓ, and BC are higher on average by 12%, 21%, 40%, and 30% respectively in MUNICH-hete-11 than in MUNICH-homo. The lower concentration difference for PM₁₀ than for the other compounds reflects the fact that non-traffic sources are more important for PM₁₀, and inversely, they are small for BC. Despite the higher concentrations in MUNICH-hete-11, BC concentrations remain strongly underestimated compared to observations, in agreement with the CFD simulations of Lin et al. (2023). For PM₁₀, NO₂, and NOₓ, the concentrations compare well to observations; e.g. the error is 33% for NO₂ and 37% for PM₁₀, with slightly better statistics using MUNICH-hete-11 than MUNICH-homo. For all pollutants, the concentrations of the intermediate level ($i = 2$) are very close to the concentrations of MUNICH.
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Monthly averaged concentrations (in µg m⁻³) of CO (a) and NO₂ (b) at JGTV monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The solid orange line represents the OSPM receptor that is close to the measurement station, and the solid red line represents the second OSPM receptor located on the other side of the street. The observations (only available for NO₂) are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** CO and NO₂ time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the JGTV street network. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for CO and NO₂ are shown in the upper- and lower-right panels (b, d) respectively for JGTV. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Monthly averaged concentrations (in µg m⁻³) of CO (a) and NO₂ (b) at JGTV monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The solid orange line represents the OSPM receptor that is close to the measurement station, and the solid red line represents the second OSPM receptor located on the other side of the street. The observations (only available for NO₂) are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 6.** CO and NO₂ time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the JGTV street network. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for CO and NO₂ are shown in the upper- and lower-right panels (b, d) respectively for JGTV. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
homo. This is due to a mixture of different parameters such as street morphology and dispersion conditions over the simulated period.

In streets, both BC and organic matter (OM) exhibit higher concentrations than in the background (Lugon et al., 2021a). OM consists of primary and secondary aerosols that are formed from the oxidation of volatile organic compounds and/or the condensation of semi-volatile organic compounds. Concentrations of both BC and OM simulated in MUNICH-homo and in the first level of MUNICH-hete are compared in Fig 8 (see Fig. E2 for NO₂ and PM₁₀). The concentrations simulated in MUNICH-hete-11 are always higher than in MUNICH-homo. In most streets of the network that are narrow and limited in terms of traffic, the increase in concentrations is limited. It is more important for more open and frequented streets, such as the Boulevard d'Alsace Lorraine.

With BC being a primary inert pollutant emitted by traffic, its concentrations are strongly influenced by the discretization as emissions are no longer artificially diluted in the whole street volume. They are now constraint at the bottom of the street, inducing an average increase of 30% compared to the homogeneous version. For OM and PM₁₀, the increase is lower (16% and 12% respectively) because of the stronger influence of non-traffic sources.

## 5 Sensitivity analysis

In this section, the influence of the aspect ratio $a_r$ on the concentrations in the heterogeneous version of MUNICH is studied. The sensitivity of the heterogeneous version of MUNICH to the presence of a street network, which influences the concentrations entering the street via the upwind intersection, is also estimated.

### 5.1 Influence of the aspect ratio

The aspect ratio $a_r$, defined as the ratio between building height $H$ and street width $W$, is used to determine general behaviours of streets with similar geometries. When $a_r$ is small, i.e. when $W$ is larger than $H$, the street is wide, such as in the case of a boulevard or an avenue. On the other hand, when $a_r$ is large, i.e. when $H$ is larger than $W$, the street is narrow and closer to a typical street canyon. In the three cases presented in this study in greater Paris and Copenhagen, streets with small $a_r$ are associated with higher traffic emissions.

Table 1 and Appendix F present the influence of the aspect ratio on concentrations in the first vertical level of the heterogeneous version of MUNICH, MUNICH-hete-11, compared to the homogeneous version. It is quantified using the normalized mean bias (NMB). The homogeneous version of MUNICH presents higher concentrations of pollutants in wide streets (interval [0, 0.5]) in greater Paris due to higher traffic emissions compared to in narrower streets. The increase in the concentration in the first level as a result of

<table><caption>**Table 1.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network in greater Paris. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5]</td><td>54.92</td><td>33.99</td><td></td></tr><tr><td>[0.5, 1]</td><td>41.13</td><td>11.95</td><td></td></tr><tr><td>[1, 1.5]</td><td>34.26</td><td>6.27</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5]</td><td>84.31</td><td>53.06</td><td></td></tr><tr><td>[0.5, 1]</td><td>59.85</td><td>22.90</td><td></td></tr><tr><td>[1, 1.5]</td><td>43.71</td><td>12.02</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">PM<sub>10</sub></td><td>[0, 0.5]</td><td>23.06</td><td>18.01</td><td></td></tr><tr><td>[0.5, 1]</td><td>20.59</td><td>6.09</td><td></td></tr><tr><td>[1, 1.5]</td><td>19.13</td><td>2.63</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">PM<sub>2.5</sub></td><td>[0, 0.5]</td><td>21.10</td><td>16.01</td><td></td></tr><tr><td>[0.5, 1]</td><td>19.01</td><td>5.32</td><td></td></tr><tr><td>[1, 1.5]</td><td>17.91</td><td>2.50</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">BC</td><td>[0, 0.5]</td><td>1.56</td><td>41.84</td><td></td></tr><tr><td>[0.5, 1]</td><td>1.19</td><td>16.28</td><td></td></tr><tr><td>[1, 1.5]</td><td>0.96</td><td>7.78</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">OM</td><td>[0, 0.5]</td><td>5.60</td><td>22.58</td><td></td></tr><tr><td>[0.5, 1]</td><td>4.91</td><td>7.98</td><td></td></tr><tr><td>[1, 1.5]</td><td>4.46</td><td>3.27</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr></tbody></table>

using the heterogeneous version of the model is more important for wide streets (up to 53% for NOₓ). Concerning particles, PM₁₀ concentrations are increased by up to 18%, while the increase is about 42% and 23% for BC and OM respectively. In streets with larger $a_r$ (interval [1, 1.5]), NOₓ concentrations are only 12% higher than in the homogeneous version, and the PM₁₀ concentrations are only 2.5% higher.

Concerning the two cases in Copenhagen, CO and NOₓ concentrations present similar behaviour in comparison to the greater Paris case. For the JGTV street network, the increase is limited compared to the HCAB network – 3% against 11.5% respectively for CO. Wide streets (interval [0, 0.5]) have higher concentrations of NO₂ compared to the homogeneous version of MUNICH, whereas for the other intervals of $a_r$, they are lower. Thus, O₃ concentrations are lower for wide streets and larger for the other intervals. This could be linked to the lack of volatile organic compounds in these two Copenhagen cases.
homo. This is due to a mixture of different parameters such as street morphology and dispersion conditions over the simulated period.

In streets, both BC and organic matter (OM) exhibit higher concentrations than in the background (Lugon et al., 2021a). OM consists of primary and secondary aerosols that are formed from the oxidation of volatile organic compounds and/or the condensation of semi-volatile organic compounds. Concentrations of both BC and OM simulated in MUNICH-homo and in the first level of MUNICH-hete are compared in Fig 8 (see Fig. E2 for NO₂ and PM₁₀). The concentrations simulated in MUNICH-hete-11 are always higher than in MUNICH-homo. In most streets of the network that are narrow and limited in terms of traffic, the increase in concentrations is limited. It is more important for more open and frequented streets, such as the Boulevard d'Alsace Lorraine.

With BC being a primary inert pollutant emitted by traffic, its concentrations are strongly influenced by the discretization as emissions are no longer artificially diluted in the whole street volume. They are now constraint at the bottom of the street, inducing an average increase of 30% compared to the homogeneous version. For OM and PM₁₀, the increase is lower (16% and 12% respectively) because of the stronger influence of non-traffic sources.

## 5 Sensitivity analysis

In this section, the influence of the aspect ratio $a_r$ on the concentrations in the heterogeneous version of MUNICH is studied. The sensitivity of the heterogeneous version of MUNICH to the presence of a street network, which influences the concentrations entering the street via the upwind intersection, is also estimated.

### 5.1 Influence of the aspect ratio

The aspect ratio $a_r$, defined as the ratio between building height $H$ and street width $W$, is used to determine general behaviours of streets with similar geometries. When $a_r$ is small, i.e. when $W$ is larger than $H$, the street is wide, such as in the case of a boulevard or an avenue. On the other hand, when $a_r$ is large, i.e. when $H$ is larger than $W$, the street is narrow and closer to a typical street canyon. In the three cases presented in this study in greater Paris and Copenhagen, streets with small $a_r$ are associated with higher traffic emissions.

Table 1 and Appendix F present the influence of the aspect ratio on concentrations in the first vertical level of the heterogeneous version of MUNICH, MUNICH-hete-11, compared to the homogeneous version. It is quantified using the normalized mean bias (NMB). The homogeneous version of MUNICH presents higher concentrations of pollutants in wide streets (interval [0, 0.5]) in greater Paris due to higher traffic emissions compared to in narrower streets. The increase in the concentration in the first level as a result of

<table><caption>**Table 1.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network in greater Paris. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5]</td><td>54.92</td><td>33.99</td><td></td></tr><tr><td>[0.5, 1]</td><td>41.13</td><td>11.95</td><td></td></tr><tr><td>[1, 1.5]</td><td>34.26</td><td>6.27</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5]</td><td>84.31</td><td>53.06</td><td></td></tr><tr><td>[0.5, 1]</td><td>59.85</td><td>22.90</td><td></td></tr><tr><td>[1, 1.5]</td><td>43.71</td><td>12.02</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">PM<sub>10</sub></td><td>[0, 0.5]</td><td>23.06</td><td>18.01</td><td></td></tr><tr><td>[0.5, 1]</td><td>20.59</td><td>6.09</td><td></td></tr><tr><td>[1, 1.5]</td><td>19.13</td><td>2.63</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">PM<sub>2.5</sub></td><td>[0, 0.5]</td><td>21.10</td><td>16.01</td><td></td></tr><tr><td>[0.5, 1]</td><td>19.01</td><td>5.32</td><td></td></tr><tr><td>[1, 1.5]</td><td>17.91</td><td>2.50</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">BC</td><td>[0, 0.5]</td><td>1.56</td><td>41.84</td><td></td></tr><tr><td>[0.5, 1]</td><td>1.19</td><td>16.28</td><td></td></tr><tr><td>[1, 1.5]</td><td>0.96</td><td>7.78</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr><tr><td rowspan="4">OM</td><td>[0, 0.5]</td><td>5.60</td><td>22.58</td><td></td></tr><tr><td>[0.5, 1]</td><td>4.91</td><td>7.98</td><td></td></tr><tr><td>[1, 1.5]</td><td>4.46</td><td>3.27</td><td></td></tr><tr><td>≥ 1.5</td><td>—</td><td>—</td><td></td></tr></tbody></table>

using the heterogeneous version of the model is more important for wide streets (up to 53% for NOₓ). Concerning particles, PM₁₀ concentrations are increased by up to 18%, while the increase is about 42% and 23% for BC and OM respectively. In streets with larger $a_r$ (interval [1, 1.5]), NOₓ concentrations are only 12% higher than in the homogeneous version, and the PM₁₀ concentrations are only 2.5% higher.

Concerning the two cases in Copenhagen, CO and NOₓ concentrations present similar behaviour in comparison to the greater Paris case. For the JGTV street network, the increase is limited compared to the HCAB network – 3% against 11.5% respectively for CO. Wide streets (interval [0, 0.5]) have higher concentrations of NO₂ compared to the homogeneous version of MUNICH, whereas for the other intervals of $a_r$, they are lower. Thus, O₃ concentrations are lower for wide streets and larger for the other intervals. This could be linked to the lack of volatile organic compounds in these two Copenhagen cases.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Average daily profile of the concentrations of NOₓ (a) and PM₁₀ (b) over the simulation period at the Boulevard d'Alsace Lorraine monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The observations are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** OM and BC time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the district of Le Perreux-sur-Marne. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for OM and BC are shown in the upper- and lower-right panels (b, d) respectively for the district of Le Perreux-sur-Marne. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Average daily profile of the concentrations of NOₓ (a) and PM₁₀ (b) over the simulation period at the Boulevard d'Alsace Lorraine monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The observations are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 8.** OM and BC time-averaged concentrations (in µg m⁻³) for MUNICH-homo are shown in the upper- and lower-left panels (a, c) respectively for the district of Le Perreux-sur-Marne. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for OM and BC are shown in the upper- and lower-right panels (b, d) respectively for the district of Le Perreux-sur-Marne. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.</figcaption></figure>
## 5.2 Sensitivity to the street network

The computation of inflow and outflow fluxes at intersections is performed by estimating the balance of fluxes entering and leaving the intersection from the different street segments attached to it. If this balance is not perfect, there are exchanges with the atmosphere above the intersection. When the total flux entering the intersection is higher than the one leaving it, the flow overload is directed to the atmosphere. When the total flux leaving the intersection is higher, a flux from the atmosphere to the intersection is considered. A further explanation is available in Kim et al. (2018, 2022).

Without a street network around the street segment of interest, the pollutant mass fluxes entering the street are determined from the background concentration. The contribution of the neighbouring streets is thus not taken into account, and concentrations of pollutants emitted in streets are expected to be lower than when there is a street network.

Table 2 and Appendix G present the influence of the neighbouring streets for the three cases of the study. It is quantified using the normalized mean error (NME) and NMB between the simulations without and with the neighbouring streets. As expected, without the street network, the concentrations are lower. The bias is between 12% and 21% for NO₂, between 18% and 27% for NOₓ, and 14% for BC. Biases are lower for OM and PM (2%) because of the stronger influence of background concentrations for those compounds.

## 6 Conclusions

The street-network model MUNICH v2.0 has been modified to introduce concentration heterogeneities in the street and to better represent population exposure. To model the vertical gradients frequently observed, the streets were discretized with three levels, thus limiting the artificial dilution of emissions and concentrations. Based on a parameterization from OSPM, a ventilation zone is considered under specific conditions to represent horizontal heterogeneities. In order to test these developments, the heterogeneous version of MUNICH (MUNICH-hete) has been applied to two cases in Copenhagen, Denmark, with comparisons to OSPM, and to one case near Paris, France. Overall, MUNICH-hete improves the comparison to observations compared to the homogeneous version. The errors compared to observations are reduced by up to 20% for NOₓ and up to 15% for BC.

As expected, in MUNICH-hete, concentrations of compounds emitted by traffic (CO, NO₂, NOₓ, PM₁₀, BC, and OM) are higher at the bottom of the street than at the top. These increases can reach up to 60% and 30% for NO₂ and PM₁₀ respectively. The intermediate level, serving as a buffer, presents concentrations higher than or similar to the homogeneous version (MUNICH-homo). Finally, concentrations in the highest level, in direct contact with the atmosphere above the street, are the lowest of the street. For the

<figure><img src="image_12.png" /><figcaption>**Figure A1.** Representation of the recirculation zone in the heterogeneous version of MUNICH.</figcaption></figure>

Danish cases, the low NO₂ concentrations observed in the lower levels could be related to the absence of volatile organic compounds in the model setup and the coarse vertical discretization limiting O₃ transport deep into the street.

A sensitivity study of the influence of the street network on concentrations in the streets shows the importance of considering neighbouring streets in MUNICH. When no network is considered, concentrations in the street are lower due to the overestimated impact of the atmosphere above. At the bottom of the street, concentrations of NO₂ and BC are reduced by up to 28% and 14% respectively without a network. PM and OM are less impacted, with a reduction of about 2% due to a strong influence of non-traffic sources.

For the next step, the ventilation zone will be fully discretized vertically to facilitate the penetration of background concentrations into the bottom of the street. The horizontal exchange fluxes between the two zones will also be modelled. Deposition and resuspension processes that were not considered in this development will be added. The fluxes that are currently assumed to be vertically homogeneous will be discretized. Finally, chemistry of volatile organic compounds could be added in the Danish cases if the background concentrations and emissions are available.

## Appendix A: Volumes of the recirculation and ventilation zones

The algorithm used in the heterogeneous version of MUNICH to compute the volumes of the recirculation and ventilation zones is based on the parameterization of OSPM (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). This algorithm is applied at the beginning of each time step as the size of the recirculation zone is dependent on wind speed and direction. Figure A1 presents the shape of the recirculation zone and the associated parameters.
## 5.2 Sensitivity to the street network

The computation of inflow and outflow fluxes at intersections is performed by estimating the balance of fluxes entering and leaving the intersection from the different street segments attached to it. If this balance is not perfect, there are exchanges with the atmosphere above the intersection. When the total flux entering the intersection is higher than the one leaving it, the flow overload is directed to the atmosphere. When the total flux leaving the intersection is higher, a flux from the atmosphere to the intersection is considered. A further explanation is available in Kim et al. (2018, 2022).

Without a street network around the street segment of interest, the pollutant mass fluxes entering the street are determined from the background concentration. The contribution of the neighbouring streets is thus not taken into account, and concentrations of pollutants emitted in streets are expected to be lower than when there is a street network.

Table 2 and Appendix G present the influence of the neighbouring streets for the three cases of the study. It is quantified using the normalized mean error (NME) and NMB between the simulations without and with the neighbouring streets. As expected, without the street network, the concentrations are lower. The bias is between 12% and 21% for NO₂, between 18% and 27% for NOₓ, and 14% for BC. Biases are lower for OM and PM (2%) because of the stronger influence of background concentrations for those compounds.

## 6 Conclusions

The street-network model MUNICH v2.0 has been modified to introduce concentration heterogeneities in the street and to better represent population exposure. To model the vertical gradients frequently observed, the streets were discretized with three levels, thus limiting the artificial dilution of emissions and concentrations. Based on a parameterization from OSPM, a ventilation zone is considered under specific conditions to represent horizontal heterogeneities. In order to test these developments, the heterogeneous version of MUNICH (MUNICH-hete) has been applied to two cases in Copenhagen, Denmark, with comparisons to OSPM, and to one case near Paris, France. Overall, MUNICH-hete improves the comparison to observations compared to the homogeneous version. The errors compared to observations are reduced by up to 20% for NOₓ and up to 15% for BC.

As expected, in MUNICH-hete, concentrations of compounds emitted by traffic (CO, NO₂, NOₓ, PM₁₀, BC, and OM) are higher at the bottom of the street than at the top. These increases can reach up to 60% and 30% for NO₂ and PM₁₀ respectively. The intermediate level, serving as a buffer, presents concentrations higher than or similar to the homogeneous version (MUNICH-homo). Finally, concentrations in the highest level, in direct contact with the atmosphere above the street, are the lowest of the street. For the

<figure><img src="image_10.png" /><figcaption>**Figure A1.** Representation of the recirculation zone in the heterogeneous version of MUNICH.</figcaption></figure>

Danish cases, the low NO₂ concentrations observed in the lower levels could be related to the absence of volatile organic compounds in the model setup and the coarse vertical discretization limiting O₃ transport deep into the street.

A sensitivity study of the influence of the street network on concentrations in the streets shows the importance of considering neighbouring streets in MUNICH. When no network is considered, concentrations in the street are lower due to the overestimated impact of the atmosphere above. At the bottom of the street, concentrations of NO₂ and BC are reduced by up to 28% and 14% respectively without a network. PM and OM are less impacted, with a reduction of about 2% due to a strong influence of non-traffic sources.

For the next step, the ventilation zone will be fully discretized vertically to facilitate the penetration of background concentrations into the bottom of the street. The horizontal exchange fluxes between the two zones will also be modelled. Deposition and resuspension processes that were not considered in this development will be added. The fluxes that are currently assumed to be vertically homogeneous will be discretized. Finally, chemistry of volatile organic compounds could be added in the Danish cases if the background concentrations and emissions are available.

## Appendix A: Volumes of the recirculation and ventilation zones

The algorithm used in the heterogeneous version of MUNICH to compute the volumes of the recirculation and ventilation zones is based on the parameterization of OSPM (Berkowicz et al., 1997; Berkowicz, 2000b; Ottosen et al., 2015). This algorithm is applied at the beginning of each time step as the size of the recirculation zone is dependent on wind speed and direction. Figure A1 presents the shape of the recirculation zone and the associated parameters.
<table><caption>**Table 2.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment of the Boulevard d'Alsace Lorraine with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th colspan="2">MUNICH-hete level</th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th><th>BC</th><th>OM</th></tr></thead><tbody><tr><td rowspan="4">11</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>58.92</td><td>104.74</td><td>24.58</td><td>22.18</td><td>1.84</td><td>6.13</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>51.50</td><td>85.46</td><td>24.15</td><td>21.97</td><td>1.58</td><td>5.98</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>23.25</td><td>23.63</td><td>17.26</td><td>17.28</td><td>20.95</td><td>20.29</td></tr><tr><td>NMB [%]</td><td>-12.58</td><td>-18.41</td><td>-1.75</td><td>-0.98</td><td>-14.02</td><td>-2.37</td></tr><tr><td rowspan="4">12</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>47.81</td><td>75.36</td><td>21.94</td><td>20.08</td><td>1.42</td><td>5.30</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>42.67</td><td>62.34</td><td>21.77</td><td>20.04</td><td>1.25</td><td>5.25</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>20.97</td><td>23.11</td><td>14.55</td><td>14.46</td><td>19.37</td><td>17.11</td></tr><tr><td>NMB [%]</td><td>-10.74</td><td>-17.28</td><td>-0.75</td><td>-0.20</td><td>-12.15</td><td>-0.96</td></tr><tr><td rowspan="4">13</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>37.97</td><td>52.63</td><td>19.91</td><td>18.41</td><td>1.09</td><td>4.69</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>35.65</td><td>47.01</td><td>19.89</td><td>18.44</td><td>1.02</td><td>4.69</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>12.83</td><td>15.21</td><td>8.11</td><td>8.05</td><td>11.78</td><td>9.47</td></tr><tr><td>NMB [%]</td><td>-6.10</td><td>-10.67</td><td>-0.08</td><td>0.15</td><td>-6.70</td><td>-0.12</td></tr></tbody></table>

The first step is to compute the length of the vortex in the direction of the wind:

$$L_{\text{vortex}} = 2Hf, \qquad (A1)$$

with *H* being the height of the street, and

$$f = \begin{cases} 1 & \text{if } u_{\text{roof}} \ge 2 \, \text{m s}^{-1} \\ \sqrt{0.5u_{\text{roof}}} & \text{if } u_{\text{roof}} < 2 \, \text{m s}^{-1}, \end{cases} \qquad (A2)$$

with *u*<sub>roof</sub> being the wind speed at roof level.

The width of the trapeze base is the projection of *L*<sub>vortex</sub> in the street:

$$W_{\text{base}} = L_{\text{vortex}} \sin(\theta), \qquad (A3)$$

with *θ* being the angle between the wind direction and the street orientation.

The width of the trapeze top is equal to half of the base:

$$W_3 = \frac{L_{\text{vortex}} \sin(\theta)}{2}. \qquad (A4)$$

Knowing these lengths and using algebraic considerations, the widths *W*<sub>1</sub> and *W*<sub>2</sub> can be calculated:

$$\begin{cases} W_1 = W_3 + \frac{(h_2 + h_3)\Delta W}{H}, \\ W_2 = W_3 + \frac{h_3\Delta W}{H} \end{cases}, \qquad (A5)$$

with $\Delta W = W_{\text{base}} - W_3$.

The horizontal surfaces for vertical exchanges between the levels and with the concentrations above the street are calculated as follows:

$$\begin{cases} S_1^h = W_1 L \\ S_2^h = W_2 L \\ S_3^h = W_3 L \end{cases}, \qquad (A6)$$

with *L* being the street length.

The vertical surfaces for advection via intersections are determined with the following:

$$\left\{ \begin{aligned} S_1^v &= \frac{W_1(h_1)^2(W_{\text{base}} - W_1)}{2} \\ S_2^v &= \frac{W_2(h_2)^2(W_1 - W_2)}{2} \\ S_3^v &= \frac{W_3(h_3)^2(W_2 - W_3)}{2} \end{aligned} \right. . \qquad (A7)$$

Finally, the volumes associated with each level of the recirculation zone are

$$\left\{ \begin{aligned} V_1 &= h_1 L \left( W_1 + \frac{(W_{\text{base}} - W_1)}{2} \right) \\ V_2 &= h_2 L \left( W_2 + \frac{(W_1 - W_2)}{2} \right) \\ V_3 &= h_3 L \left( W_3 + \frac{(W_2 - W_3)}{2} \right) \end{aligned} \right. . \qquad (A8)$$

In the current version of MUNICH-hete, it is assumed that traffic emissions are all affected to the recirculation zone. Therefore, its base *W*<sub>base</sub> has to be superior or equal to the street width *W*. If this is not the case, the recirculation zone is considered to fill the whole street volume.

When considered, the other widths are also limited by the street width:

$$\begin{cases} W_1 = \min(W, W_1) \\ W_2 = \min(W, W_2) \\ W_3 = \min(W, W_3) \end{cases}. \qquad (A9)$$

## Appendix B: Statistical indicators

For evaluation of the simulations compared to observations, the following statistical indicators are used. *o* and *s* represent
<table><caption>**Table 2.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment of the Boulevard d'Alsace Lorraine with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th colspan="2">MUNICH-hete level</th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th><th>BC</th><th>OM</th></tr></thead><tbody><tr><td rowspan="4">11</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>58.92</td><td>104.74</td><td>24.58</td><td>22.18</td><td>1.84</td><td>6.13</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>51.50</td><td>85.46</td><td>24.15</td><td>21.97</td><td>1.58</td><td>5.98</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>23.25</td><td>23.63</td><td>17.26</td><td>17.28</td><td>20.95</td><td>20.29</td></tr><tr><td>NMB [%]</td><td>-12.58</td><td>-18.41</td><td>-1.75</td><td>-0.98</td><td>-14.02</td><td>-2.37</td></tr><tr><td rowspan="4">12</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>47.81</td><td>75.36</td><td>21.94</td><td>20.08</td><td>1.42</td><td>5.30</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>42.67</td><td>62.34</td><td>21.77</td><td>20.04</td><td>1.25</td><td>5.25</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>20.97</td><td>23.11</td><td>14.55</td><td>14.46</td><td>19.37</td><td>17.11</td></tr><tr><td>NMB [%]</td><td>-10.74</td><td>-17.28</td><td>-0.75</td><td>-0.20</td><td>-12.15</td><td>-0.96</td></tr><tr><td rowspan="4">13</td><td rowspan="2">With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>37.97</td><td>52.63</td><td>19.91</td><td>18.41</td><td>1.09</td><td>4.69</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>35.65</td><td>47.01</td><td>19.89</td><td>18.44</td><td>1.02</td><td>4.69</td></tr><tr><td rowspan="2">Without network</td><td>NME [%]</td><td>12.83</td><td>15.21</td><td>8.11</td><td>8.05</td><td>11.78</td><td>9.47</td></tr><tr><td>NMB [%]</td><td>-6.10</td><td>-10.67</td><td>-0.08</td><td>0.15</td><td>-6.70</td><td>-0.12</td></tr></tbody></table>

The first step is to compute the length of the vortex in the direction of the wind:

$$L_{\text{vortex}} = 2Hf, \qquad (A1)$$

with *H* being the height of the street, and

$$f = \begin{cases} 1 & \text{if } u_{\text{roof}} \ge 2 \, \text{m s}^{-1} \\ \sqrt{0.5u_{\text{roof}}} & \text{if } u_{\text{roof}} < 2 \, \text{m s}^{-1}, \end{cases} \qquad (A2)$$

with *u*<sub>roof</sub> being the wind speed at roof level.

The width of the trapeze base is the projection of *L*<sub>vortex</sub> in the street:

$$W_{\text{base}} = L_{\text{vortex}} \sin(\theta), \qquad (A3)$$

with *θ* being the angle between the wind direction and the street orientation.

The width of the trapeze top is equal to half of the base:

$$W_3 = \frac{L_{\text{vortex}} \sin(\theta)}{2}. \qquad (A4)$$

Knowing these lengths and using algebraic considerations, the widths *W*<sub>1</sub> and *W*<sub>2</sub> can be calculated:

$$\begin{cases} W_1 = W_3 + \frac{(h_2 + h_3)\Delta W}{H} \\ W_2 = W_3 + \frac{h_3\Delta W}{H} \end{cases}, \qquad (A5)$$

with $\Delta W = W_{\text{base}} - W_3$.

The horizontal surfaces for vertical exchanges between the levels and with the concentrations above the street are calculated as follows:

$$\begin{cases} S_1^h = W_1 L \\ S_2^h = W_2 L \\ S_3^h = W_3 L \end{cases}, \qquad (A6)$$

with *L* being the street length.

The vertical surfaces for advection via intersections are determined with the following:

$$\left\{ \begin{aligned} S_1^v &= \frac{W_1(h_1)^2(W_{\text{base}} - W_1)}{2} \\ S_2^v &= \frac{W_2(h_2)^2(W_1 - W_2)}{2} \\ S_3^v &= \frac{W_3(h_3)^2(W_2 - W_3)}{2} \end{aligned} \right. . \qquad (A7)$$

Finally, the volumes associated with each level of the recirculation zone are

$$\left\{ \begin{aligned} V_1 &= h_1 L \left( W_1 + \frac{(W_{\text{base}} - W_1)}{2} \right) \\ V_2 &= h_2 L \left( W_2 + \frac{(W_1 - W_2)}{2} \right) \\ V_3 &= h_3 L \left( W_3 + \frac{(W_2 - W_3)}{2} \right) \end{aligned} \right. . \qquad (A8)$$

In the current version of MUNICH-hete, it is assumed that traffic emissions are all affected to the recirculation zone. Therefore, its base *W*<sub>base</sub> has to be superior or equal to the street width *W*. If this is not the case, the recirculation zone is considered to fill the whole street volume.

When considered, the other widths are also limited by the street width:

$$\left\{ \begin{aligned} W_1 &= \min(W, W_1) \\ W_2 &= \min(W, W_2) \\ W_3 &= \min(W, W_3) \end{aligned} \right. . \qquad (A9)$$

## Appendix B: Statistical indicators

For evaluation of the simulations compared to observations, the following statistical indicators are used. *o* and *s* represent
the observed and the simulated concentrations respectively. The overbar represents the average.

- Mean fractional error (MFE):

$$MFE = 2 \times \left( \frac{|s - o|}{s + o} \right).$$

- Mean fractional bias (MFB):

$$MFB = 2 \times \left( \frac{s - o}{s + o} \right).$$

- Factor of 2 (FAC2): fraction of data that satisfy $0.5 \le \frac{s}{o} \le 2.0$.

For comparison of the simulations, the normalized mean bias (NMB) and normalized mean error (NME) are used. X represents concentrations, with $X_0$ being the reference simulation and $X_i$ being the compared simulation. The overbar represents the average.

- Normalized mean bias:

$$NMB = \frac{\overline{(X_i - X_0)}}{\overline{X_0}}.$$

- Normalized mean error:

$$NME = \frac{|X_i - X_0|}{\overline{X_0}}.$$
the observed and the simulated concentrations respectively.
The overbar represents the average.

- Mean fractional error (MFE):

$$MFE = 2 \times \left( \frac{|s - o|}{s + o} \right).$$

- Mean fractional bias (MFB):

$$MFB = 2 \times \left( \frac{s - o}{s + o} \right).$$

- Factor of 2 (FAC2): fraction of data that satisfy $0.5 \le \frac{s}{o} \le 2.0$.

For comparison of the simulations, the normalized mean bias (NMB) and normalized mean error (NME) are used. X represents concentrations, with $X_0$ being the reference simulation and $X_i$ being the compared simulation. The overbar represents the average.

- Normalized mean bias:

$$NMB = \frac{\overline{(X_i - X_0)}}{\overline{X_0}}.$$

- Normalized mean error:

$$NME = \frac{|X_i - X_0|}{\overline{X_0}}.$$
Appendix C: Additional information for HCAB

<figure><img src="image_13.png" /></figure>

**Figure C1.** Monthly averaged concentrations (in µg m⁻³) of NOₓ (a) and O₃ (b) at HCAB monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The solid red line represents the OSPM receptor that is close to the measurement station, and the solid orange line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.

<table><caption>**Table C1.** Statistical indicators of the evaluation of the hourly simulated concentrations compared to observations at HCAB monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="2">Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>285.35</td><td>34.64</td><td>72.30</td><td>44.33</td></tr><tr><td>MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>210.88</td><td>25.34</td><td>43.77</td><td>50.44</td></tr><tr><td rowspan="4">MUNICH-hete-11</td><td>MFE</td><td>0.36</td><td>0.48</td><td>0.55</td><td>0.38</td></tr><tr><td>MFB</td><td>-0.29</td><td>-0.29</td><td>-0.42</td><td>0.14</td></tr><tr><td>FAC2</td><td>0.91</td><td>0.73</td><td>0.63</td><td>0.83</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>249.29</td><td>32.37</td><td>74.95</td><td>47.34</td></tr><tr><td rowspan="4">MUNICH-hete-12</td><td>MFE</td><td>0.28</td><td>0.35</td><td>0.35</td><td>0.37</td></tr><tr><td>MFB</td><td>-0.14</td><td>-0.06</td><td>0.05</td><td>0.07</td></tr><tr><td>FAC2</td><td>0.96</td><td>0.87</td><td>0.86</td><td>0.84</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>223.39</td><td>27.69</td><td>53.86</td><td>49.36</td></tr><tr><td rowspan="4">MUNICH-hete-13</td><td>MFE</td><td>0.32</td><td>0.42</td><td>0.43</td><td>0.37</td></tr><tr><td>MFB</td><td>-0.24</td><td>-0.21</td><td>-0.24</td><td>0.12</td></tr><tr><td>FAC2</td><td>0.93</td><td>0.80</td><td>0.78</td><td>0.83</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>197.59</td><td>21.87</td><td>32.90</td><td>52.59</td></tr><tr><td rowspan="4">OSPM-R1</td><td>MFE</td><td>0.39</td><td>0.57</td><td>0.73</td><td>0.39</td></tr><tr><td>MFB</td><td>-0.34</td><td>-0.43</td><td>-0.65</td><td>0.20</td></tr><tr><td>FAC2</td><td>0.88</td><td>0.61</td><td>0.42</td><td>0.82</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>272.86</td><td>40.59</td><td>86.60</td><td>39.50</td></tr><tr><td rowspan="4">OSPM-R2</td><td>MFE</td><td>0.26</td><td>0.34</td><td>0.40</td><td>0.35</td></tr><tr><td>MFB</td><td>-0.05</td><td>0.19</td><td>0.21</td><td>-0.11</td></tr><tr><td>FAC2</td><td>0.96</td><td>0.87</td><td>0.82</td><td>0.86</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>227.84</td><td>30.53</td><td>55.43</td><td>45.77</td></tr><tr><td rowspan="4"></td><td>MFE</td><td>0.33</td><td>0.46</td><td>0.61</td><td>0.36</td></tr><tr><td>MFB</td><td>-0.23</td><td>-0.16</td><td>-0.33</td><td>0.05</td></tr><tr><td>FAC2</td><td>0.93</td><td>0.74</td><td>0.61</td><td>0.84</td></tr></tbody></table>
Appendix C: Additional information for HCAB

<figure><img src="image_11.png" /></figure>

**Figure C1.** Monthly averaged concentrations (in µg m⁻³) of NOₓ (a) and O₃ (b) at HCAB monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (l1) is shown with square markers, the intermediate level (l2) is shown with crossed markers, and the top level (l3) is shown with diamond markers. The solid red line represents the OSPM receptor that is close to the measurement station, and the solid orange line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.

<table><caption>**Table C1.** Statistical indicators of the evaluation of the hourly simulated concentrations compared to observations at HCAB monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="2">Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>285.35</td><td>34.64</td><td>72.30</td><td>44.33</td></tr><tr><td>MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>210.88</td><td>25.34</td><td>43.77</td><td>50.44</td></tr><tr><td rowspan="4">MUNICH-hete-11</td><td>MFE</td><td>0.36</td><td>0.48</td><td>0.55</td><td>0.38</td></tr><tr><td>MFB</td><td>-0.29</td><td>-0.29</td><td>-0.42</td><td>0.14</td></tr><tr><td>FAC2</td><td>0.91</td><td>0.73</td><td>0.63</td><td>0.83</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>249.29</td><td>32.37</td><td>74.95</td><td>47.34</td></tr><tr><td rowspan="4">MUNICH-hete-12</td><td>MFE</td><td>0.28</td><td>0.35</td><td>0.35</td><td>0.37</td></tr><tr><td>MFB</td><td>-0.14</td><td>-0.06</td><td>0.05</td><td>0.07</td></tr><tr><td>FAC2</td><td>0.96</td><td>0.87</td><td>0.86</td><td>0.84</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>223.39</td><td>27.69</td><td>53.86</td><td>49.36</td></tr><tr><td rowspan="4">MUNICH-hete-13</td><td>MFE</td><td>0.32</td><td>0.42</td><td>0.43</td><td>0.37</td></tr><tr><td>MFB</td><td>-0.24</td><td>-0.21</td><td>-0.24</td><td>0.12</td></tr><tr><td>FAC2</td><td>0.93</td><td>0.80</td><td>0.78</td><td>0.83</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>197.59</td><td>21.87</td><td>32.90</td><td>52.59</td></tr><tr><td rowspan="4">OSPM-R1</td><td>MFE</td><td>0.39</td><td>0.57</td><td>0.73</td><td>0.39</td></tr><tr><td>MFB</td><td>-0.34</td><td>-0.43</td><td>-0.65</td><td>0.20</td></tr><tr><td>FAC2</td><td>0.88</td><td>0.61</td><td>0.42</td><td>0.82</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>272.86</td><td>40.59</td><td>86.60</td><td>39.50</td></tr><tr><td rowspan="4">OSPM-R2</td><td>MFE</td><td>0.26</td><td>0.34</td><td>0.40</td><td>0.35</td></tr><tr><td>MFB</td><td>-0.05</td><td>0.19</td><td>0.21</td><td>-0.11</td></tr><tr><td>FAC2</td><td>0.96</td><td>0.87</td><td>0.82</td><td>0.86</td></tr><tr><td>Mean concentration [µg m<sup>-3</sup>]</td><td>227.84</td><td>30.53</td><td>55.43</td><td>45.77</td></tr><tr><td rowspan="4"></td><td>MFE</td><td>0.33</td><td>0.46</td><td>0.61</td><td>0.36</td></tr><tr><td>MFB</td><td>-0.23</td><td>-0.16</td><td>-0.33</td><td>0.05</td></tr><tr><td>FAC2</td><td>0.93</td><td>0.74</td><td>0.61</td><td>0.84</td></tr></tbody></table>
## Appendix D: Additional information for JGTV

<figure><img src="image_14.png" /></figure>

**Figure D1.** Monthly averaged concentrations (in µg m⁻³) of NOₓ (a) and O₃ (b) at JGTV monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (l1) is shown with square markers, the intermediate level (l2) is shown with crossed markers, and the top level (l3) is shown with diamond markers. The solid orange line represents the OSPM receptor that is close to the measurement station, and the solid red line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.

<table><caption>**Table D1.** Statistical indicators of the evaluation of the hourly simulated concentrations compared to observations at JGTV monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th></tr></thead><tbody><tr><td>Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>24.50</td><td>46.94</td></tr><tr><td rowspan="4">MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>21.69</td><td>32.10</td></tr><tr><td>MFE</td><td>0.46</td><td>0.53</td></tr><tr><td>MFB</td><td>-0.03</td><td>-0.19</td></tr><tr><td>FAC2</td><td>0.75</td><td>0.66</td></tr><tr><td rowspan="4">MUNICH-hete-l1</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>19.77</td><td>35.99</td></tr><tr><td>MFE</td><td>0.49</td><td>0.51</td></tr><tr><td>MFB</td><td>-0.13</td><td>-0.11</td></tr><tr><td>FAC2</td><td>0.75</td><td>0.68</td></tr><tr><td rowspan="4">MUNICH-hete-l2</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>18.37</td><td>27.78</td></tr><tr><td>MFE</td><td>0.52</td><td>0.58</td></tr><tr><td>MFB</td><td>-0.20</td><td>-0.32</td></tr><tr><td>FAC2</td><td>0.71</td><td>0.61</td></tr><tr><td rowspan="4">MUNICH-hete-l3</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>16.81</td><td>20.10</td></tr><tr><td>MFE</td><td>0.57</td><td>0.74</td></tr><tr><td>MFB</td><td>-0.29</td><td>-0.58</td></tr><tr><td>FAC2</td><td>0.68</td><td>0.46</td></tr><tr><td rowspan="4">OSPM-R1</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>25.60</td><td>40.46</td></tr><tr><td>MFE</td><td>0.53</td><td>0.63</td></tr><tr><td>MFB</td><td>0.15</td><td>0.03</td></tr><tr><td>FAC2</td><td>0.67</td><td>0.56</td></tr><tr><td rowspan="4">OSPM-R2</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>28.41</td><td>47.21</td></tr><tr><td>MFE</td><td>0.41</td><td>0.42</td></tr><tr><td>MFB</td><td>0.25</td><td>0.18</td></tr><tr><td>FAC2</td><td>0.80</td><td>0.79</td></tr></tbody></table>
## Appendix D: Additional information for JGTV

<figure><img src="image_12.png" /><figcaption>**Figure D1.** Monthly averaged concentrations (in µg m⁻³) of NOₓ (a) and O₃ (b) at JGTV monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (l1) is shown with square markers, the intermediate level (l2) is shown with crossed markers, and the top level (l3) is shown with diamond markers. The solid orange line represents the OSPM receptor that is close to the measurement station, and the solid red line represents the second OSPM receptor located on the other side of the street. The observations are in solid black, and the background concentrations are in solid grey.</figcaption></figure>

<table><caption>**Table D1.** Statistical indicators of the evaluation of the hourly simulated concentrations compared to observations at JGTV monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th></tr></thead><tbody><tr><td>Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>24.50</td><td>46.94</td></tr><tr><td rowspan="4">MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>21.69</td><td>32.10</td></tr><tr><td>MFE</td><td>0.46</td><td>0.53</td></tr><tr><td>MFB</td><td>-0.03</td><td>-0.19</td></tr><tr><td>FAC2</td><td>0.75</td><td>0.66</td></tr><tr><td rowspan="4">MUNICH-hete-l1</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>19.77</td><td>35.99</td></tr><tr><td>MFE</td><td>0.49</td><td>0.51</td></tr><tr><td>MFB</td><td>-0.13</td><td>-0.11</td></tr><tr><td>FAC2</td><td>0.75</td><td>0.68</td></tr><tr><td rowspan="4">MUNICH-hete-l2</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>18.37</td><td>27.78</td></tr><tr><td>MFE</td><td>0.52</td><td>0.58</td></tr><tr><td>MFB</td><td>-0.20</td><td>-0.32</td></tr><tr><td>FAC2</td><td>0.71</td><td>0.61</td></tr><tr><td rowspan="4">MUNICH-hete-l3</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>16.81</td><td>20.10</td></tr><tr><td>MFE</td><td>0.57</td><td>0.74</td></tr><tr><td>MFB</td><td>-0.29</td><td>-0.58</td></tr><tr><td>FAC2</td><td>0.68</td><td>0.46</td></tr><tr><td rowspan="4">OSPM-R1</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>25.60</td><td>40.46</td></tr><tr><td>MFE</td><td>0.53</td><td>0.63</td></tr><tr><td>MFB</td><td>0.15</td><td>0.03</td></tr><tr><td>FAC2</td><td>0.67</td><td>0.56</td></tr><tr><td rowspan="4">OSPM-R2</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>28.41</td><td>47.21</td></tr><tr><td>MFE</td><td>0.41</td><td>0.42</td></tr><tr><td>MFB</td><td>0.25</td><td>0.18</td></tr><tr><td>FAC2</td><td>0.80</td><td>0.79</td></tr></tbody></table>
## Appendix E: Additional information for the district of Le Perreux-sur-Marne

<figure><img src="image_15.png" /></figure>

**Figure E1.** Average daily profile of the concentrations of NO (a) and NO₂ (b) over the simulation period at the Boulevard d'Alsace Lorraine monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The observations are in solid black, and the background concentrations are in solid grey.

<table><caption>**Table E1.** Statistical indicators of the evaluation of the hourly simulated concentrations to observations at the Boulevard d'Alsace Lorraine monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th><th>BC</th></tr></thead><tbody><tr><td>Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>52.16</td><td>147.19</td><td>23.42</td><td>12.72</td><td>6.10</td></tr><tr><td rowspan="4">MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>48.49</td><td>74.62</td><td>21.89</td><td>20.06</td><td>1.41</td></tr><tr><td>MFE</td><td>0.36</td><td>0.65</td><td>0.38</td><td>0.53</td><td>1.17</td></tr><tr><td>MFB</td><td>-0.11</td><td>-0.60</td><td>-0.08</td><td>0.38</td><td>-1.16</td></tr><tr><td>FAC2</td><td>0.87</td><td>0.52</td><td>0.84</td><td>0.67</td><td>0.09</td></tr><tr><td rowspan="4">MUNICH-hete-11</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>58.92</td><td>104.74</td><td>24.58</td><td>22.18</td><td>1.84</td></tr><tr><td>MFE</td><td>0.33</td><td>0.41</td><td>0.37</td><td>0.57</td><td>1.01</td></tr><tr><td>MFB</td><td>0.08</td><td>-0.30</td><td>0.03</td><td>0.47</td><td>-1.00</td></tr><tr><td>FAC2</td><td>0.91</td><td>0.80</td><td>0.86</td><td>0.61</td><td>0.16</td></tr><tr><td rowspan="4">MUNICH-hete-12</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>47.81</td><td>75.36</td><td>21.94</td><td>20.08</td><td>1.42</td></tr><tr><td>MFE</td><td>0.36</td><td>0.64</td><td>0.38</td><td>0.53</td><td>1.17</td></tr><tr><td>MFB</td><td>-0.12</td><td>-0.59</td><td>-0.08</td><td>0.38</td><td>-1.16</td></tr><tr><td>FAC2</td><td>0.86</td><td>0.53</td><td>0.84</td><td>0.67</td><td>0.08</td></tr><tr><td rowspan="4">MUNICH-hete-13</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>37.97</td><td>52.63</td><td>19.91</td><td>18.41</td><td>1.09</td></tr><tr><td>MFE</td><td>0.49</td><td>0.90</td><td>0.42</td><td>0.50</td><td>1.30</td></tr><tr><td>MFB</td><td>-0.35</td><td>-0.88</td><td>-0.18</td><td>0.30</td><td>-1.30</td></tr><tr><td>FAC2</td><td>0.71</td><td>0.25</td><td>0.80</td><td>0.70</td><td>0.06</td></tr></tbody></table>
## Appendix E: Additional information for the district of Le Perreux-sur-Marne

<figure><img src="image_13.png" /></figure>

**Figure E1.** Average daily profile of the concentrations of NO (a) and NO₂ (b) over the simulation period at the Boulevard d'Alsace Lorraine monitoring station. The solid blue line represents the homogeneous version of MUNICH. The three dashed green lines represent the three levels of the heterogeneous version of MUNICH – the lowest level (11) is shown with square markers, the intermediate level (12) is shown with crossed markers, and the top level (13) is shown with diamond markers. The observations are in solid black, and the background concentrations are in solid grey.

<table><caption>**Table E1.** Statistical indicators of the evaluation of the hourly simulated concentrations to observations at the Boulevard d'Alsace Lorraine monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th></th><th></th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>PM<sub>10</sub></th><th>PM<sub>2.5</sub></th><th>BC</th></tr></thead><tbody><tr><td rowspan="2">Observations</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>52.16</td><td>147.19</td><td>23.42</td><td>12.72</td><td>6.10</td></tr><tr><td>MUNICH-homo</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>48.49</td><td>74.62</td><td>21.89</td><td>20.06</td><td>1.41</td></tr><tr><td rowspan="4">MUNICH-hete-11</td><td>MFE</td><td>0.36</td><td>0.65</td><td>0.38</td><td>0.53</td><td>1.17</td></tr><tr><td>MFB</td><td>-0.11</td><td>-0.60</td><td>-0.08</td><td>0.38</td><td>-1.16</td></tr><tr><td>FAC2</td><td>0.87</td><td>0.52</td><td>0.84</td><td>0.67</td><td>0.09</td></tr><tr><td>MUNICH-hete-12</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>58.92</td><td>104.74</td><td>24.58</td><td>22.18</td><td>1.84</td></tr><tr><td rowspan="4">MUNICH-hete-13</td><td>MFE</td><td>0.33</td><td>0.41</td><td>0.37</td><td>0.57</td><td>1.01</td></tr><tr><td>MFB</td><td>0.08</td><td>-0.30</td><td>0.03</td><td>0.47</td><td>-1.00</td></tr><tr><td>FAC2</td><td>0.91</td><td>0.80</td><td>0.86</td><td>0.61</td><td>0.16</td></tr><tr><td>MUNICH-hete-14</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>47.81</td><td>75.36</td><td>21.94</td><td>20.08</td><td>1.42</td></tr><tr><td rowspan="4">MUNICH-hete-15</td><td>MFE</td><td>0.36</td><td>0.64</td><td>0.38</td><td>0.53</td><td>1.17</td></tr><tr><td>MFB</td><td>-0.12</td><td>-0.59</td><td>-0.08</td><td>0.38</td><td>-1.16</td></tr><tr><td>FAC2</td><td>0.86</td><td>0.53</td><td>0.84</td><td>0.67</td><td>0.08</td></tr><tr><td>MUNICH-hete-16</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>37.97</td><td>52.63</td><td>19.91</td><td>18.41</td><td>1.09</td></tr><tr><td rowspan="4">MUNICH-hete-17</td><td>MFE</td><td>0.49</td><td>0.90</td><td>0.42</td><td>0.50</td><td>1.30</td></tr><tr><td>MFB</td><td>-0.35</td><td>-0.88</td><td>-0.18</td><td>0.30</td><td>-1.30</td></tr><tr><td>FAC2</td><td>0.71</td><td>0.25</td><td>0.80</td><td>0.70</td><td>0.06</td></tr></tbody></table>
<figure><img src="image_16.png" /></figure>

**Figure E2.** NO₂ and PM₁₀ time-averaged concentrations (in µg m⁻³) for MUNICH-homo in the left panel (a, c) for the district of Le Perreux-sur-Marne. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for NO₂ and PM₁₀ in the right panel (b, d) for the district of Le Perreux-sur-Marne. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.
<figure><img src="image_14.png" /></figure>

**Figure E2.** NO₂ and PM₁₀ time-averaged concentrations (in µg m⁻³) for MUNICH-homo in the left panel (a, c) for the district of Le Perreux-sur-Marne. Relative differences (in %) between the first level of MUNICH-hete and MUNICH-homo for NO₂ and PM₁₀ in the right panel (b, d) for the district of Le Perreux-sur-Marne. A positive relative difference indicates higher concentrations for MUNICH-hete. The brown cross in panels (a) and (c) represents the position of the measurement station.
## Appendix F: Additional information for the influence of the aspect ratio

<table><caption>**Table F1.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network around HCAB. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">CO</td><td>[0, 0.5[</td><td>201.26</td><td>11.61</td><td></td></tr><tr><td>[0.5, 1[</td><td>196.24</td><td>4.71</td><td></td></tr><tr><td>[1, 1.5[</td><td>192.05</td><td>3.12</td><td></td></tr><tr><td>≥ 1.5</td><td>193.23</td><td>2.12</td><td></td></tr><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5[</td><td>23.18</td><td>14.46</td><td></td></tr><tr><td>[0.5, 1[</td><td>21.69</td><td>-0.09</td><td></td></tr><tr><td>[1, 1.5[</td><td>20.62</td><td>-1.24</td><td></td></tr><tr><td>≥ 1.5</td><td>20.56</td><td>-3.27</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5[</td><td>34.84</td><td>51.16</td><td></td></tr><tr><td>[0.5, 1[</td><td>29.67</td><td>21.34</td><td></td></tr><tr><td>[1, 1.5[</td><td>26.40</td><td>15.84</td><td></td></tr><tr><td>≥ 1.5</td><td>25.87</td><td>9.29</td><td></td></tr><tr><td rowspan="4">O<sub>3</sub></td><td>[0, 0.5[</td><td>51.48</td><td>-2.08</td><td></td></tr><tr><td>[0.5, 1[</td><td>52.33</td><td>1.68</td><td></td></tr><tr><td>[1, 1.5[</td><td>53.01</td><td>1.57</td><td></td></tr><tr><td>≥ 1.5</td><td>53.00</td><td>1.94</td><td></td></tr></tbody></table>

<table><caption>**Table F2.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network around JGTV. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">CO</td><td>[0, 0.5[</td><td>180.72</td><td>2.82</td><td></td></tr><tr><td>[0.5, 1[</td><td>181.49</td><td>2.06</td><td></td></tr><tr><td>[1, 1.5[</td><td>179.85</td><td>1.28</td><td></td></tr><tr><td>≥ 1.5</td><td>176.26</td><td>0.26</td><td></td></tr><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5[</td><td>17.06</td><td>3.76</td><td></td></tr><tr><td>[0.5, 1[</td><td>17.62</td><td>-0.75</td><td></td></tr><tr><td>[1, 1.5[</td><td>17.31</td><td>-2.23</td><td></td></tr><tr><td>≥ 1.5</td><td>16.14</td><td>-1.69</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5[</td><td>21.59</td><td>16.75</td><td></td></tr><tr><td>[0.5, 1[</td><td>22.02</td><td>11.23</td><td></td></tr><tr><td>[1, 1.5[</td><td>20.67</td><td>7.64</td><td></td></tr><tr><td>≥ 1.5</td><td>18.30</td><td>1.91</td><td></td></tr><tr><td rowspan="4">O<sub>3</sub></td><td>[0, 0.5[</td><td>56.06</td><td>-0.31</td><td></td></tr><tr><td>[0.5, 1[</td><td>55.54</td><td>0.85</td><td></td></tr><tr><td>[1, 1.5[</td><td>55.67</td><td>1.11</td><td></td></tr><tr><td>≥ 1.5</td><td>56.57</td><td>0.58</td><td></td></tr></tbody></table>
## Appendix F: Additional information for the influence of the aspect ratio

<table><caption>**Table F1.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network around HCAB. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">CO</td><td>[0, 0.5[</td><td>201.26</td><td>11.61</td><td></td></tr><tr><td>[0.5, 1[</td><td>196.24</td><td>4.71</td><td></td></tr><tr><td>[1, 1.5[</td><td>192.05</td><td>3.12</td><td></td></tr><tr><td>≥ 1.5</td><td>193.23</td><td>2.12</td><td></td></tr><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5[</td><td>23.18</td><td>14.46</td><td></td></tr><tr><td>[0.5, 1[</td><td>21.69</td><td>-0.09</td><td></td></tr><tr><td>[1, 1.5[</td><td>20.62</td><td>-1.24</td><td></td></tr><tr><td>≥ 1.5</td><td>20.56</td><td>-3.27</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5[</td><td>34.84</td><td>51.16</td><td></td></tr><tr><td>[0.5, 1[</td><td>29.67</td><td>21.34</td><td></td></tr><tr><td>[1, 1.5[</td><td>26.40</td><td>15.84</td><td></td></tr><tr><td>≥ 1.5</td><td>25.87</td><td>9.29</td><td></td></tr><tr><td rowspan="4">O<sub>3</sub></td><td>[0, 0.5[</td><td>51.48</td><td>-2.08</td><td></td></tr><tr><td>[0.5, 1[</td><td>52.33</td><td>1.68</td><td></td></tr><tr><td>[1, 1.5[</td><td>53.01</td><td>1.57</td><td></td></tr><tr><td>≥ 1.5</td><td>53.00</td><td>1.94</td><td></td></tr></tbody></table>

<table><caption>**Table F2.** Statistical indicators of the influence of the aspect ratio on concentrations in MUNICH-hete-11 for the street network around JGTV. Indicators are presented in Appendix B.</caption><thead><tr><th rowspan="2">Pollutant</th><th rowspan="2">a<sub>r</sub></th><th colspan="2">Mean concentration<br/>[µg m<sup>-3</sup>]</th><th rowspan="2">NMB<br/>[%]</th></tr><tr><th>MUNICH-homo</th><th>MUNICH-hete-11</th></tr></thead><tbody><tr><td rowspan="4">CO</td><td>[0, 0.5[</td><td>180.72</td><td>2.82</td><td></td></tr><tr><td>[0.5, 1[</td><td>181.49</td><td>2.06</td><td></td></tr><tr><td>[1, 1.5[</td><td>179.85</td><td>1.28</td><td></td></tr><tr><td>≥ 1.5</td><td>176.26</td><td>0.26</td><td></td></tr><tr><td rowspan="4">NO<sub>2</sub></td><td>[0, 0.5[</td><td>17.06</td><td>3.76</td><td></td></tr><tr><td>[0.5, 1[</td><td>17.62</td><td>-0.75</td><td></td></tr><tr><td>[1, 1.5[</td><td>17.31</td><td>-2.23</td><td></td></tr><tr><td>≥ 1.5</td><td>16.14</td><td>-1.69</td><td></td></tr><tr><td rowspan="4">NO<sub>x</sub></td><td>[0, 0.5[</td><td>21.59</td><td>16.75</td><td></td></tr><tr><td>[0.5, 1[</td><td>22.02</td><td>11.23</td><td></td></tr><tr><td>[1, 1.5[</td><td>20.67</td><td>7.64</td><td></td></tr><tr><td>≥ 1.5</td><td>18.30</td><td>1.91</td><td></td></tr><tr><td rowspan="4">O<sub>3</sub></td><td>[0, 0.5[</td><td>56.06</td><td>-0.31</td><td></td></tr><tr><td>[0.5, 1[</td><td>55.54</td><td>0.85</td><td></td></tr><tr><td>[1, 1.5[</td><td>55.67</td><td>1.11</td><td></td></tr><tr><td>≥ 1.5</td><td>56.57</td><td>0.58</td><td></td></tr></tbody></table>
## Appendix G: Additional information for the sensitivity to street network

CC were responsible for funding acquisition. TS and KS

<table><caption>**Table G1.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment HCAB with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th>MUNICH-hete level</th><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="4">11</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>249.29</td><td>32.37</td><td>74.95</td><td>47.34</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>214.27</td><td>23.28</td><td>46.69</td><td>52.98</td></tr><tr><td>NME [%]</td><td>14.05</td><td>28.06</td><td>37.70</td><td>11.90</td></tr><tr><td>NMB [%]</td><td>-14.05</td><td>-28.06</td><td>-37.70</td><td>11.90</td></tr><tr><td rowspan="4">12</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>223.39</td><td>27.69</td><td>53.86</td><td>49.36</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>198.12</td><td>20.78</td><td>33.47</td><td>53.79</td></tr><tr><td>NME [%]</td><td>11.32</td><td>24.95</td><td>37.86</td><td>8.99</td></tr><tr><td>NMB [%]</td><td>-11.32</td><td>-24.95</td><td>-37.86</td><td>8.99</td></tr><tr><td rowspan="4">13</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>197.59</td><td>21.87</td><td>32.90</td><td>52.59</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>186.96</td><td>18.65</td><td>24.33</td><td>54.78</td></tr><tr><td>NME [%]</td><td>5.38</td><td>14.73</td><td>26.05</td><td>4.17</td></tr><tr><td>NMB [%]</td><td>-5.38</td><td>-14.73</td><td>-26.05</td><td>4.17</td></tr></tbody></table>

<table><caption>**Table G2.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment JGTV with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th>MUNICH-hete level</th><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="4">11</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>201.21</td><td>19.77</td><td>35.99</td><td>55.19</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>189.24</td><td>17.46</td><td>27.56</td><td>56.45</td></tr><tr><td>NME [%]</td><td>6.09</td><td>12.38</td><td>23.63</td><td>3.13</td></tr><tr><td>NMB [%]</td><td>-5.95</td><td>-11.70</td><td>-23.43</td><td>2.29</td></tr><tr><td rowspan="4">12</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>189.49</td><td>18.37</td><td>27.78</td><td>55.54</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>181.01</td><td>16.52</td><td>21.80</td><td>56.64</td></tr><tr><td>NME [%]</td><td>4.71</td><td>10.87</td><td>21.85</td><td>2.85</td></tr><tr><td>NMB [%]</td><td>-4.47</td><td>-10.03</td><td>-21.53</td><td>1.99</td></tr><tr><td rowspan="4">13</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>178.59</td><td>16.81</td><td>20.10</td><td>56.12</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>176.03</td><td>15.92</td><td>18.30</td><td>56.80</td></tr><tr><td>NME [%]</td><td>1.84</td><td>6.16</td><td>9.57</td><td>1.83</td></tr><tr><td>NMB [%]</td><td>-1.43</td><td>-5.30</td><td>-8.92</td><td>1.21</td></tr></tbody></table>

*Code and data availability.* MUNICH-hete is available at https://doi.org/10.5281/zenodo.7778271 (Sarica et al., 2023a). The configuration files, the input data, and the scripts to generate the figures and statistics are available at Sarica et al. (2023a).

*Author contributions.* TS, AM, KS, and YR conceptualized the development of the study. TS, KS, and YK developed the software, and TS performed the simulations. MK and SSJ provided help for the simulations with OSPM and the data for the Danish cases. TS, KS, and YR performed the analysis. KS, YR, MK, and SSJ supervised the study. KS and

wrote the original draft, and all the authors reviewed and edited the final draft.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.
## Appendix G: Additional information for the sensitivity to street network

CC were responsible for funding acquisition. TS and KS

<table><caption>**Table G1.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment HCAB with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th>MUNICH-hete level</th><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="4">11</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>249.29</td><td>32.37</td><td>74.95</td><td>47.34</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>214.27</td><td>23.28</td><td>46.69</td><td>52.98</td></tr><tr><td>NME [%]</td><td>14.05</td><td>28.06</td><td>37.70</td><td>11.90</td></tr><tr><td>NMB [%]</td><td>-14.05</td><td>-28.06</td><td>-37.70</td><td>11.90</td></tr><tr><td rowspan="4">12</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>223.39</td><td>27.69</td><td>53.86</td><td>49.36</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>198.12</td><td>20.78</td><td>33.47</td><td>53.79</td></tr><tr><td>NME [%]</td><td>11.32</td><td>24.95</td><td>37.86</td><td>8.99</td></tr><tr><td>NMB [%]</td><td>-11.32</td><td>-24.95</td><td>-37.86</td><td>8.99</td></tr><tr><td rowspan="4">13</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>197.59</td><td>21.87</td><td>32.90</td><td>52.59</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>186.96</td><td>18.65</td><td>24.33</td><td>54.78</td></tr><tr><td>NME [%]</td><td>5.38</td><td>14.73</td><td>26.05</td><td>4.17</td></tr><tr><td>NMB [%]</td><td>-5.38</td><td>-14.73</td><td>-26.05</td><td>4.17</td></tr></tbody></table>

<table><caption>**Table G2.** Statistical indicators of the influence of the street network on concentrations simulated with MUNICH-hete for the street segment JGTV with the monitoring station. Indicators are presented in Appendix B.</caption><thead><tr><th>MUNICH-hete level</th><th></th><th></th><th>CO</th><th>NO<sub>2</sub></th><th>NO<sub>x</sub></th><th>O<sub>3</sub></th></tr></thead><tbody><tr><td rowspan="4">11</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>201.21</td><td>19.77</td><td>35.99</td><td>55.19</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>189.24</td><td>17.46</td><td>27.56</td><td>56.45</td></tr><tr><td>NME [%]</td><td>6.09</td><td>12.38</td><td>23.63</td><td>3.13</td></tr><tr><td>NMB [%]</td><td>-5.95</td><td>-11.70</td><td>-23.43</td><td>2.29</td></tr><tr><td rowspan="4">12</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>189.49</td><td>18.37</td><td>27.78</td><td>55.54</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>181.01</td><td>16.52</td><td>21.80</td><td>56.64</td></tr><tr><td>NME [%]</td><td>4.71</td><td>10.87</td><td>21.85</td><td>2.85</td></tr><tr><td>NMB [%]</td><td>-4.47</td><td>-10.03</td><td>-21.53</td><td>1.99</td></tr><tr><td rowspan="4">13</td><td>With network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>178.59</td><td>16.81</td><td>20.10</td><td>56.12</td></tr><tr><td rowspan="2">Without network</td><td>Mean concentration [µg m<sup>-3</sup>]</td><td>176.03</td><td>15.92</td><td>18.30</td><td>56.80</td></tr><tr><td>NME [%]</td><td>1.84</td><td>6.16</td><td>9.57</td><td>1.83</td></tr><tr><td>NMB [%]</td><td>-1.43</td><td>-5.30</td><td>-8.92</td><td>1.21</td></tr></tbody></table>

*Code and data availability.* MUNICH-hete is available at https://doi.org/10.5281/zenodo.7778271 (Sarica et al., 2023a). The configuration files, the input data, and the scripts to generate the figures and statistics are available at Sarica et al. (2023a).

*Author contributions.* TS, AM, KS, and YR conceptualized the development of the study. TS, KS, and YK developed the software, and TS performed the simulations. MK and SSJ provided help for the simulations with OSPM and the data for the Danish cases. TS, KS, and YR performed the analysis. KS, YR, MK, and SSJ supervised the study. KS and

wrote the original draft, and all the authors reviewed and edited the final draft.

*Competing interests.* The contact author has declared that none of the authors has any competing interests.
*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Special issue statement.* This article is part of the special issue "Air quality research at street level – Part II (ACP/GMD inter-journal SI)". It is not associated with a conference.

*Review statement.* This paper was edited by Christoph Knote and reviewed by two anonymous referees.

## References

Amato, F., Pérez, N., López, M., Ripoll, A., Alastuey, A., Pandolfi, M., Karanasiou, A., Salmatonidis, A., Padoan, E., Frasca, D., Marcoccia, M., Viana, M., Moreno, T., Reche, C., Martins, V., Brines, M., Minguillón, M., Ealo, M., Rivas, I., van Drooge, B., Benavides, J., Craviotto, J., and Querol, X.: Vertical and horizontal fall-off of black carbon and $NO_2$ within urban blocks, Sci. Total Environ., 686, 236–245, https://doi.org/10.1016/j.scitotenv.2019.05.434, 2019.

Archambeau, F., Méchitoua, N., and Sakiz, M.: Code Saturne: A Finite Volume Code for the computation of turbulent incompressible flows – Industrial Applications, International Journal on Finite Volumes, 1, https://hal.science/hal-01115371 (last access: 12 September 2023), 2004.

Atkinson, R.: Atmospheric chemistry of VOCs and $NO_x$, Atmos. Environ., 34, 2063–2101, https://doi.org/10.1016/S1352-2310(99)00460-4, 2000.

Benavides, J., Snyder, M., Guevara, M., Soret, A., Pérez García-Pando, C., Amato, F., Querol, X., and Jorba, O.: CALIOPE-Urban v1.0: coupling R-LINE with a mesoscale air quality modelling system for urban air quality forecasts over Barcelona city (Spain), Geosci. Model Dev., 12, 2811–2835, https://doi.org/10.5194/gmd-12-2811-2019, 2019.

Berkowicz, R.: A simple model for urban background pollution, Environ. Monit. Assess., 65, 259–267, https://doi.org/10.1023/a:1006466025186, 2000a.

Berkowicz, R.: OSPM – a parameterised street pollution model, Environ. Monit. Assess., 65, 323–331, https://doi.org/10.1023/A:1006448321977, 2000b.

Berkowicz, R., Hertel, O., Larsen, S., Sørensen, N., and Nielsen, M.: Modelling traffic pollution in streets, https://backend.orbit.dtu.dk/ws/portalfiles/portal/128001317/Modelling_traffic_pollution_in_streets.pdf (last access: 13 February 2023), 1997.

Cimorelli, A., Perry, S., Venkatram, A., Weil, J., Paine, R., Wilson, R., Lee, R., Peters, W., Brode, R., and Paumier, J.: AERMOD: description of model formulation, US Environmental Protection Agency, Tech. rep., EPA-454/R-03-004, https://gaftp.epa.gov/Air/aqmg/SCRAM/models/preferred/aermod/aermod_mfd_454-R-03-004.pdf (last access: 12 September 2023), 2004.

Coceal, O. and Belcher, S.: A canopy model of mean winds through urban areas, Q. J. Roy. Meteor. Soc., 130, 1349–1372, https://doi.org/10.1256/qj.03.40, 2004.

Collett, R. S. and Oduyemi, K.: Air quality modelling: a technical review of mathematical approaches, Meteorol. Appl., 4, 235–246, https://doi.org/10.1017/S1350482797000455, 1997.

Conti, G. O., Heibati, B., Kloog, I., Fiore, M., and Ferrante, M.: A review of AirQ Models and their applications for forecasting the air pollution health outcomes, Environ. Sci. Pollut. R., 24, 6426–6445, https://doi.org/10.1007/s11356-016-8180-1, 2017.

Dai, Y., Cai, X., Zhong, J., and MacKenzie, A. R.: Modelling chemistry and transport in urban street canyons: Comparing offline multi-box models with large-eddy simulation, Atmos. Environ., 264, 118709, https://doi.org/10.1016/j.atmosenv.2021.118709, 2021.

de la Paz, D., Borge, R., Vedrenne, M., Lumbreras, J., Amato, F., Karanasiou, A., Boldo, E., and Moreno, T.: Implementation of road dust resuspension in air quality simulations of particulate matter in Madrid (Spain), Front. Environ. Sci., 3, https://doi.org/10.3389/fenvs.2015.00072, 2015.

El-Harbawi, M.: Air quality modelling, simulation, and computational methods: a review, Environ. Rev., 21, 149–179, https://doi.org/10.1139/er-2012-0056, 2013.

Falasca, S. and Curci, G.: High-resolution air quality modelling: Sensitivity tests to horizontal resolution and urban canopy with WRF-CHIMERE, Atmos. Environ., 187, 241–254, https://doi.org/10.1016/j.atmosenv.2018.05.048, 2018.

Fuller, R., Landrigan, P. J., Balakrishnan, K., Bathan, G., Bose-O'Reilly, S., Brauer, M., Caravanos, J., Chiles, T., Cohen, A., Corra, L., Cropper, M., Ferraro, G., Hanna, J., Hanrahan, D., Hu, H., Hunter, D., Janata, G., Kupka, R., Lanphear, B., Lichtveld, M., Martin, K., Mustapha, A., Sanchez-Triana, E., Sandilya, K., Schaefli, L., Shaw, J., Seddon, J., Suk, W., Téllez-Rojo, M. M., and Yan, C.: Pollution and health: a progress update, Lancet Planet. Health, 6, e535–e547, https://doi.org/10.1016/S2542-5196(22)00090-0, 2022.

Gao, Z., Bresson, R., Qu, Y., Milliez, M., de Munck, C., and Carissimo, B.: High resolution unsteady RANS simulation of wind, thermal effects and pollution dispersion for studying urban renewal scenarios in a neighborhood of Toulouse, Urban Clim., 23, 114–130, https://doi.org/10.1016/j.uclim.2016.11.002, 2018.

Hood, C., Stocker, J., Seaton, M., Johnson, K., O'Neill, J., Thorne, L., and Carruthers, D.: Comprehensive evaluation of an advanced street canyon air pollution model, J. Air Waste Manage., 71, 247–267, https://doi.org/10.1080/10962247.2020.1803158, 2021.

Jensen, S. S., Ketzel, M., Brandt, J., Becker, T., Plejdrup, M., Winther, M., Ellermann, T., Christensen, J. H., Nielsen, O.-K., Hertel, O., and Fuglsang, M. W.: 04 – Air Quality at Your Street – Public Digital Map of Air Quality in Denmark, in: Proceedings, COST Association, Brussels, Belgium, Sixth Scientific Meeting EuNetAir 5–7 October 2016, Academy of Sciences, Prague, Czech Republic, 14–17, https://doi.org/10.5162/6eunetair2016/04, 2016.

Karl, M., Walker, S.-E., Solberg, S., and Ramacher, M. O. P.: The Eulerian urban dispersion model EPISODE – Part 2: Extensions to the source dispersion and photochemistry for EPISODE-CityChem v1.2 and its application to the city of Hamburg, Geosci. Model Dev., 12, 3357–3399, https://doi.org/10.5194/gmd-12-3357-2019, 2019.

Khan, S. and Quamrul, H.: Review of developments in air quality modelling and air quality dispersion models, J. Environ. Eng. Sci., 16, 1–10, https://doi.org/10.1680/jenes.20.00004, 2021.
*Disclaimer.* Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

*Special issue statement.* This article is part of the special issue "Air quality research at street level – Part II (ACP/GMD inter-journal SI)". It is not associated with a conference.

*Review statement.* This paper was edited by Christoph Knote and reviewed by two anonymous referees.

## References

Amato, F., Pérez, N., López, M., Ripoll, A., Alastuey, A., Pandolfi, M., Karanasiou, A., Salmatonidis, A., Padoan, E., Frasca, D., Marcoccia, M., Viana, M., Moreno, T., Reche, C., Martins, V., Brines, M., Minguillón, M., Ealo, M., Rivas, I., van Drooge, B., Benavides, J., Craviotto, J., and Querol, X.: Vertical and horizontal fall-off of black carbon and $NO_2$ within urban blocks, Sci. Total Environ., 686, 236–245, https://doi.org/10.1016/j.scitotenv.2019.05.434, 2019.

Archambeau, F., Méchitoua, N., and Sakiz, M.: Code Saturne: A Finite Volume Code for the computation of turbulent incompressible flows – Industrial Applications, International Journal on Finite Volumes, 1, https://hal.science/hal-01115371 (last access: 12 September 2023), 2004.

Atkinson, R.: Atmospheric chemistry of VOCs and $NO_x$, Atmos. Environ., 34, 2063–2101, https://doi.org/10.1016/S1352-2310(99)00460-4, 2000.

Benavides, J., Snyder, M., Guevara, M., Soret, A., Pérez García-Pando, C., Amato, F., Querol, X., and Jorba, O.: CALIOPE-Urban v1.0: coupling R-LINE with a mesoscale air quality modelling system for urban air quality forecasts over Barcelona city (Spain), Geosci. Model Dev., 12, 2811–2835, https://doi.org/10.5194/gmd-12-2811-2019, 2019.

Berkowicz, R.: A simple model for urban background pollution, Environ. Monit. Assess., 65, 259–267, https://doi.org/10.1023/a:1006466025186, 2000a.

Berkowicz, R.: OSPM – a parameterised street pollution model, Environ. Monit. Assess., 65, 323–331, https://doi.org/10.1023/A:1006448321977, 2000b.

Berkowicz, R., Hertel, O., Larsen, S., Sørensen, N., and Nielsen, M.: Modelling traffic pollution in streets, https://backend.orbit.dtu.dk/ws/portalfiles/portal/128001317/Modelling_traffic_pollution_in_streets.pdf (last access: 13 February 2023), 1997.

Cimorelli, A., Perry, S., Venkatram, A., Weil, J., Paine, R., Wilson, R., Lee, R., Peters, W., Brode, R., and Paumier, J.: AERMOD: description of model formulation, US Environmental Protection Agency, Tech. rep., EPA-454/R-03-004, https://gaftp.epa.gov/Air/aqmg/SCRAM/models/preferred/aermod/aermod_mfd_454-R-03-004.pdf (last access: 12 September 2023), 2004.

Coceal, O. and Belcher, S.: A canopy model of mean winds through urban areas, Q. J. Roy. Meteor. Soc., 130, 1349–1372, https://doi.org/10.1256/qj.03.40, 2004.

Collett, R. S. and Oduyemi, K.: Air quality modelling: a technical review of mathematical approaches, Meteorol. Appl., 4, 235–246, https://doi.org/10.1017/S1350482797000455, 1997.

Conti, G. O., Heibati, B., Kloog, I., Fiore, M., and Ferrante, M.: A review of AirQ Models and their applications for forecasting the air pollution health outcomes, Environ. Sci. Pollut. R., 24, 6426–6445, https://doi.org/10.1007/s11356-016-8180-1, 2017.

Dai, Y., Cai, X., Zhong, J., and MacKenzie, A. R.: Modelling chemistry and transport in urban street canyons: Comparing offline multi-box models with large-eddy simulation, Atmos. Environ., 264, 118709, https://doi.org/10.1016/j.atmosenv.2021.118709, 2021.

de la Paz, D., Borge, R., Vedrenne, M., Lumbreras, J., Amato, F., Karanasiou, A., Boldo, E., and Moreno, T.: Implementation of road dust resuspension in air quality simulations of particulate matter in Madrid (Spain), Front. Environ. Sci., 3, https://doi.org/10.3389/fenvs.2015.00072, 2015.

El-Harbawi, M.: Air quality modelling, simulation, and computational methods: a review, Environ. Rev., 21, 149–179, https://doi.org/10.1139/er-2012-0056, 2013.

Falasca, S. and Curci, G.: High-resolution air quality modeling: Sensitivity tests to horizontal resolution and urban canopy with WRF-CHIMERE, Atmos. Environ., 187, 241–254, https://doi.org/10.1016/j.atmosenv.2018.05.048, 2018.

Fuller, R., Landrigan, P. J., Balakrishnan, K., Bathan, G., Bose-O'Reilly, S., Brauer, M., Caravanos, J., Chiles, T., Cohen, A., Corra, L., Cropper, M., Ferraro, G., Hanna, J., Hanrahan, D., Hu, H., Hunter, D., Janata, G., Kupka, R., Lanphear, B., Lichtveld, M., Martin, K., Mustapha, A., Sanchez-Triana, E., Sandilya, K., Schaefli, L., Shaw, J., Seddon, J., Suk, W., Téllez-Rojo, M. M., and Yan, C.: Pollution and health: a progress update, Lancet Planet. Health, 6, e535–e547, https://doi.org/10.1016/S2542-5196(22)00090-0, 2022.

Gao, Z., Bresson, R., Qu, Y., Milliez, M., de Munck, C., and Carissimo, B.: High resolution unsteady RANS simulation of wind, thermal effects and pollution dispersion for studying urban renewal scenarios in a neighborhood of Toulouse, Urban Clim., 23, 114–130, https://doi.org/10.1016/j.uclim.2016.11.002, 2018.

Hood, C., Stocker, J., Seaton, M., Johnson, K., O'Neill, J., Thorne, L., and Carruthers, D.: Comprehensive evaluation of an advanced street canyon air pollution model, J. Air Waste Manage., 71, 247–267, https://doi.org/10.1080/10962247.2020.1803158, 2021.

Jensen, S. S., Ketzel, M., Brandt, J., Becker, T., Plejdrup, M., Winther, M., Ellermann, T., Christensen, J. H., Nielsen, O.-K., Hertel, O., and Fuglsang, M. W.: 04 – Air Quality at Your Street – Public Digital Map of Air Quality in Denmark, in: Proceedings, COST Association, Brussels, Belgium, Sixth Scientific Meeting EuNetAir 5–7 October 2016, Academy of Sciences, Prague, Czech Republic, 14–17, https://doi.org/10.5162/6eunetair2016/04, 2016.

Karl, M., Walker, S.-E., Solberg, S., and Ramacher, M. O. P.: The Eulerian urban dispersion model EPISODE – Part 2: Extensions to the source dispersion and photochemistry for EPISODE-CityChem v1.2 and its application to the city of Hamburg, Geosci. Model Dev., 12, 3357–3399, https://doi.org/10.5194/gmd-12-3357-2019, 2019.

Khan, S. and Quamrul, H.: Review of developments in air quality modelling and air quality dispersion models, J. Environ. Eng. Sci., 16, 1–10, https://doi.org/10.1680/jenes.20.00004, 2021.
Kim, Y., Wu, Y., Seigneur, C., and Roustan, Y.: Multi-scale modeling of urban air pollution: development and application of a Street-in-Grid model (v1.0) by coupling MUNICH (v1.0) and Polair3D (v1.8.1), Geosci. Model Dev., 11, 611–629, https://doi.org/10.5194/gmd-11-611-2018, 2018.

Kim, Y., Lugon, L., Maison, A., Sarica, T., Roustan, Y., Valari, M., Zhang, Y., André, M., and Sartelet, K.: MUNICH v2.0: a street-network model coupled with SSH-aerosol (v1.2) for multi-pollutant modelling, Geosci. Model Dev., 15, 7371–7396, https://doi.org/10.5194/gmd-15-7371-2022, 2022.

Kwak, K.-H. and Baik, J.-J.: A CFD modeling study of the impacts of NOₓ and VOC emissions on reactive pollutant dispersion in and above a street canyon, Atmos. Environ., 46, 71–80, https://doi.org/10.1016/j.atmosenv.2011.10.024, 2012.

Landrigan, P. J., Fuller, R., Acosta, N. J. R., Adeyi, O., Arnold, R., Basu, N. N., Baldé, A. B., Bertollini, R., Bose-O'Reilly, S., Boufford, J. I., Breysse, P. N., Chiles, T., Mahidol, C., Coll-Seck, A. M., Cropper, M. L., Fobil, J., Fuster, V., Greenstone, M., Haines, A., Hanrahan, D., Hunter, D., Khare, M., Krupnick, A., Lanphear, B., Lohani, B., Martin, K., Mathiasen, K. V., McTeer, M. A., Murray, C. J. L., Ndahimananjara, J. D., Perera, F., Potočnik, J., Preker, A. S., Ramesh, J., Rockström, J., Salinas, C., Samson, L. D., Sandilya, K., Sly, P. D., Smith, K. R., Steiner, A., Stewart, R. B., Suk, W. A., van Schayck, O. C. P., Yadama, G. N., Yumkella, K., and Zhong, M.: The Lancet Commission on pollution and health, Lancet, 391, 462–512, https://doi.org/10.1016/S0140-6736(17)32345-0, 2018.

Lateb, M., Meroney, R., Yataghene, M., Fellouah, H., Saleh, F., and Boufadel, M.: On the use of numerical modelling for near-field pollutant dispersion in urban environments – A review, Environ. Pollut., 208, 271–283, https://doi.org/10.1016/j.envpol.2015.07.039, 2016.

Leighton, P. A.: Photochemistry of air pollution, Academic Press, New York, ISBN 9780323156455, 1961.

Liang, M., Chao, Y., Tu, Y., and Xu, T.: Vehicle Pollutant Dispersion in the Urban Atmospheric Environment: A Review of Mechanism, Modeling, and Application, Atmosphere, 14, 279, https://doi.org/10.3390/atmos14020279, 2023.

Lin, C., Wang, Y., Ooka, R., Flageul, C., Kim, Y., Kikumoto, H., Wang, Z., and Sartelet, K.: Modeling of street-scale pollutant dispersion by coupled simulation of chemical reaction, aerosol dynamics, and CFD, Atmos. Chem. Phys., 23, 1421–1436, https://doi.org/10.5194/acp-23-1421-2023, 2023.

Lugon, L., Sartelet, K., Kim, Y., Vigneron, J., and Chrétien, O.: Nonstationary modeling of NO₂, NO and NOₓ in Paris using the Street-in-Grid model: coupling local and regional scales with a two-way dynamic approach, Atmos. Chem. Phys., 20, 7717–7740, https://doi.org/10.5194/acp-20-7717-2020, 2020.

Lugon, L., Sartelet, K., Kim, Y., Vigneron, J., and Chrétien, O.: Simulation of primary and secondary particles in the streets of Paris using MUNICH, Faraday Discuss., 226, 432–456, https://doi.org/10.1039/D0FD00092B, 2021a.

Lugon, L., Vigneron, J., Debert, C., Chrétien, O., and Sartelet, K.: Black carbon modeling in urban areas: investigating the influence of resuspension and non-exhaust emissions in streets using the Street-in-Grid model for inert particles (SinG-inert), Geosci. Model Dev., 14, 7001–7019, https://doi.org/10.5194/gmd-14-7001-2021, 2021b.

Maison, A., Flageul, C., Carissimo, B., Tuzet, A., and Sartelet, K.: Parametrization of Horizontal and Vertical Transfers for the Street-Network Model MUNICH Using the CFD Model Code_Saturne, Atmosphere, 13, 527, https://doi.org/10.3390/atmos13040527, 2022.

Mallet, V., Quélo, D., Sportisse, B., Ahmed de Biasi, M., Debry, É., Korsakissok, I., Wu, L., Roustan, Y., Sartelet, K., Tombette, M., and Foudhil, H.: Technical Note: The air quality modeling system Polyphemus, Atmos. Chem. Phys., 7, 5479–5487, https://doi.org/10.5194/acp-7-5479-2007, 2007.

McHugh, C., Carruthers, D., and Edmunds, H.: ADMS and ADMS-Urban, Int. J. Environ. Pollut., 8, 438–440, https://doi.org/10.1504/IJEP.1997.028193, 1997.

Menut, L., Bessagnet, B., Briant, R., Cholakian, A., Couvidat, F., Mailler, S., Pennel, R., Siour, G., Tuccella, P., Turquety, S., and Valari, M.: The CHIMERE v2020r1 online chemistry-transport model, Geosci. Model Dev., 14, 6781–6811, https://doi.org/10.5194/gmd-14-6781-2021, 2021.

Olesen, H., Ketzel, M., Jensen, S., Løfstrøm, P., Im, U., and Becker, T.: User Guide to OML-Highway. A tool for air pollution assessments along highways., Tech. rep., Aarhus University, DCE – Danish Centre for Environment and Energy, https://dce2.au.dk/pub/TR59.pdf (last access: 14 February 2023), 2015.

Otosen, T.-B., Kakosimos, K. E., Johansson, C., Hertel, O., Brandt, J., Skov, H., Berkowicz, R., Ellermann, T., Jensen, S. S., and Ketzel, M.: Analysis of the impact of inhomogeneous emissions in the Operational Street Pollution Model (OSPM), Geosci. Model Dev., 8, 3231–3245, https://doi.org/10.5194/gmd-8-3231-2015, 2015.

Rood, A. S.: Performance evaluation of AERMOD, CALPUFF, and legacy air dispersion models using the Winter Validation Tracer Study dataset, Atmos. Environ., 89, 707–720, https://doi.org/10.1016/j.atmosenv.2014.02.054, 2014.

Sanchez, B., Santiago, J.-L., Martilli, A., Palacios, M., and Kirchner, F.: CFD modeling of reactive pollutant dispersion in simplified urban configurations with different chemical mechanisms, Atmos. Chem. Phys., 16, 12143–12157, https://doi.org/10.5194/acp-16-12143-2016, 2016.

Sarica, T., Sartelet, K., Roustan, Y., Kim, Y., Lugon, L., André, M., Marques, B., D'Anna, B., Chaillou, C., and Larrieu, C.: Modelling Pollutant Concentrations in Streets: A Sensitivity Analysis to Asphalt and Traffic Related Emissions, in: Air Pollution Modeling and its Application XXVIII, edited by: Mensink, C. and Jorba, O., Springer International Publishing, Cham, 287–293, https://doi.org/10.1007/978-3-031-12786-1_39, 2022.

Sarica, T., Maison, A., Roustan, Y., Ketzel, M., Jensen, S. S., Kim, Y., and Sartelet, K.: The Model of Urban Network of Intersecting Canyons and Highways (MUNICH), Zenodo [code and data set], https://doi.org/10.5281/zenodo.7778271, 2023a.

Sarica, T., Sartelet, K., Roustan, Y., Kim, Y., Lugon, L., Marques, B., D'Anna, B., Chaillou, C., and Larrieu, C.: Sensitivity of pollutant concentrations in urban streets to asphalt and traffic-related emissions, Environ. Pollut., 332, 121955, https://doi.org/10.1016/j.envpol.2023.121955, 2023b.

Sartelet, K., Zhu, S., Moukhtar, S., André, M., André, J., Gros, V., Favez, O., Brasseur, A., and Redaelli, M.: Emission of intermediate, semi and low volatile organic compounds from traffic and their impact on secondary organic aerosol concen-
Kim, Y., Wu, Y., Seigneur, C., and Roustan, Y.: Multi-scale modeling of urban air pollution: development and application of a Street-in-Grid model (v1.0) by coupling MUNICH (v1.0) and Polair3D (v1.8.1), Geosci. Model Dev., 11, 611–629, https://doi.org/10.5194/gmd-11-611-2018, 2018.

Kim, Y., Lugon, L., Maison, A., Sarica, T., Roustan, Y., Valari, M., Zhang, Y., André, M., and Sartelet, K.: MUNICH v2.0: a street-network model coupled with SSH-aerosol (v1.2) for multi-pollutant modelling, Geosci. Model Dev., 15, 7371–7396, https://doi.org/10.5194/gmd-15-7371-2022, 2022.

Kwak, K.-H. and Baik, J.-J.: A CFD modeling study of the impacts of NOₓ and VOC emissions on reactive pollutant dispersion in and above a street canyon, Atmos. Environ., 46, 71–80, https://doi.org/10.1016/j.atmosenv.2011.10.024, 2012.

Landrigan, P. J., Fuller, R., Acosta, N. J. R., Adeyi, O., Arnold, R., Basu, N. N., Baldé, A. B., Bertollini, R., Bose-O'Reilly, S., Boufford, J. I., Breysse, P. N., Chiles, T., Mahidol, C., Coll-Seck, A. M., Cropper, M. L., Fobil, J., Fuster, V., Greenstone, M., Haines, A., Hanrahan, D., Hunter, D., Khare, M., Krupnick, A., Lanphear, B., Lohani, B., Martin, K., Mathiasen, K. V., McTeer, M. A., Murray, C. J. L., Ndahimananjara, J. D., Perera, F., Potočnik, J., Preker, A. S., Ramesh, J., Rockström, J., Salinas, C., Samson, L. D., Sandilya, K., Sly, P. D., Smith, K. R., Steiner, A., Stewart, R. B., Suk, W. A., van Schayck, O. C. P., Yadama, G. N., Yumkella, K., and Zhong, M.: The Lancet Commission on pollution and health, Lancet, 391, 462–512, https://doi.org/10.1016/S0140-6736(17)32345-0, 2018.

Lateb, M., Meroney, R., Yataghene, M., Fellouah, H., Saleh, F., and Boufadel, M.: On the use of numerical modelling for near-field pollutant dispersion in urban environments – A review, Environ. Pollut., 208, 271–283, https://doi.org/10.1016/j.envpol.2015.07.039, 2016.

Leighton, P. A.: Photochemistry of air pollution, Academic Press, New York, ISBN 9780323156455, 1961.

Liang, M., Chao, Y., Tu, Y., and Xu, T.: Vehicle Pollutant Dispersion in the Urban Atmospheric Environment: A Review of Mechanism, Modeling, and Application, Atmosphere, 14, 279, https://doi.org/10.3390/atmos14020279, 2023.

Lin, C., Wang, Y., Ooka, R., Flageul, C., Kim, Y., Kikumoto, H., Wang, Z., and Sartelet, K.: Modeling of street-scale pollutant dispersion by coupled simulation of chemical reaction, aerosol dynamics, and CFD, Atmos. Chem. Phys., 23, 1421–1436, https://doi.org/10.5194/acp-23-1421-2023, 2023.

Lugon, L., Sartelet, K., Kim, Y., Vigneron, J., and Chrétien, O.: Nonstationary modeling of NO₂, NO and NOₓ in Paris using the Street-in-Grid model: coupling local and regional scales with a two-way dynamic approach, Atmos. Chem. Phys., 20, 7717–7740, https://doi.org/10.5194/acp-20-7717-2020, 2020.

Lugon, L., Sartelet, K., Kim, Y., Vigneron, J., and Chrétien, O.: Simulation of primary and secondary particles in the streets of Paris using MUNICH, Faraday Discuss., 226, 432–456, https://doi.org/10.1039/D0FD00092B, 2021a.

Lugon, L., Vigneron, J., Debert, C., Chrétien, O., and Sartelet, K.: Black carbon modeling in urban areas: investigating the influence of resuspension and non-exhaust emissions in streets using the Street-in-Grid model for inert particles (SinG-inert), Geosci. Model Dev., 14, 7001–7019, https://doi.org/10.5194/gmd-14-7001-2021, 2021b.

Maison, A., Flageul, C., Carissimo, B., Tuzet, A., and Sartelet, K.: Parametrization of Horizontal and Vertical Transfers for the Street-Network Model MUNICH Using the CFD Model Code_Saturne, Atmosphere, 13, 527, https://doi.org/10.3390/atmos13040527, 2022.

Mallet, V., Quélo, D., Sportisse, B., Ahmed de Biasi, M., Debry, É., Korsakissok, I., Wu, L., Roustan, Y., Sartelet, K., Tombette, M., and Foudhil, H.: Technical Note: The air quality modeling system Polyphemus, Atmos. Chem. Phys., 7, 5479–5487, https://doi.org/10.5194/acp-7-5479-2007, 2007.

McHugh, C., Carruthers, D., and Edmunds, H.: ADMS and ADMS-Urban, Int. J. Environ. Pollut., 8, 438–440, https://doi.org/10.1504/IJEP.1997.028193, 1997.

Menut, L., Bessagnet, B., Briant, R., Cholakian, A., Couvidat, F., Mailler, S., Pennel, R., Siour, G., Tuccella, P., Turquety, S., and Valari, M.: The CHIMERE v2020r1 online chemistry-transport model, Geosci. Model Dev., 14, 6781–6811, https://doi.org/10.5194/gmd-14-6781-2021, 2021.

Olesen, H., Ketzel, M., Jensen, S., Løfstrøm, P., Im, U., and Becker, T.: User Guide to OML-Highway. A tool for air pollution assessments along highways., Tech. rep., Aarhus University, DCE – Danish Centre for Environment and Energy, https://dce2.au.dk/pub/TR59.pdf (last access: 14 February 2023), 2015.

Otosen, T.-B., Kakosimos, K. E., Johansson, C., Hertel, O., Brandt, J., Skov, H., Berkowicz, R., Ellermann, T., Jensen, S. S., and Ketzel, M.: Analysis of the impact of inhomogeneous emissions in the Operational Street Pollution Model (OSPM), Geosci. Model Dev., 8, 3231–3245, https://doi.org/10.5194/gmd-8-3231-2015, 2015.

Rood, A. S.: Performance evaluation of AERMOD, CALPUFF, and legacy air dispersion models using the Winter Validation Tracer Study dataset, Atmos. Environ., 89, 707–720, https://doi.org/10.1016/j.atmosenv.2014.02.054, 2014.

Sanchez, B., Santiago, J.-L., Martilli, A., Palacios, M., and Kirchner, F.: CFD modeling of reactive pollutant dispersion in simplified urban configurations with different chemical mechanisms, Atmos. Chem. Phys., 16, 12143–12157, https://doi.org/10.5194/acp-16-12143-2016, 2016.

Sarica, T., Sartelet, K., Roustan, Y., Kim, Y., Lugon, L., André, M., Marques, B., D'Anna, B., Chaillou, C., and Larrieu, C.: Modelling Pollutant Concentrations in Streets: A Sensitivity Analysis to Asphalt and Traffic Related Emissions, in: Air Pollution Modeling and its Application XXVIII, edited by: Mensink, C. and Jorba, O., Springer International Publishing, Cham, 287–293, https://doi.org/10.1007/978-3-031-12786-1_39, 2022.

Sarica, T., Maison, A., Roustan, Y., Ketzel, M., Jensen, S. S., Kim, Y., and Sartelet, K.: The Model of Urban Network of Intersecting Canyons and Highways (MUNICH), Zenodo [code and data set], https://doi.org/10.5281/zenodo.7778271, 2023a.

Sarica, T., Sartelet, K., Roustan, Y., Kim, Y., Lugon, L., Marques, B., D'Anna, B., Chaillou, C., and Larrieu, C.: Sensitivity of pollutant concentrations in urban streets to asphalt and traffic-related emissions, Environ. Pollut., 332, 121955, https://doi.org/10.1016/j.envpol.2023.121955, 2023b.

Sartelet, K., Zhu, S., Moukhtar, S., André, M., André, J., Gros, V., Favez, O., Brasseur, A., and Redaelli, M.: Emission of intermediate, semi and low volatile organic compounds from traffic and their impact on secondary organic aerosol concen-
trations over Greater Paris, Atmos. Environ., 180, 126–137,
https://doi.org/10.1016/j.atmosenv.2018.02.031, 2018.

Sartelet, K., Couvidat, F., Wang, Z., Flageul, C., and Kim, Y.: SSH-Aerosol v1.1: A Modular Box Model to Simulate the Evolution of Primary and Secondary Aerosols, Atmosphere, 11, 525, https://doi.org/10.3390/atmos11050525, 2020.

Skamarock, W., Klemp, J., Dudhia, J., Gill, D., Barker, D., Wang, W., Huang, X.-Y., and Duda, M.: A Description of the Advanced Research WRF Version 3, Tech. rep., UCAR/NCAR, https://doi.org/10.5065/D68S4MVH, 2008.

Solazzo, E., Cai, X., and Vardoulakis, S.: Modelling wind flow and vehicle-induced turbulence in urban streets, Atmos. Environ., 42, 4918–4931, https://doi.org/10.1016/j.atmosenv.2008.02.032, 2008.

Soulhac, L., Garbero, V., Salizzoni, P., Mejean, P., and Perkins, R.: Flow and dispersion in street intersections, Atmos. Environ., 43, 2981–2996, https://doi.org/10.1016/j.atmosenv.2009.02.061, 2009.

Soulhac, L., Salizzoni, P., Cierco, F.-X., and Perkins, R.: The model SIRANE for atmospheric urban pollutant dispersion; part I, presentation of the model, Atmos. Environ., 45, 7379–7395, https://doi.org/10.1016/j.atmosenv.2011.07.008, 2011.

Soulhac, L., Fellini, S., Nguyen, C. V., and Salizzoni, P.: Evaluation of Photostationary and Non-Photostationary Operational Models for NOX Pollution in a Street Canyon, Atmos. Environ., 297, 119589, https://doi.org/10.1016/j.atmosenv.2023.119589, 2023.

United Nations: World Urbanization Prospects: The 2018 Revision (ST/ESA/SER.A/420), Tech. rep., United Nations, https://population.un.org/wup/Publications/Files/WUP2018-Report.pdf (last access: 12 September 2023), 2019.

Vardoulakis, S., Fisher, B., Pericleous, K., and Gonzalez-Flesca, N.: Modelling air quality in street canyons: a review, Atmos. Environ., 37, 155–182, https://doi.org/10.1016/S1352-2310(02)00857-9, 2003.

Vardoulakis, S., Solazzo, E., and Lumbreras, J.: Intra-urban and street scale variability of BTEX, NO₂ and O₃ in Birmingham, UK: Implications for exposure assessment, Atmos. Environ., 45, 5069–5078, https://doi.org/10.1016/j.atmosenv.2011.06.038, 2011.

Wang, W.: Analytically modelling mean wind and stress profiles in canopies, Bound.-Lay. Meteorol., 151, 239–256, https://doi.org/10.1007/s10546-013-9899-6, 2014.

Wolf, T., Pettersson, L. H., and Esau, I.: A very high-resolution assessment and modelling of urban air quality, Atmos. Chem. Phys., 20, 625–647, https://doi.org/10.5194/acp-20-625-2020, 2020.

Wong, D. C., Pleim, J., Mathur, R., Binkowski, F., Otte, T., Gilliam, R., Pouliot, G., Xiu, A., Young, J. O., and Kang, D.: WRF-CMAQ two-way coupled system with aerosol feedback: software development and preliminary results, Geosci. Model Dev., 5, 299–312, https://doi.org/10.5194/gmd-5-299-2012, 2012.

Xie, S., Zhang, Y., Qi, L., and Tang, X.: Spatial distribution of traffic-related pollutant concentrations in street canyons, Atmos. Environ., 37, 3213–3224, https://doi.org/10.1016/S1352-2310(03)00321-2, 2003.

Zhong, J., Cai, X.-M., and Bloss, W. J.: Large eddy simulation of reactive pollutants in a deep urban street canyon: Coupling dynamics with O₃-NOₓ-VOC chemistry, Environ. Pollut., 224, 171–184, https://doi.org/10.1016/j.envpol.2017.01.076, 2017.
trations over Greater Paris, Atmos. Environ., 180, 126–137,
https://doi.org/10.1016/j.atmosenv.2018.02.031, 2018.

Sartelet, K., Couvidat, F., Wang, Z., Flageul, C., and Kim, Y.: SSH-Aerosol v1.1: A Modular Box Model to Simulate the Evolution of Primary and Secondary Aerosols, Atmosphere, 11, 525, https://doi.org/10.3390/atmos11050525, 2020.

Skamarock, W., Klemp, J., Dudhia, J., Gill, D., Barker, D., Wang, W., Huang, X.-Y., and Duda, M.: A Description of the Advanced Research WRF Version 3, Tech. rep., UCAR/NCAR, https://doi.org/10.5065/D68S4MVH, 2008.

Solazzo, E., Cai, X., and Vardoulakis, S.: Modelling wind flow and vehicle-induced turbulence in urban streets, Atmos. Environ., 42, 4918–4931, https://doi.org/10.1016/j.atmosenv.2008.02.032, 2008.

Soulhac, L., Garbero, V., Salizzoni, P., Mejean, P., and Perkins, R.: Flow and dispersion in street intersections, Atmos. Environ., 43, 2981–2996, https://doi.org/10.1016/j.atmosenv.2009.02.061, 2009.

Soulhac, L., Salizzoni, P., Cierco, F.-X., and Perkins, R.: The model SIRANE for atmospheric urban pollutant dispersion; part I, presentation of the model, Atmos. Environ., 45, 7379–7395, https://doi.org/10.1016/j.atmosenv.2011.07.008, 2011.

Soulhac, L., Fellini, S., Nguyen, C. V., and Salizzoni, P.: Evaluation of Photostationary and Non-Photostationary Operational Models for NOX Pollution in a Street Canyon, Atmos. Environ., 297, 119589, https://doi.org/10.1016/j.atmosenv.2023.119589, 2023.

United Nations: World Urbanization Prospects: The 2018 Revision (ST/ESA/SER.A/420), Tech. rep., United Nations, https://population.un.org/wup/Publications/Files/WUP2018-Report.pdf (last access: 12 September 2023), 2019.

Vardoulakis, S., Fisher, B., Pericleous, K., and Gonzalez-Flesca, N.: Modelling air quality in street canyons: a review, Atmos. Environ., 37, 155–182, https://doi.org/10.1016/S1352-2310(02)00857-9, 2003.

Vardoulakis, S., Solazzo, E., and Lumbreras, J.: Intra-urban and street scale variability of BTEX, NO₂ and O₃ in Birmingham, UK: Implications for exposure assessment, Atmos. Environ., 45, 5069–5078, https://doi.org/10.1016/j.atmosenv.2011.06.038, 2011.

Wang, W.: Analytically modelling mean wind and stress profiles in canopies, Bound.-Lay. Meteorol., 151, 239–256, https://doi.org/10.1007/s10546-013-9899-6, 2014.

Wolf, T., Pettersson, L. H., and Esau, I.: A very high-resolution assessment and modelling of urban air quality, Atmos. Chem. Phys., 20, 625–647, https://doi.org/10.5194/acp-20-625-2020, 2020.

Wong, D. C., Pleim, J., Mathur, R., Binkowski, F., Otte, T., Gilliam, R., Pouliot, G., Xiu, A., Young, J. O., and Kang, D.: WRF-CMAQ two-way coupled system with aerosol feedback: software development and preliminary results, Geosci. Model Dev., 5, 299–312, https://doi.org/10.5194/gmd-5-299-2012, 2012.

Xie, S., Zhang, Y., Qi, L., and Tang, X.: Spatial distribution of traffic-related pollutant concentrations in street canyons, Atmos. Environ., 37, 3213–3224, https://doi.org/10.1016/S1352-2310(03)00321-2, 2003.

Zhong, J., Cai, X.-M., and Bloss, W. J.: Large eddy simulation of reactive pollutants in a deep urban street canyon: Coupling dynamics with O₃-NOₓ-VOC chemistry, Environ. Pollut., 224, 171–184, https://doi.org/10.1016/j.envpol.2017.01.076, 2017.