# Parameterizing the aerodynamic effect of trees in street canyons for the street network model MUNICH using the CFD model Code_Saturne

Alice Maison, Cédric Flageul, Bertrand Carissimo, Yunyi Wang, Andrée

Tuzet, Karine Sartelet

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

## To cite this version:

Alice Maison, Cédric Flageul, Bertrand Carissimo, Yunyi Wang, Andrée Tuzet, et al., Parameterizing the aerodynamic effect of trees in street canyons for the street network model MUNICH using the CFD model Code_Saturne. Atmospheric Chemistry and Physics, 2022, 22 (14), pp.9369-9388. 10.5194/acp-22-9369-2022. hal-03819560

HAL Id: hal-03819560

https://hal.science/hal-03819560

Submitted on 14 Sep 2023

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
<figure><img src="image_3.png" /></figure>

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

# Parameterizing the aerodynamic effect of trees in street canyons for the street network model MUNICH using the CFD model Code_Saturne

Alice Maison¹,², Cédric Flageul³, Bertrand Carissimo¹, Yunyi Wang¹, Andrée Tuzet², and Karine Sartelet¹

¹CEREA, École des Ponts, EDF R&D, Marne-la-Vallée, France

²Université Paris-Saclay, INRAE, AgroParisTech, UMR EcoSys, 78850 Thiverval-Grignon, France

³PPRIME Institute, Curiosity Group, Université de Poitiers, CNRS, ISAE-ENSMA, Poitiers, France

**Correspondence:** Alice Maison (alice.maison@enpc.fr) and Karine Sartelet (karine.sartelet@enpc.fr)

Received: 15 April 2022 – Discussion started: 29 April 2022

Revised: 24 June 2022 – Accepted: 1 July 2022 – Published: 20 July 2022

**Abstract.** Trees provide many ecosystem services in cities such as urban heat island reduction, water runoff limitation, and carbon storage. However, the presence of trees in street canyons reduces the wind velocity in the street and limits pollutant dispersion. Thus, to obtain accurate simulations of pollutant concentrations, the aerodynamic effect of trees should be taken into account in air quality models at the street level.

The Model of Urban Network of Intersecting Canyons and Highways (MUNICH) simulates the pollutant concentrations in a street network, considering dispersion and physico-chemical processes. It can be coupled to a regional-scale chemical transport model to simulate air quality over districts or cities. The aerodynamic effect of the tree crown is parameterized here through its impact on the average wind velocity in the street direction and the vertical transfer coefficient associated with the dispersion of a tracer. The parameterization is built using local-scale simulations performed with the computational fluid dynamics (CFDs) code Code_Saturne. The two-dimensional CFD simulations in an infinite street canyon are used to quantify the effect of trees, depending on the tree characteristics (leaf area index, crown volume fraction, and tree height to street height ratio) using a drag porosity approach. The tree crown slows down the flow and produces turbulent kinetic energy in the street, thus impacting the tracer dispersion. This effect increases with the leaf area index and the crown volume fraction of the trees, and the average horizontal velocity in the street is reduced by up to 68 %, while the vertical transfer coefficient by up to 23 % in the simulations performed here.

A parameterization of these effects on horizontal and vertical transfers for the street model MUNICH is proposed. Existing parameterizations in MUNICH are modified based on Code_Saturne simulations to account for both building and tree effects on vertical and horizontal transfers. The parameterization is built to obtain similar tree effects (quantified by a relative deviation between the cases without and with trees) between Code_Saturne and MUNICH. The vertical wind profile and mixing length depend on leaf area index, crown radius, and tree height to street height ratio. The interaction between the trees and the street aspect ratio is also considered.
# 1 Introduction

Cities are, by definition, areas with high densities of people, infrastructure, and activities, and this urbanization causes many issues. First, air quality is poor because of the numerous air pollutants emitted by anthropic activities such as traffic, industries, or residential activities, and the reduction in air flow by high buildings limits the dispersion of these pollutants (Faiz, 1993; Akimoto, 2003; Yuan et al., 2014; Zhang et al., 2020). In addition to air pollution issues, radiative and water budgets are strongly modified in cities compared to the countryside (Bozonnet et al., 2015). Temperatures are, on average, higher than in the countryside because of the urban heat island created by additional anthropogenic energy released, storage of radiative energy in dark materials, radiation multi-reflection, and lack of vegetation and associated evaporation (Oke, 1982; Pigeon et al., 2007; Stewart, 2011; Hebbert and Jankovic, 2013). Impervious soils also decrease water infiltration and intensify runoff (Leopold, 1968). In addition, growing urbanization and increasing extreme events (due to climate change) such as pollution peaks, heat waves, and floods have negative consequences on the environment and human health (Robine et al., 2007; Angel et al., 2011; Pascal et al., 2013; West et al., 2016; IPCC, 2021).

One nature-based solution is to green the city by planting vegetation as lawns, planting trees in streets or in parks, and growing green walls and roofs (Livesley et al., 2016; Revelli and Porporato, 2018). Vegetation and especially trees contribute to improve human thermal comfort by creating a favorable micro-climate with lower air temperature (through solar radiation interception and creation of shade) and higher evaporation (through ground and vegetation evapotranspiration; Taha et al., 1991; Bowler et al., 2010; Gillner et al., 2015; Klemm et al., 2015; Lobaccaro and Acero, 2015;GNANWARDENA et al., 2017). This positive effect of trees is significant in particular during heat wave episodes, which will be more frequent in the future due to climate change (IPCC, 2021). Trees and vegetated areas also favor infiltration in soils that contributes to offset water runoff induced by soil artificialization (Armson et al., 2013; Berland et al., 2017). Besides, vegetation is known to store carbon (Nowak and Crane, 2002; Svirejeva-Hopkins et al., 2004) and to enhance human well-being (van Dillen et al., 2012; Bertram and Reddanz, 2015; Krekel et al., 2015). For all these ecosystem services in urban areas, city greening is often promoted and, for example, the city of Paris has about 205 000 trees of which 52 % are roadside trees (Direction des Espaces Verts et de l'Environnement - Mairie de Paris, 2021).

Many studies have tried to figure out the impact of trees on air pollution in a street canyon, and they have shown that trees are an important parameter to take into account if we want to understand and accurately simulate the pollutant concentrations in the streets (Beckett et al., 1998; Nowak et al., 1998; Jayasooriya et al., 2017). Vegetation, and especially trees, represent surfaces available for pollutant dry de-

position and, hence, can contribute to reducing air pollutant concentrations. However, pollutant removal and its impact on air quality vary greatly, depending on tree characteristics, tree species, and pollutant type (Nowak et al., 2006; Hwang et al., 2011; Selmi et al., 2016; Xue and Li, 2017; Ozdemir, 2019). Trees may also affect atmospheric chemistry by emitting biogenic volatile organic compounds (BVOCs), which may lead to the formation of ozone and secondary organic aerosols (Calfapietra et al., 2013; Préndez et al., 2019; Gu et al., 2021). Furthermore, trees may alter air quality by influencing aerodynamic processes and limiting the pollutant dispersion (Buccolieri et al., 2011; Wania et al., 2012; Vos et al., 2013; Gromke and Ruck, 2007, 2009, 2012; Gromke and Blocken, 2015). The aerodynamic effect of trees is defined as the drag force resulting from the friction between the air and the leaves. Since the tree crown can be seen as a porous medium as air passes through it but is slowed down, and the drag force increases with the leaf surface.

For wind perpendicular to a street, the air recirculates inside the street canyon (Harman et al., 2004). Pollutants emitted at the bottom of the street (by traffic) accumulate on the leeward side of the street, inducing higher local concentrations (Vardoulakis et al., 2003; Cai et al., 2008; Huang et al., 2019). Obstacles in the street, such as trees, can strongly impact the air flow and intensify the pollutant accumulation (Vos et al., 2013). This effect has been studied using computational fluid dynamics (CFD) models where the aerodynamic influence of trees on the flow is represented by a porosity model (Buccolieri et al., 2009; Zaïdi et al., 2013; Wei et al., 2016; Jeanjean et al., 2017; Santiago et al., 2017). Vos et al. (2013) showed that pollutant concentrations may increase by 20 % in a street because of the presence of two rows of trees in the street. The effects of vegetation depend on the height, width, and density (leaf area index) of trees (Vos et al., 2013; Janhäll, 2015), as well as the height to width ratio of streets (Wania et al., 2012). Finally, it is necessary to accurately assess the effect of trees on pollutant dispersion in order to find what configurations are more effective in reducing air pollution in street canyon and to guide urban development policy (Janhäll, 2015).

The effect of trees on aerodynamic processes should also be considered at the street level in air quality models. As discussed previously, CFD models including trees are used to study wind fields and pollutant transport in street canyons (Li et al., 2006). However, as the street is discretized with a fine mesh, the computational cost is high, and simulations at the city scale are too expensive today. Fast-running codes, such as simplified street network or street-in-grid models, are developed to simulate street pollutant concentrations over neighborhoods or cities, but they do not take into account the effect of trees in the streets. The objective of this study is to parameterize the effect of trees on air flow in the Model of Urban Network of Intersecting Canyons and Highways (MUNICH; Kim et al., 2018; Lugon et al., 2020, http://cerea.enpc.fr/munich/, last access: 17 December 2021). To build this pa-
rameterization, simulations in street canyons are performed with Code_Saturne (Archambeau et al., 2004, https://www. code-saturne.org/, last access: 17 December 2021), a CFD code, which can represent the tree aerodynamic effect with a drag porosity approach (Katul et al., 2004) and has previ- ously been compared with field measurements (Zaïdi et al., 2013). MUNICH parameterizations have already been com- pared with Code_Saturne results in a treeless canyon, and a new parameterization for horizontal and vertical transfers has been developed in MUNICH based on Wang (2012, 2014) and Code_Saturne simulations (Maison et al., 2022). In the present study, the tree aerodynamic effect is added to this parameterization, and Code_Saturne (version 6.0) is used as a reference to parameterize the aerodynamic effect of trees in the street network model MUNICH. CFD simulations are performed in three streets of aspect ratios varying from 0.3 to 1.0, and a large range of tree leaf area index, crown radius, and heights is tested. The tree aerodynamic effect quantified with Code_Saturne is analyzed depending on these street and tree characteristics.

The approach to modeling the dispersion of pollutants in MUNICH and Code_Saturne is fundamentally different due to their physical modeling and discretization in space and time. MUNICH is a street network model that simulates air pollutant concentration in an urban canopy. Street dimensions and pollutant concentrations are assumed to be homogeneous in each street segment. Air flow is divided into a horizontal flux from one street to another and a vertical flux between the street and the background (Maison et al., 2022). Background concentrations above the street can be computed by 3D chemistry transport models (CTMs), such as Polair3D (Sartelet et al., 2018; Lugon et al., 2020, 2021). To build the tree parameterization, the CFD simulation setup is adapted for the comparison with MUNICH, and several simulations are performed with Code_Saturne, considering a range of street and tree characteristics.

The structure of the paper is as follows. The MUNICH and Code_Saturne models are presented in Sect. 2. Then, the tree effect on horizontal wind speed and vertical transfer coefficient is quantified with Code_Saturne simulations in Sect. 3 and parameterized in MUNICH in Sect. 4. Conclusions are presented in Sect. 5.

## 2 Materials and methods

### 2.1 Description of MUNICH

In MUNICH, each street segment is assumed to be homogeneous, i.e., with uniform building height $H$ and street width $W$ and of length $L$ (m). The street is characterized by its height-to-width ratio called the aspect ratio, $a_r = H/W$ (-). Only the average pollutant concentrations over the street are considered. Pollutants are transported by the horizontal wind speed (advection) in the street network and by a vertical transfer coefficient between the streets and the background.

Several parameterizations of the horizontal wind speed and of the vertical transfer coefficient exist in MUNICH. The ones recently developed in Maison et al. (2022) and based on Code_Saturne simulations are used and detailed here. The vertical profile of the wind speed in the street direction is calculated as an attenuation of the wind speed in the street direction and at the roof level, $U_{H,φ}$ (m s⁻¹), as follows (Maison et al., 2022):

$$ U(z) = U_{H,φ} [C_1 I_0(g(z)) + C_2 K_0(g(z))] \\ \text{with } U_{H,φ} = U_H |\cos(\varphi)|, \qquad (1) $$

where $C_1$ and $C_2$ are integration coefficients, and $I_0$ and $K_0$ are the first and second type modified Bessel functions of order 0. Besides, the wind speed at the roof level $U_H$ has to be multiplied by $|\cos(\varphi)|$ to select the component of the wind speed in the street direction, where $\varphi$ is the angle between the wind direction and the street orientation (°). This vertical wind profile is then integrated between the soil roughness $z_{0s}$ (m) and $H$ to compute the average horizontal wind speed in the street direction. The function $g$ is calculated as follows (Wang, 2012, 2014):

$$ g(z) = 2\sqrt{\alpha \frac{z}{H}}, \qquad (2) $$

and

$$ C_1 = \frac{1}{I_0(g(H)) - I_0(g(z_{0s})) K_0(g(H))/K_0(g(z_{0s}))} \\ \text{and } C_2 = -\frac{C_1 I_0(g(z_{0s}))}{K_0(g(z_{0s}))}, \qquad (3) $$

where $\alpha$ is a dimensionless coefficient expressing the effects of wind angle on wind attenuation in the street. It is computed as follows:

$$ \alpha = \frac{C_B a_r}{\kappa s_H} \quad \text{with } C_B = 0.31 [1 - \exp(-1.6 a_r)] f_{\varphi} \qquad (4) $$

$$ \text{and } f_{\varphi} = \begin{cases} |\cos(2\varphi)|^3 & \text{if } \varphi \in [0, 45^{\circ}] \cup [135, 225^{\circ}] \\ & \cup [315, 360^{\circ}] \\ 0 & \text{if } \varphi \in [45, 135^{\circ}] \cup [225, 315^{\circ}], \end{cases} \qquad (5) $$

where $s_H = s(z = H)$ is a dimensionless factor describing the effect of canopy on the mixing length $l_m$ (Wang, 2012, 2014). The mixing length is calculated as follows:

$$ \frac{1}{l_m} = \frac{1}{\kappa z} + \frac{1}{l_{c_b}} \Rightarrow l_m = \kappa z \frac{l_{c_b}}{l_{c_b} + \kappa z} = \kappa z s(z), \qquad (6) $$

where $\kappa z$ corresponds to the mixing length over a rough bare soil (without canopy), and $l_{c_b}$ is the characteristic length (m) in the street canyon, corresponding to the mixing length of the urban canopy alone ($l_{c_b} = 0.5$ W).

The vertical transfer coefficient that drives pollutant exchange between the street and the background zone is calculated at $z = H$ as follows:

$$ q_{\text{vert}} = \sigma_W l_m = \sigma_W \kappa H s_H, \qquad (7) $$
where $\sigma_W$ (in m s⁻¹) is a velocity scale equal to the standard deviation of the vertical wind speed. It depends on both the friction velocity above the urban canopy, $u_*$, and on atmospheric stability (Soulhac et al., 2011).

In Sect. 4, the tree effect is parameterized by taking into account the characteristic length of the trees in the mixing length $l_m$ (Eq. 19, which leads to modifications in the coefficient $\alpha$ Eq. 23).

## 2.2 Description of Code_Saturne

### 2.2.1 Street and tree modeling setup

In the present study, the tree effect is studied in three street canyons of the following different street aspect ratios: a wide street canyon (WC), an intermediate canyon (IC), and a narrow canyon (NC). Their characteristics are presented in Table 1.

A color code is introduced to simplify the reading of the figures, with red symbols for the wide canyon (WC), purple ones for the intermediate canyon (IC), and blue ones for the narrow canyon (NC). Abbreviations, parameters, and variables used are listed in Appendix A.

The $k-\epsilon$ linear production turbulence model is used in Code_Saturne. Stationary simulations are performed with a thermally neutral atmosphere. A 2D infinite street canyon is modeled with periodic condition on the y axis, but the flow and the wind speed vector are 3D. The mesh is composed of hexahedral cells of 1 m in the y-axis direction and 0.5 m on the x and z axes. The vertical profiles of U, k, and $\epsilon$ are set in the inlet (top left border of the domain). The complete description of Code_Saturne simulation setup in treeless canyons can be found in Maison et al. (2022).

The tree geometry and its representation in Code_Saturne are shown in Fig. 1. Here, two rows of trees ($n=2$) are considered, with one on each side of the street, and whose positions on the x axis are $x=34.5$ m and $x=48.0$ m (the wall positions are $x=27.5$ and $x=55.0$ m). The tree crown centers are located in a position on the x axis so that the largest crowns (CVF ≈ 25%) do not reach the street walls for the three street canyons studied. r is the tree radius (m), $h$, $h_{min}$, and $h_{max}$ correspond, respectively, to the middle, minimum, and maximum heights of the tree crown (m). Note that all the simulations verify $h_{max} \le H$, i.e., the top of the trees do not exceed the top of the street.

In the cells containing the trees, an additional drag term is added to the Navier-Stokes equations as follows (Katul et al., 2004; Zaïdi et al., 2013):

$$S_{u,i} = -\rho \text{LAD} C_{D_t} |U| U_i, \quad (8)$$

where $|U|$ is the modulus of U, $U_i$ is the velocity in the i-direction, LAD is the leaf area density in square meters of leaf surface per cubic meter of crown volume (m² leaves m⁻³ crown), $\rho$ is the air density, and $C_{D_t}$ is the tree drag coefficient set to 0.2, which is a representative value for

trees (Katul and Albertson, 1998). The leaf area index (LAI) is the one-sided green leaf area per unit ground surface area, and the leaf area density (LAD) is the one-sided green leaf area per unit volume. They are calculated as follows:

$$ \text{LAI} = \frac{\text{surface of leaves}}{\text{soil projected surface area}} \left( m_{\text{leaf}}^2 m_{\text{soil}}^{-2} \right) \\ \text{and LAD} = \frac{\text{surface of leaves}}{\text{volume of tree crown}} \left( m_{\text{leaf}}^2 m_{\text{crown}}^{-3} \right). \quad (9) $$

The source terms for trees are also implemented in the $k-\epsilon$ equations as follows:

$$S_k = \rho \text{LAD} C_{D_t} (\beta_p |U|^3 - \beta_d |U|k) \quad (10)$$

$$S_{\epsilon} = \rho \text{LAD} C_{D_t} \left( C_{4\epsilon} \beta_p |U|^3 \frac{\epsilon}{k} - C_{5\epsilon} \beta_d |U| \epsilon \right), \quad (11)$$

where $C_{4\epsilon} = C_{5\epsilon} = 0.9$, $\beta_p = 1.0$, and $\beta_d = 5.03$ are constants of the model (Zaïdi et al., 2013). Note that this type of tree aerodynamic effect modeling is commonly used and evaluated by comparison with experimental results (Buccolieri et al., 2018).

Only the impact of the tree leaves is considered, and the impact of the tree trunk and branches on the flow is not modeled. As the streets are modeled in 2D and an infinite length in the y direction is assumed (see Sect. 2), the LAI considered here is the equivalent cylindrical LAI₂D. Usually, studies consider the tree LAI, which is well defined for a 3D tree. The relation between the tree LAI₃D and the equivalent cylindrical LAI₂D is given in Appendix B. To check that this effect of trees is well modeled using this simplified 2D setup, 3D simulations were performed and compared to some of the 2D simulations presented here, showing very good agreement between the 2D and 3D simulations in terms of the influence of the trees on the flow (Appendix B).

To gain a complete vision of the effect of trees, 45 simulations are performed per canyon, and the impact of three independent tree parameters are studied, including the following:

- The 2D-equivalent LAI, i.e., LAI₂D = 0.5, 1, 2, 3, and 4 (see Appendix B, for the conversion between LAI₃D and equivalent LAI₂D).

- The crown volume fraction (CVF), calculated as the tree (2D cylindrical) crown volume divided by the street volume, i.e., CVF = $\frac{n\pi r^2}{HW}$. This ratio is already used in studies quantifying the effect of trees on street pollution such as Gromke and Blocken (2015). The relation between 2D and 3D CVF is shown in Appendix B. Note that, apart from Appendix B, all CVF mentioned in the figures and in the text refer to the 2D CVF. In total, the following three CVF ranges are simulated: CVF ≈ 5, 10, and 25%.

- The tree-to-street height ratio calculated as the crown middle height divided by the building height ($h/H$).
<table><caption>**Table 1.** Characteristics of the three canyons studied.</caption><thead><tr><th>Canyon</th><th>Building height<br>H (m)</th><th>Street width<br>W (m)</th><th>Street aspect ratio<br>a<sub>r</sub> (-)</th><th>Maximum height of the<br>domain (3H) (m)</th></tr></thead><tbody><tr><td>WC</td><td>8.5</td><td>27.5</td><td>0.3</td><td>25.5</td></tr><tr><td>IC</td><td>14.0</td><td>27.5</td><td>0.5</td><td>42.0</td></tr><tr><td>NC</td><td>27.5</td><td>27.5</td><td>1.0</td><td>82.5</td></tr></tbody></table>

<figure><img src="image_5.png" /><figcaption>**Figure 1.** The 3D scheme of the street canyon with cylindrical tree crown dimensions.</figcaption></figure>

In total, the following three $h/H$ ranges are simulated: $h/H \approx 1/3$, $1/2$, and $2/3$.

Note that these three tree characteristics are normalized by street characteristics, and similar values are chosen in the three canyons to be able to compare them and to quantify if the street aspect ratio influences the tree effect (i.e., if there is an interaction between tree and building effects). However, the normalized tree characteristics are not exactly equal in the three canyons due to the street sizes and the 0.5 m mesh cells that limit the possible range of tree sizes. The detailed list of the tree characteristics and the different values used for the tree parameters are presented in Table A4.

### 2.2.2 Calculation of vertical and horizontal transfers in Code_Saturne for comparison to MUNICH

To evaluate the vertical transfer between the street and the background zone in Code_Saturne, a passive tracer is emitted in each mesh cell of the street, with an arbitrary stationary emission rate $e = 1000 \mu g s^{-1}$ for a street canyon of length $L = 1$ m. For the tracer to be dispersed only by vertical transfers and not by horizontal winds within the street, simulations are performed with the wind perpendicular to the street at the inlet ($\varphi = 90^{\circ}$). The initial street and background tracer

concentrations are zero. At the end of the simulation (stationary state reached), the tracer concentration is averaged in the street and in the background zone (denoted as $C_{street}$ and $C_{bg}$ in $\mu g m^{-3}$) to reproduce a MUNICH homogeneous street assumption. The background zone corresponds to the area of the same volume as the street but located just above the street (see Maison et al., 2022, for more details). In the absence of other processes (horizontal transport, deposition, and chemical reactions), the tracer mass balance in the street yields the following:

$$Q_{\text{vert}} = e \qquad (12)$$

$$\Rightarrow q_{\text{vert}} \text{WL} \frac{C_{\text{street}} - C_{\text{bg}}}{H} = e \qquad (13)$$

$$\Rightarrow q_{\text{vert}} = \frac{eH}{\text{WL}(C_{\text{street}} - C_{\text{bg}})}, \qquad (14)$$

where $Q_{\text{vert}}$ is the vertical flux of pollutant at the roof level for the whole street ($\mu g s^{-1}$), $q_{\text{vert}}$ is the vertical transfer coefficient ($m^2 s^{-1}$), and WL is the exchange surface ($m^2$). Thus, the vertical transfer coefficient can be compared between both models. In Code_Saturne, $q_{\text{vert}}$ is calculated from the emission rate $e$ and the concentration gradient ($\frac{C_{\text{street}} - C_{\text{bg}}}{H}$),
following Eq. (14), and in MUNICH, it is calculated from Eq. (7).

In MUNICH, the horizontal transfer velocity is equal to the street average horizontal wind speed in the street direction $U_{street}$, which is calculated by integrating Eq. (1) between $z = z_0$ and $z = H$. In Code_Saturne, $U_{street}$ is estimated from the wind speed in the y direction averaged over the street mesh cells. In addition to $U_{street}$, the MUNICH vertical profile of the wind speed (Eq. 1) can also be compared to Code_Saturne by averaging the wind speed in the y direction over the street width (x axis).

As the objective of the study is to parameterize the aerodynamic effect of trees, the chemistry and deposition on built and vegetated surfaces are not considered here. Further details on the simulation setup of MUNICH and Code_Saturne and on the comparison of vertical and horizontal transfers in a treeless canyon are presented in Maison et al. (2022).

## 3 Quantification of the tree crown effect on horizontal and vertical transfers by Code_Saturne simulations

To quantify tree effect on the horizontal wind speed along the street and on the vertical transfer coefficient, Code_Saturne simulations are performed for a large range of tree characteristics (LAI2D, CVF, and height ratio $h/H$), as summarized in Table A4. The tree effect is expressed as a relative deviation between the simulations without and with trees.

### 3.1 Tree effect on horizontal transfer

To quantify and compare the effect of tree crowns in Code_Saturne and MUNICH, and thus to overcome eventual differences between the two models observed in a treeless canyon for the horizontal velocity $U_{street}$, a relative deviation (RD) of $U_{street}$, between the simulations with and without trees (%), is computed as follows:

$$RD_{U_{street}} = 100 \times \frac{U_{street} - U_{street0}}{U_{street0}}, \qquad (15)$$

where $U_{street0}$ stands for the average wind velocity in a treeless street, and $U_{street}$ is the average wind velocity in a street with trees, as computed with Code_Saturne. This RD between the Code_Saturne simulations with and without trees is shown in Fig. 2a for WC, 2b for IC, and 2c for NC. It shows that RD$U_{street}$ becomes increasingly negative, meaning that $U_{street}$ decreases as the LAI2D and the CVF increase. On the opposite, the tree height has either no impact or a small impact compared to LAI2D and CVF. When the effect of tree height is noticeable, an increase in the tree height induces a decrease in RD.

In the range of the tree characteristics studied, $U_{street}$ is attenuated from 7.3% to 62.3%. The tree effect on $U_{street}$ can be compared between the three canyons. The tree effect

on $U_{street}$ increases as the canyon is deeper, highlighting the complex interaction between the street dimensions and the tree effect on the velocity. This observation is consistent with the study of Wania et al. (2012).

### 3.2 Tree effect on vertical transfer

The relative deviation of the vertical transfer coefficient (RD$_{q_{vert}}$) between simulations with and without trees is introduced to quantify the tree effect on the vertical transfer coefficient. Similar to the relative deviation of $U_{street}$ (Eq. 15), the relative deviation RD$_{q_{vert}}$ (%) is expressed as follows:

$$RD_{q_{vert}} = 100 \times \frac{q_{vert} - q_{vert0}}{q_{vert0}}, \qquad (16)$$

where $q_{vert0}$ stands for the vertical transfer coefficient in a treeless street, and $q_{vert}$ is the vertical transfer coefficient in a street with trees. Code_Saturne RD$_{q_{vert}}$ is plotted for different tree parameters LAI2D, CVF, and height ratio $h/H$ in Fig. 3a for WC, 3b for IC, and 3c for NC.

For WC, RD$_{q_{vert}}$ increases with tree LAI2D, CVF, and height ratio (Fig. 3a). For IC, RD$_{q_{vert}}$ also increases with CVF and height ratio, but RD$_{q_{vert}}$ tends to slightly decrease as the LAI increase when the LAI is high and the ratio $h/H$ is small (differences between LAI = 3 and 4 for $h/H = 0.36$ and 0.50; Fig. 3b). For NC, a small increase in RD$_{q_{vert}}$ with LAI2D and CVF is observed for $h/H = 0.65$. However, for the two other smaller height ratios, the tree effect is very low ($-1.1 \le RD_{q_{vert}} \le 0.5$; Fig. 3c).

This reduction in the tree effect on $q_{vert}$ when LAI or CVF increases for small $h/H$ ratios in IC and NC can be explained by micro-scale effects (modified air flow path and turbulent viscosity) and is left out of the scope of the present study due to the corresponding low amplitude. Besides, investigation of such micro-scale effects would probably require more advanced turbulence models, for example, switching from a first-order model, $k-\epsilon$, to a second-order one (Rij-SSG; Speziale et al., 1991) or a 3D large eddy simulation (LES).

Note that, unlike $U_{street}$ (Sect. 3.1), the tree effect on $q_{vert}$ decreases when the canyon becomes deeper. For example, the tree effect, as quantified by RD$_{q_{vert}}$, ranges between $-1.0\%$ and $-20.3\%$ for WC, $-1.0\%$ and $-18.7\%$ for IC, and $0.5\%$ and $-2.7\%$ for NC. The averaged tree effect is relatively less strong when the canyon is deeper. In other words, regarding vertical transfers, the street effect dominates over the tree effect.

For vertical transfer, a wind perpendicular to the street ($\varphi = 90^\circ$) is used to focus on the effects on vertical transfers only. In this case of perpendicular wind, the air flow occurring in street canyons is complex and leads to heterogeneous tracer concentration in the street. In general, for the three canyons, the presence of the two tree crowns tends to increase the tracer concentration on the leeward side of the street and to decrease it on the windward side (Buccolieri et al., 2009; Gromke and Ruck, 2012). Depending on the
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Relative deviation (RD) of $U_{\text{street}}$ computed from Code_Saturne simulations for different tree LAI<sub>2D</sub>, CVF, and height ratio for WC (a), IC (b), and NC (c). The graphic is divided into three columns corresponding to the three height ratios, and higher CVF correspond to darker colors. $|RD_{U_{\text{street}}}|$ values are specified with data labels.</figcaption></figure>

street aspect ratio and, therefore, on the air flow regime (Oke, 1988; Harman et al., 2004) and on the tree characteristics, the tree crown effect on local tracer concentration (leeward versus windward side) is more or less important. For example, for NC, the flow regime is skimming, and on average, the variation in $C_{\text{street}}$ due to the presence of trees compensates between the two sides of the street, explaining why $RD_{q_{\text{vert}}}$ is very low for NC.

This section demonstrated that the tree effect parameterizations should depend on tree characteristics and also on building characteristics to account for the building–tree interactions. The next section aims to parameterize in MUNICH the tree effect on $U_{\text{street}}$ and $q_{\text{vert}}$ observed in Code_Saturne simulations.

# 4 Parameterization of the aerodynamic effect of tree crowns in MUNICH

## 4.1 Model description

The MUNICH parameterizations of horizontal and vertical transfers detailed in Maison et al. (2022) are modified to take into account the tree effects. These parameterizations are based on Wang (2012, 2014) equations, which were originally developed for homogeneous vegetated cover. To remain consistent with this hypothesis, the parameterization will depend on the homogeneous leaf area index in the street, denoted as LAI<sub>street</sub>. As for the conversion from LAI<sub>3D</sub> to LAI<sub>2D</sub>, LAI<sub>street</sub> is estimated from LAI<sub>2D</sub> conserving the leaf surface. LAI<sub>street</sub> is calculated by spreading the tree crown
<figure><img src="image_7.png" /><figcaption>Figure 3. Relative deviation (RD) of $q_{vert}$ computed from Code_Saturne simulations for different tree LAI2D, CVF, and height ratio for WC (a), IC (b), and NC (c). The graphic is divided into three columns corresponding to the three height ratios, and higher CVF correspond to darker colors. $|RD_{q_{vert}}|$ values are specified with data labels.</figcaption></figure>

cylindrical LAI2D over the whole street width as follows:

$$ LAI_{street} \times S_{street} = LAI_{2D} \times S_{2D} \qquad (17) $$

$$ \begin{aligned} \Rightarrow LAI_{street} &= LAI_{2D} \times \frac{S_{2D}}{S_{street}} = \frac{2rLnLAI_{2D}}{WL} \\ &= \frac{2rLnLAI_{2D}}{W}, \end{aligned} \qquad (18) $$

where $S_{street}$ is the soil-projected area of the street homogeneous tree crown, and $S_{2D}$ is the soil-projected area of the Code_Saturne 2D cylindrical tree crown (in m²). Note that, regardless of the tree crown geometry, LAI$_{street}$ is always equal to the street total leaf surface divided by the street ground area (WL).

To account for the tree effect on the mixing length, the characteristic length of the trees ($l_{ct}$ in m) is added into the

equation describing the mixing length $l_m$ as follows:

$$ \frac{1}{l_m} = \frac{1}{\kappa H} + \frac{1}{l_{cb}} + \frac{1}{l_{ct} f_{b \times t}}, \qquad (19) $$

$$ \text{with } l_{ct} = \frac{E_t H}{C_D \frac{1}{2} LAI_{street}}, \qquad (20) $$

where $\frac{1}{2} LAI_{street}$ corresponds to the leaf frontal area density, assuming a random leaf orientation distribution, $E_t$ is a proportionality constant taken equal to 0.054 for vegetated cover, as suggested by Wang (2014), and $f_{b \times t}$ is a function parameterized based on Code_Saturne simulations in Sect. 4.2 and representing the interaction between buildings, trees, and the tree crown height.

The dimensionless factor $s_H$ expressing the effects of the tree canopy on the mixing length $l_m$ is calculated using its definition in Eq. (6) and the expression of $l_m$ in Eq. (19) as
follows:

$$l_m = \kappa H s_H \quad \text{at roof level} \qquad (21)$$

$$\text{with } s_H = \begin{cases} \frac{l_{cb}}{l_{cb} + \kappa H} & \text{without tree} \\ \frac{l_{cb} l_{ct} f_{b \times t}}{\kappa H (l_{cb} + l_{ct} f_{b \times t}) + l_{cb} l_{ct} f_{b \times t}} & \text{with trees.} \end{cases} \qquad (22)$$

The simulations with and without trees have to be distinguished to avoid any convergence issue since $l_{ct}$ tends to $+\infty$ when LAI<sub>street</sub> tends to 0. This $s_H$ factor now includes rough soil, buildings, and tree effects on the mixing length and is then used in the calculation of the attenuation coefficient $\alpha$ (Eq. 4). In the numerator, the expressions for buildings and trees are added as follows:

$$\alpha = \frac{C_B a_r + C_{D_t} C_u \frac{1}{2} \text{LAI}_{\text{street}}}{\kappa s_H}, \qquad (23)$$

where $C_{D_t}$ is the tree drag coefficient (dimensionless) taken equal to the one, as used in Code_Saturne ($C_{D_t} = 0.2$). Wang (2014) presents $C_u$ as a dimensionless coefficient homogeneous on the vertical axis, but that can depend on canopy features. This coefficient has to be determined based on experimental or simulated observation data and will be parameterized in Sect. 4.2 with CFD simulations.

## 4.2 Parameter determination based on Code_Saturne simulations

There are two parameters introduced in $l_m$ and $\alpha$ equations, i.e., $f_{b \times t}$ and $C_u$ (Eqs. 19 and 23), that have to be determined based on Code_Saturne RD<sub>qvert</sub> and RD<sub>Ustreet</sub>. The function $f_{b \times t}$ depends on $a_r$ to account for the interaction between trees and buildings observed in Fig. 3 and also on $h_{max}$ to account for the effect of the tree height. Note that, since tree crowns are assumed to be homogeneous within the canopy in the original Wang formulation, the tree crown height $h_{max}$ is not taken into account (see Eq. 20), and therefore, it needs to be included in the function $f_{b \times t}$. The $f_{b \times t}$ expression is determined by maximizing the fit between Code_Saturne RD<sub>qvert</sub> and MUNICH RD<sub>qven</sub> as follows:

$$f_{b \times t} = \frac{a_0 + a_1 \exp(a_2 a_r)}{(h_{\max}/H)^2} \qquad (24)$$
with $a_0 = 3.26$, $a_1 = 0.0256$ and $a_2 = 6.70.

A comparison of Code_Saturne RD<sub>qvert</sub> and the RD<sub>qvert</sub> parameterized with Eqs. (7), (19), (20), and (24) is presented for the three canyons in Fig. 4a, c, e. The normalized mean absolute error (NMAE) and bias (NMB) are calculated to compare Code_Saturne and MUNICH results (see Appendix C for the definition of the statistical indicators). Figure 4a, c, e show a good agreement between Code_Saturne RD<sub>qvert</sub> and the parameterized RD<sub>qven</sub> because the $f_{b \times t}$ function was determined to minimize bias. The NMAE can be high for RD<sub>qvert</sub> (up to 63% in NC) because the parameterization does

<table><caption>**Table 2.** Statistical indicators (normalized mean absolute error, NMAE, and bias, NMB, in percent) for the comparison of Code_Saturne and MUNICH $q_{vert}$ and $U_{street}$.</caption><thead><tr><th rowspan="2">Canyon</th><th colspan="2">q<sub>vert</sub></th><th colspan="2">U<sub>street</sub></th></tr><tr><th>NMAE</th><th>NMB</th><th>NMAE</th><th>NMB</th></tr></thead><tbody><tr><td>WC</td><td>2.2</td><td>-2.1</td><td>3.4</td><td>-1.1</td></tr><tr><td>IC</td><td>3.5</td><td>2.8</td><td>4.4</td><td>-1.7</td></tr><tr><td>NC</td><td>4.1</td><td>-4.1</td><td>6.8</td><td>-2.8</td></tr></tbody></table>

not reproduce the slight decrease in RD<sub>qvert</sub> when LAI<sub>2D</sub> increases from LAI<sub>2D</sub> = 3 to 4 when $h/H \approx 1/3$ and $1/2$ (see Sect. 3.2 and Fig. 3b and c). But this is not an issue for NC since the $q_{vert}$ values are low.

For $U_{street}$, the parameter $C_u$ also needs to be determined. A constant value of $C_u = 6.7$ is sufficient to obtain a good fit between Code_Saturne RD<sub>Ustreet</sub> and MUNICH RD<sub>Ustreet</sub>. No dependency of $C_u$ on building or tree features is needed, as shown in Fig. 4b, d, f, which compares Code_Saturne RD<sub>Ustreet</sub> and the parameterized RD<sub>Ustreet</sub>. Note that the parameterizations of $l_{ct}$ and $f_{b \times t}$ impact not only the vertical transfers but also the horizontal wind speed because they are involved in the calculation of the $s_H$ factor and, hence, of the $\alpha$ coefficient. The good comparisons between Code_Saturne RD<sub>Ustreet</sub> and MUNICH RD<sub>Ustreet</sub> also show that the parameterized $s_H$ reproduces the horizontal wind speed well and, therefore, also the interactions between trees and building effects and the influence of $h_{max}$.

## 4.3 Comparison of $q_{vert}$, $U_{street}$, and wind profiles

Figure 5a and b present, respectively, a comparison of Code_Saturne and parameterized $q_{vert}$ and $U_{street}$. The statistical indicators are presented in Table 2. As MUNICH was parameterized to reproduce the tree effect observed in Code_Saturne well, and as the two models agree well in a treeless canyon (Maison et al., 2022), the vertical transfer coefficient and wind speed with trees are close between the two models, as expected (Fig. 5 and Table 2). The parameterized vertical transfer coefficients with trees agree well with the Code_Saturne ones, with normalized mean absolute errors ranging from 2.2% to 4.1% and bias from -4.1% to 2.8%. The parameterized average wind speed with trees agree well with Code_Saturne ones, with normalized mean absolute errors ranging from 3.4% to 6.8% and a normalized mean bias from -2.8% to -1.1%.

For each street canyon, Fig. 6 compares the Code_Saturne and parameterized vertical wind profiles for fixed CVF and $h/H$ ratio but for five different LAI<sub>2D</sub>.

Figure 6 shows that an increase in the LAI<sub>2D</sub> induces a decrease in the wind velocity. Concerning the vertical profile shape, trees induce a lower wind velocity on the entire street and not only in the tree crown. They even slightly
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Comparison of (a, c, e) RD<sub>qvert</sub> and (b, d, f) RD<sub>Ustreet</sub> computed from Code_Saturne simulations and parameterized in MUNICH for different tree LAI<sub>2D</sub>, CVF, and height ratio and for (a, b) WC, (c, d) IC, and (e, f) NC.</figcaption></figure>

impact the velocity just above the street. The maximum of attenuation of the wind is located in the middle of the tree crown and the wind velocity is re-accelerated under the tree crown. In the middle of the tree crown, the parameterized wind speed is close to the one of Code_Saturne. In the lower

part and under the tree crown, the parameterized wind speed is underestimated. In fact, the re-acceleration under the tree crown is complex to consider in parameterized models. Besides, in real streets the tree height is not homogeneous, so this re-acceleration under the tree crown might be unrealized.
<figure><img src="image_9.png" /><figcaption>Figure 5. Normalized $U_{street}$ and $q_{vert}$ computed from Code_Saturne simulations for different tree LAI<sub>2D</sub>, CVF, and height ratios.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Figure 6. Comparison of Code_Saturne (solid lines) and parameterized (dotted lines) vertical wind profiles. The profiles in a treeless canyon are in black and darker colors correspond to increasing LAI<sub>2D</sub>. For the three canyons, $h/H \approx 0.5$ and CVF $\approx 25\%$.</figcaption></figure>

In Code_Saturne simulations. In the parameterized wind profile with trees, the reduction in the wind speed just above the street due to the presence of trees was neglected. So the parameterized wind speed above and in the upper part of the tree crown is overestimated compared to Code_Saturne. Note that Kent et al. (2017) proposed a method to include the vegetation effect in aerodynamic roughness parameters and so to account for vegetation in the above-urban canopy wind profile of Macdonald et al. (1998). This method considers vegetation at the city scale and is not applicable here, since we work at the street scale. For now, as MUNICH assumes a homogeneous street canyon, the profile is averaged to compute $U_{street}$, and the shape of the profile is not used.

Moreover, to illustrate the impact of the developed parametrization on pollutant concentrations, MUNICH simulations are performed with carbon monoxide emissions in IC without and with trees of various characteristics. The results are presented in Appendix D.

# 5 Conclusions

Although the discretizations and the physical modeling used in MUNICH and Code_Saturne are fundamentally different, the setup of the CFD simulations was adapted to compare the two models through the average horizontal wind speed along the street and the calculation of the vertical transfer coefficient at the roof level by averaging a passive tracer concentration in the street and the background (above the streets). To build a parameterization suited to the street network model MUNICH, the CFD simulations were simplified with several hypothesis such as an infinite street canyon with 2D CFD simulations and homogeneous emissions over the street. Furthermore, quantities of interest were averaged over the street in the CFD simulations.

The aerodynamic effect of trees in street canyons is quantified with Code_Saturne simulations. Simulations were performed with two rows of trees and a large range of tree
leaf area index, crown radius, and heights. In the range of the simulations performed, the average horizontal velocity in the street is reduced by -7.3% to -62.3% and the vertical transfer coefficient by +0.5% to -20.3%. This highlights the necessity of adding the tree aerodynamic effects into street models such as MUNICH.

Maison et al. (2022) proposed a parameterization of horizontal and vertical transfers in a treeless canyon, and the present study adds the tree effect to the building effect in this parameterization. The differences in wind speed and vertical transfer coefficient in a treeless canyon are very low between the two models (Maison et al., 2022). However, to overcome any model-specific difference and to build a parameterization for MUNICH of the tree effect on the aerodynamic parameters, relative deviation ratios between cases without and with trees are defined.

First, the Maison et al. (2022) vertical transfer coefficient parameterization is modified by adding a term representing the tree effect and the interaction between the trees and the buildings in the mixing-length expression. This term is a function of two dimensionless parameters to characterize the trees, street leaf area index (LAI<sub>street</sub>) and tree-to-street height ratio (h<sub>max</sub>/H), and also of street characteristics, building height (H), and street width (W). A parameterization of the tree mixing length was defined from Code_Saturne simulations to obtain a tree effect as close as possible in the two models. The parameterized vertical transfer coefficients with trees agree well with those of Code_Saturne, with the normalized mean absolute error ranging from 2.2% to 4.1% and bias from -4.1% to 2.8%.

Second, this new mixing length expression is also used to compute the vertical wind profile and average wind speed along the street. Only one constant is fixed in the wind attenuation coefficient to maximize the fit between Code_Saturne and MUNICH tree effect on average wind speed. The comparison of the average wind speed gives a normalized mean absolute error ranging from 3.4% to 6.8% and a normalized mean bias from -2.8% to -1.1%.

MUNICH now includes a relatively simple parameterization of the tree effect on both horizontal and vertical aerodynamic processes, based on commonly used tree parameters that can be easily computed from urban databases and can reproduce the main effects obtained in much more detailed and costly CFD simulations. This parameterization can also be used in urban climate models to compute water and heat transfer in tree-lined streets. The perspectives of this study are to quantify the effect of street trees on air quality from the street level to the scale of the city of Paris. Dry deposition of gaseous pollutants and aerosols on tree leaves and emission of biogenic organic volatile compounds related to tree water stress will be considered in MUNICH. The contribution of biogenic and anthropic precursors to the formation of organic aerosols over an entire city will be compared.

## Appendix A: Lists of abbreviations, variable parameters, and tree dimensions

Each line of the Table A4 corresponds to the parameters of a simulation. For each case, simulations with different LAI<sub>2D</sub> were performed as follows: 0.5, 1, 2, 3, and 4.

<table><caption>Table A1. List of abbreviations.</caption><thead><tr><td>Acronym</td><td>Definition</td></tr></thead><tbody><tr><td>WC</td><td>Wide canyon</td></tr><tr><td>IC</td><td>Intermediate canyon</td></tr><tr><td>NC</td><td>Narrow canyon</td></tr><tr><td>CTM</td><td>Chemistry transport model</td></tr><tr><td>CFD</td><td>Computational fluid dynamics</td></tr><tr><td>RD</td><td>Relative deviation</td></tr><tr><td>BVOC</td><td>Biogenic volatile organic compound</td></tr></tbody></table>
<table><caption>**Table A2.** List of parameters.</caption><thead><tr><th>Symbol</th><th>Definition</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>κ</td><td>Von Kàrmàn constant</td><td>0.42</td><td>-</td></tr><tr><td>PBLH</td><td>Planetary boundary layer height</td><td>1000</td><td>m</td></tr><tr><td>z<sub>0s</sub></td><td>Code_Saturne inside street walls roughness length</td><td>0.10</td><td>m</td></tr><tr><td>u<sub>*</sub></td><td>Friction velocity</td><td>0.727</td><td>m s<sup>-1</sup></td></tr><tr><td>E</td><td>Parameter in modified Wang (2014) parameterization</td><td>0.5</td><td>-</td></tr><tr><td>n</td><td>Number of tree rows</td><td>2</td><td>-</td></tr><tr><td>C<sub>Dt</sub></td><td>Tree crown drag coefficient (Katul and Albertson, 1998)</td><td>0.2</td><td>-</td></tr></tbody></table>

<table><caption>**Table A3.** List of variables.</caption><thead><tr><th>Group of variables</th><th>Symbol</th><th>Definition</th><th>Unit</th></tr></thead><tbody><tr><td rowspan="5">Street characteristics</td><td>H</td><td>Buildings height</td><td>m</td></tr><tr><td>W</td><td>Street width</td><td>m</td></tr><tr><td>L</td><td>Street length</td><td>m</td></tr><tr><td>a<sub>r</sub></td><td>Street aspect ratio</td><td>-</td></tr><tr><td>l<sub>cb</sub></td><td>Characteristic length in the street</td><td>m</td></tr><tr><td rowspan="14">Tree characteristics</td><td>LAD</td><td>Tree leaf area density</td><td>m<sup>2</sup><sub>leaf</sub> m<sup>-3</sup><sub>tree crown</sub></td></tr><tr><td>LAI</td><td>Tree leaf area index</td><td>m<sup>2</sup><sub>leaf</sub> m<sup>-2</sup><sub>soil</sub></td></tr><tr><td>LAI<sub>3D</sub></td><td>3D leaf area index of the spherical tree crown</td><td>m<sup>2</sup><sub>leaf</sub> m<sup>-2</sup><sub>soil</sub></td></tr><tr><td>LAI<sub>2D</sub></td><td>2D equivalent leaf area index of the cylindrical tree crown</td><td>m<sup>2</sup><sub>leaf</sub> m<sup>-2</sup><sub>soil</sub></td></tr><tr><td>LAI<sub>street</sub></td><td>Leaf area index of the homogeneous street tree crown</td><td>m<sup>2</sup><sub>leaf</sub> m<sup>-2</sup><sub>soil</sub></td></tr><tr><td>S<sub>3D</sub></td><td>Soil projected area of the 3D spherical tree crown</td><td>m<sup>2</sup></td></tr><tr><td>S<sub>2D</sub></td><td>Soil projected area of the 2D cylindrical tree crown</td><td>m<sup>2</sup></td></tr><tr><td>S<sub>street</sub></td><td>Soil projected area of the street homogeneous tree crown</td><td>m<sup>2</sup></td></tr><tr><td>r</td><td>Tree radius</td><td>m</td></tr><tr><td>CVF</td><td>Crown volume fraction</td><td>- or %</td></tr><tr><td>h</td><td>Middle crown height</td><td>m</td></tr><tr><td>h<sub>max</sub></td><td>Maximum tree crown height</td><td>m</td></tr><tr><td>h<sub>min</sub></td><td>Minimum tree crown height</td><td>m</td></tr><tr><td>δ</td><td>Spacing between two trees within a row</td><td>m</td></tr><tr><td>l<sub>ct</sub></td><td>Tree characteristic length</td><td>m</td></tr><tr><td rowspan="10">Horizontal wind speed</td><td>U<sub>street</sub></td><td>Average street horizontal wind speed</td><td>m s<sup>-1</sup></td></tr><tr><td>U</td><td>Norm of the horizontal wind speed</td><td>m s<sup>-1</sup></td></tr><tr><td>U<sub>x</sub></td><td>Horizontal wind speed in the x direction</td><td>m s<sup>-1</sup></td></tr><tr><td>U<sub>y</sub></td><td>Horizontal wind speed in the y direction</td><td>m s<sup>-1</sup></td></tr><tr><td>U<sub>H</sub></td><td>Average horizontal wind speed at roof level</td><td>m s<sup>-1</sup></td></tr><tr><td>φ</td><td>Angle of the wind direction</td><td>rad or °</td></tr><tr><td>α</td><td>Wind attenuation coefficient</td><td>-</td></tr><tr><td>s<sub>H</sub></td><td>Characteristic length factor</td><td>-</td></tr><tr><td>C<sub>u</sub></td><td>Empiric coefficient in α equation</td><td>-</td></tr><tr><td>C<sub>B</sub></td><td>Function of a<sub>r</sub> and φ</td><td>-</td></tr><tr><td rowspan="7">Vertical transfer</td><td>Q<sub>vert</sub></td><td>Vertical flux of pollutant</td><td>µg s<sup>-1</sup></td></tr><tr><td>q<sub>vert</sub></td><td>Vertical transfer coefficient</td><td>m<sup>2</sup> s<sup>-1</sup></td></tr><tr><td>σ<sub>W</sub></td><td>Standard deviation of the vertical wind speed at z = H (Soulhac et al., 2011)</td><td>m s<sup>-1</sup></td></tr><tr><td>l<sub>m</sub></td><td>Mixing length in the street</td><td>m</td></tr><tr><td>e</td><td>Passive tracer emission rate for the street</td><td>µg s<sup>-1</sup></td></tr><tr><td>C<sub>street</sub></td><td>Street concentration</td><td>µg m<sup>-3</sup></td></tr><tr><td>C<sub>bg</sub></td><td>Background concentration</td><td>µg m<sup>-3</sup></td></tr></tbody></table>
<table><caption>**Table A4.** List of tree dimensions simulated with Code_Saturne (r is the radius, h is the middle crown height, $h_{min}$ is the crown bottom height, and $h_{max}$ is the crown top height) and calculated tree parameters (CVF is the crown volume fraction, and $h/H$ is the height ratio).</caption><thead><tr><th>Canyon</th><th>r (m)</th><th>h (m)</th><th>h<sub>min</sub> (m)</th><th>h<sub>max</sub> (m)</th><th>CVF (%)</th><th>h/H (-)</th></tr></thead><tbody><tr><td rowspan="10">WC</td><td>1.5</td><td>3.5</td><td>2.0</td><td>5.0</td><td>6.0</td><td>0.41</td></tr><tr><td>1.5</td><td>4.5</td><td>3.0</td><td>6.0</td><td>6.0</td><td>0.53</td></tr><tr><td>1.5</td><td>5.5</td><td>4.0</td><td>7.0</td><td>6.0</td><td>0.65</td></tr><tr><td>2.0</td><td>3.5</td><td>1.5</td><td>5.5</td><td>10.8</td><td>0.41</td></tr><tr><td>2.0</td><td>4.5</td><td>2.5</td><td>6.5</td><td>10.8</td><td>0.53</td></tr><tr><td>2.0</td><td>5.5</td><td>3.5</td><td>7.5</td><td>10.8</td><td>0.65</td></tr><tr><td>3.0</td><td>3.5</td><td>0.5</td><td>6.5</td><td>24.2</td><td>0.41</td></tr><tr><td>3.0</td><td>4.5</td><td>1.5</td><td>7.5</td><td>24.2</td><td>0.53</td></tr><tr><td>3.0</td><td>5.5</td><td>2.5</td><td>8.5</td><td>24.2</td><td>0.65</td></tr><tr><td rowspan="10">IC</td><td>2.0</td><td>5.0</td><td>3.0</td><td>7.0</td><td>6.5</td><td>0.36</td></tr><tr><td>2.0</td><td>7.0</td><td>5.0</td><td>9.0</td><td>6.5</td><td>0.50</td></tr><tr><td>2.0</td><td>9.0</td><td>7.0</td><td>11.0</td><td>6.5</td><td>0.64</td></tr><tr><td>2.5</td><td>5.0</td><td>2.5</td><td>7.5</td><td>10.2</td><td>0.36</td></tr><tr><td>2.5</td><td>7.0</td><td>4.5</td><td>9.5</td><td>10.2</td><td>0.50</td></tr><tr><td>2.5</td><td>9.0</td><td>6.5</td><td>11.5</td><td>10.2</td><td>0.64</td></tr><tr><td>4.0</td><td>5.0</td><td>1.0</td><td>9.0</td><td>26.1</td><td>0.36</td></tr><tr><td>4.0</td><td>7.0</td><td>3.0</td><td>10.0</td><td>26.1</td><td>0.50</td></tr><tr><td>4.0</td><td>9.0</td><td>5.0</td><td>13.0</td><td>26.1</td><td>0.64</td></tr><tr><td rowspan="10">NC</td><td>2.5</td><td>9.0</td><td>6.5</td><td>11.5</td><td>5.2</td><td>0.33</td></tr><tr><td>2.5</td><td>14.0</td><td>11.5</td><td>16.5</td><td>5.2</td><td>0.51</td></tr><tr><td>2.5</td><td>18.0</td><td>15.5</td><td>20.5</td><td>5.2</td><td>0.65</td></tr><tr><td>3.5</td><td>9.0</td><td>5.5</td><td>12.5</td><td>10.2</td><td>0.33</td></tr><tr><td>3.5</td><td>14.0</td><td>10.5</td><td>17.5</td><td>10.2</td><td>0.51</td></tr><tr><td>3.5</td><td>18.0</td><td>14.5</td><td>21.5</td><td>10.2</td><td>0.65</td></tr><tr><td>5.5</td><td>9.0</td><td>3.5</td><td>14.5</td><td>25.1</td><td>0.33</td></tr><tr><td>5.5</td><td>14.0</td><td>8.5</td><td>19.5</td><td>25.1</td><td>0.51</td></tr><tr><td>5.5</td><td>18.0</td><td>12.5</td><td>23.5</td><td>25.1</td><td>0.65</td></tr></tbody></table>

## Appendix B: Comparison of 2D and 3D Code_Saturne simulations with trees

Code_Saturne simulations were performed in a periodic 2D canyon of length $L = 1$ m, and the street was therefore considered infinite due to this periodicity. In this case, the tree crown is represented as an infinite cylinder of radius $r$. To take into account the fact that, in reality, most of the tree crowns are spherical and are spaced from each other, an equivalent cylindrical LAI<sub>2D</sub> is calculated in the 2D simulations. For the tree effect to be similar, the surface of the leaves is kept constant between the 2D and 3D simulations, as follows:

$$LAI_{3D} \times S_{3D} = LAI_{2D} \times S_{2D} \qquad (B1)$$

$$LAI_{2D} = LAI_{3D} \times \frac{S_{3D}}{S_{2D}} = LAI_{3D} \times \frac{\pi r^2}{2r(2r + \delta)}, \qquad (B2)$$

where $\delta$ is the spacing between two tree crowns within the row, and $2r + \delta$ is the length of the street section for one spherical tree (m). $S_{3D}$ and $S_{2D}$ are, respectively, the soil-projected area of the 3D and 2D tree crowns (m$^2$). In addition to the leaf surface and tree height, the tree radius is also conserved between 2D and 3D simulations; however, the value of the crown volume fraction (CVF) varies, as follows:

$$\begin{gathered} \text{CVF}_{2D} = \frac{V_{2D}}{V_{\text{street}}} = \frac{n\pi r^2 L}{HWL} \\ \text{and } \text{CVF}_{3D} = \frac{V_{3D}}{V_{\text{street}}} = \frac{n\frac{4}{3}\pi r^3}{HW(2r + \delta)}. \end{gathered} \qquad (B3)$$

$V_{\text{street}}$ is the street volume (m$^3$), and $V_{3D}$ and $V_{2D}$ are, respectively, the 3D and 2D tree crown volumes (m$^3$). The relation between 2D and 3D CVF is therefore as follows:

$$\text{CVF}_{2D} = \frac{3}{4} \times \frac{2r + \delta}{r} \times \text{CVF}_{3D}. \qquad (B4)$$
2D and 3D simulations with different LAI are performed to check the validity of the 2D infinite street canyon hypothesis. 3D simulations are realized with a periodic street with a 10 m length and cell meshes of 0.5 × 0.5 × 0.5 m in the x, y, and z directions. In both 2D and 3D simulations, $C_{\text{street}}$ and $C_{\text{bg}}$ are computed with a normal wind to the street ($\varphi = 90^{\circ}$) and $U_{\text{street}}$ with $\varphi = 45^{\circ}$. The results are presented in Figs. B1a–d for WC, B1e–h for IC, and B1i–l for NC.

<figure><img src="image_11.png" /><figcaption>**Figure B1.** Comparison of $C_{\text{street}}$ (a, e, i), $C_{\text{bg}}$ (b, f, j), $q_{\text{vert}}$ (c, g, k), and $U_{\text{street}}$ (d, h, l) for 2D and 3D simulations in the three canyons without and with trees of LAI$_{2D}$ = 1 and 4. CVF$_{2D}$ ≈ 10% and h/H ≈ 0.5 are fixed. The percentages written in the middle of the histogram bars correspond to the RD between the 2D and 3D cases.</figcaption></figure>
For each canyon, the cases without tree (LAI = 0) are equivalent. C<sub>street</sub>, C<sub>bg</sub>, q<sub>vert</sub>, and U<sub>street</sub> are similar between the 2D and 3D cases, with a maximum relative difference of 0.1 % for U<sub>street</sub>. For the cases with trees, concerning C<sub>street</sub>, C<sub>bg</sub>, and q<sub>vert</sub>, the differences between the 2D and 3D cases are relatively low over the range of simulations performed with a maximum of 3.3 % for C<sub>street</sub> and C<sub>bg</sub> in WC. For U<sub>street</sub>, slightly larger differences are observed (up to -5.4 % in WC for LAI = 4(2D)/12.7(3D)). These differences can be explained by the distribution of the turbulence around the crown. For each variable and canyon, the relative deviation between the 2D and 3D cases increases with the tree LAI. Given the low differences observed between the 2D and 3D simulations, the hypothesis of a 2D canyon with a cylindrical tree crown associated with an equivalent LAI<sub>2D</sub> to represent a spherical tree crown spaced of δ m in a street of length L is reasonable.

## Appendix C: Definition of the statistical indicators

Code_Saturne and MUNICH simulations are denoted $cs_i$ and $m_i$, respectively. In this section, $n$ is the total number of simulations, which is equal to 45 per street canyon.

- Normalized mean absolute error (%):

$$NMAE = 100 \times \frac{\sum_{i=1}^{n} |m_i - cs_i|}{\left| \sum_{i=1}^{n} cs_i \right|} \quad (C1)$$

- Normalized mean bias (%):

$$NMB = 100 \times \frac{\sum_{i=1}^{n} (m_i - cs_i)}{\sum_{i=1}^{n} cs_i} \quad (C2)$$

## Appendix D: Aerodynamic tree effect on concentration

To illustrate the impact of the developed parameterization on pollutant concentrations, MUNICH simulations are performed in the Intermediate Canyon (IC) with carbon monoxide (CO) emissions in the street. The input parameters are detailed in Table D1.

In the simulations, only the tree aerodynamic effect is considered, and there is no chemistry, no deposition of pollutants on urban or tree surfaces, and no biogenic volatile organic compound (BVOC) emission. One street is modeled without any inflow of pollutant from nearby streets. The following three processes are taken into account: pollutant emission in the street, pollutant outflow by horizontal transfer, $Q_{outflow}$ (function of the parameterized U<sub>street</sub>), and vertical transfer

<table><caption>Table D1. List of parameters fixed in MUNICH simulation.</caption><thead><tr><th>Fixed parameter</th><th>Symbol</th><th>Value</th><th>Unit</th></tr></thead><tbody><tr><td>Building height</td><td>H</td><td>14</td><td>m</td></tr><tr><td>Street width</td><td>W</td><td>27.5</td><td>m</td></tr><tr><td>Street length</td><td>L</td><td>200</td><td>m</td></tr><tr><td>Wind angle</td><td>φ</td><td>45</td><td>°</td></tr><tr><td>Friction velocity</td><td>u*</td><td>0.7</td><td>ms<sup>-1</sup></td></tr><tr><td>CO emission rate</td><td>e</td><td>1000</td><td>µg s<sup>-1</sup> m<sup>-1</sup></td></tr><tr><td>CO background concentration</td><td>C<sub>bg</sub></td><td>100</td><td>µg m<sup>-3</sup></td></tr></tbody></table>

between the street and the background, $Q_{vert}$ (function of the parameterized $q_{vert}$).

The street average CO concentrations ($C_{street}$) are compared in Fig. D1 for a canyon without trees (Fig. D1a) and for a canyon with trees of various LAI, CVF, and height ratios $h/H$ (Fig. D1b, c).

Figure D1 shows that the increase in tree LAI and CVF induces higher average street CO concentrations. The effect of the tree height ratio $h/H$ is negligible for low LAI and increases when the LAI increases. The variations in the concentrations are due to both the effect of trees on the horizontal transfer velocity, U<sub>street</sub>, and on the vertical transfer coefficient, $q_{vert}$.

In this example, the aerodynamical effect of trees leads to a large increase in the CO concentration from 4 % to 45 %. However, to simulate the global effect of trees on air quality at street level, other processes such as dry deposition on leaves and BVOC emissions need to be considered.
<figure><img src="image_12.png" /><figcaption>**Figure D1.** Comparison of CO concentration ($C_{street}$) simulated with MUNICH in IC for different tree LAI$_{2D}$, CVF, and tree height ratios $h/H$. The graphic is divided into three columns corresponding to the treeless canyon (a) and to canyons with trees of height ratios $h/H = 0.36$ (b) and $0.64$ (c). The concentrations $C_{street}$ (in µg m$^{-3}$) are specified with data labels.</figcaption></figure>

**Code availability.** The last version of the MUNICH source code is available online at https://doi.org/10.5281/zenodo.4168984 (Kim et al., 2022) and https://github.com/cerea-lab/munich (last access: 5 January 2022).

**Data availability.** For the three 2D canyons considered in the present study, the mesh, the source code, and the XML setup file allowing the reproduction of the CFD results using Code_Saturne version 6.0 are available online at https://gitlab.enpc.fr/alice.maison/tree_parametrization (last access: 3 December 2021) and at https://doi.org/10.17632/fzfrjsz3mv.2 (Maison and Flageul, 2021) under the GNU GPL2.0 license.

**Author contributions.** KS, CF, BC, and AM were responsible for the conceptualization. AM developed the MUNICH model code and performed the Code_Saturne simulations. AM, KS, CF, and BC performed the formal analysis. AM conducted the visualization. AM, KS, and CF were responsible for writing the original draft, and BC and AT reviewed it. CF, BC, and YW provided support for the Code_Saturne computing resources. KS and AT were responsible for the funding acquisition.

**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 (ACP/GMD inter-journal SI)". It is not associated with a conference.

**Acknowledgements.** The authors thank Youngseob Kim and Lya Lugon, for their support in the development of the MUNICH model, and Martin Ferrand, for his support in the understanding of the Code_Saturne model.

**Financial support.** This work has partially been funded by the sTREEt ANR project (grant no. ANR-19-CE22-0012), the DIM Qi2 (Air Quality Research Network on air quality in the Île-de-France region), and Paris Île-de-France Region.

**Review statement.** This paper was edited by Yang Zhang and reviewed by two anonymous referees.

## References

* Akimoto, H.: Global Air Quality and Pollution, Science, 302, 1716–1719, https://doi.org/10.1126/science.1092666, 2003.

* Angel, S., Parent, J., Civco, D. L., Blei, A., and Potere, D.: The dimensions of global urban expansion: Estimates and projections for all countries, 2000–2050, Prog. Plann., 75, 53–107, https://doi.org/10.1016/j progress.2011.04.001, 2011.

* 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, 1–62, https://hal.archives-ouvertes.fr/hal-01115371
(last access: 18 July 2022), 2004.

Armson, D., Stringer, P., and Ennos, A.: The effect of street trees and amenity grass on urban surface water runoff in Manchester, UK, Urban For. Urban Gree., 12, 282–286, https://doi.org/10.1016/j.ufug.2013.04.001, 2013.

Beckett, K., Freer-Smith, P., and Taylor, G.: Urban woodlands: their role in reducing the effects of particulate pollution, Environ. Pollut., 99, 347–360, https://doi.org/10.1016/S0269-7491(98)00016-5, 1998.

Berland, A., Shiflett, S. A., Shuster, W. D., Garmes-
tani, A. S., Goddard, H. C., Herrmann, D. L., and
Hopton, M. E.: The role of trees in urban stormwa-
ter management, Landscape Urban Plan., 162, 167–177,
https://doi.org/10.1016/j.landurbplan.2017.02.017, 2017.

Bertram, C. and Rehdanz, K.: The role of urban green space for human well-being, Ecol. Econ., 120, 139–152, https://doi.org/10.1016/j.ecolecon.2015.10.013, 2015.

Bowler, D. E., Buyung-Ali, L., Knight, T. M., and Pullin, A. S.: Urban greening to cool towns and cities: A systematic review of the empirical evidence, Landscape Urban Plan., 97, 147–155, https://doi.org/10.1016/j.landurbplan.2010.05.006, 2010.

Bozonnet, E., Musy, M., Calmet, I., and Rodriguez, F.: Modeling methods to assess urban fluxes and heat island mitigation measures from street to city scale, International Journal of Low-Carbon Technologies, 10, 62–77, https://doi.org/10.1093/ijlct/ctt049, 2015.

Buccolieri, R., Gromke, C., Di Sabatino, S., and Ruck, B.: Aerodynamic effects of trees on pollutant concentration in street canyons, Sci. Total Environ., 407, 5247–5256, https://doi.org/10.1016/j.scitotenv.2009.06.016, 2009.

Buccolieri, R., Salim, M., Leo, L. S., Di Sabatino, S., Chan, A., Ielpo, P., de Gennaro, G., and Gromke, C.: Analysis of local scale tree–atmosphere interaction on pollutant concentration in idealized street canyons and application to a real urban junction, Atmos. Environ., 45, 1702–1713, https://doi.org/10.1016/j.atmosenv.2010.12.058, 2011.

Buccolieri, R., Santiago, J.-L., Rivas, E., and Sanchez, B.: Review on urban tree modelling in CFD simulations: Aerodynamic, deposition and thermal effects, Urban For. Urban Gree., 31, 212–220, https://doi.org/10.1016/j.ufug.2018.03.003, 2018.

Cai, X.-M., Barlow, J., and Belcher, S.: Dispersion and transfer of passive scalars in and above street canyons – Large-eddy simulations, Atmos. Environ., 42, 5885–5895, https://doi.org/10.1016/j.atmosenv.2008.03.040, 2008.

Calfapietra, C., Fares, S., Manes, F., Morani, A., Sgrigna, G., and Loreto, F.: Role of Biogenic Volatile Organic Compounds (BVOC) emitted by urban trees on ozone concentration in cities: A review, Environ. Pollut., 183, 71–80, https://doi.org/10.1016/j.envpol.2013.03.012, 2013.

Faiz, A.: Automotive emissions in developing countries-
relative implications for global warming, acidification and
urban air quality, Transportation Res. A-Pol., 27, 167–186,
https://doi.org/10.1016/0965-8564(93)90057-R, 1993.

Gillner, S., Vogt, J., Tharang, A., Dettmann, S., and Roloff, A.: Role of street trees in mitigating effects of heat and drought at highly sealed urban sites, Landscape Urban Plan., 143, 33–42, https://doi.org/10.1016/j.landurbplan.2015.06.005, 2015.

Gromke, C. and Blocken, B.: Influence of avenue-trees on air qual-
ity at the urban neighborhood scale. Part II: Traffic pollutant con-
centrations at pedestrian level, Environ. Pollut., 196, 176–184,
https://doi.org/10.1016/j.envpol.2014.10.015, 2015.

Gromke, C. and Ruck, B.: Influence of trees on the dispersion of pollutants in an urban street canyon – Experimental investigation of the flow and concentration field., Atmos. Environ., 41, 3287–3302, https://doi.org/10.1016/j.atmosenv.2006.12.043, 2007.

Gromke, C. and Ruck, B.: On the Impact of Trees on Dispers-
sion Processes of Traffic Emissions in Street Canyons, Bound.-
Lay. Meteorol., 131, 19–34, https://doi.org/10.1007/s10546-008-
9301-2, 2009.

Gromke, C. and Ruck, B.: Pollutant Concentrations in Street Canyons of Different Aspect Ratio with Avenues of Trees for Various Wind Directions, Bound.-Lay. Meteorol., 144, 41–64, https://doi.org/10.1007/s10546-012-9703-z, 2012.

Gu, S., Guenther, A., and Faiola, C.: Effects of Anthropogenic and Biogenic Volatile Organic Compounds on Los Angeles Air Quality, Environ. Sci. Technol., 55, 12191–12201, https://doi.org/10.1021/acs.est.1c01481, 2021.

Gunawardena, K., Wells, M., and Kershaw, T.: Utilis-
ing green and bluespace to mitigate urban heat is-
land intensity, Sci. Total Environ., 584–585, 1040–1055,
https://doi.org/10.1016/j.scitotenv.2017.01.158, 2017.

Harman, I. N., Barlow, J. F., and Belcher, S. E.: Scalar fluxes from urban street canyons. Part II Model, Bound.-Lay. Meteorol., 113, 387–409, https://doi.org/10.1007/s10546-004-6205-7, 2004.

Hebbert, M. and Jankovic, V.: Cities and Climate Change: The Precedents and Why They Matter, Urban Stud., 50, 1332–1347, https://doi.org/10.1177/0042098013480970, 2013.

Huang, Y.-D., Hou, R.-W., Liu, Z.-Y., Song, Y., Cui, P.-Y., and Kim, C.-N.: Effects of Wind Direction on the Airflow and Pollutant Dispersion inside a Long Street Canyon, Aerosol Air Qual. Res., 19, 1152–1171, https://doi.org/10.4209/aaqr.2018.09.0344, 2019.

Hwang, H.-J., Yook, S.-J., and Ahn, K.-H.: Experimental in-
vestigation of submicron and ultrafine soot particle re-
moval by tree leaves, Atmos. Environ., 45, 6987–6994,
https://doi.org/10.1016/j.atmosenv.2011.09.019, 2011.

IPCC: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Report, Intergovernmental Panel on Climate Change, United Nations, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R., Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, https://www.ipcc.ch/report/ar6/wg1/ (last access: 18 July 2022), 2021.

Janhäll, S.: Review on urban vegetation and particle air pollution – Deposition and dispersion, Atmos. Environ., 105, 130–137, https://doi.org/10.1016/j.atmosenv.2015.01.052, 2015.

Jayasooriya, V., Ng, A., Muthukumaran, S., and Perera, B.: Green infrastructure practices for improvement of urban air quality, Urban For. Urban Gree., 21, 34–47, https://doi.org/10.1016/j.ufug.2016.11.007, 2017.

Jeanjean, A., Buccolierib, R., Eddy, J., Monks, P., and Leigh, R.: Air quality affected by trees in real street canyons: The case of Marylebone neighbourhood in
central London., Urban For. Urban Gree., 22, 41–43,
https://doi.org/10.1016/j.ufug.2017.01.009, 2017.

Katul, G. and Albertson, J.: An Investigation of Higher-Order Closure Models for a Forested Canopy, Bound.-Lay. Meteorol., 89, 47–74, https://doi.org/10.1023/A:1001509106381, 1998.

Katul, G. G., Mahrt, L., Poggi, D., and Sanz, C.: ONE- and TWO-Equation Models for Canopy Turbulence, Bound.-Lay. Meteorol., 113, 81–109, https://doi.org/10.1023/B:BOUN.0000037333.48760.e5, 2004.

Kent, C. W., Grimmond, S., and Gatey, D.: Aerodynamic roughness parameters in cities: Inclusion of vegetation, J. Wind Eng. Ind. Aerod., 169, 168–176, https://doi.org/10.1016/j.jweia.2017.07.016, 2017.

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., Sartelet, K., Lugon, L., Roustan, Y., Sarica, T., Maison, A., Valari, M., Zhang, Y., and André, M.: The Model of Urban Network of Intersecting Canyons and Highways (MUNICH), Zenodo [code], https://doi.org/10.5281/zenodo.6167477, 2022.

Klemm, W., Heusinkveld, B. G., Lenzholzer, S., and van Hove, B.: Street greenery and its physical and psychological impact on thermal comfort, Landscape Urban Plan., 138, 87–98, https://doi.org/10.1016/j.landurbplan.2015.02.009, 2015.

Krekel, K., Kolbe, J., and Wüstemann, H.: The Greener, The Happier? The Effects of Urban Green and Abandoned Areas on Residential Well-Being, The German Socio-Economic Panel study at DIW Berlin, 728, 65, https://doi.org/10.2139/ssrn.2554477, 2015.

Leopold, L. B.: Hydrology for urban land planning – A guidebook on the hydrologic effects of urban land use, US. Geological Survey, 554, 18, https://doi.org/10.3133/cir554, 1968.

Li, X., Liu, C., Leung, D., and Lam, K.: Recent progress in CFD modelling of wind field and pollutant transport in street canyons, Atmos. Environ., 40, 5640–5658, https://doi.org/10.1016/j.atmosenv.2006.04.055, 2006.

Livesley, S. J., McPherson, E. G., and Calfapietra, C.: The Urban Forest and Ecosystem Services: Impacts on Urban Water, Heat, and Pollution Cycles at the Tree, Street, and City Scale, J. Environ. Qual., 45, 119–124, https://doi.org/10.2134/jeq2015.11.0567, 2016.

Lobaccaro, G. and Acero, J.: Comparative analysis of green actions to improve outdoor thermal comfort inside typical urban street canyons, Urban Climate, 14, 251–267, https://doi.org/10.1016/j.uclim.2015.10.002, 2015.

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., https://doi.org/10.1039/D0FD00092B, 2021.

Macdonald, R., Griffiths, R., and Hall, D.: An improved method for the estimation of surface roughness of obstacle arrays, At-

mos. Environ., 32, 1857–1864, https://doi.org/10.1016/S1352-2310(97)00403-2, 1998.

Maison, A. and Flageul, C.: Parameterizing the aerodynamic effect of trees in street canyons for the street-network model MUNICH using the CFD model Code_Saturne – Code_Saturne simulation dataset, V2, Mendeley Data [code], https://doi.org/10.17632/fzfrjsz3mv.2, 2021.

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.

Direction des Espaces Verts et de l'Environnement – Mairie de Paris: Les arbres – OpenDataParis, https://opendata.paris.fr/ explore/dataset/les-arbres/, last access: 17 December 2021.

Nowak, D. J., McHale, P. J., Ibarra, M., Crane, D., Stevens, J. C., and Luley, C. J.: Modeling the Effects of Urban Vegetation on Air Pollution, in: Air Pollution Modeling and Its Application XII, edited by: Gryning, S. E. and Chaumerliac, N., NATO, Challenges of Modern Society, vol 22., Springer, Boston, MA, 1998.

Nowak, D. J. and Crane, D. E.: Carbon storage and sequestration by urban trees in the USA, Environ. Pollut., 116, 381–389, https://doi.org/10.1016/S0269-7491(01)00214-7, 2002.

Nowak, D. J., Crane, D. E., and Stevens, J. C.: Air pollution removal by urban trees and shrubs in the United States, Urban For. Urban Gree., 4, 115–123, https://doi.org/10.1016/j.ufug.2006.01.007, 2006.

Oke, T.: Street design and urban canopy layer climate, Energ. Buildings, 11, 103–113, https://doi.org/10.1016/0378-7788(88)90026-6, 1988.

Oke, T. R.: The energetic basis of the urban heat island, Q. J. Roy. Meteor. Soc., 108, 1–24, https://doi.org/10.1002/qj.49710845502, 1982.

Ozdemir, H.: Mitigation impact of roadside trees on fine particle pollution, Sci. Total Environ., 659, 1176–1185, https://doi.org/10.1016/j.scitotenv.2018.12.262, 2019.

Pascal, M., Corso, M., Chanel, O., Declercq, C., Badaloni, C., Cesaroni, G., Henschel, S., Meister, K., Haluza, D., Martin-Olmedo, P., and Medina, S.: Assessing the public health impacts of urban air pollution in 25 European cities: Results of the Aphekom project, Sci. Total Environ., 449, 390–400, https://doi.org/10.1016/j.scitotenv.2013.01.077, 2013.

Pigeon, G., Legain, D., Durand, P., and Masson, V.: Anthropogenic heat release in an old European agglomeration (Toulouse, France), Int. J. Climatol., 27, 1969–1981, https://doi.org/10.1002/joc.1530, 2007.

Préndez, M., Araya, M., Criollo, C., Egas, C., Farías, I., Fuentealba, R., and González, E.: Urban Trees and Their Relationship with Air Pollution by Particulate Matter and Ozone in Santiago, Chile, in: Urban Climates in Latin America, edited by: Henríquez, C. and Romero, H., Springer International Publishing, pp. 167–206, https://doi.org/10.1007/978-3-319-97013-4_8, 2019.

Revelli, R. and Porporato, A.: Ecohydrological model for the quantification of ecosystem services provided by urban street trees, Urban Ecosyst., 21, 489–504, https://doi.org/10.1007/s11252-018-0741-2, 2018.

Robine, J., Cheung, S., and Roy, S. L.: Report on excess mortality in Europe during summer 2003, Tech. Rep., EU Community Action Programme for Public Health, http://ec.europa.eu/health/
phProjects/2005/action1/docs/action1_2005_a2_15_en.pdf
(last access: 18 July 2022), 2007.

Santiago, J.-L., Rivas, E., Sanchez, B., Buccolieri, R., and Martin, F.: The Impact of Planting Trees on NOₓ Concentrations: The Case of the Plaza de la Cruz Neighborhood in Pamplona (Spain), Atmosphere, 8, 131, https://doi.org/10.3390/atmos8070131, 2017.

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 concentrations over Greater Paris, Atmos. Environ., 180, 126–137, 2018.

Selmi, W., Weber, C., Rivière, E., Blond, N., Mehdi, L., and Nowak, D.: Air pollution removal by trees in public green spaces in Strasbourg city, France, Urban For. Urban Gree., 17, 192–201, https://doi.org/10.1016/j.ufug.2016.04.010, 2016.

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.

Speziale, C., Sarkar, S., and Gatski, T.: Modelling the pressure-strain correlation of turbulence – An invariant dynamical systems approach, J. Fluid Mech., 227, 245–272, https://doi.org/10.1017/S0022112091000101, 1991.

Stewart, I.: A systematic review and scientific critique of methodology in modern urban heat island literature, Int. J. Climatol., 31, 200–217, https://doi.org/10.1002/joc.2141, 2011.

Svirejeva-Hopkins, A., Schellnhuber, H., and Pomaz, V.: Urbanised territories as a specific component of the Global Carbon Cycle, Ecol. Model., 173, 295–312, https://doi.org/10.1016/j.ecolmodel.2003.09.022, 2004.

Taha, H., Akbari, H., and Rosenfeld, A.: Heat island and oasis effects of vegetative canopies: Micro-meteorological field-measurements, Theor. Appl. Climatol., 44, 123–138, https://doi.org/10.1007/BF00867999, 1991.

van Dillen, S. M. E., de Vries, S., Groenewegen, P. P., and Spreeuwenberg, P.: Greenspace in urban neighbourhoods and residents' health: adding quality to quantity, J. Epidemiol. Commun. H., 66, e8, https://doi.org/10.1136/jech.2009.104695, 2012.

Vardoulakis, S., Fisher, B. E., 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.

Vos, P. E., Maiheu, B., Vankerkom, J., and Janssen, S.: Improving local air quality in cities: To tree or not to tree?, Environ. Pollut., 183, 113–122, https://doi.org/10.1016/j.envpol.2012.10.021, 2013.

Wang, W.: An Analytical Model for Mean Wind Profiles in Sparse Canopies, Bound.-Lay. Meteorol., 142, 383–399, https://doi.org/10.1007/s10546-011-9687-0, 2012.

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.

Wania, A., Bruse, M., Blond, N., and Weber, C.: Analysing the influence of different street vegetation on traffic-induced particle dispersion using microscale simulations., J. Env. Manag., 94, 91–101, https://doi.org/10.1016/j.jenvman.2011.06.036, 2012.

Wei, X., Dupont, E., Gilbert, E., Musson-Genon, L., and Carissimo, B.: Experimental and Numerical Study of Wind and Turbulence in a Near-Field Dispersion Campaign at an Inhomogeneous Site, Bound.-Lay. Meteorol., 160, 475–499, https://doi.org/10.1007/s10546-016-0148-7, 2016.

West, J. J., Cohen, A., Dentener, F., Brunekreef, B., Zhu, T., Armstrong, B., Bell, M. L., Brauer, M., Carmichael, G., Costa, D. L., Dockery, D. W., Kleeman, M., Krzyzanowski, M., Künzli, N., Liousse, C., Lung, S.-C. C., Martin, R. V., Pöschl, U., Pope, C. A., Roberts, J. M., Russell, A. G., and Wiedinmyer, C.: What We Breathe Impacts Our Health: Improving Understanding of the Link between Air Pollution and Health, Environ. Sci. Technol., 50, 4895–4904, https://doi.org/10.1021/acs.est.5b03827, 2016.

Xue, F. and Li, X.: The impact of roadside trees on traffic released PM₁₀ in urban street canyon: Aerodynamic and deposition effects, Sustain. Cities Soc., 30, 195–204, https://doi.org/10.1016/j.scs.2017.02.001, 2017.

Yuan, C., Ng, E., and Norford, L. K.: Improving air quality in high-density cities by understanding the relationship between air pollutant dispersion and urban morphologies, Build. Environ., 71, 245–258, https://doi.org/10.1016/j.buildenv.2013.10.008, 2014.

Zaïdi, H., Dupont, E., Milliez, M., Musson-Genon, L., and Carissimo, B.: Numerical Simulations of the Microscale Heterogeneities of Turbulence Observed on a Complex Site, Bound.-Lay. Meteorol., 147, 237–259, https://doi.org/10.1007/s10546-012-9783-9, 2013.

Zhang, Y., Gu, Z., and Yu, C. W.: Impact Factors on Airflow and Pollutant Dispersion in Urban Street Canyons and Comprehensive Simulations: a Review, Current Pollution Report, 6, 425–439, https://doi.org/10.1007/s40726-020-00166-0, 2020.