# Significant impact of urban tree biogenic emissions on air quality estimated by a bottom-up inventory and chemistry transport modeling

Alice Maison, Lya Lugon, Soo-Jin Park, Alexia Baudic, Christopher Cantrell, Florian Couvidat, Barbara d'Anna, Claudia Di Biagio, Aline Gratien, Valérie Gros, et al.

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

## To cite this version:

Alice Maison, Lya Lugon, Soo-Jin Park, Alexia Baudic, Christopher Cantrell, et al. Significant impact of urban tree biogenic emissions on air quality estimated by a bottom-up inventory and chemistry transport modeling. Atmospheric Chemistry and Physics, 2024, 24 (10), pp.6011-6046. 10.5194/acp-24-6011-2024. hal-04589996

HAL Id: hal-04589996

https://hal.science/hal-04589996

Submitted on 28 May 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>

# Significant impact of urban tree biogenic emissions on air quality estimated by a bottom-up inventory and chemistry transport modeling

Alice Maison<fnref n="1" />,<fnref n="2" />, Lya Lugon<fnref n="1" />, Soo-Jin Park<fnref n="1" />, Alexia Baudic<fnref n="3" />, Christopher Cantrell<fnref n="4" />, Florian Couvidat<fnref n="5" />,
Barbara D'Anna<fnref n="6" />, Claudia Di Biagio<fnref n="7" />, Aline Gratien<fnref n="7" />, Valérie Gros<fnref n="8" />, Carmen Kalalian<fnref n="8" />,a,
Julien Kammer<fnref n="6" />, Vincent Michoud<fnref n="7" />, Jean-Eudes Petit<fnref n="8" />, Marwa Shahin<fnref n="6" />, Leila Simon<fnref n="8" />,b, Myrto Valari<fnref n="9" />,
Jérémy Vigneron<fnref n="3" />, Andrée Tuzet<fnref n="2" />, and Karine Sartelet<fnref n="1" />

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

Received: 22 November 2023 – Discussion started: 20 December 2023

Revised: 13 March 2024 – Accepted: 27 March 2024 – Published: 24 May 2024

**Abstract.** Biogenic volatile organic compounds (BVOCs) are emitted by vegetation and react with other compounds to form ozone and secondary organic matter (OM). In regional air quality models, biogenic emissions are often calculated using a plant functional type approach, which depends on the land use category. However, over cities, the land use is urban, so trees and their emissions are not represented. Here, we develop a bottom-up inventory of urban tree biogenic emissions in which the location of trees and their characteristics are derived from the tree database of the Paris city combined with allometric equations. Biogenic emissions are then computed for each tree based on their leaf dry biomass, tree-species-dependent emission factors, and activity factors representing the effects of light and temperature. Emissions are integrated in WRF-CHIMERE air quality simulations performed over June–July 2022. Over Paris city, the urban tree emissions have a significant impact on OM, inducing an average increase in the OM of about 5%, reaching 14% locally during the heatwaves. Ozone concentrations increase by 1.0% on average and by 2.4% during heatwaves, with a local increase of up to 6%. The concentration increase remains spatially localized over Paris, extending to the Paris suburbs in the case of ozone during heatwaves. The inclusion of urban tree emissions improves the estimation of OM concentrations compared to in situ measurements, but they are still underestimated as trees are still missing from the inventory. OM concentrations are sensitive to terpene emissions, highlighting the importance of favoring urban tree species with low-terpene emissions.

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

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

<footnote n="3">³Airparif, Association Agrée pour la Surveillance de la Qualité de l'Air en région Île-de-France,
7 rue Crillon, 75004 Paris, France</footnote>

<footnote n="4">⁴Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, 94010 Créteil, France</footnote>

<footnote n="5">⁵Institut National de l'Environnement Industriel et des Risques, 60550 Verneuil-en-Halatte, France</footnote>

<footnote n="6">⁶Aix Marseille Univ, CNRS, LCE, 13331 Marseille, France</footnote>

<footnote n="7">⁷Université Paris Cité and Univ Paris Est Creteil, CNRS, LISA, 75013 Paris, France</footnote>

<footnote n="8">⁸Laboratoire des Sciences du Climat et l'Environnement, CEA/Orme des Merisiers,
91191 Gif-sur-Yvette, France</footnote>

<footnote n="9">⁹Laboratoire de Météorologie Dynamique-IPSL, Sorbonne Université/CNRS/École Normale Supérieure-PSL
Université/École Polytechnique-Institut Polytechnique de Paris, 75005 Paris, France</footnote>

<footnote>a now at: Université Paris-Saclay, INRAE, AgroParisTech, UMR EcoSys, 91120 Palaiseau, France</footnote>

<footnote>b now at: Atmospheric Composition Research, Finnish Meteorological Institute, 00101 Helsinki, Finland</footnote>
<figure><img src="image_2.png" /></figure>

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

# Significant impact of urban tree biogenic emissions on air quality estimated by a bottom-up inventory and chemistry transport modeling

Alice Maison<fnref n="1" />,<fnref n="2" />, Lya Lugon<fnref n="1" />, Soo-Jin Park<fnref n="1" />, Alexia Baudic<fnref n="3" />, Christopher Cantrell<fnref n="4" />, Florian Couvidat<fnref n="5" />,
Barbara D'Anna<fnref n="6" />, Claudia Di Biagio<fnref n="7" />, Aline Gratien<fnref n="7" />, Valérie Gros<fnref n="8" />, Carmen Kalalian<fnref n="8" />,a,
Julien Kammer<fnref n="6" />, Vincent Michoud<fnref n="7" />, Jean-Eudes Petit<fnref n="8" />, Marwa Shahin<fnref n="6" />, Leila Simon<fnref n="8" />,b, Myrto Valari<fnref n="9" />,
Jérémy Vigneron<fnref n="3" />, Andrée Tuzet<fnref n="2" />, and Karine Sartelet<fnref n="1" />

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

Received: 22 November 2023 – Discussion started: 20 December 2023

Revised: 13 March 2024 – Accepted: 27 March 2024 – Published: 24 May 2024

**Abstract.** Biogenic volatile organic compounds (BVOCs) are emitted by vegetation and react with other compounds to form ozone and secondary organic matter (OM). In regional air quality models, biogenic emissions are often calculated using a plant functional type approach, which depends on the land use category. However, over cities, the land use is urban, so trees and their emissions are not represented. Here, we develop a bottom-up inventory of urban tree biogenic emissions in which the location of trees and their characteristics are derived from the tree database of the Paris city combined with allometric equations. Biogenic emissions are then computed for each tree based on their leaf dry biomass, tree-species-dependent emission factors, and activity factors representing the effects of light and temperature. Emissions are integrated in WRF-CHIMERE air quality simulations performed over June–July 2022. Over Paris city, the urban tree emissions have a significant impact on OM, inducing an average increase in the OM of about 5%, reaching 14% locally during the heatwaves. Ozone concentrations increase by 1.0% on average and by 2.4% during heatwaves, with a local increase of up to 6%. The concentration increase remains spatially localized over Paris, extending to the Paris suburbs in the case of ozone during heatwaves. The inclusion of urban tree emissions improves the estimation of OM concentrations compared to in situ measurements, but they are still underestimated as trees are still missing from the inventory. OM concentrations are sensitive to terpene emissions, highlighting the importance of favoring urban tree species with low-terpene emissions.

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

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

<footnote n="3">³Airparif, Association Agrée pour la Surveillance de la Qualité de l'Air en région Île-de-France,
7 rue Crillon, 75004 Paris, France</footnote>

<footnote n="4">⁴Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, 94010 Créteil, France</footnote>

<footnote n="5">⁵Institut National de l'Environnement Industriel et des Risques, 60550 Verneuil-en-Halatte, France</footnote>

<footnote n="6">⁶Aix Marseille Univ, CNRS, LCE, 13331 Marseille, France</footnote>

<footnote n="7">⁷Université Paris Cité and Univ Paris Est Creteil, CNRS, LISA, 75013 Paris, France</footnote>

<footnote n="8">⁸Laboratoire des Sciences du Climat et l'Environnement, CEA/Orme des Merisiers,
91191 Gif-sur-Yvette, France</footnote>

<footnote n="9">⁹Laboratoire de Météorologie Dynamique-IPSL, Sorbonne Université/CNRS/École Normale Supérieure-PSL
Université/École Polytechnique-Institut Polytechnique de Paris, 75005 Paris, France</footnote>

<footnote>a now at: Université Paris-Saclay, INRAE, AgroParisTech, UMR EcoSys, 91120 Palaiseau, France</footnote>

<footnote>b now at: Atmospheric Composition Research, Finnish Meteorological Institute, 00101 Helsinki, Finland</footnote>
# 1 Introduction

With an increasing number of people living in cities, urban areas are experiencing continuous expansion (Angel et al., 2011; United Nations, 2018). Artificial surfaces with darker, impermeable materials and high buildings, as well as release of anthropogenic heat, strongly modify the energy budget of the urban area (Taha et al., 1988; Taha, 1997; Pigeon et al., 2007a; Kuttler, 2008; Oke et al., 2017; Masson et al., 2020). An increase in temperature in the city compared to the surrounding countryside is often observed and is called the urban heat island (UHI) effect (Oke, 1982; Kim, 1992). Due to the high local emission sources (traffic) and the modification of airflows by tall buildings which limits the pollutant dispersion, concentrations of several pollutants, such as NO₂ and particles, are higher in cities than the surrounding areas (Lyons et al., 1990; Fenger, 1999; Thunis, 2018; Li et al., 2019; Yang et al., 2020).

To mitigate the negative effects of urbanization, urban vegetation and trees in particular are now widely promoted (Livesley et al., 2016; Chang et al., 2017; Roeland et al., 2019). Trees can locally reduce air and surface temperatures by creating shade and by evaporating water through transpiration (Jamei et al., 2016; Taleghani, 2018; Lai et al., 2019; Hami et al., 2019; Nasrollahi et al., 2020). Trees can also remove gaseous and particulate pollutants from the atmosphere by dry deposition, although this effect is quantitatively questionable due to the large variability and uncertainties (Nowak et al., 2006; Escobedo and Nowak, 2009; Setälä et al., 2013; Selmi et al., 2016; Nemitz et al., 2020; Lindén et al., 2023).

Trees are known to naturally emit biogenic volatile organic compounds (BVOCs). The term BVOC includes gaseous non-methane hydrocarbons and includes many families of molecules: isoprene, terpenes, alkanes, alkenes, carbonyls, alcohols, esters, ethers, and acids. BVOC emissions are involved in stress resistance mechanisms (due to heat, water shortage, oxidation, and herbivore or pathogen attack) and communication (plant–plant and plant–insect interactions) (Kesselmeier and Staudt, 1999). Emission rates depend on abiotic factors such as temperature and light, and biotic factors include things such as tree species, leaf age, and stress level (Niinemets et al., 2004; Loreto and Schnitzler, 2010; Niinemets, 2010). BVOC emissions are therefore highly variable in space and time. Unlike specific anthropogenic volatile organic compounds (AVOCs) such as benzene, emitted BVOCs may not be directly harmful to human health. However, BVOCs may form secondary pollutants, such as ozone (Calfapietra, 2013; Atkinson and Arey, 2003a; Churkina et al., 2017) and secondary organic aerosols (Salvador et al., 2020; Minguillón et al., 2016; Churkina et al., 2017; Lehtipalo et al., 2018). BVOCs emitted in the gaseous phase are oxidized in the atmosphere, forming more functionalized compounds that are semi-volatile and may be

absorbed into aerosols. In the urban VOC-limited environment with high nitrogen oxides (NOₓ) (emitted by traffic), ozone formation strongly depends on the VOC concentrations. There are also feedbacks between the urban environment, which is stressful for trees (higher temperatures and concentrations of oxidizing pollutants and difficult access to water) (Lüttge and Buckeridge, 2023; Czaja et al., 2020), and BVOC emissions.

To understand processes and forecast the evolution of pollutant concentrations, numerical models are widely used. Air quality models of various types and resolutions exist, depending on the scale and processes studied. Chemistry transport models (CTMs) are Eulerian models that represent the chemistry and transport of pollutants in three-dimensional grid cells, e.g., CHIMERE (Menut et al., 2021), Polair3D (Boutahar et al., 2004), WRF-Chem (NOAA/ESRL, 2023), CMAQ (Byun and Schere, 2006; Appel et al., 2021), and MOCAGE (CNRM, 2023). Their typical horizontal resolution varies between 1 and 10² km, and they are used from the global to the regional scales (Mailler et al., 2017). Many input data are necessary: surface characterization, spatiotemporal emissions of each pollutant, and boundary and initial conditions. Due to the coarse resolution, the surface type is characterized by land use categories such as open water, urban, forest, and crop. Forest trees are usually divided into one to five categories based on general characteristics (evergreen or deciduous; broadleaf or needleleaf). Most of the CTMs compute the BVOC emissions from forest and crops based on plant functional type (PFT) and the MEGAN empirical model (Model of Emissions of Gases and Aerosols from Nature) (Guenther et al., 2006, 2012; Matthias et al., 2018). The heterogeneity of the vegetation species is not explicitly modeled, but the model gives a rough estimate of the BVOC emission rates in the grid cells containing vegetation. CTMs can be used to compute air quality over large urban areas, but at this spatial resolution, the land use is urban, and no biogenic emission from urban vegetation is taken into account. In parallel, tree inventories are being developed in many cities (Bennett, 2023) and give us a much more accurate description of the urban forest. They can contain the precise locations, species, and sizes of trees, allowing the study of the impact of urban vegetation on air quality (Mircea et al., 2023).

Based on the tree inventory implemented in Paris by the municipality (Municipality of Paris, 2023) and a series of allometric equations (McPherson et al., 2016), a method is developed to estimate the BVOC emissions by urban trees in Paris. This “bottom-up” inventory of BVOC emissions by urban trees is used to estimate emissions from Paris trees over June and July 2022. This period is chosen because biogenic emissions are expected to be higher in summer, especially during heatwaves, and also because of the numerous in situ measurements that have been performed in the Paris region.
# 1 Introduction

With an increasing number of people living in cities, urban areas are experiencing continuous expansion (Angel et al., 2011; United Nations, 2018). Artificial surfaces with darker, impermeable materials and high buildings, as well as release of anthropogenic heat, strongly modify the energy budget of the urban area (Taha et al., 1988; Taha, 1997; Pigeon et al., 2007a; Kuttler, 2008; Oke et al., 2017; Masson et al., 2020). An increase in temperature in the city compared to the surrounding countryside is often observed and is called the urban heat island (UHI) effect (Oke, 1982; Kim, 1992). Due to the high local emission sources (traffic) and the modification of airflows by tall buildings which limits the pollutant dispersion, concentrations of several pollutants, such as NO₂ and particles, are higher in cities than the surrounding areas (Lyons et al., 1990; Fenger, 1999; Thunis, 2018; Li et al., 2019; Yang et al., 2020).

To mitigate the negative effects of urbanization, urban vegetation and trees in particular are now widely promoted (Livesley et al., 2016; Chang et al., 2017; Roeland et al., 2019). Trees can locally reduce air and surface temperatures by creating shade and by evaporating water through transpiration (Jamei et al., 2016; Taleghani, 2018; Lai et al., 2019; Hami et al., 2019; Nasrollahi et al., 2020). Trees can also remove gaseous and particulate pollutants from the atmosphere by dry deposition, although this effect is quantitatively questionable due to the large variability and uncertainties (Nowak et al., 2006; Escobedo and Nowak, 2009; Setälä et al., 2013; Selmi et al., 2016; Nemitz et al., 2020; Lindén et al., 2023).

Trees are known to naturally emit biogenic volatile organic compounds (BVOCs). The term BVOC includes gaseous non-methane hydrocarbons and includes many families of molecules: isoprene, terpenes, alkanes, alkenes, carbonyls, alcohols, esters, ethers, and acids. BVOC emissions are involved in stress resistance mechanisms (due to heat, water shortage, oxidation, and herbivore or pathogen attack) and communication (plant–plant and plant–insect interactions) (Kesselmeier and Staudt, 1999). Emission rates depend on abiotic factors such as temperature and light, and biotic factors include things such as tree species, leaf age, and stress level (Niinemets et al., 2004; Loreto and Schnitzler, 2010; Niinemets, 2010). BVOC emissions are therefore highly variable in space and time. Unlike specific anthropogenic volatile organic compounds (AVOCs) such as benzene, emitted BVOCs may not be directly harmful to human health. However, BVOCs may form secondary pollutants, such as ozone (Calfapietra, 2013; Atkinson and Arey, 2003a; Churkina et al., 2017) and secondary organic aerosols (Salvador et al., 2020; Minguillón et al., 2016; Churkina et al., 2017; Lehtipalo et al., 2018). BVOCs emitted in the gaseous phase are oxidized in the atmosphere, forming more functionalized compounds that are semi-volatile and may be

absorbed into aerosols. In the urban VOC-limited environment with high nitrogen oxides (NOₓ) (emitted by traffic), ozone formation strongly depends on the VOC concentrations. There are also feedbacks between the urban environment, which is stressful for trees (higher temperatures and concentrations of oxidizing pollutants and difficult access to water) (Lüttge and Buckeridge, 2023; Czaja et al., 2020), and BVOC emissions.

To understand processes and forecast the evolution of pollutant concentrations, numerical models are widely used. Air quality models of various types and resolutions exist, depending on the scale and processes studied. Chemistry transport models (CTMs) are Eulerian models that represent the chemistry and transport of pollutants in three-dimensional grid cells, e.g., CHIMERE (Menut et al., 2021), Polair3D (Boutahar et al., 2004), WRF-Chem (NOAA/ESRL, 2023), CMAQ (Byun and Schere, 2006; Appel et al., 2021), and MOCAGE (CNRM, 2023). Their typical horizontal resolution varies between 1 and 10² km, and they are used from the global to the regional scales (Mailler et al., 2017). Many input data are necessary: surface characterization, spatiotemporal emissions of each pollutant, and boundary and initial conditions. Due to the coarse resolution, the surface type is characterized by land use categories such as open water, urban, forest, and crop. Forest trees are usually divided into one to five categories based on general characteristics (evergreen or deciduous; broadleaf or needleleaf). Most of the CTMs compute the BVOC emissions from forest and crops based on plant functional type (PFT) and the MEGAN empirical model (Model of Emissions of Gases and Aerosols from Nature) (Guenther et al., 2006, 2012; Matthias et al., 2018). The heterogeneity of the vegetation species is not explicitly modeled, but the model gives a rough estimate of the BVOC emission rates in the grid cells containing vegetation. CTMs can be used to compute air quality over large urban areas, but at this spatial resolution, the land use is urban, and no biogenic emission from urban vegetation is taken into account. In parallel, tree inventories are being developed in many cities (Bennett, 2023) and give us a much more accurate description of the urban forest. They can contain the precise locations, species, and sizes of trees, allowing the study of the impact of urban vegetation on air quality (Mircea et al., 2023).

Based on the tree inventory implemented in Paris by the municipality (Municipality of Paris, 2023) and a series of allometric equations (McPherson et al., 2016), a method is developed to estimate the BVOC emissions by urban trees in Paris. This “bottom-up” inventory of BVOC emissions by urban trees is used to estimate emissions from Paris trees over June and July 2022. This period is chosen because biogenic emissions are expected to be higher in summer, especially during heatwaves, and also because of the numerous in situ measurements that have been performed in the Paris region.
The effect on isoprene (C5H8), monoterpane, ozone (O3), organic matter (OM), and particulate matter (PM) concentrations over the Île-de-France (IDF) region is quantified using the CTM CHIMERE.

# 2 Materials and methods

## 2.1 Tree-based BVOC emission model

To compute BVOC emissions in CTMs, an empirical approach is usually used. The emission rate of each BVOC species is computed as the product of several factors: the amount of vegetation (surface of the land use category and leaf area index or mass), an emission factor at the standard conditions (leaf temperature of 30 °C and photosynthetic photon flux density (PPFD) of 1000 µmol photons m⁻² s⁻¹), and activity factors representing physiological or meteorological effects. One emission factor is associated with each PFT. The development of BVOC rate measurements at the leaf or branch scale with chambers and tree inventories allows the estimation of BVOC emission rates at the tree level (Owen et al., 2001, 2003; Stewart et al., 2003; Karl et al., 2009; Steinbrecher et al., 2009). The emission rate of a BVOC k for a tree t (µg h⁻¹ per tree) can be estimated as

$$ER_{k,t} = DB_t \cdot EF_{k,t} \cdot \gamma_k, \quad (1)$$

where

- $EF_{k,t}$ (µg g_DW⁻¹ h⁻¹) is the emission factor (or potential) at standard conditions;

- $DB_t$ is the dry leaf biomass ($g_{DW}$, where DW stands for dry weight); and

- $\gamma_k$ combines the different dimensionless emission activity factors.

## 2.2 Tree inventory and characteristics

The Paris tree database (https://opendata.paris/fr/explore/dataset/les-arbres/map/, last access: 3 March 2023) regroups an inventory of the public trees. Much mapped information is available for each tree, including the precise location (coordinates), address, type (roadside, garden, cemetery, etc.), tree species, height, trunk circumference, and development stage. It is regularly updated, and the version used in this study was downloaded in March 2023. A map of trees around Avenue des Champs-Élysées taken from the database is shown in Fig. B1. The proportion (P) of the tree genera found in Paris is presented in Fig. B2, and the distributions of their trunk circumferences and crown heights are shown in Fig. B3. The municipality of Paris estimates that around one-third of the Parisian trees are missing from their database, mainly trees located in private areas. Without further information on these private trees, they are not taken into account in this study.

To compute the BVOC emissions (Eq. 1) of each individual tree, an estimation of the leaf dry biomass is necessary. Dry biomass such as leaf area and crown dimensions can be estimated using allometric equations. These allometric relationships are statistical models based on a sample of measurements predicting tree size as a function of parameters such as trunk diameter or age since planting. Many studies propose equations for forest trees (Burton et al., 1991; Bartelink, 1997; Karlik and McKay, 2002), but studies on urban trees are more scarce (Nowak, 1996). The open database of McPherson et al. (2016) is chosen in this study because it was developed specifically for urban trees and includes many genera found in Paris (84% of the trees in the Paris inventory) (365 growth equations for 174 tree species). For missing tree genera, an equation from another tree genus in the same family is selected, as described in Sect. 3. It includes allometric tree measurements for different climates in the United States, so assumptions are necessary to select the climate for each tree species that is the closest to that of the Paris region (see Sect. 3).

## 2.3 Description of regional-scale air quality simulations

To quantify the impact of the Parisian tree emissions on air quality, regional-scale simulations using the 3D CTM CHIMERE v2020_r3 (Menut et al., 2021) coupled to the chemical module SSH-aerosol v1.3 (Sartelet et al., 2020) are performed. The gas-phase chemical scheme is MELCHIOR2, modified to represent secondary organic aerosol formation, as described in Sartelet et al. (2020). Biogenic emissions are estimated using the MEGANv2.1 algorithm implemented in CHIMERE (Couvidat et al., 2018), which corresponds to a land use approach. The following section describes the simulation setup.

Simulations are performed during summer 2022, between 6 June and 31 July 2022, with a 5 d spin-up period (1–5 June). Summertime is chosen as biogenic emissions are the highest during this period due to meteorological conditions. In France, summer 2022 was exceptionally hot and sunny, with little precipitation (on average 1 to 3 °C above seasonal values over most of France) (Meteo France, 2023). The domain of study corresponds to the Île-de-France region, with a 1 km × 1 km spatial resolution (IDF1). Initial and boundary conditions are taken from two additional nesting simulations: one over France with a 9 km × 9 km spatial resolution (FRA9), and one over the northwest of France with a 3 km × 3 km spatial resolution (IDF3), as shown in Fig. 1. For the FRA9 domain, boundary and initial conditions are obtained from the CAMS platform (Inness et al., 2019), with a 10 km × 10 km spatial resolution.

Meteorological data for all domains are computed using the Weather Research and Forecasting (WRF) model v3.7.1 available in CHIMERE (Powers et al., 2017; Menut et al., 2021). Even if CHIMERE and WRF simulations are performed simultaneously, one-way coupling is used, and then,
The effect on isoprene (C₅H₈), monoterpane, ozone (O₃), organic matter (OM), and particulate matter (PM) concentrations over the Île-de-France (IDF) region is quantified using the CTM CHIMERE.

# 2 Materials and methods

## 2.1 Tree-based BVOC emission model

To compute BVOC emissions in CTMs, an empirical approach is usually used. The emission rate of each BVOC species is computed as the product of several factors: the amount of vegetation (surface of the land use category and leaf area index or mass), an emission factor at the standard conditions (leaf temperature of 30 °C and photosynthetic photon flux density (PPFD) of 1000 µmol photons m⁻² s⁻¹), and activity factors representing physiological or meteorological effects. One emission factor is associated with each PFT. The development of BVOC rate measurements at the leaf or branch scale with chambers and tree inventories allows the estimation of BVOC emission rates at the tree level (Owen et al., 2001, 2003; Stewart et al., 2003; Karl et al., 2009; Steinbrecher et al., 2009). The emission rate of a BVOC k for a tree t (µg h⁻¹ per tree) can be estimated as

$$ER_{k,t} = DB_t \cdot EF_{k,t} \cdot \gamma_k, \quad (1)$$

where

- $EF_{k,t}$ (µg g_DW⁻¹ h⁻¹) is the emission factor (or potential) at standard conditions;

- $DB_t$ is the dry leaf biomass ($g_{DW}$, where DW stands for dry weight); and

- $\gamma_k$ combines the different dimensionless emission activity factors.

## 2.2 Tree inventory and characteristics

The Paris tree database (https://opendata.paris/fr/explore/dataset/les-arbres/map/, last access: 3 March 2023) regroups an inventory of the public trees. Much mapped information is available for each tree, including the precise location (coordinates), address, type (roadside, garden, cemetery, etc.), tree species, height, trunk circumference, and development stage. It is regularly updated, and the version used in this study was downloaded in March 2023. A map of trees around Avenue des Champs-Élysées taken from the database is shown in Fig. B1. The proportion (P) of the tree genera found in Paris is presented in Fig. B2, and the distributions of their trunk circumferences and crown heights are shown in Fig. B3. The municipality of Paris estimates that around one-third of the Parisian trees are missing from their database, mainly trees located in private areas. Without further information on these private trees, they are not taken into account in this study.

To compute the BVOC emissions (Eq. 1) of each individual tree, an estimation of the leaf dry biomass is necessary. Dry biomass such as leaf area and crown dimensions can be estimated using allometric equations. These allometric relationships are statistical models based on a sample of measurements predicting tree size as a function of parameters such as trunk diameter or age since planting. Many studies propose equations for forest trees (Burton et al., 1991; Bartelink, 1997; Karlik and McKay, 2002), but studies on urban trees are more scarce (Nowak, 1996). The open database of McPherson et al. (2016) is chosen in this study because it was developed specifically for urban trees and includes many genera found in Paris (84% of the trees in the Paris inventory) (365 growth equations for 174 tree species). For missing tree genera, an equation from another tree genus in the same family is selected, as described in Sect. 3. It includes allometric tree measurements for different climates in the United States, so assumptions are necessary to select the climate for each tree species that is the closest to that of the Paris region (see Sect. 3).

## 2.3 Description of regional-scale air quality simulations

To quantify the impact of the Parisian tree emissions on air quality, regional-scale simulations using the 3D CTM CHIMERE v2020_r3 (Menut et al., 2021) coupled to the chemical module SSH-aerosol v1.3 (Sartelet et al., 2020) are performed. The gas-phase chemical scheme is MELCHIOR2, modified to represent secondary organic aerosol formation, as described in Sartelet et al. (2020). Biogenic emissions are estimated using the MEGANv2.1 algorithm implemented in CHIMERE (Couvidat et al., 2018), which corresponds to a land use approach. The following section describes the simulation setup.

Simulations are performed during summer 2022, between 6 June and 31 July 2022, with a 5 d spin-up period (1–5 June). Summertime is chosen as biogenic emissions are the highest during this period due to meteorological conditions. In France, summer 2022 was exceptionally hot and sunny, with little precipitation (on average 1 to 3 °C above seasonal values over most of France) (Meteo France, 2023). The domain of study corresponds to the Île-de-France region, with a 1 km × 1 km spatial resolution (IDF1). Initial and boundary conditions are taken from two additional nesting simulations: one over France with a 9 km × 9 km spatial resolution (FRA9), and one over the northwest of France with a 3 km × 3 km spatial resolution (IDF3), as shown in Fig. 1. For the FRA9 domain, boundary and initial conditions are obtained from the CAMS platform (Inness et al., 2019), with a 10 km × 10 km spatial resolution.

Meteorological data for all domains are computed using the Weather Research and Forecasting (WRF) model v3.7.1 available in CHIMERE (Powers et al., 2017; Menut et al., 2021). Even if CHIMERE and WRF simulations are performed simultaneously, one-way coupling is used, and then,
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Representation of the simulated domains. The blue rectangles represent the location of the different nested domains.</figcaption></figure>

the concentrations computed in CHIMERE are assumed to have no influence on the meteorological fields computed by WRF. WRF simulations are performed with 33 vertical levels from 0 to 20 km altitude. A more refined vertical discretization is employed in the first four vertical levels (average heights of 0.12, 25, 50, and 83 m, respectively), which contains almost all buildings in Paris region.

The spatial distribution of each land use category used in WRF simulations is based on CORINE Land Cover (European Environment Agency, 2020). It was chosen for its very fine resolution (250 m) since, for 1 km resolution simulations over a city, a detailed description of the land use is required to correctly describe the urban fabric. To ensure the compatibility with the Noah Land Surface Model (LSM), we converted the classification from CORINE Land Cover into MODIS (Moderate Resolution Imaging Spectroradiometer) categories, following Vogel and Afshari (2020). Three urban categories are employed to differentiate street and building dimensions, as well as heat transfer parameters in commercial, high-, and low-intensity residential areas.

In order to represent more precisely the meteorological fields in urban areas, the single-layer Urban Canopy Model (UCM) (Kusaka et al., 2001) is used in the IDF3 and IDF1 domains. The UCM was chosen in WRF because it allows the input of anthropogenic heat (AH) fluxes for different urban categories; AH is assumed to be 45 W m⁻² for commercial areas, 10 W m⁻² for high-intensity residential areas, and 5 W m⁻² for low-intensity residential areas, based on Pigeon et al. (2007b) and Sailor et al. (2015). AH is crucial to correctly model the heat island effect, as well as the friction velocities above buildings. Table 1 summarizes the other physical options employed in the WRF simulations.

As expected, the WRF model simulates higher temperatures in urban areas than in rural areas (fields or forests), as shown in Fig. C1, which presents the 2-month average air temperature at 2 m simulated by WRF. The simulated downwards shortwave (SW) radiation at ground surface is also used to compute tree biogenic emissions. It is quite homogeneous over the Paris region, with a 2-month average of 500 W m⁻² (daytime) and spatial variations within 5 % of the mean.

Anthropogenic emissions in the domains FRA9 and IDF3 are from the latest (2020) EMEP emission inventory (EMEP, 2019) (0.1° × 0.1° horizontal resolution), and in IDF1, they are from the latest (2019) regional emission inventory of the Air Quality Monitoring Network (AQMN) Airparif for the greater Paris area (https://www.airparif.asso.fr/, last access: 5 May 2024) (1 km × 1 km spatial resolution). Traffic emissions correspond to those of the summer 2022, calculated using the bottom-up traffic emissions model HEAVEN (https://www.airparif.asso.fr/heaven-emissions-du-trafic-en-temps-reel, last access: 5 May 2024), while non-traffic anthropogenic emissions correspond to the 2019 Airparif inventory.

## 2.4 Description of the air quality experimental measurements

The results of the simulations are compared to experimental measurements performed at different sites in the Paris region. In Sect. 4.1.2 and 4.2, temporal variations in the observed and simulated concentrations are presented in three main sites: the Halles site, a permanent air quality monitoring station located in the city center and operated by Airparif; the PRG (Paris Rive Gauche) site, located on the seventh floor of the Lamarck B building of Université Paris Cité (30 m above ground level (a.g.l.)) in the southeastern side of the city and set up as part of the ACROSS campaign (Cantrell and Michoud, 2022); and the site of SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), an atmospheric observatory located 20 km southwest of Paris which is integrated in the ACTRIS European Research Infrastructure Consortium (https://www.actris.eu, last access: 5 May 2024) (Haeffelin et al., 2005). The Halles and PRG stations are both urban background sites, while SIRTA is a suburban background site. The three sites and the measurements performed are described in Table 2, and a more complete description of the measurements and their associated uncertainties is provided in Appendix A.

The reference simulations (without urban trees biogenic emissions) are validated in Sect. 4.1 with the observation sites of the Airparif network. These sites correspond to 21 permanent air quality monitoring stations included within a
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Representation of the simulated domains. The blue rectangles represent the location of the different nested domains.</figcaption></figure>

the concentrations computed in CHIMERE are assumed to have no influence on the meteorological fields computed by WRF. WRF simulations are performed with 33 vertical levels from 0 to 20 km altitude. A more refined vertical discretization is employed in the first four vertical levels (average heights of 0.12, 25, 50, and 83 m, respectively), which contains almost all buildings in Paris region.

The spatial distribution of each land use category used in WRF simulations is based on CORINE Land Cover (European Environment Agency, 2020). It was chosen for its very fine resolution (250 m) since, for 1 km resolution simulations over a city, a detailed description of the land use is required to correctly describe the urban fabric. To ensure the compatibility with the Noah Land Surface Model (LSM), we converted the classification from CORINE Land Cover into MODIS (Moderate Resolution Imaging Spectroradiometer) categories, following Vogel and Afshari (2020). Three urban categories are employed to differentiate street and building dimensions, as well as heat transfer parameters in commercial, high-, and low-intensity residential areas.

In order to represent more precisely the meteorological fields in urban areas, the single-layer Urban Canopy Model (UCM) (Kusaka et al., 2001) is used in the IDF3 and IDF1 domains. The UCM was chosen in WRF because it allows the input of anthropogenic heat (AH) fluxes for different urban categories; AH is assumed to be 45 W m⁻² for commercial areas, 10 W m⁻² for high-intensity residential areas, and 5 W m⁻² for low-intensity residential areas, based on Pigeon et al. (2007b) and Sailor et al. (2015). AH is crucial to correctly model the heat island effect, as well as the friction velocities above buildings. Table 1 summarizes the other physical options employed in the WRF simulations.

As expected, the WRF model simulates higher temperatures in urban areas than in rural areas (fields or forests), as shown in Fig. C1, which presents the 2-month average air temperature at 2 m simulated by WRF. The simulated downwards shortwave (SW) radiation at ground surface is also used to compute tree biogenic emissions. It is quite homogeneous over the Paris region, with a 2-month average of 500 W m⁻² (daytime) and spatial variations within 5 % of the mean.

Anthropogenic emissions in the domains FRA9 and IDF3 are from the latest (2020) EMEP emission inventory (EMEP, 2019) (0.1° × 0.1° horizontal resolution), and in IDF1, they are from the latest (2019) regional emission inventory of the Air Quality Monitoring Network (AQMN) Airparif for the greater Paris area (https://www.airparif.asso.fr/, last access: 5 May 2024) (1 km × 1 km spatial resolution). Traffic emissions correspond to those of the summer 2022, calculated using the bottom-up traffic emissions model HEAVEN (https://www.airparif.asso.fr/heaven-emissions-du-trafic-en-temps-reel, last access: 5 May 2024), while non-traffic anthropogenic emissions correspond to the 2019 Airparif inventory.

## 2.4 Description of the air quality experimental measurements

The results of the simulations are compared to experimental measurements performed at different sites in the Paris region. In Sect. 4.1.2 and 4.2, temporal variations in the observed and simulated concentrations are presented in three main sites: the Halles site, a permanent air quality monitoring station located in the city center and operated by Airparif; the PRG (Paris Rive Gauche) site, located on the seventh floor of the Lamarck B building of Université Paris Cité (30 m above ground level (a.g.l.)) in the southeastern side of the city and set up as part of the ACROSS campaign (Cantrell and Michoud, 2022); and the site of SIRTA (Site Instrumental de Recherche par Télédétection Atmosphérique), an atmospheric observatory located 20 km southwest of Paris which is integrated in the ACTRIS European Research Infrastructure Consortium (https://www.actris.eu, last access: 5 May 2024) (Haeffelin et al., 2005). The Halles and PRG stations are both urban background sites, while SIRTA is a suburban background site. The three sites and the measurements performed are described in Table 2, and a more complete description of the measurements and their associated uncertainties is provided in Appendix A.

The reference simulations (without urban trees biogenic emissions) are validated in Sect. 4.1 with the observation sites of the Airparif network. These sites correspond to 21 permanent air quality monitoring stations included within a
<table><caption>**Table 1.** Main physical options employed in WRF simulations.</caption><thead><tr><td>Option in WRF namelist</td><td>Option complete name</td><td>Option selected</td></tr></thead><tbody><tr><td>mpphysics</td><td>Microphysics</td><td>Thompson graupel scheme</td></tr><tr><td>cuphysics</td><td>Cumulus</td><td>Grell–Devenyi ensemble scheme</td></tr><tr><td>ra_lwphysics</td><td>Longwave radiation</td><td>rrtmg scheme</td></tr><tr><td>ra_swphysics</td><td>Shortwave radiation</td><td>rrtmg scheme</td></tr><tr><td>bl_pblphysics</td><td>Boundary layer</td><td>YSU scheme</td></tr><tr><td>sfSFclayPhysics</td><td>Surface layer</td><td>Monin–Obukhov similarity scheme</td></tr><tr><td>sf물 surfacePhysics</td><td>Land surface</td><td>Noah Land Surface Model</td></tr><tr><td>sf urbanPhysics</td><td>Urban Canopy Model</td><td>Single-layer (only in IDF3 and IDF1)</td></tr></tbody></table>

<table><caption>**Table 2.** Description of the experimental measurements performed at different sites and used in this study. ACSM is the aerosol chemical speciation monitor (Petit et al., 2015), PTR-MS is the proton transfer reaction mass spectrometry (Simon et al., 2023), and GC-FID is the gas chromatograph with a flame ionization detector (Gros et al., 2011).</caption><thead><tr><td>Site</td><td>Location</td><td>Typology</td><td>Species measured</td><td>Instrument</td></tr></thead><tbody><tr><td rowspan="4">Halles</td><td rowspan="4">1st district of Paris city (48.862128° N, 2.344622° E)</td><td rowspan="4">Urban background</td><td>NO<sub>2</sub></td><td>AC32M</td></tr><tr><td>O<sub>3</sub></td><td>O3 42e</td></tr><tr><td>PM</td><td>FIDAS 200</td></tr><tr><td>OM</td><td>ACSM</td></tr><tr><td rowspan="3">PRG</td><td rowspan="3">13th district of Paris city (48.827778° N, 2.380562° E)</td><td rowspan="3">Urban background</td><td>C<sub>5</sub>H<sub>8</sub></td><td>PTR-MS</td></tr><tr><td>Monoterpenes</td><td>PTR-MS</td></tr><tr><td>OM</td><td>ACSM</td></tr><tr><td rowspan="3">SIRTA</td><td rowspan="3">20 km southwest of Paris (48.709890° N, 2.147938° E)</td><td rowspan="3">Suburban</td><td>C<sub>5</sub>H<sub>8</sub></td><td>GC-FID</td></tr><tr><td>Monoterpenes</td><td>PTR-MS</td></tr><tr><td>OM</td><td>ACSM</td></tr></tbody></table>

large operational stations network operated by Airparif (see Table A1 and https://www.airparif.asso.fr/carte-des-stations, last access: 5 May 2024). The map in Fig. 2 shows the location of all measurement stations which are used to evaluate the simulations (Sect. 4.1 and 4.2). It also presents the land use from GLOBCOVER (Bontemps et al., 2011) used in the CHIMERE simulation over IDF1. It is mainly composed of agricultural lands, forests of varying sizes, and a large urban area including Paris and its suburbs.

# 3 Bottom-up inventory of tree BVOC emissions and comparison to the land use approach

## 3.1 Calculation of BVOC emissions at the tree level

### 3.1.1 Estimation of the tree dry biomass

The total tree leaf dry biomass (in grams of dry weight, g<sub>DW</sub>) is computed based on the McPherson et al. (2016) allometric equation database. Tree data were collected in 17 reference cities representative of the different US climate zones and analyzed to obtain growth equations. The database contains equations to estimate the tree characteristics from the tree species, climate, and the trunk diameter at breast height (DBH) (at 1.3 m). To find the correspondence between Parisian trees and this database, the US climates were first ranked from closest to farthest from the Parisian climate based on a qualitative comparison of annual rainfall and temperatures (see Table D1). For each tree species in the Paris tree inventory, the allometric equations are obtained from the US database by selecting the closest tree category in terms of tree species and climate, following the decision tree shown in Fig. D1. The default species is the plane tree (Planatus × hispanica), which is the predominant species in Paris.

Then, the trunk diameter at breast height, DBH (cm), is computed from the trunk circumference, CIRC (cm), available in the Paris tree inventory for each tree, assuming a cylindrical tree trunk, as DBH = CIRC / π. The total tree leaf area (LA in m²) is then computed from each Parisian tree using the selected equation form and coefficient and the computed DBH. For example, the function LA = f(DBH) is shown for three tree species which have different allometric equation forms in Eqs. (D1), (D2), and (D3). Finally, the dry biomass (DB in g<sub>DW</sub>) is the product of the leaf area and the dry weight density (DWD in g<sub>DW</sub> m⁻²), calculated as DB<sub>t</sub> = LA<sub>t</sub> × DWD<sub>t</sub>. The dry weight density depends on tree species, and it is also given in the McPherson et al. (2016) database. For instance, DWD = 500, 520, and 560 g<sub>DW</sub> m⁻²
<table><caption>**Table 1.** Main physical options employed in WRF simulations.</caption><thead><tr><th>Option in WRF namelist</th><th>Option complete name</th><th>Option selected</th></tr></thead><tbody><tr><td>mpphysics</td><td>Microphysics</td><td>Thompson graupel scheme</td></tr><tr><td>cuphysics</td><td>Cumulus</td><td>Grell–Devenyi ensemble scheme</td></tr><tr><td>ra_lwphysics</td><td>Longwave radiation</td><td>rrtmg scheme</td></tr><tr><td>ra_swphysics</td><td>Shortwave radiation</td><td>rrtmg scheme</td></tr><tr><td>bl_pblphysics</td><td>Boundary layer</td><td>YSU scheme</td></tr><tr><td>sfSFclayPhysics</td><td>Surface layer</td><td>Monin–Obukhov similarity scheme</td></tr><tr><td>sf물 surfacePhysics</td><td>Land surface</td><td>Noah Land Surface Model</td></tr><tr><td>sf urbanPhysics</td><td>Urban Canopy Model</td><td>Single-layer (only in IDF3 and IDF1)</td></tr></tbody></table>

<table><caption>**Table 2.** Description of the experimental measurements performed at different sites and used in this study. ACSM is the aerosol chemical speciation monitor (Petit et al., 2015), PTR-MS is the proton transfer reaction mass spectrometry (Simon et al., 2023), and GC-FID is the gas chromatograph with a flame ionization detector (Gros et al., 2011).</caption><thead><tr><th>Site</th><th>Location</th><th>Typology</th><th>Species measured</th><th>Instrument</th></tr></thead><tbody><tr><td rowspan="4">Halles</td><td rowspan="4">1st district of Paris city<br>(48.862128° N, 2.344622° E)</td><td rowspan="4">Urban background</td><td>NO<sub>2</sub></td><td>AC32M</td></tr><tr><td>O<sub>3</sub></td><td>O3 42e</td></tr><tr><td>PM</td><td>FIDAS 200</td></tr><tr><td>OM</td><td>ACSM</td></tr><tr><td rowspan="3">PRG</td><td rowspan="3">13th district of Paris city<br>(48.827778° N, 2.380562° E)</td><td rowspan="3">Urban background</td><td>C<sub>5</sub>H<sub>8</sub></td><td>PTR-MS</td></tr><tr><td>Monoterpenes</td><td>PTR-MS</td></tr><tr><td>OM</td><td>ACSM</td></tr><tr><td rowspan="3">SIRTA</td><td rowspan="3">20 km southwest of Paris<br>(48.709890° N, 2.147938° E)</td><td rowspan="3">Suburban</td><td>C<sub>5</sub>H<sub>8</sub></td><td>GC-FID</td></tr><tr><td>Monoterpenes</td><td>PTR-MS</td></tr><tr><td>OM</td><td>ACSM</td></tr></tbody></table>

large operational stations network operated by Airparif (see Table A1 and https://www.airparif.asso.fr/carte-des-stations, last access: 5 May 2024). The map in Fig. 2 shows the location of all measurement stations which are used to evaluate the simulations (Sect. 4.1 and 4.2). It also presents the land use from GLOBCOVER (Bontemps et al., 2011) used in the CHIMERE simulation over IDF1. It is mainly composed of agricultural lands, forests of varying sizes, and a large urban area including Paris and its suburbs.

# 3 Bottom-up inventory of tree BVOC emissions and comparison to the land use approach

## 3.1 Calculation of BVOC emissions at the tree level

### 3.1.1 Estimation of the tree dry biomass

The total tree leaf dry biomass (in grams of dry weight, g<sub>DW</sub>) is computed based on the McPherson et al. (2016) allometric equation database. Tree data were collected in 17 reference cities representative of the different US climate zones and analyzed to obtain growth equations. The database contains equations to estimate the tree characteristics from the tree species, climate, and the trunk diameter at breast height (DBH) (at 1.3 m). To find the correspondence between Parisian trees and this database, the US climates were first ranked from closest to farthest from the Parisian climate based on a qualitative comparison of annual rainfall and temperatures (see Table D1). For each tree species in the Paris tree inventory, the allometric equations are obtained from the US database by selecting the closest tree category in terms of tree species and climate, following the decision tree shown in Fig. D1. The default species is the plane tree (Planatus × hispanica), which is the predominant species in Paris.

Then, the trunk diameter at breast height, DBH (cm), is computed from the trunk circumference, CIRC (cm), available in the Paris tree inventory for each tree, assuming a cylindrical tree trunk, as DBH = CIRC / π. The total tree leaf area (LA in m²) is then computed from each Parisian tree using the selected equation form and coefficient and the computed DBH. For example, the function LA = f(DBH) is shown for three tree species which have different allometric equation forms in Eqs. (D1), (D2), and (D3). Finally, the dry biomass (DB in g<sub>DW</sub>) is the product of the leaf area and the dry weight density (DWD in g<sub>DW</sub> m⁻²), calculated as DB<sub>t</sub> = LA<sub>t</sub> × DWD<sub>t</sub>. The dry weight density depends on tree species, and it is also given in the McPherson et al. (2016) database. For instance, DWD = 500, 520, and 560 g<sub>DW</sub> m⁻²
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Map of the GLOBCOVER major land use in each grid cell used in IDF1 CHIMERE simulations. The crosses represent the locations of the different measurement stations.</figcaption></figure>

for *Planatus × hispanica*, *Acer platanoides*, and *Prunus serrulata*, respectively. The computed LA and DB are shown in Fig. 3 for the predominant tree species ($P > 1\%$) as a function of the DBH. It shows that LA and DB increase with DBH, but there is a large variability depending on the tree species. For example, for a tree of DBH = 100 cm, the estimation of the leaf surface is equal to 1151.5 m² for *Planatus × hispanica*, 582.5 m² for *Acer platanoides*, and 1147.7 m² for *Prunus serrulata*.

As the simulation is performed during the late spring and summer periods, the tree foliage is assumed to be fully developed, so that the leaf area and dry biomass are constant over time. For longer simulation periods, the temporal evolution of leaf area and dry biomass should be introduced.

### 3.1.2 Emission factors by tree species and computation of activity factors

The emission factors by tree species are taken from MEGANv3.2 code, downloaded at https://bai.ess.uci.edu/megan/data-and-code/megan32 (last access: 10 July 2023). The EF presented by tree species is assumed to be identical within the same tree genus. Therefore, EF by tree genus is used for all trees except for the *Quercus* genus (oak), whose species are known to have very different BVOC emission profiles (Loreto, 2002). The EF for the *Quercus* species is taken from Ciccioli et al. (2023) for isoprene and monoterpenes. For the *Quercus* species missing in Ciccioli et al. (2023) but with a known emission type (Loreto, 2002), the EF values are taken from MEGANv3.2. For unknown emission types, the EF value is set by default to oak isoprene emitters in MEGANv3.2. For all tree species, the EF values for sesquiterpenes and oxygenated VOC are taken from MEGANv3.2. The emission factors of nitric oxide (NO)

and carbon monoxide (CO) are fixed for all tree species and are equal to $EF_{NO} = 0.05$ and $EF_{CO} = 1.0 \mu g \text{ g}_{\text{DW}}^{-1} \text{ h}^{-1}$, as suggested in MEGANv3.2. The emission factors of isoprene (ISOP), total monoterpenes (MTS), total sesquiterpenes (SQTs), and total other VOCs (OVOCs) are shown in Table D2 for the predominant tree genera and oak species.

Emission factors, which are measured at standard conditions, are then multiplied by dimensionless factors representing variations in the emissions as a function of biotic and abiotic processes. A detailed description of the calculation of activity factors for light (photosynthetic photon flux density, PPFD), $\gamma_P$, and for temperature, $\gamma_T$, is given in Appendix D3. To illustrate the variation in the activity factors with light and temperature, Fig. 4 shows averaged $\gamma_T$ and $\gamma_P$ for isoprene, $\alpha$-pinene, and $\beta$-pinene. Figure 4a shows that BVOC emissions increase with temperature. At high temperatures, isoprene emissions are capped, while monoterpene emissions rise sharply. BVOC emissions also increase with light (Fig. 4b), and the activity factor reaches its maximum value of 1 after PPFD = 1000 $\mu mol$ photons m$^{-2}$ s$^{-1}$.

Other activity factors could be added to represent the effect of leaf age and water stress. In this study, emissions are calculated per measure of leaf biomass, considering an average emission for all leaves in the canopy. In addition, we assume that in June and July, tree foliage is fully developed, and leaf area and dry biomass are constant. Therefore, no activity factor is added to modulate emissions according to the fraction of new, growing, mature, and old foliage (Guenther et al., 2012).

Several studies also introduce an activity factor to represent the impact of soil moisture and water stress on isoprene emissions (Guenther et al., 2012; Jiang et al., 2018; Bonn et al., 2019; Otu-Larbi et al., 2020; Wang et al., 2022). Al-
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Map of the GLOBCOVER major land use in each grid cell used in IDF1 CHIMERE simulations. The crosses represent the locations of the different measurement stations.</figcaption></figure>

for *Planatus × hispanica*, *Acer platanoides*, and *Prunus serrulata*, respectively. The computed LA and DB are shown in Fig. 3 for the predominant tree species ($P > 1\%$) as a function of the DBH. It shows that LA and DB increase with DBH, but there is a large variability depending on the tree species. For example, for a tree of DBH = 100 cm, the estimation of the leaf surface is equal to 1151.5 m² for *Planatus × hispanica*, 582.5 m² for *Acer platanoides*, and 1147.7 m² for *Prunus serrulata*.

As the simulation is performed during the late spring and summer periods, the tree foliage is assumed to be fully developed, so that the leaf area and dry biomass are constant over time. For longer simulation periods, the temporal evolution of leaf area and dry biomass should be introduced.

### 3.1.2 Emission factors by tree species and computation of activity factors

The emission factors by tree species are taken from MEGANv3.2 code, downloaded at https://bai.ess.uci.edu/megan/data-and-code/megan32 (last access: 10 July 2023). The EF presented by tree species is assumed to be identical within the same tree genus. Therefore, EF by tree genus is used for all trees except for the *Quercus* genus (oak), whose species are known to have very different BVOC emission profiles (Loreto, 2002). The EF for the *Quercus* species is taken from Ciccioli et al. (2023) for isoprene and monoterpenes. For the *Quercus* species missing in Ciccioli et al. (2023) but with a known emission type (Loreto, 2002), the EF values are taken from MEGANv3.2. For unknown emission types, the EF value is set by default to oak isoprene emitters in MEGANv3.2. For all tree species, the EF values for sesquiterpenes and oxygenated VOC are taken from MEGANv3.2. The emission factors of nitric oxide (NO)

and carbon monoxide (CO) are fixed for all tree species and are equal to $EF_{NO} = 0.05$ and $EF_{CO} = 1.0 \mu g \text{ g}_{\text{DW}}^{-1} \text{ h}^{-1}$, as suggested in MEGANv3.2. The emission factors of isoprene (ISOP), total monoterpenes (MTs), total sesquiterpenes (SQTs), and total other VOCs (OVOCs) are shown in Table D2 for the predominant tree genera and oak species.

Emission factors, which are measured at standard conditions, are then multiplied by dimensionless factors representing variations in the emissions as a function of biotic and abiotic processes. A detailed description of the calculation of activity factors for light (photosynthetic photon flux density, PPFD), $\gamma_P$, and for temperature, $\gamma_T$, is given in Appendix D3. To illustrate the variation in the activity factors with light and temperature, Fig. 4 shows averaged $\gamma_T$ and $\gamma_P$ for isoprene, $\alpha$-pinene, and $\beta$-pinene. Figure 4a shows that BVOC emissions increase with temperature. At high temperatures, isoprene emissions are capped, while monoterpene emissions rise sharply. BVOC emissions also increase with light (Fig. 4b), and the activity factor reaches its maximum value of 1 after PPFD = 1000 $\mu mol$ photons m$^{-2}$ s$^{-1}$.

Other activity factors could be added to represent the effect of leaf age and water stress. In this study, emissions are calculated per measure of leaf biomass, considering an average emission for all leaves in the canopy. In addition, we assume that in June and July, tree foliage is fully developed, and leaf area and dry biomass are constant. Therefore, no activity factor is added to modulate emissions according to the fraction of new, growing, mature, and old foliage (Guenther et al., 2012).

Several studies also introduce an activity factor to represent the impact of soil moisture and water stress on isoprene emissions (Guenther et al., 2012; Jiang et al., 2018; Bonn et al., 2019; Otu-Larbi et al., 2020; Wang et al., 2022). Al-
<figure><img src="image_4.png" /><figcaption>Figure 3. (a) Leaf area (LA) and (b) dry biomass (DB) computed for the predominant tree species (P > 1%) in the Paris city inventory as a function of DBH.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>Figure 4. Dependence of activity factors on (a) temperature and (b) light variations for three BVOCs (equations from Guenther et al., 1995, 2012; T24 and T240 are fixed to 294 K in this figure).</figcaption></figure>

though urban trees planted in reduced soil volumes may be subject to water stress (Lüttge and Buckeridge, 2023), the resolution of the CTM does not allow us to accurately simulate the soil water content in an urban environment, so no activity factor modulating isoprene emissions as a function of water content is taken into account here.

## 3.2 Integration of individual tree BVOC emissions in CHIMERE

After estimating the biogenic emissions of each tree in the city of Paris, these emissions are integrated into the CHIMERE CTM. To do this, they must be spatialized and specified, as detailed in this section.

### 3.2.1 Integration of individual tree BVOC emissions on the CTM grid

First, each tree is located within the CTM grid using its precise position given in the Paris tree inventory and the coordinates of the CTM grid. The product of the dry biomass and the emission factor (DBtEFt,k) is then summed for all trees belonging to the same cell to compute the BVOC emission rates (ERi,j,k in µg m⁻² h⁻¹) as

$$ER_{i,j,k} = \frac{1}{\Delta x_{i,j} \Delta y_{i,j}} \sum_{t \in i,j} (DB_t EF_{t,k}) \gamma_{T_{i,j,k}} \gamma_{P_{i,j,k}} \quad (2)$$

where Δx_{i,j} and Δy_{i,j} are the size of the cell (i, j) in the x and y directions (m), here both equal to 1000 m. A map of the dry biomass integrated on the CTM grid cells is shown in Fig. 5.
<figure><img src="image_6.png" /><figcaption>Figure 3. (a) Leaf area (LA) and (b) dry biomass (DB) computed for the predominant tree species (P > 1%) in the Paris city inventory as a function of DBH.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Figure 4. Dependence of activity factors on (a) temperature and (b) light variations for three BVOCs (equations from Guenther et al., 1995, 2012; T24 and T240 are fixed to 294 K in this figure).</figcaption></figure>

though urban trees planted in reduced soil volumes may be subject to water stress (Lüttge and Buckeridge, 2023), the resolution of the CTM does not allow us to accurately simulate the soil water content in an urban environment, so no activity factor modulating isoprene emissions as a function of water content is taken into account here.

## 3.2 Integration of individual tree BVOC emissions in CHIMERE

After estimating the biogenic emissions of each tree in the city of Paris, these emissions are integrated into the CHIMERE CTM. To do this, they must be spatialized and specified, as detailed in this section.

### 3.2.1 Integration of individual tree BVOC emissions on the CTM grid

First, each tree is located within the CTM grid using its precise position given in the Paris tree inventory and the coordinates of the CTM grid. The product of the dry biomass and the emission factor (DBtEFt,k) is then summed for all trees belonging to the same cell to compute the BVOC emission rates (ERi,j,k in µg m⁻² h⁻¹) as

$$ER_{i,j,k} = \frac{1}{\Delta x_{i,j} \Delta y_{i,j}} \sum_{t \in i,j} (DB_t EF_{t,k}) \gamma_{T_{i,j,k}} \cdot \gamma_{P_{i,j,k}}, \quad (2)$$

where Δx_{i,j} and Δy_{i,j} are the size of the cell (i, j) in the x and y directions (m), here both equal to 1000 m. A map of the dry biomass integrated on the CTM grid cells is shown in Fig. 5.
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Tree leaf dry biomass computed over Paris from the Paris tree database and McPherson et al. (2016) with a spatial resolution of 1 km × 1 km.</figcaption></figure>

The average cell dry biomass over Paris is 130 $g m^{-2}$ and can reach 390 $g m^{-2}$ in cells containing large parks or cemeteries. The Paris tree inventory does not include all the trees in the Vincennes and Boulogne woods; however, the emissions of these large woods are already modeled at the regional scale using the land use approach. Thus, they are not considered in the Paris tree inventory added in the simulation bioparis in order to avoid overlapping of emissions (Fig. 6). Note that except for Fig. 2, all the maps presented in this study represent the average value of the variable in the 1 km × 1 km grid cells without any post-processing.

### 3.2.2 Speciation and aggregation of BVOC species

The emission factors of MEGANv3.2 are estimated for different categories of BVOCs, which are presented in the rows of Table D3. These BVOC categories need to be disaggregated into model species to be used in the CTM. The chemical scheme used in CHIMERE corresponds to MELCHIOR2, and the model species are shown in the columns of Table D3. To disaggregate the BVOC categories into model species, the BVOC categories are first speciated into detailed real species, which are then aggregated into the model species. The speciation in real BVOC species is done with a speciation matrix available in MEGANv3.2 code (downloaded at https://bai.ess.uci.edu/megan/data-and-code/megan32, last access: 10 July 2023). Then, the real BVOC species used in CHIMERE are speciated and aggregated into MELCHIOR2 species. The product of the two matrices gives the speciation/aggregation matrix, as described in Table D3. Note that no specific speciation is applied to sesquiterpenes, which are all included in the model species humulene (HUMULE). Monoterpenes are speciated as α-pinene (APINEN), β-pinene (BPINEN), limonene (LIMONE), and ocimene (OCIMEN); other VOCs (OVOCs) represent ethylene (C₂H₄) and oxygenated VOCs (CH₃OH, CH₃CHO, CH₃COE, and MEMALD).

Then, the BVOC emissions from urban trees in Paris are added to the regional-scale BVOC emissions to estimate the BVOC emissions over the Île-de-France region. Section 3.3 below details the complementarity between the bottom-up inventory for urban trees and the regional-scale PFT-based emissions.

### 3.3 Complementarity of the emissions computed by the bottom-up and the land use approaches

At the regional-scale, biogenic emissions are estimated using a land use approach with emission factors that depend on the land use and PFT, as described in Guenther et al. (2012). As the land use is urban over Paris, vegetation is not considered, and there are no biogenic emissions, as shown in Fig. 6a, which represents the 2-month-averaged isoprene emissions computed with the bottom-up inventory, and in Fig. 6b, which represents isoprene emissions computed with the land use approach in CHIMERE.

The bottom-up inventory allows accounting of local biogenic emissions in the city, but Fig. 6 shows that tree inventory and emissions are probably still missing in the Paris suburban area because there is currently no tree inventory for most of the urban areas outside Paris city. The order of magnitude of isoprene emissions computed by the bottom-up inventory seems coherent compared to regional-scale emissions. Emission rates in Paris (0.12 $µg m^{-2} s^{-1}$ on average for isoprene) are lower than those simulated over the large Île-de-France forests. This is also the case for other BVOC species, as shown in Appendix D5. The relative distribution of monoterpenes emitted is different between the urban and the regional scales, as shown in Fig. D3. In particular, there is relatively more β-pinene in the regional-scale emissions. This is due to the different vegetation species between the city and the regional scale and to the speciation of monoterpenes, which may be different.

The temporal variation in the spatially averaged emissions of different biogenic compounds is shown in Fig. 7. For all compounds, emissions are strongly correlated with temperature and sunlight. Over the 2-month periods, there are three emission peaks corresponding to periods of heatwaves with clear-sky conditions and air temperature reaching 35 °C. The impact of BVOC emissions on air quality is expected to be higher during these periods. Therefore, the effect of emissions on pollutant concentrations will be calculated both on the 2-month period and on the heatwave periods, which correspond to the following days: 15 to 18 June, 11 to 14 July, and 17 to 19 July. In terms of emitted compounds, isoprene is the most emitted biogenic species, followed by OVOCs. Monoterpenes and CO are emitted to a lesser extent, followed by sesquiterpenes and NO. This distribution of emissions is fairly typical of emissions calculated using the MEGAN model (Guenther et al., 2012; Ciaccioli et al., 2023). In terms of emission intensity, some recent studies computing the BVOC emissions over Europe
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Tree leaf dry biomass computed over Paris from the Paris tree database and McPherson et al. (2016) with a spatial resolution of 1 km × 1 km.</figcaption></figure>

The average cell dry biomass over Paris is 130 $g m^{-2}$ and can reach 390 $g m^{-2}$ in cells containing large parks or cemeteries. The Paris tree inventory does not include all the trees in the Vincennes and Boulogne woods; however, the emissions of these large woods are already modeled at the regional scale using the land use approach. Thus, they are not considered in the Paris tree inventory added in the simulation bioparis in order to avoid overlapping of emissions (Fig. 6). Note that except for Fig. 2, all the maps presented in this study represent the average value of the variable in the 1 km × 1 km grid cells without any post-processing.

### 3.2.2 Speciation and aggregation of BVOC species

The emission factors of MEGANv3.2 are estimated for different categories of BVOCs, which are presented in the rows of Table D3. These BVOC categories need to be disaggregated into model species to be used in the CTM. The chemical scheme used in CHIMERE corresponds to MELCHIOR2, and the model species are shown in the columns of Table D3. To disaggregate the BVOC categories into model species, the BVOC categories are first speciated into detailed real species, which are then aggregated into the model species. The speciation in real BVOC species is done with a speciation matrix available in MEGANv3.2 code (downloaded at https://bai.ess.uci.edu/megan/data-and-code/megan32, last access: 10 July 2023). Then, the real BVOC species used in CHIMERE are speciated and aggregated into MELCHIOR2 species. The product of the two matrices gives the speciation/aggregation matrix, as described in Table D3. Note that no specific speciation is applied to sesquiterpenes, which are all included in the model species humulene (HUMULE). Monoterpenes are speciated as α-pinene (APINEN), β-pinene (BPINEN), limonene (LIMONE), and ocimene (OCIMEN); other VOCs (OVOCs) represent ethylene (C₂H₄) and oxygenated VOCs (CH₃OH, CH₃CHO, CH₃COE, and MEMALD).

Then, the BVOC emissions from urban trees in Paris are added to the regional-scale BVOC emissions to estimate the BVOC emissions over the Île-de-France region. Section 3.3 below details the complementarity between the bottom-up inventory for urban trees and the regional-scale PFT-based emissions.

### 3.3 Complementarity of the emissions computed by the bottom-up and the land use approaches

At the regional-scale, biogenic emissions are estimated using a land use approach with emission factors that depend on the land use and PFT, as described in Guenther et al. (2012). As the land use is urban over Paris, vegetation is not considered, and there are no biogenic emissions, as shown in Fig. 6a, which represents the 2-month-averaged isoprene emissions computed with the bottom-up inventory, and in Fig. 6b, which represents isoprene emissions computed with the land use approach in CHIMERE.

The bottom-up inventory allows accounting of local biogenic emissions in the city, but Fig. 6 shows that tree inventory and emissions are probably still missing in the Paris suburban area because there is currently no tree inventory for most of the urban areas outside Paris city. The order of magnitude of isoprene emissions computed by the bottom-up inventory seems coherent compared to regional-scale emissions. Emission rates in Paris (0.12 $µg m^{-2} s^{-1}$ on average for isoprene) are lower than those simulated over the large Île-de-France forests. This is also the case for other BVOC species, as shown in Appendix D5. The relative distribution of monoterpenes emitted is different between the urban and the regional scales, as shown in Fig. D3. In particular, there is relatively more β-pinene in the regional-scale emissions. This is due to the different vegetation species between the city and the regional scale and to the speciation of monoterpenes, which may be different.

The temporal variation in the spatially averaged emissions of different biogenic compounds is shown in Fig. 7. For all compounds, emissions are strongly correlated with temperature and sunlight. Over the 2-month periods, there are three emission peaks corresponding to periods of heatwaves with clear-sky conditions and air temperature reaching 35 °C. The impact of BVOC emissions on air quality is expected to be higher during these periods. Therefore, the effect of emissions on pollutant concentrations will be calculated both on the 2-month period and on the heatwave periods, which correspond to the following days: 15 to 18 June, 11 to 14 July, and 17 to 19 July. In terms of emitted compounds, isoprene is the most emitted biogenic species, followed by OVOCs. Monoterpenes and CO are emitted to a lesser extent, followed by sesquiterpenes and NO. This distribution of emissions is fairly typical of emissions calculated using the MEGAN model (Guenther et al., 2012; Ciaccioli et al., 2023). In terms of emission intensity, some recent studies computing the BVOC emissions over Europe
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Comparison of the 2-month-averaged isoprene emissions with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.</figcaption></figure>

with plant-emission-specific models instead of using the PFT approach of MEGAN have reported that isoprene emissions may be overestimated by a factor of 3 in MEGANv2.1, while monoterpene and sesquiterpene emissions may be underestimated by a factor of 3, especially in summer (Jiang et al., 2019; Ciccioli et al., 2023). These discrepancies were attributed to the different vegetation classifications and emission factors at standard conditions. Using plant-emission-specific models, Jiang et al. (2019) found a better comparison to observations for isoprene and organic aerosol concentrations at the European scale, while around the Paris basin in summer, differences in emissions mainly concern monoterpenes and sesquiterpenes. In order to take these emission uncertainties into account in our study, sensitivity simulations are carried out by multiplying monoterpene and sesquiterpene emissions by a factor of 2 or 3.

## 4 Quantification of the impact of BVOC emissions from urban trees on air quality at the regional scale

Before studying the impacts of the bottom-up inventory, comparisons of simulated and observed key variables are performed to evaluate the simulation performance. Meteorological variables are first analyzed in Sect. 4.1.1, as biogenic

emissions are strongly related to them. Then, Sect. 4.1.2 presents comparisons of modeled and observed pollutant concentrations at different air quality stations in Île-de-France. Simulations are performed with the emission factors presented above (REF) and with monoterpene and sesquiterpene emissions multiplied by a factor of 2 (REF-TX2) and 3 (REF-TX3).

Then, to quantify the impacts of urban trees on air quality, simulations with the biogenic emissions from urban trees are performed and compared to the simulations without trees in Sect. 4.2. Three simulations with urban trees are performed, with one for each monoterpene and sesquiterpene emission scenario, which are referred to as bioparis, bioparis-TX2, and bioparis-TX3. All the simulations performed and the corresponding emissions are presented in Table 3.

### 4.1 Validation of the reference simulations

#### 4.1.1 Meteorology

The surface meteorological fields simulated by WRF-CHIMERE are compared to measurements performed at SIRTA (with 10 min time steps) in Fig. E1 and Table E1 and at seven weather stations operated by Météo-France (MF) (with hourly time steps) in Table E2. The comparison shows
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Comparison of the 2-month-averaged isoprene emissions with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.</figcaption></figure>

with plant-emission-specific models instead of using the PFT approach of MEGAN have reported that isoprene emissions may be overestimated by a factor of 3 in MEGANv2.1, while monoterpane and sesquiterpene emissions may be underestimated by a factor of 3, especially in summer (Jiang et al., 2019; Ciccioli et al., 2023). These discrepancies were attributed to the different vegetation classifications and emission factors at standard conditions. Using plant-emission-specific models, Jiang et al. (2019) found a better comparison to observations for isoprene and organic aerosol concentrations at the European scale, while around the Paris basin in summer, differences in emissions mainly concern monoterpenes and sesquiterpenes. In order to take these emission uncertainties into account in our study, sensitivity simulations are carried out by multiplying monoterpane and sesquiterpene emissions by a factor of 2 or 3.

## 4 Quantification of the impact of BVOC emissions from urban trees on air quality at the regional scale

Before studying the impacts of the bottom-up inventory, comparisons of simulated and observed key variables are performed to evaluate the simulation performance. Meteorological variables are first analyzed in Sect. 4.1.1, as biogenic

emissions are strongly related to them. Then, Sect. 4.1.2 presents comparisons of modeled and observed pollutant concentrations at different air quality stations in Île-de-France. Simulations are performed with the emission factors presented above (REF) and with monoterpane and sesquiterpene emissions multiplied by a factor of 2 (REF-TX2) and 3 (REF-TX3).

Then, to quantify the impacts of urban trees on air quality, simulations with the biogenic emissions from urban trees are performed and compared to the simulations without trees in Sect. 4.2. Three simulations with urban trees are performed, with one for each monoterpane and sesquiterpene emission scenario, which are referred to as bioparis, bioparis-TX2, and bioparis-TX3. All the simulations performed and the corresponding emissions are presented in Table 3.

### 4.1 Validation of the reference simulations

#### 4.1.1 Meteorology

The surface meteorological fields simulated by WRF-ChIMERE are compared to measurements performed at SIRTA (with 10 min time steps) in Fig. E1 and Table E1 and at seven weather stations operated by Météo-France (MF) (with hourly time steps) in Table E2. The comparison shows
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Temporal variation in the spatially averaged biogenic emissions computed over Paris with the bottom-up inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

<table><caption>**Table 3.** Simulation list with corresponding emission scenarios. ER stands for emission rate, SQT is for sesquiterpene, MT is for monoterpene, LUA is for land use approach, and BUI is for bottom-up inventory.</caption><thead><tr><th>Simulation name</th><th>Emissions computed with the land use approach over IDF</th><th>Emissions computed with the bottom-up inventory over Paris</th></tr></thead><tbody><tr><td>REF</td><td>Yes</td><td>No</td></tr><tr><td>REF-TX2</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 2</td><td>No</td></tr><tr><td>REF-TX3</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 3</td><td>No</td></tr><tr><td>bioparis</td><td>Yes</td><td>Yes</td></tr><tr><td>bioparis-TX2</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 2</td><td>Yes, ER<sup>BUI</sup><sub>MT&SQT</sub> × 2</td></tr><tr><td>bioparis-TX3</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 3</td><td>Yes, ER<sup>BUI</sup><sub>MT&SQT</sub> × 3</td></tr></tbody></table>

that the model simulates the 2 m air temperature with slight bias of +0.8 to 0.9 °C and the global shortwave radiation with a bias around +50 W m⁻² well. The other variables are also rather well simulated, with a slight overestimation of the relative humidity (bias of -6.3 at SIRTA and -9.4 % at MF stations) and of the surface atmospheric pressure (bias of +2.7 at SIRTA and +1.0 hPa at MF stations). The wind speed is overestimated (bias of +1.9 at SIRTA and +1.2 m s⁻¹ at MF stations) mainly because of the difference in representativeness between the average wind speed simulated in the first 24 m high vertical mesh and the wind measured at 10 m. A detailed description of the validation is provided in Appendix E.

### 4.1.2 Model-data comparisons of gas and particle concentrations

NO₂, O₃, OM, PM₂.₅, isoprene, and monoterpene concentrations simulated by CHIMERE for the three emission scenarios REF, REF-TX2, and REF-TX3 are compared to observations performed at different measurement stations over Île-de-France in Appendix F. As shown in Table F1, NO₂ and O₃

concentrations compare well to observations (the most strict performance criteria defined by Hanna and Chang, 2012, are respected for all statistical indicators), and the impact of emission scenarios is low. PM₂.₅ and OM concentrations are more sensitive to the emission scenarios, and the comparisons to observations are better at suburban and urban stations with the REF-TX3 scenario. Nevertheless, the less strict performance criteria by Hanna and Chang (2012) are respected for all statistical indicators in REF-TX2 and REF-TX3 simulations. The simulations tend to underestimate isoprene and monoterpene concentrations compared to observations. This could be partly explained by the absence of biogenic emissions inside Paris in the reference simulations. A more detailed description of comparisons between simulations and observations is given in Appendix F.

## 4.2 Impact of biogenic emissions from urban trees on isoprene, monoterpene, ozone, organic matter, and PM₂.₅ concentrations

### 4.2.1 Impact of urban tree biogenic emissions on isoprene and monoterpene concentrations

Comparisons of the hourly concentrations of isoprene and monoterpenes observed and simulated in the reference case (REF-TX2) and with the urban tree biogenic emissions (bioparis-TX2) at PRG presented in Fig. 8.

Figure 8a shows that isoprene concentrations simulated at PRG are underestimated in the reference simulation compared to the measurements. The inclusion of the urban tree biogenic emissions allows a better representation of the isoprene concentrations (decrease in the normalized absolute difference (NAD) from 0.57 (REF-TX2) to 0.38 (bioparis-TX2) and increase in the correlation from 0.38 (REF-TX2) to 0.42 (bioparis-TX2)). However, in the bioparis-TX2 simulation, the daytime concentrations are overestimated on 16, 17, 19, and 20 June and between 2 and 10 July by about a factor of 1.5, but the concentration peak around 18 June is underestimated. At night, non-zero isoprene concentrations
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Temporal variation in the spatially averaged biogenic emissions computed over Paris with the bottom-up inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

<table><caption>**Table 3.** Simulation list with corresponding emission scenarios. ER stands for emission rate, SQT is for sesquiterpene, MT is for monoterpene, LUA is for land use approach, and BUI is for bottom-up inventory.</caption><thead><tr><th>Simulation name</th><th>Emissions computed with the land use approach over IDF</th><th>Emissions computed with the bottom-up inventory over Paris</th></tr></thead><tbody><tr><td>REF</td><td>Yes</td><td>No</td></tr><tr><td>REF-TX2</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 2</td><td>No</td></tr><tr><td>REF-TX3</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 3</td><td>No</td></tr><tr><td>bioparis</td><td>Yes</td><td>Yes</td></tr><tr><td>bioparis-TX2</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 2</td><td>Yes, ER<sup>BUI</sup><sub>MT&SQT</sub> × 2</td></tr><tr><td>bioparis-TX3</td><td>Yes, ER<sup>LUA</sup><sub>MT&SQT</sub> × 3</td><td>Yes, ER<sup>BUI</sup><sub>MT&SQT</sub> × 3</td></tr></tbody></table>

that the model simulates the 2 m air temperature with slight bias of +0.8 to 0.9 °C and the global shortwave radiation with a bias around +50 W m⁻² well. The other variables are also rather well simulated, with a slight overestimation of the relative humidity (bias of -6.3 at SIRTA and -9.4 % at MF stations) and of the surface atmospheric pressure (bias of +2.7 at SIRTA and +1.0 hPa at MF stations). The wind speed is overestimated (bias of +1.9 at SIRTA and +1.2 m s⁻¹ at MF stations) mainly because of the difference in representativeness between the average wind speed simulated in the first 24 m high vertical mesh and the wind measured at 10 m. A detailed description of the validation is provided in Appendix E.

### 4.1.2 Model-data comparisons of gas and particle concentrations

NO₂, O₃, OM, PM₂.₅, isoprene, and monoterpene concentrations simulated by CHIMERE for the three emission scenarios REF, REF-TX2, and REF-TX3 are compared to observations performed at different measurement stations over Île-de-France in Appendix F. As shown in Table F1, NO₂ and O₃

concentrations compare well to observations (the most strict performance criteria defined by Hanna and Chang, 2012, are respected for all statistical indicators), and the impact of emission scenarios is low. PM₂.₅ and OM concentrations are more sensitive to the emission scenarios, and the comparisons to observations are better at suburban and urban stations with the REF-TX3 scenario. Nevertheless, the less strict performance criteria by Hanna and Chang (2012) are respected for all statistical indicators in REF-TX2 and REF-TX3 simulations. The simulations tend to underestimate isoprene and monoterpene concentrations compared to observations. This could be partly explained by the absence of biogenic emissions inside Paris in the reference simulations. A more detailed description of comparisons between simulations and observations is given in Appendix F.

## 4.2 Impact of biogenic emissions from urban trees on isoprene, monoterpene, ozone, organic matter, and PM₂.₅ concentrations

### 4.2.1 Impact of urban tree biogenic emissions on isoprene and monoterpene concentrations

Comparisons of the hourly concentrations of isoprene and monoterpenes observed and simulated in the reference case (REF-TX2) and with the urban tree biogenic emissions (bioparis-TX2) at PRG presented in Fig. 8.

Figure 8a shows that isoprene concentrations simulated at PRG are underestimated in the reference simulation compared to the measurements. The inclusion of the urban tree biogenic emissions allows a better representation of the isoprene concentrations (decrease in the normalized absolute difference (NAD) from 0.57 (REF-TX2) to 0.38 (bioparis-TX2) and increase in the correlation from 0.38 (REF-TX2) to 0.42 (bioparis-TX2)). However, in the bioparis-TX2 simulation, the daytime concentrations are overestimated on 16, 17, 19, and 20 June and between 2 and 10 July by about a factor of 1.5, but the concentration peak around 18 June is underestimated. At night, non-zero isoprene concentrations
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Observed and simulated hourly concentrations of (a) isoprene and (b) monoterpenes at the PRG station with (bioparis-TX2) and without (REF-TX2) the bottom-up biogenic emission inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

are measured, and the simulated concentrations are almost zero because isoprene is emitted only during the day by biogenic emissions and has a short lifetime ($\tau_{OH} \approx 1.5$ h with $[OH] = 10^6$ molec. cm$^{-3}$, Atkinson and Arey, 2003b). Isoprene is also emitted by road traffic, according to the VOC speciation used (Theloke and Friedrich, 2007; Baudic et al., 2016), but in the model, traffic emissions are too low at night to represent the measured concentrations. This model-measurement discrepancy could be due to a measurement artifact or to missing anthropogenic sources of isoprene at night. In view of the uncertainties in the measurements, the model provides a satisfactory representation of the order of magnitude of isoprene concentrations.

Table 4 presents the averaged isoprene and monoterpane concentrations and the relative impact of bioparis during the 2-month period and heatwaves. Note that the relative difference in the concentrations is calculated on an hourly time step and then averaged over the 2-month or heatwave periods. The min and max columns in the tables correspond to the minimum and maximum values of the time-averaged concentration or relative difference. As seen previously, biogenic emissions are driven by environmental variables, in particular temperature and solar radiation. To determine whether the effect of local biogenic emissions is greater during heatwaves, isoprene concentrations are also compared during these periods. It is especially relevant to quantify this effect because the frequency of these episodes is expected to increase in the future due to climate change (IPCC, 2021). The

heatwave period refers to the averaged concentrations on the following days: 15, 16, 17, and 18 June and 11, 12, 13, 14, 17, 18, and 19 July. During these periods, high air temperatures and clear-sky conditions were observed as shown in Fig. E1. Table 4 shows that at the scale of the city of Paris, local isoprene emissions significantly increase isoprene concentrations (+1100% on average). The effect of bioparis during the heatwave periods is higher (+2400% on average) because emissions during that period are higher. As the TX2 and TX3 scenarios do not modify isoprene emissions, there is no impact on isoprene concentrations.

The comparison of monoterpane concentrations presented in Fig. 8b shows that monoterpenes are underestimated at PRG in the reference simulation, and the addition of the urban tree emissions strongly increases the monoterpane concentrations. However, the simulated concentrations still underestimate the observations, even with the bioparis-TX2 scenario, probably because of missing anthropogenic sources (Jo et al., 2023). Like isoprene, Table 4 shows that the addition of monoterpane emissions greatly increases the monoterpane concentrations by $6.4 \times 10^{12}$% on average over the 2-month period and by $1.4 \times 10^8$% during the heatwave periods. Monoterpane concentrations logically increase when their emissions are multiplied by 2 (TX2) or 3 (TX3), but the simulated concentrations underestimate the measurements. This discrepancy raises the question of potentially missing local sources of vegetation that emit monoterpenes in the area of measurement.
<figure><img src="image_11.png" /><figcaption>**Figure 8.** Observed and simulated hourly concentrations of (a) isoprene and (b) monoterpenes at the PRG station with (bioparis-TX2) and without (REF-TX2) the bottom-up biogenic emission inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

are measured, and the simulated concentrations are almost zero because isoprene is emitted only during the day by biogenic emissions and has a short lifetime ($\tau_{OH} \approx 1.5$ h with $[OH] = 10^6$ molec. cm$^{-3}$, Atkinson and Arey, 2003b). Isoprene is also emitted by road traffic, according to the VOC speciation used (Theloke and Friedrich, 2007; Baudic et al., 2016), but in the model, traffic emissions are too low at night to represent the measured concentrations. This model-measurement discrepancy could be due to a measurement artifact or to missing anthropogenic sources of isoprene at night. In view of the uncertainties in the measurements, the model provides a satisfactory representation of the order of magnitude of isoprene concentrations.

Table 4 presents the averaged isoprene and monoterpane concentrations and the relative impact of bioparis during the 2-month period and heatwaves. Note that the relative difference in the concentrations is calculated on an hourly time step and then averaged over the 2-month or heatwave periods. The min and max columns in the tables correspond to the minimum and maximum values of the time-averaged concentration or relative difference. As seen previously, biogenic emissions are driven by environmental variables, in particular temperature and solar radiation. To determine whether the effect of local biogenic emissions is greater during heatwaves, isoprene concentrations are also compared during these periods. It is especially relevant to quantify this effect because the frequency of these episodes is expected to increase in the future due to climate change (IPCC, 2021). The

heatwave period refers to the averaged concentrations on the following days: 15, 16, 17, and 18 June and 11, 12, 13, 14, 17, 18, and 19 July. During these periods, high air temperatures and clear-sky conditions were observed as shown in Fig. E1. Table 4 shows that at the scale of the city of Paris, local isoprene emissions significantly increase isoprene concentrations (+1100% on average). The effect of bioparis during the heatwave periods is higher (+2400% on average) because emissions during that period are higher. As the TX2 and TX3 scenarios do not modify isoprene emissions, there is no impact on isoprene concentrations.

The comparison of monoterpane concentrations presented in Fig. 8b shows that monoterpenes are underestimated at PRG in the reference simulation, and the addition of the urban tree emissions strongly increases the monoterpane concentrations. However, the simulated concentrations still underestimate the observations, even with the bioparis-TX2 scenario, probably because of missing anthropogenic sources (Jo et al., 2023). Like isoprene, Table 4 shows that the addition of monoterpane emissions greatly increases the monoterpane concentrations by $6.4 \times 10^{12}$% on average over the 2-month period and by $1.4 \times 10^8$% during the heatwave periods. Monoterpane concentrations logically increase when their emissions are multiplied by 2 (TX2) or 3 (TX3), but the simulated concentrations underestimate the measurements. This discrepancy raises the question of potentially missing local sources of vegetation that emit monoterpenes in the area of measurement.
<table><caption>**Table 4.** Comparison of minimum, mean, and maximum isoprene and monoterpane concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td rowspan="9">Isoprene</td><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.08</td><td>0.11</td><td>0.23</td></tr><tr><td>REF_TX2</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.08</td><td>0.11</td><td>0.23</td></tr><tr><td>REF_TX3</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.09</td><td>0.12</td><td>0.22</td></tr><tr><td>bioparis</td><td>0.04</td><td>0.28</td><td>0.70</td><td>0.12</td><td>0.61</td><td>1.51</td></tr><tr><td>bioparis_TX2</td><td>0.04</td><td>0.28</td><td>0.68</td><td>0.12</td><td>0.61</td><td>1.48</td></tr><tr><td>bioparis_TX3</td><td>0.04</td><td>0.27</td><td>0.67</td><td>0.12</td><td>0.60</td><td>1.45</td></tr><tr><td>Relative difference (%) between</td><td>Min</td><td>Mean</td><td>Max</td><td>Min</td><td>Mean</td><td>Max</td></tr><tr><td>bioparis and REF</td><td>40</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>58</td><td>2.4 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>38</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>53</td><td>2.3 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>36</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>48</td><td>2.3 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td rowspan="9">Monoterpenes</td><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>0.005</td><td>0.009</td><td>0.12</td><td>0.009</td><td>0.016</td><td>0.19</td></tr><tr><td>REF_TX2</td><td>0.01</td><td>0.02</td><td>0.23</td><td>0.02</td><td>0.03</td><td>0.38</td></tr><tr><td>REF_TX3</td><td>0.02</td><td>0.03</td><td>0.35</td><td>0.04</td><td>0.07</td><td>0.57</td></tr><tr><td>bioparis</td><td>0.007</td><td>0.02</td><td>0.12</td><td>0.01</td><td>0.03</td><td>0.21</td></tr><tr><td>bioparis_TX2</td><td>0.01</td><td>0.04</td><td>0.24</td><td>0.03</td><td>0.08</td><td>0.41</td></tr><tr><td>bioparis_TX3</td><td>0.03</td><td>0.07</td><td>0.36</td><td>0.05</td><td>0.13</td><td>0.61</td></tr><tr><td>Relative difference (%) between</td><td>Min</td><td>Mean</td><td>Max</td><td>Min</td><td>Mean</td><td>Max</td></tr><tr><td>bioparis and REF</td><td>3.6</td><td>6.4 × 10<sup>12</sup></td><td>9.0 × 10<sup>13</sup></td><td>6.3</td><td>1.4 × 10<sup>8</sup></td><td>1.0 × 10<sup>9</sup></td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>3.5</td><td>1.1 × 10<sup>13</sup></td><td>1.5 × 10<sup>14</sup></td><td>6.2</td><td>1.5 × 10<sup>8</sup></td><td>1.2 × 10<sup>9</sup></td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>3.4</td><td>1.7 × 10<sup>13</sup></td><td>2.0 × 10<sup>14</sup></td><td>6.1</td><td>1.6 × 10<sup>8</sup></td><td>1.6 × 10<sup>9</sup></td></tr></tbody></table>

### 4.2.2 Impact of urban tree biogenic emissions on organic matter and particles concentrations

Figure 9 compares the observed and simulated OM concentrations at the (Fig. 9a) PRG and (Fig. 9b) Halles stations. It shows that the impact of the urban biogenic emissions on the PRG site is smaller for OM concentrations than for isoprene and monoterpane concentrations. The increase in OM concentrations with urban tree biogenic emissions in the Halles site is mainly visible during the heatwaves (Fig. 9b). The urban biogenic emissions lead to an increase in OM concentrations on average over Paris of 4.6% during the 2-month period, as shown in Table 5. The increase in the OM concentrations is slightly larger when terpene emissions are doubled (+5.6%) and tripled (+6.1%). Due to larger biogenic emissions, the increase in OM concentrations is also larger during the heatwave (+5.4%).

The impact of the urban biogenic emissions is also less visible on hourly concentrations of PM<sub>2.5</sub>, so the relative differences in OM and PM<sub>2.5</sub> concentrations are mapped in Figs. 10 and 11. The two top panels present the REF-TX2 concentrations and the relative difference between bioparis-TX2 and REF-TX2 concentrations averaged for the 2-month period. The same maps are presented in the two lower panels but with the concentrations averaged for the heatwave periods.

Figure 10 shows the spatial variability in the local biogenic emission effect and that the increase in OM concentrations is due to emissions from urban trees remaining localized over Paris. It is greater in cells with a large tree biomass (Fig. 5), where biogenic emissions are also larger, in particular monoterpenes (Fig. D2) and sesquiterpenes (Fig. D4). This correlation shows that biogenic emissions from urban trees contribute strongly to local OM formation. The increase in OM concentrations is slightly larger during the heatwave periods, as shown in Table 5. The impact of urban emissions extends a little further during these periods.

The increase in PM<sub>2.5</sub> concentrations is lower than for OM, but the spatial distribution is similar. The impact remains localized over Paris (+0.6% on average) and is strongest during heatwaves (+1.3%), as shown by the maps in Fig. 11 and Table 5. The increase in PM<sub>2.5</sub> is larger when monoterpane and sesquiterpene emissions are doubled (TX2) and tripled (TX3) (Table 5). This underlines the importance of terpenes in the formation of particulate matter.

### 4.2.3 Impact on ozone concentrations

Ozone concentrations also increase with the urban tree biogenic emissions (+1% on average), especially during the heatwave periods (+2.4%). The increase in O<sub>3</sub> concentrations is also mostly localized in the Paris city and extends to
<table><caption>**Table 4.** Comparison of minimum, mean, and maximum isoprene and monoterpane concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td rowspan="10">Isoprene</td><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.08</td><td>0.11</td><td>0.23</td></tr><tr><td>REF_TX2</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.08</td><td>0.11</td><td>0.23</td></tr><tr><td>REF_TX3</td><td>0.03</td><td>0.05</td><td>0.13</td><td>0.09</td><td>0.12</td><td>0.22</td></tr><tr><td>bioparis</td><td>0.04</td><td>0.28</td><td>0.70</td><td>0.12</td><td>0.61</td><td>1.51</td></tr><tr><td>bioparis_TX2</td><td>0.04</td><td>0.28</td><td>0.68</td><td>0.12</td><td>0.61</td><td>1.48</td></tr><tr><td>bioparis_TX3</td><td>0.04</td><td>0.27</td><td>0.67</td><td>0.12</td><td>0.60</td><td>1.45</td></tr><tr><td>Relative difference (%) between bioparis and REF</td><td>40</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>58</td><td>2.4 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>38</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>53</td><td>2.3 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>36</td><td>1.1 × 10<sup>3</sup></td><td>7.1 × 10<sup>3</sup></td><td>48</td><td>2.3 × 10<sup>3</sup></td><td>1.5 × 10<sup>4</sup></td></tr><tr><td rowspan="10">Monoterpenes</td><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>0.005</td><td>0.009</td><td>0.12</td><td>0.009</td><td>0.016</td><td>0.19</td></tr><tr><td>REF_TX2</td><td>0.01</td><td>0.02</td><td>0.23</td><td>0.02</td><td>0.03</td><td>0.38</td></tr><tr><td>REF_TX3</td><td>0.02</td><td>0.03</td><td>0.35</td><td>0.04</td><td>0.07</td><td>0.57</td></tr><tr><td>bioparis</td><td>0.007</td><td>0.02</td><td>0.12</td><td>0.01</td><td>0.03</td><td>0.21</td></tr><tr><td>bioparis_TX2</td><td>0.01</td><td>0.04</td><td>0.24</td><td>0.03</td><td>0.08</td><td>0.41</td></tr><tr><td>bioparis_TX3</td><td>0.03</td><td>0.07</td><td>0.36</td><td>0.05</td><td>0.13</td><td>0.61</td></tr><tr><td>Relative difference (%) between bioparis and REF</td><td>3.6</td><td>6.4 × 10<sup>12</sup></td><td>9.0 × 10<sup>13</sup></td><td>6.3</td><td>1.4 × 10<sup>8</sup></td><td>1.0 × 10<sup>9</sup></td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>3.5</td><td>1.1 × 10<sup>13</sup></td><td>1.5 × 10<sup>14</sup></td><td>6.2</td><td>1.5 × 10<sup>8</sup></td><td>1.2 × 10<sup>9</sup></td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>3.4</td><td>1.7 × 10<sup>13</sup></td><td>2.0 × 10<sup>14</sup></td><td>6.1</td><td>1.6 × 10<sup>8</sup></td><td>1.6 × 10<sup>9</sup></td></tr></tbody></table>

### 4.2.2 Impact of urban tree biogenic emissions on organic matter and particles concentrations

Figure 9 compares the observed and simulated OM concentrations at the (Fig. 9a) PRG and (Fig. 9b) Halles stations. It shows that the impact of the urban biogenic emissions on the PRG site is smaller for OM concentrations than for isoprene and monoterpane concentrations. The increase in OM concentrations with urban tree biogenic emissions in the Halles site is mainly visible during the heatwaves (Fig. 9b). The urban biogenic emissions lead to an increase in OM concentrations on average over Paris of 4.6% during the 2-month period, as shown in Table 5. The increase in the OM concentrations is slightly larger when terpene emissions are doubled (+5.6%) and tripled (+6.1%). Due to larger biogenic emissions, the increase in OM concentrations is also larger during the heatwave (+5.4%).

The impact of the urban biogenic emissions is also less visible on hourly concentrations of PM<sub>2.5</sub>, so the relative differences in OM and PM<sub>2.5</sub> concentrations are mapped in Figs. 10 and 11. The two top panels present the REF-TX2 concentrations and the relative difference between bioparis-TX2 and REF-TX2 concentrations averaged for the 2-month period. The same maps are presented in the two lower panels but with the concentrations averaged for the heatwave periods.

Figure 10 shows the spatial variability in the local biogenic emission effect and that the increase in OM concentrations is due to emissions from urban trees remaining localized over Paris. It is greater in cells with a large tree biomass (Fig. 5), where biogenic emissions are also larger, in particular monoterpenes (Fig. D2) and sesquiterpenes (Fig. D4). This correlation shows that biogenic emissions from urban trees contribute strongly to local OM formation. The increase in OM concentrations is slightly larger during the heatwave periods, as shown in Table 5. The impact of urban emissions extends a little further during these periods.

The increase in PM<sub>2.5</sub> concentrations is lower than for OM, but the spatial distribution is similar. The impact remains localized over Paris (+0.6% on average) and is strongest during heatwaves (+1.3%), as shown by the maps in Fig. 11 and Table 5. The increase in PM<sub>2.5</sub> is larger when monoterpane and sesquiterpene emissions are doubled (TX2) and tripled (TX3) (Table 5). This underlines the importance of terpenes in the formation of particulate matter.

### 4.2.3 Impact on ozone concentrations

Ozone concentrations also increase with the urban tree biogenic emissions (+1% on average), especially during the heatwave periods (+2.4%). The increase in O<sub>3</sub> concentrations is also mostly localized in the Paris city and extends to
<figure><img src="image_10.png" /><figcaption>**Figure 9.** Observed and simulated hourly concentrations of organic matter (OM) at the (a) PRG and (b) Halles stations with (bioparis-TX2) and without (REF-TX2) the urban biogenic emission inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

<table><caption>**Table 5.** Comparison of minimum, mean, and maximum organic matter (OM) and PM<sub>2.5</sub> concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td rowspan="10">OM</td><td>Concentration (µg m<sup>-3</sup>)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>1.27</td><td>1.35</td><td>1.57</td><td>2.91</td><td>3.01</td><td>3.34</td></tr><tr><td>REF_TX2</td><td>2.34</td><td>2.45</td><td>2.79</td><td>5.45</td><td>5.62</td><td>6.31</td></tr><tr><td>REF_TX3</td><td>3.45</td><td>3.58</td><td>4.14</td><td>7.94</td><td>8.16</td><td>9.27</td></tr><tr><td>bioparis</td><td>1.32</td><td>1.40</td><td>1.63</td><td>2.95</td><td>3.13</td><td>3.48</td></tr><tr><td>bioparis_TX2</td><td>2.43</td><td>2.55</td><td>2.87</td><td>5.58</td><td>5.87</td><td>6.59</td></tr><tr><td>bioparis_TX3</td><td>3.57</td><td>3.74</td><td>4.26</td><td>8.15</td><td>8.55</td><td>9.69</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>bioparis and REF</td><td>0.6</td><td>4.60</td><td>11.51</td><td>0.2</td><td>5.44</td><td>14.37</td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>0.6</td><td>5.58</td><td>15.86</td><td>0.3</td><td>6.08</td><td>18.00</td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>0.6</td><td>6.12</td><td>18.42</td><td>0.3</td><td>6.45</td><td>20.58</td></tr><tr><td rowspan="10">PM<sub>2.5</sub></td><td>Concentration (µg m<sup>-3</sup>)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>8.48</td><td>9.25</td><td>10.85</td><td>10.83</td><td>11.64</td><td>13.36</td></tr><tr><td>REF_TX2</td><td>9.85</td><td>10.65</td><td>12.26</td><td>14.13</td><td>14.97</td><td>16.76</td></tr><tr><td>REF_TX3</td><td>11.32</td><td>12.15</td><td>13.84</td><td>17.65</td><td>18.52</td><td>20.45</td></tr><tr><td>bioparis</td><td>8.54</td><td>9.31</td><td>10.91</td><td>10.85</td><td>11.78</td><td>13.53</td></tr><tr><td>bioparis_TX2</td><td>9.98</td><td>10.77</td><td>12.38</td><td>14.18</td><td>15.25</td><td>17.16</td></tr><tr><td>bioparis_TX3</td><td>11.53</td><td>12.34</td><td>14.06</td><td>17.72</td><td>18.96</td><td>21.07</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF and bioparis</td><td>0.12</td><td>0.64</td><td>1.60</td><td>0.09</td><td>1.25</td><td>3.09</td></tr><tr><td>REF_TX2 and bioparis_TX2</td><td>0.20</td><td>1.12</td><td>3.06</td><td>0.14</td><td>2.06</td><td>5.78</td></tr><tr><td>REF_TX3 and bioparis_TX3</td><td>0.25</td><td>1.55</td><td>4.52</td><td>0.17</td><td>2.69</td><td>8.29</td></tr></tbody></table>
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Observed and simulated hourly concentrations of organic matter (OM) at the (a) PRG and (b) Halles stations with (bioparis-TX2) and without (REF-TX2) the urban biogenic emission inventory. Heatwave periods are indicated by shaded orange areas.</figcaption></figure>

<table><caption>**Table 5.** Comparison of minimum, mean, and maximum organic matter (OM) and PM<sub>2.5</sub> concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Species</th><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td rowspan="10">OM</td><td>Concentration (µg m<sup>-3</sup>)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>1.27</td><td>1.35</td><td>1.57</td><td>2.91</td><td>3.01</td><td>3.34</td></tr><tr><td>REF_TX2</td><td>2.34</td><td>2.45</td><td>2.79</td><td>5.45</td><td>5.62</td><td>6.31</td></tr><tr><td>REF_TX3</td><td>3.45</td><td>3.58</td><td>4.14</td><td>7.94</td><td>8.16</td><td>9.27</td></tr><tr><td>bioparis</td><td>1.32</td><td>1.40</td><td>1.63</td><td>2.95</td><td>3.13</td><td>3.48</td></tr><tr><td>bioparis_TX2</td><td>2.43</td><td>2.55</td><td>2.87</td><td>5.58</td><td>5.87</td><td>6.59</td></tr><tr><td>bioparis_TX3</td><td>3.57</td><td>3.74</td><td>4.26</td><td>8.15</td><td>8.55</td><td>9.69</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>bioparis and REF</td><td>0.6</td><td>4.60</td><td>11.51</td><td>0.2</td><td>5.44</td><td>14.37</td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>0.6</td><td>5.58</td><td>15.86</td><td>0.3</td><td>6.08</td><td>18.00</td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>0.6</td><td>6.12</td><td>18.42</td><td>0.3</td><td>6.45</td><td>20.58</td></tr><tr><td rowspan="10">PM<sub>2.5</sub></td><td>Concentration (µg m<sup>-3</sup>)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>8.48</td><td>9.25</td><td>10.85</td><td>10.83</td><td>11.64</td><td>13.36</td></tr><tr><td>REF_TX2</td><td>9.85</td><td>10.65</td><td>12.26</td><td>14.13</td><td>14.97</td><td>16.76</td></tr><tr><td>REF_TX3</td><td>11.32</td><td>12.15</td><td>13.84</td><td>17.65</td><td>18.52</td><td>20.45</td></tr><tr><td>bioparis</td><td>8.54</td><td>9.31</td><td>10.91</td><td>10.85</td><td>11.78</td><td>13.53</td></tr><tr><td>bioparis_TX2</td><td>9.98</td><td>10.77</td><td>12.38</td><td>14.18</td><td>15.25</td><td>17.16</td></tr><tr><td>bioparis_TX3</td><td>11.53</td><td>12.34</td><td>14.06</td><td>17.72</td><td>18.96</td><td>21.07</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF and bioparis</td><td>0.12</td><td>0.64</td><td>1.60</td><td>0.09</td><td>1.25</td><td>3.09</td></tr><tr><td>REF_TX2 and bioparis_TX2</td><td>0.20</td><td>1.12</td><td>3.06</td><td>0.14</td><td>2.06</td><td>5.78</td></tr><tr><td>REF_TX3 and bioparis_TX3</td><td>0.25</td><td>1.55</td><td>4.52</td><td>0.17</td><td>2.69</td><td>8.29</td></tr></tbody></table>
<figure><img src="image_11.png" /><figcaption>**Figure 10.** Average OM concentrations (µg m⁻³) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the OM concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

the Paris suburbs during heatwaves (Fig. 12). Table 6 shows that doubled or tripled monoterpene and sesquiterpene emissions increase ozone concentrations, but the increase is relatively lower than for OM and PM₂.₅. This suggests that ozone formation is less sensitive to monoterpene and sesquiterpene emissions, which mostly impact the formation of organic matter.

The urban biogenic emissions mainly increase O₃ concentrations during the day, as the concentrations of biogenic species are higher, and O₃ is formed during daytime influenced by solar radiation and temperature. The impact of the biogenic bottom-up inventory on maximal daily ozone concentrations (8 h moving average) is also evaluated, as this value is used in the French air quality standards (LCSQA, 2016). The bottom-up inventory increases the ozone maximal 8 h concentrations on average over the 2 months by around 0.6% and by 1.2% during heatwaves in all scenarios. The maximal impact goes from 4.0% to 4.8%, according to the biogenic emission factor scenario, on average over the 2-month period and from 7.6% to 8.5% during heatwaves.

## 5 Conclusions

This study presents the development of an inventory of biogenic emissions from urban trees using a bottom-up approach and based on the city tree inventory, tree allometric relations and empirical emission equations. The emissions are computed for individual urban trees and integrated into CHIMERE-WRF simulations to quantify the impact of urban trees on pollutant concentrations.

Air quality simulations over the Paris region during the 2 months of June and July 2022 lead to simulated NO₂, O₃, and PM₂.₅ concentrations that are globally consistent with measurements. OM, isoprene, and monoterpene concentrations are underestimated, but they increase when emissions from urban trees are taken into account. Over Paris city, urban trees induce a significant increase in OM concentrations of 4.6% on average over the 2 months and of 5.4% during the heatwave periods. This increase can reach 11.5% locally on average over the 2 months and 14.4% during the heatwave period. The increase in OM concentrations is sensitive to monoterpene and sesquiterpene emissions and remains localized over Paris city where the urban trees are located. O₃ concentrations also slightly increase due to the urban tree emissions by 1.0% on average over the 2 months
<figure><img src="image_13.png" /><figcaption>**Figure 10.** Average OM concentrations (µg m⁻³) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the OM concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

the Paris suburbs during heatwaves (Fig. 12). Table 6 shows that doubled or tripled monoterpene and sesquiterpene emissions increase ozone concentrations, but the increase is relatively lower than for OM and PM₂.₅. This suggests that ozone formation is less sensitive to monoterpene and sesquiterpene emissions, which mostly impact the formation of organic matter.

The urban biogenic emissions mainly increase O₃ concentrations during the day, as the concentrations of biogenic species are higher, and O₃ is formed during daytime influenced by solar radiation and temperature. The impact of the biogenic bottom-up inventory on maximal daily ozone concentrations (8 h moving average) is also evaluated, as this value is used in the French air quality standards (LCSQA, 2016). The bottom-up inventory increases the ozone maximal 8 h concentrations on average over the 2 months by around 0.6% and by 1.2% during heatwaves in all scenarios. The maximal impact goes from 4.0% to 4.8%, according to the biogenic emission factor scenario, on average over the 2-month period and from 7.6% to 8.5% during heatwaves.

## 5 Conclusions

This study presents the development of an inventory of biogenic emissions from urban trees using a bottom-up approach and based on the city tree inventory, tree allometric relations and empirical emission equations. The emissions are computed for individual urban trees and integrated into CHIMERE-WRF simulations to quantify the impact of urban trees on pollutant concentrations.

Air quality simulations over the Paris region during the 2 months of June and July 2022 lead to simulated NO₂, O₃, and PM₂.₅ concentrations that are globally consistent with measurements. OM, isoprene, and monoterpene concentrations are underestimated, but they increase when emissions from urban trees are taken into account. Over Paris city, urban trees induce a significant increase in OM concentrations of 4.6% on average over the 2 months and of 5.4% during the heatwave periods. This increase can reach 11.5% locally on average over the 2 months and 14.4% during the heatwave period. The increase in OM concentrations is sensitive to monoterpene and sesquiterpene emissions and remains localized over Paris city where the urban trees are located. O₃ concentrations also slightly increase due to the urban tree emissions by 1.0% on average over the 2 months
<figure><img src="image_12.png" /><figcaption>**Figure 11.** Average PM<sub>2.5</sub> concentrations (µg m<sup>-3</sup>) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the PM<sub>2.5</sub> concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

<table><caption>**Table 6.** Comparison of minimum, mean, and maximum ozone concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>34.92</td><td>41.54</td><td>43.81</td><td>40.96</td><td>46.94</td><td>49.53</td></tr><tr><td>REF_TX2</td><td>35.42</td><td>42.03</td><td>44.30</td><td>41.65</td><td>47.68</td><td>50.28</td></tr><tr><td>REF_TX3</td><td>35.87</td><td>42.46</td><td>44.73</td><td>42.27</td><td>48.34</td><td>50.95</td></tr><tr><td>bioparis</td><td>35.48</td><td>42.05</td><td>44.47</td><td>41.63</td><td>48.26</td><td>51.15</td></tr><tr><td>bioparis_TX2</td><td>36.06</td><td>42.62</td><td>45.05</td><td>42.36</td><td>49.14</td><td>52.10</td></tr><tr><td>bioparis_TX3</td><td>36.59</td><td>43.12</td><td>45.57</td><td>43.01</td><td>49.94</td><td>52.95</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>bioparis and REF</td><td>0.28</td><td>1.03</td><td>2.38</td><td>0.40</td><td>2.42</td><td>5.72</td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>0.32</td><td>1.17</td><td>2.67</td><td>0.48</td><td>2.65</td><td>6.21</td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>0.35</td><td>1.30</td><td>2.94</td><td>0.45</td><td>2.87</td><td>6.66</td></tr></tbody></table>

and by 2.4% during the heatwaves. This increase can locally reach 2.4% on average over the 2 months and 5.7% during the heatwaves. The increase in O<sub>3</sub> concentrations during the heatwave periods extends to the Paris suburbs, which is further than for OM. These values correspond to temporal averages, but the effect of urban emissions on OM, PM<sub>2.5</sub>, and O<sub>3</sub> is higher during the daytime when biogenic emissions and

photolysis occur, aggravating O<sub>3</sub> peaks during heatwaves. This shows that urban tree emissions have a large impact on air quality, and low-emitting tree species should be favored in cities. In particular, we recommend choosing to plant tree species that emit few terpenes.

OM concentrations are particularly sensitive to terpene emissions. It is essential to better estimate terpene emission
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Average PM<sub>2.5</sub> concentrations (µg m<sup>-3</sup>) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the PM<sub>2.5</sub> concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

<table><caption>**Table 6.** Comparison of minimum, mean, and maximum ozone concentrations averaged in Paris for each simulation and relative difference between the bioparis and the reference simulations during the 2-month and heatwave periods.</caption><thead><tr><th rowspan="2">Simulation</th><th colspan="3">2-month period</th><th colspan="3">Heatwave periods</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td>Concentration (ppbv)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>REF</td><td>34.92</td><td>41.54</td><td>43.81</td><td>40.96</td><td>46.94</td><td>49.53</td></tr><tr><td>REF_TX2</td><td>35.42</td><td>42.03</td><td>44.30</td><td>41.65</td><td>47.68</td><td>50.28</td></tr><tr><td>REF_TX3</td><td>35.87</td><td>42.46</td><td>44.73</td><td>42.27</td><td>48.34</td><td>50.95</td></tr><tr><td>bioparis</td><td>35.48</td><td>42.05</td><td>44.47</td><td>41.63</td><td>48.26</td><td>51.15</td></tr><tr><td>bioparis_TX2</td><td>36.06</td><td>42.62</td><td>45.05</td><td>42.36</td><td>49.14</td><td>52.10</td></tr><tr><td>bioparis_TX3</td><td>36.59</td><td>43.12</td><td>45.57</td><td>43.01</td><td>49.94</td><td>52.95</td></tr><tr><td>Relative difference (%) between</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>bioparis and REF</td><td>0.28</td><td>1.03</td><td>2.38</td><td>0.40</td><td>2.42</td><td>5.72</td></tr><tr><td>bioparis_TX2 and REF_TX2</td><td>0.32</td><td>1.17</td><td>2.67</td><td>0.48</td><td>2.65</td><td>6.21</td></tr><tr><td>bioparis_TX3 and REF_TX3</td><td>0.35</td><td>1.30</td><td>2.94</td><td>0.45</td><td>2.87</td><td>6.66</td></tr></tbody></table>

and by 2.4% during the heatwaves. This increase can locally reach 2.4% on average over the 2 months and 5.7% during the heatwaves. The increase in O<sub>3</sub> concentrations during the heatwave periods extends to the Paris suburbs, which is further than for OM. These values correspond to temporal averages, but the effect of urban emissions on OM, PM<sub>2.5</sub>, and O<sub>3</sub> is higher during the daytime when biogenic emissions and

photolysis occur, aggravating O<sub>3</sub> peaks during heatwaves. This shows that urban tree emissions have a large impact on air quality, and low-emitting tree species should be favored in cities. In particular, we recommend choosing to plant tree species that emit few terpenes.

OM concentrations are particularly sensitive to terpene emissions. It is essential to better estimate terpene emission
<figure><img src="image_13.png" /><figcaption>**Figure 12.** Average O₃ concentrations (ppbv, parts per billion by volume) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the O₃ concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

factors of urban and suburban trees. Furthermore, it should be noted that part of the urban vegetation (in private areas) and part of the suburban vegetation were not taken into account in this study, as the tree inventory is only available for the public trees of Paris city. The effect of urban and suburban trees on air quality is therefore probably underestimated, and this may contribute to the underestimation of monoterpane and OM concentrations observed. Tree inventories should be set up systematically in more cities and their suburbs. This could be completed with methods for characterizing urban vegetation using aerial images, for example. This methodology for building a BVOC bottom-up inventory could be easily applied to other cities that have a tree inventory.

Further work would involve improving the estimation of the tree-scale biogenic emissions by improving the spatial resolution of the meteorological fields. Speciation of monoterpenes and oxygenated VOCs emitted into model chemical species is assumed to be identical for each tree species. A speciation of monoterpenes according to the tree species, as done in Steinbrecher et al. (2009), could be introduced. However, as this speciation does not include all the tree species found in Paris, the speciation should be enriched with other data. Finally, urbanization may induce very local modifications of climate in streets, with potentially higher

temperatures, modified solar radiation due to building shading, and water stress if trees are planted in limited soil volume. These processes are not taken into account in this study, where the spatial resolution is typical of regional urban studies, i.e., 1 km × 1 km.

# Appendix A: Detailed description of the experimental measurements

## A1 Measurements performed at SIRTA

Isoprene was measured with a gas chromatograph equipped with a flame ionization detector (GC-FID), airmoVOC C2–C6 (Chromatotec, Saint-Antoine, France). The instrument is described in detail in Gros et al. (2011). Calibration was performed using a NPL (National Physical Laboratory, Teddington, UK) standard. Uncertainty is estimated to be less than 15%. Monoterpenes were measured at SIRTA using a proton transfer reaction quadrupole mass spectrometer (PTR-Q-MS) from Ionicon Analytik (Innsbruck, Austria) with a time resolution of 5 min. This instrument was implemented at SIRTA for long-term measurements in early 2020, and its operating conditions are described in Simon et al. (2023). The ambient air was sampled at 15 m, 1 h blank measurements
<figure><img src="image_15.png" /><figcaption>**Figure 12.** Average O₃ concentrations (ppbv, parts per billion by volume) simulated by CHIMERE (REF-TX2) during (a) the whole period and (c) the heatwave, and the relative difference in the O₃ concentrations with the urban tree biogenic emissions (bioparis-TX2) during (b) the whole period and (d) the heatwave.</figcaption></figure>

factors of urban and suburban trees. Furthermore, it should be noted that part of the urban vegetation (in private areas) and part of the suburban vegetation were not taken into account in this study, as the tree inventory is only available for the public trees of Paris city. The effect of urban and suburban trees on air quality is therefore probably underestimated, and this may contribute to the underestimation of monoterpane and OM concentrations observed. Tree inventories should be set up systematically in more cities and their suburbs. This could be completed with methods for characterizing urban vegetation using aerial images, for example. This methodology for building a BVOC bottom-up inventory could be easily applied to other cities that have a tree inventory.

Further work would involve improving the estimation of the tree-scale biogenic emissions by improving the spatial resolution of the meteorological fields. Speciation of monoterpenes and oxygenated VOCs emitted into model chemical species is assumed to be identical for each tree species. A speciation of monoterpenes according to the tree species, as done in Steinbrecher et al. (2009), could be introduced. However, as this speciation does not include all the tree species found in Paris, the speciation should be enriched with other data. Finally, urbanization may induce very local modifications of climate in streets, with potentially higher

temperatures, modified solar radiation due to building shading, and water stress if trees are planted in limited soil volume. These processes are not taken into account in this study, where the spatial resolution is typical of regional urban studies, i.e., 1 km × 1 km.

# Appendix A: Detailed description of the experimental measurements

## A1 Measurements performed at SIRTA

Isoprene was measured with a gas chromatograph equipped with a flame ionization detector (GC-FID), airmoVOC C2–C6 (Chromatotec, Saint-Antoine, France). The instrument is described in detail in Gros et al. (2011). Calibration was performed using a NPL (National Physical Laboratory, Teddington, UK) standard. Uncertainty is estimated to be less than 15 %. Monoterpenes were measured at SIRTA using a proton transfer reaction quadrupole mass spectrometer (PTR-Q-MS) from Ionicon Analytik (Innsbruck, Austria) with a time resolution of 5 min. This instrument was implemented at SIRTA for long-term measurements in early 2020, and its operating conditions are described in Simon et al. (2023). The ambient air was sampled at 15 m, 1 h blank measurements
were performed every 13 h, and calibrations were done every month with a NPL standard containing α-pinene. Monoterpenes were measured at the mass-to-charge ratio (m/z) 137, and the associated uncertainties for the period of June–July 2022 were 32%, while the mean detection limit was of 25 ppt (parts per trillion).

## A2 Measurements performed at PRG

Gas and aerosol sampling at the PRG site are performed at 30 m a.g.l. VOCs were measured at the PRG site using a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS; PTR 4000 X2, Ionicon Analytik, Austria) equipped with a CHARON inlet, which is already extensively described elsewhere (Jordan et al., 2009; Eichler et al., 2015; Müller et al., 2017; Leglise et al., 2019). The instrument has been programmed to automatically switch between gas and particle phases and was working at 2.6 mbar and at $E/N = 120$ Td. Gas was sampled at the top of a seventh floor building through a 12 m long Teflon tubing with a 17.5 mm inner diameter. The flow in this main line was fixed at $40 \text{ L min}^{-1}$ until a glass manifold where all gas-phase instruments sampled ambient air. Sensitivity and background have been regularly controlled during the course of the experiment using a pure nitrogen cylinder (99.99999% purity; Linde) and a certified gas standard (containing 10 VOCs at 100 ppb; NPL), thus providing quantitative measurement with an uncertainty typically of the order of 10 ppt.

## A3 Measurements performed at the Halles

In the Halles station, NO₂ concentrations are measured by chemiluminescence detection with an AC32M analyzer from ENVEA (formerly Environnement SA) with a measurement uncertainty of 10%. O₃ concentrations are measured by ultraviolet (UV) photometry with an O342e analyzer from ENVEA with a measurement uncertainty of 11%. PM₂.₅ is measured with a Fidas 200 analyzer from Palas and certified technically compliant by the Laboratoire Central de la Surveillance de la Qualité de l'Air (LCSQA) for continuous, real-time regulatory monitoring of PM₁₀ and PM₂.₅ fractions based on the optical detection of light scattered by aerosols (Lorenz-Mie solution). The uncertainties associated with the measurement are estimated to 9%. More information on the certified devices for regulatory air quality measurement is available (in French) at https://www.lcsqa.org/system/files/media/documents/Liste%20appareils%20conforme%20mesure%20 qualit% C3%A9%20air%20M%C3%A0J_13-05-20_v2_0.pdf (last access: 5 May 2024).

<table><caption>**Table A1.** List of Airparif stations with the species measured and used in this study.</caption><thead><tr><th colspan="4">Measured</th></tr></thead><tbody><tr><td>Paris 1er Les Halles</td><td>48.862128° N, 2.3446227° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub>, OM</td></tr><tr><td>Paris 7eme</td><td>48.8571944° N, 2.2932778° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 12eme</td><td>48.8371944° N, 2.3938056° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 13eme</td><td>48.8284722° N, 2.3595583° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Paris 15eme</td><td>48.8303889° N, 2.2698861° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 18eme</td><td>48.8917278° N, 2.345575° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Aubervilliers</td><td>48.9039444° N, 2.3847222° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Argenteuil</td><td>48.8278324° N, 2.3805391° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Bobigny</td><td>48.9024111° N, 2.4526167° E</td><td>Urban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Champigny-sur-Marne</td><td>48.816692° N, 2.516669° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Evry</td><td>48.8276389° N, 2.3267111° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Lognes</td><td>48.8403167° N, 2.6346611° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Montgeron</td><td>48.7065833° N, 2.4570833° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Neuilly-sur-Seine</td><td>48.8813333° N, 2.2773167° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Gennevilliers</td><td>48.9298219° N, 2.291413° E</td><td>Urban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Vitry-sur-Seine</td><td>48.7756628° N, 2.374005° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Gonesse</td><td>48.9908583° N, 2.4447722° E</td><td>Suburban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Mantes-la-Jolie</td><td>48.996225° N, 1.7032972° E</td><td>Suburban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Melun</td><td>48.5281028° N, 2.6539472° E</td><td>Suburban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Fontainebleau forest</td><td>48.4562391° N, 2.6793973° E</td><td>Rural</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Saint-Martin-du-Tertre</td><td>49.1082856° N, 2.1531876° E</td><td>Rural</td><td>O<sub>3</sub>, PM<sub>2.5</sub></td></tr></tbody></table>
were performed every 13 h, and calibrations were done every month with a NPL standard containing α-pinene. Monoterpenes were measured at the mass-to-charge ratio (m/z) 137, and the associated uncertainties for the period of June–July 2022 were 32%, while the mean detection limit was of 25 ppt (parts per trillion).

## A2 Measurements performed at PRG

Gas and aerosol sampling at the PRG site are performed at 30 m a.g.l. VOCs were measured at the PRG site using a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS; PTR 4000 X2, Ionicon Analytik, Austria) equipped with a CHARON inlet, which is already extensively described elsewhere (Jordan et al., 2009; Eichler et al., 2015; Müller et al., 2017; Leglise et al., 2019). The instrument has been programmed to automatically switch between gas and particle phases and was working at 2.6 mbar and at $E/N = 120$ Td. Gas was sampled at the top of a seventh floor building through a 12 m long Teflon tubing with a 17.5 mm inner diameter. The flow in this main line was fixed at $40 \text{ L min}^{-1}$ until a glass manifold where all gas-phase instruments sampled ambient air. Sensitivity and background have been regularly controlled during the course of the experiment using a pure nitrogen cylinder (99.99999% purity; Linde) and a certified gas standard (containing 10 VOCs at 100 ppb; NPL), thus providing quantitative measurement with an uncertainty typically of the order of 10 ppt.

## A3 Measurements performed at the Halles

In the Halles station, $NO_2$ concentrations are measured by chemiluminescence detection with an AC32M analyzer from ENVEA (formerly Environnement SA) with a measurement uncertainty of 10%. $O_3$ concentrations are measured by ultraviolet (UV) photometry with an O342e analyzer from ENVEA with a measurement uncertainty of 11%. $PM_{2.5}$ is measured with a Fidas 200 analyzer from Palas and certified technically compliant by the Laboratoire Central de la Surveillance de la Qualité de l'Air (LCSQA) for continuous, real-time regulatory monitoring of $PM_{10}$ and $PM_{2.5}$ fractions based on the optical detection of light scattered by aerosols (Lorenz-Mie solution). The uncertainties associated with the measurement are estimated to 9%. More information on the certified devices for regulatory air quality measurement is available (in French) at https://www.lcsqa.org/system/files/media/documents/Liste%20appareils%20conforme%20mesure%20 qualit% C3%A9%20air%20M%C3%A0J_13-05-20_v2_0.pdf (last access: 5 May 2024).

<table><caption>**Table A1.** List of Airparif stations with the species measured and used in this study.</caption><thead><tr><th colspan="4">Measured</th></tr></thead><tbody><tr><td>Paris 1er Les Halles</td><td>48.862128° N, 2.3446227° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub>, OM</td></tr><tr><td>Paris 7eme</td><td>48.8571944° N, 2.2932778° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 12eme</td><td>48.8371944° N, 2.3938056° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 13eme</td><td>48.8284722° N, 2.3595583° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Paris 15eme</td><td>48.8303889° N, 2.2698861° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Paris 18eme</td><td>48.8917278° N, 2.345575° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Aubervilliers</td><td>48.9039444° N, 2.3847222° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Argenteuil</td><td>48.8278324° N, 2.3805391° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Bobigny</td><td>48.9024111° N, 2.4526167° E</td><td>Urban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Champigny-sur-Marne</td><td>48.816692° N, 2.516669° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Evry</td><td>48.8276389° N, 2.3267111° E</td><td>Urban background</td><td>NO<sub>2</sub></td></tr><tr><td>Lognes</td><td>48.8403167° N, 2.6346611° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Montgeron</td><td>48.7065833° N, 2.4570833° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Neuilly-sur-Seine</td><td>48.8813333° N, 2.2773167° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Gennevilliers</td><td>48.9298219° N, 2.291413° E</td><td>Urban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Vitry-sur-Seine</td><td>48.7756628° N, 2.374005° E</td><td>Urban background</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Gonesse</td><td>48.9908583° N, 2.4447722° E</td><td>Suburban background</td><td>NO<sub>2</sub>, PM<sub>2.5</sub></td></tr><tr><td>Mantes-la-Jolie</td><td>48.996225° N, 1.7032972° E</td><td>Suburban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Melun</td><td>48.5281028° N, 2.6539472° E</td><td>Suburban background</td><td>NO<sub>2</sub>, O<sub>3</sub></td></tr><tr><td>Fontainebleau forest</td><td>48.4562391° N, 2.6793973° E</td><td>Rural</td><td>NO<sub>2</sub>, O<sub>3</sub>, PM<sub>2.5</sub></td></tr><tr><td>Saint-Martin-du-Tertre</td><td>49.1082856° N, 2.1531876° E</td><td>Rural</td><td>O<sub>3</sub>, PM<sub>2.5</sub></td></tr></tbody></table>
Appendix B: Characteristics of trees in the Paris tree database

<figure><img src="image_14.png" /><figcaption>Figure B1. Screenshot of the Paris tree database near Avenue des Champs-Élysées (Municipality of Paris, 2023).</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure B2. Proportion (%) of each tree genus in Paris (only genera with *P* > 1 % are shown; the rest of the trees are in the “other” category).</figcaption></figure>
Appendix B: Characteristics of trees in the Paris tree database

<figure><img src="image_16.png" /><figcaption>Figure B1. Screenshot of the Paris tree database near Avenue des Champs-Élysées (Municipality of Paris, 2023).</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>Figure B2. Proportion (%) of each tree genus in Paris (only genera with *P* > 1% are shown; the rest of the trees are in the “other” category).</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Figure B3.** Box plot of the (a) trunk circumference and (b) tree height of the most dominant tree genera.</figcaption></figure>

**Appendix C: Average air temperature simulated by WRF**

<figure><img src="image_17.png" /><figcaption>**Figure C1.** The 2-month-averaged air temperature at 2 m simulated by WRF in the IDF1 domain.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure B3.** Box plot of the (a) trunk circumference and (b) tree height of the most dominant tree genera.</figcaption></figure>

## Appendix C: Average air temperature simulated by WRF

<figure><img src="image_19.png" /><figcaption>**Figure C1.** The 2-month-averaged air temperature at 2 m simulated by WRF in the IDF1 domain.</figcaption></figure>
## Appendix D: BVOC emissions

### D1 Estimation of the tree dry biomass

<table><caption>**Table D1.** US reference cities and climates used in the McPherson et al. (2016) study ranked from the closest to the farthest from the Parisian climate. The last column refers to Köppen climate classification (Paris region is Cfb).</caption><thead><tr><th>Rank</th><th>Region Code</th><th>Region Name</th><th>City</th><th>State</th><th>Climate class</th></tr></thead><tbody><tr><td>1</td><td>NoEast</td><td>Northeast</td><td>Queens</td><td>New York</td><td>Cfa</td></tr><tr><td>2</td><td>Piedmt</td><td>South</td><td>Charlotte</td><td>North Carolina</td><td>Cfa</td></tr><tr><td>3</td><td>LoMidW</td><td>Lower Midwest</td><td>Indianapolis</td><td>Indiana</td><td>Cfa</td></tr><tr><td>4</td><td>GulfCo</td><td>Coastal Plain</td><td>Charleston</td><td>South Carolina</td><td>Cfa</td></tr><tr><td>5</td><td>CenFla</td><td>Central Florida</td><td>Orlando</td><td>Florida</td><td>Cfa</td></tr><tr><td>6</td><td>PacfNW</td><td>Pacific northwest</td><td>Longview</td><td>Oregon</td><td>Csb</td></tr><tr><td>7</td><td>TpIntW</td><td>Temperate interior west</td><td>Boise</td><td>Idaho</td><td>Csa</td></tr><tr><td>8</td><td>NoCalC</td><td>Northern California coast</td><td>Berkeley</td><td>California</td><td>Csb</td></tr><tr><td>9</td><td>InlEmp</td><td>Inland empire</td><td>Claremont</td><td>California</td><td>Csb</td></tr><tr><td>10</td><td>SoCalC</td><td>Southern California coast</td><td>Santa Monica</td><td>California</td><td>Csb</td></tr><tr><td>11</td><td>SacVal</td><td>Sacramento Valley</td><td>Sacramento</td><td>California</td><td>Csa</td></tr><tr><td>12</td><td>NMtnPr</td><td>North</td><td>Fort Collins</td><td>Colorado</td><td>Dfb</td></tr><tr><td>13</td><td>InterW</td><td>Interior west</td><td>Albuquerque</td><td>New Mexico</td><td>Bsk</td></tr><tr><td>14</td><td>MidWst</td><td>Midwest</td><td>Minneapolis</td><td>Minnesota</td><td>Dfa</td></tr><tr><td>15</td><td>InlVal</td><td>Inland valleys</td><td>Modesto</td><td>California</td><td>Bsk</td></tr><tr><td>16</td><td>SWDsrt</td><td>Southwest desert</td><td>Glendale</td><td>Arizona</td><td>Bwh</td></tr><tr><td>17</td><td>Tropic</td><td>Tropical</td><td>Honolulu</td><td>Hawaii</td><td>As</td></tr></tbody></table>

<figure><img src="image_18.png" /><figcaption>**Figure D1.** Decision tree to select the tree category to be used for each Paris tree. The tree category and corresponding allometric database refers to McPherson et al. (2016).</figcaption></figure>

For *Planatus × hispanica* (London plane),

$$LA = \exp \left[ a + b \ln(\ln(\text{DBH} + 1)) + \frac{\text{MSE}}{2} \right], \quad (\text{D1})$$

with $a = -2.06877$, $b = 5.77886$, and MSE (mean squared error) = 0.27978.

For *Acer platanoides* (Norway maple),

$$LA = \exp \left[ a + b \ln(\ln(\text{DBH} + 1)) + \left( \sqrt{\text{DBH}} \times \frac{\text{MSE}}{2} \right) \right], \quad (\text{D2})$$

with $a = -0.55184$, $b = 4.27852$, and MSE = 0.07518.

For *Prunus serrulata* (Japanese cherry),

$$LA = a + b \text{DBH} + c \text{DBH}^2 + d \text{DBH}^3, \quad (\text{D3})$$

with $a = -18.045$, $b = 4.6553$, $c = -0.12798$, and $d = 0.00198$, where $a, b, c, d$, and MSE are dimensionless model coefficients.
## Appendix D: BVOC emissions

### D1 Estimation of the tree dry biomass

<table><caption>**Table D1.** US reference cities and climates used in the McPherson et al. (2016) study ranked from the closest to the farthest from the Parisian climate. The last column refers to Köppen climate classification (Paris region is Cfb).</caption><thead><tr><th>Rank</th><th>Region Code</th><th>Region Name</th><th>City</th><th>State</th><th>Climate class</th></tr></thead><tbody><tr><td>1</td><td>NoEast</td><td>Northeast</td><td>Queens</td><td>New York</td><td>Cfa</td></tr><tr><td>2</td><td>Piedmt</td><td>South</td><td>Charlotte</td><td>North Carolina</td><td>Cfa</td></tr><tr><td>3</td><td>LoMidW</td><td>Lower Midwest</td><td>Indianapolis</td><td>Indiana</td><td>Cfa</td></tr><tr><td>4</td><td>GulfCo</td><td>Coastal Plain</td><td>Charleston</td><td>South Carolina</td><td>Cfa</td></tr><tr><td>5</td><td>CenFla</td><td>Central Florida</td><td>Orlando</td><td>Florida</td><td>Cfa</td></tr><tr><td>6</td><td>PacfNW</td><td>Pacific northwest</td><td>Longview</td><td>Oregon</td><td>Csb</td></tr><tr><td>7</td><td>TpIntW</td><td>Temperate interior west</td><td>Boise</td><td>Idaho</td><td>Csa</td></tr><tr><td>8</td><td>NoCalC</td><td>Northern California coast</td><td>Berkeley</td><td>California</td><td>Csb</td></tr><tr><td>9</td><td>InlEmp</td><td>Inland empire</td><td>Claremont</td><td>California</td><td>Csb</td></tr><tr><td>10</td><td>SoCalC</td><td>Southern California coast</td><td>Santa Monica</td><td>California</td><td>Csb</td></tr><tr><td>11</td><td>SacVal</td><td>Sacramento Valley</td><td>Sacramento</td><td>California</td><td>Csa</td></tr><tr><td>12</td><td>NMtnPr</td><td>North</td><td>Fort Collins</td><td>Colorado</td><td>Dfb</td></tr><tr><td>13</td><td>InterW</td><td>Interior west</td><td>Albuquerque</td><td>New Mexico</td><td>Bsk</td></tr><tr><td>14</td><td>MidWst</td><td>Midwest</td><td>Minneapolis</td><td>Minnesota</td><td>Dfa</td></tr><tr><td>15</td><td>InlVal</td><td>Inland valleys</td><td>Modesto</td><td>California</td><td>Bsk</td></tr><tr><td>16</td><td>SWDsrt</td><td>Southwest desert</td><td>Glendale</td><td>Arizona</td><td>Bwh</td></tr><tr><td>17</td><td>Tropic</td><td>Tropical</td><td>Honolulu</td><td>Hawaii</td><td>As</td></tr></tbody></table>

<figure><img src="image_20.png" /><figcaption>**Figure D1.** Decision tree to select the tree category to be used for each Paris tree. The tree category and corresponding allometric database refers to McPherson et al. (2016).</figcaption></figure>

For *Planatus × hispanica* (London plane),

$$LA = \exp \left[ a + b \ln(\ln(\text{DBH} + 1)) + \frac{\text{MSE}}{2} \right], \quad (\text{D1})$$

with $a = -2.06877$, $b = 5.77886$, and MSE (mean squared error) = 0.27978.

For *Acer platanoides* (Norway maple),

$$LA = \exp \left[ a + b \ln(\ln(\text{DBH} + 1)) + \left( \sqrt{\text{DBH}} \times \frac{\text{MSE}}{2} \right) \right], \quad (\text{D2})$$

with $a = -0.55184$, $b = 4.27852$, and MSE = 0.07518.

For *Prunus serrulata* (Japanese cherry),

$$LA = a + b \text{DBH} + c \text{DBH}^2 + d \text{DBH}^3, \quad (\text{D3})$$

with $a = -18.045$, $b = 4.6553$, $c = -0.12798$, and $d = 0.00198$, where $a, b, c, d$, and MSE are dimensionless model coefficients.
## D2 Emission factors

<table><caption>Table D2. Emission factors (EFs in µg g<sub>DW</sub><sup>-1</sup> h<sup>-1</sup>) of BVOCs for the predominant tree genera found in Paris (P > 1 %) and for the predominant *Quercus* species. ISOP is for isoprene, MT is for monoter- pene, SQT is for sesquiterpene, and OVOC is for other VOC.</caption><thead><tr><th>Genus</th><th>Species</th><th>% of trees in Paris</th><th>EF<sub>ISOP</sub></th><th>EF<sub>MT</sub></th><th>EF<sub>SQT</sub></th><th>EFOVOC</th></tr></thead><tbody><tr><td>Platanus</td><td>All</td><td>22.7</td><td>24</td><td>0.51</td><td>0.10</td><td>4.64</td></tr><tr><td>Aesculus</td><td>All</td><td>11.9</td><td>0</td><td>0.58</td><td>0.10</td><td>4.64</td></tr><tr><td>Tilia</td><td>All</td><td>10.3</td><td>0</td><td>0.53</td><td>0.10</td><td>4.64</td></tr><tr><td>Acer</td><td>All</td><td>7.7</td><td>0</td><td>0.51</td><td>0.10</td><td>4.64</td></tr><tr><td>Sophora</td><td>All</td><td>6.3</td><td>5.0</td><td>0.53</td><td>0.10</td><td>4.64</td></tr><tr><td>Prunus</td><td>All</td><td>3.9</td><td>0</td><td>1.18</td><td>0.10</td><td>4.64</td></tr><tr><td>Fraxinus</td><td>All</td><td>2.6</td><td>0</td><td>0.26</td><td>0.10</td><td>4.64</td></tr><tr><td>Pyrus</td><td>All</td><td>2.6</td><td>0</td><td>0.68</td><td>0.10</td><td>4.64</td></tr><tr><td>Celtis</td><td>All</td><td>2.3</td><td>0</td><td>0.33</td><td>0.10</td><td>4.64</td></tr><tr><td>Pinus</td><td>All</td><td>2.2</td><td>0</td><td>1.43</td><td>0.15</td><td>6.94</td></tr><tr><td>Carpinus</td><td>All</td><td>1.5</td><td>0</td><td>1.07</td><td>0.10</td><td>4.64</td></tr><tr><td>Populus</td><td>All</td><td>1.5</td><td>37</td><td>0.44</td><td>0.10</td><td>4.64</td></tr><tr><td>Malus</td><td>All</td><td>1.5</td><td>0</td><td>0.44</td><td>0.10</td><td>4.64</td></tr><tr><td>Corylus</td><td>All</td><td>1.4</td><td>1.0</td><td>1.81</td><td>0.10</td><td>4.64</td></tr><tr><td>Robinia</td><td>All</td><td>1.2</td><td>20</td><td>0.23</td><td>0.10</td><td>4.64</td></tr><tr><td>Ulmus</td><td>All</td><td>1.1</td><td>0</td><td>0.62</td><td>0.10</td><td>4.64</td></tr><tr><td>Taxus</td><td>All</td><td>1.1</td><td>0</td><td>0.58</td><td>0.15</td><td>4.64</td></tr><tr><td>Betula</td><td>All</td><td>1.0</td><td>0</td><td>0.66</td><td>0.10</td><td>4.64</td></tr><tr><td>Gleditsia</td><td>All</td><td>1.0</td><td>0</td><td>0.56</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>ilex</td><td>0.485</td><td>0.1</td><td>43</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>robur</td><td>0.365</td><td>70</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>rubra</td><td>0.272</td><td>35</td><td>0.1</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>cerris</td><td>0.257</td><td>0.1</td><td>0.6</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>petraea</td><td>0.045</td><td>45</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>pubescens</td><td>0.044</td><td>70</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>frainetto</td><td>0.036</td><td>85</td><td>0.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>palustris</td><td>0.035</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>coccinea</td><td>0.025</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>suber</td><td>0.018</td><td>0.2</td><td>20</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>coccifera</td><td>0.016</td><td>0.1</td><td>25</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>phellos</td><td>0.013</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>imbricaria</td><td>0.011</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr></tbody></table>

## D3 Computation of activity factors

For each BVOC, the activity factors for light (PPFD), $\gamma_{P_{k,t \in (i,j)}}$, and for temperature, $\gamma_{T_{k,t \in (i,j)}}$, are computed as the weighted average of a light-dependent (LDF$_k$) and light-independent fraction ($LIF_k = 1 - LDF_k$):

$$ \gamma_{P_{k,t \in (i,j)}} = (1 - \text{LDF}_k) + \text{LDF}_k \gamma_{P \text{ LDF}_{k,t \in (i,j)}}, \quad (\text{D4}) $$

$$ \gamma_{T_{k,t \in (i,j)}} = (1 - \text{LDF}_k) \gamma_{T \text{ LIF}_{k,t \in (i,j)}} + \text{LDF}_k \gamma_{T \text{ LDF}_{k,t \in (i,j)}}. \quad (\text{D5}) $$

The LDF$_k$ factor depends on the BVOC compound and can be found in Table 4 of Guenther et al. (2012).

### Light effect $\gamma_{P_{k,t \in (i,j)}}$

PPFD is the flux of photons in the 400–700 nm spectral range of solar radiation that is used for photosynthesis. It is ex-

pressed in µmol m<sup>-2</sup> s<sup>-1</sup> and is calculated from the simulated solar radiation in the grid cell where the tree is located (i, j) as

$$ \text{PPFD}_{t \in (i, j)} = 4.5 \times 0.5 \times \text{SWg}_{t \in (i, j)}, \quad (\text{D6}) $$

where SWg is the global solar radiation (shortwave), 4.5 is a factor to convert the W m<sup>-2</sup> into µmol m<sup>-2</sup> s<sup>-1</sup>, and 0.5 is an approximation of the fraction of the solar radiation energy that is in the 400–700 nm spectral range (Meek et al., 1984).

As no canopy model is used to consider the shadow effects inside the canopy, no distinction between the sunlit and shaded leaves can be done. Therefore, the dependency with the past PPFD that requires this distinction is not included, and the light activity factor is computed as (Guenther et al., 1995)

$$ \gamma_{P \text{ LDF}_{k,t \in (i,j)}} = \frac{C_P \alpha \text{PPFD}_{t \in (i,j)}}{\sqrt{1 + \alpha^2 \text{PPFD}_{t \in (i,j)}^2}}, \quad (\text{D7}) $$

with $\alpha = 0.004$, and $C_P = 1.03$.

### Temperature effect $\gamma_{T_{k,t \in (i,j)}}$

$$ \gamma_{T \text{ LDF}_{k,t \in (i,j)}} = \frac{E_{\text{opt}_{k,t \in (i,j)}} C_{T_2} \exp(C_{T_1 k} x_{t \in (i,j)})}{C_{T_2} - (C_{T_1 k} [1 - \exp(C_{T_2} x_{t \in (i,j)})])}, \quad (\text{D8}) $$

$$ \gamma_{T \text{ LIF}_{k,t \in (i,j)}} = \exp[\beta_k (T_{t \in (i,j)} - T_s)], \quad (\text{D9}) $$

$$ \text{with } x_{t \in (i,j)} = \frac{1}{R} \left( \frac{1}{T_{\text{opt}_{t \in (i,j)}}} - \frac{1}{T_{t \in (i,j)}} \right), \quad (\text{D10}) $$

and $R = 0.00831$. $T_{t \in (i,j)}$ is the leaf surface temperature (K), which is approximated here by the air temperature at 2 m a.g.l. in the horizontal grid cell (i, j) to which the tree t belongs.

$C_{T_2} = 230, T_{\text{opt}_{t \in (i,j)}}, E_{\text{opt}_{k,t \in (i,j)}}$ are empirical coefficients

$$ T_{\text{opt}_{t \in (i,j)}} = 313 + 0.6 (T_{240_{t \in (i,j)}} - T_s), \quad (\text{D11}) $$

$$ \begin{aligned} E_{\text{opt}_{k,t \in (i,j)}} = & C_{\text{eo}_k} \exp[0.05 (T_{24_{t \in (i,j)}} - T_s)] \\ & \exp[0.05 (T_{240_{t \in (i,j)}} - T_s)], \end{aligned} \quad (\text{D12}) $$

where $T_{24_{t \in (i,j)}}$ and $T_{240_{t \in (i,j)}}$ are the temperature averages over the past 24 and 240 h, and $T_s = 297$ K is the standard condition for leaf temperature. $C_{T_1 k}$, $C_{\text{eo}_k}$, and $\beta_k$ are BVOC-dependent empirical coefficients that can be found in Table 4 of Guenther et al. (2012).
## D2 Emission factors

<table><caption>**Table D2.** Emission factors (EFs in µg g<sub>DW</sub><sup>-1</sup> h<sup>-1</sup>) of BVOCs for the predominant tree genera found in Paris (P > 1 %) and for the predominant *Quercus* species. ISOP is for isoprene, MT is for monoter- pene, SQT is for sesquiterpene, and OVOC is for other VOC.</caption><thead><tr><th>Genus</th><th>Species</th><th>% of trees in Paris</th><th>EF<sub>ISOP</sub></th><th>EF<sub>MT</sub></th><th>EF<sub>SQT</sub></th><th>EFOVOC</th></tr></thead><tbody><tr><td>Platanus</td><td>All</td><td>22.7</td><td>24</td><td>0.51</td><td>0.10</td><td>4.64</td></tr><tr><td>Aesculus</td><td>All</td><td>11.9</td><td>0</td><td>0.58</td><td>0.10</td><td>4.64</td></tr><tr><td>Tilia</td><td>All</td><td>10.3</td><td>0</td><td>0.53</td><td>0.10</td><td>4.64</td></tr><tr><td>Acer</td><td>All</td><td>7.7</td><td>0</td><td>0.51</td><td>0.10</td><td>4.64</td></tr><tr><td>Sophora</td><td>All</td><td>6.3</td><td>5.0</td><td>0.53</td><td>0.10</td><td>4.64</td></tr><tr><td>Prunus</td><td>All</td><td>3.9</td><td>0</td><td>1.18</td><td>0.10</td><td>4.64</td></tr><tr><td>Fraxinus</td><td>All</td><td>2.6</td><td>0</td><td>0.26</td><td>0.10</td><td>4.64</td></tr><tr><td>Pyrus</td><td>All</td><td>2.6</td><td>0</td><td>0.68</td><td>0.10</td><td>4.64</td></tr><tr><td>Celtis</td><td>All</td><td>2.3</td><td>0</td><td>0.33</td><td>0.10</td><td>4.64</td></tr><tr><td>Pinus</td><td>All</td><td>2.2</td><td>0</td><td>1.43</td><td>0.15</td><td>6.94</td></tr><tr><td>Carpinus</td><td>All</td><td>1.5</td><td>0</td><td>1.07</td><td>0.10</td><td>4.64</td></tr><tr><td>Populus</td><td>All</td><td>1.5</td><td>37</td><td>0.44</td><td>0.10</td><td>4.64</td></tr><tr><td>Malus</td><td>All</td><td>1.5</td><td>0</td><td>0.44</td><td>0.10</td><td>4.64</td></tr><tr><td>Corylus</td><td>All</td><td>1.4</td><td>1.0</td><td>1.81</td><td>0.10</td><td>4.64</td></tr><tr><td>Robinia</td><td>All</td><td>1.2</td><td>20</td><td>0.23</td><td>0.10</td><td>4.64</td></tr><tr><td>Ulmus</td><td>All</td><td>1.1</td><td>0</td><td>0.62</td><td>0.10</td><td>4.64</td></tr><tr><td>Taxus</td><td>All</td><td>1.1</td><td>0</td><td>0.58</td><td>0.15</td><td>4.64</td></tr><tr><td>Betula</td><td>All</td><td>1.0</td><td>0</td><td>0.66</td><td>0.10</td><td>4.64</td></tr><tr><td>Gleditsia</td><td>All</td><td>1.0</td><td>0</td><td>0.56</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>ilex</td><td>0.485</td><td>0.1</td><td>43</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>robur</td><td>0.365</td><td>70</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>rubra</td><td>0.272</td><td>35</td><td>0.1</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>cerris</td><td>0.257</td><td>0.1</td><td>0.6</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>petraea</td><td>0.045</td><td>45</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>pubescens</td><td>0.044</td><td>70</td><td>0.3</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>frainetto</td><td>0.036</td><td>85</td><td>0.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>palustris</td><td>0.035</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>coccinea</td><td>0.025</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>suber</td><td>0.018</td><td>0.2</td><td>20</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>coccifera</td><td>0.016</td><td>0.1</td><td>25</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>phellos</td><td>0.013</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr><tr><td>Quercus</td><td>imbricaria</td><td>0.011</td><td>34</td><td>1.0</td><td>0.10</td><td>4.64</td></tr></tbody></table>

## D3 Computation of activity factors

For each BVOC, the activity factors for light (PPFD), $\gamma_{P_{k,t \in (i,j)}}$, and for temperature, $\gamma_{T_{k,t \in (i,j)}}$, are computed as the weighted average of a light-dependent (LDF$_k$) and light-independent fraction ($LIF_k = 1 - LDF_k$):

$$ \gamma_{P_{k,t \in (i,j)}} = (1 - \text{LDF}_k) + \text{LDF}_k \gamma_{P \text{ LDF}_{k,t \in (i,j)}}, \quad (\text{D4}) $$

$$ \gamma_{T_{k,t \in (i,j)}} = (1 - \text{LDF}_k) \gamma_{T \text{ LIF}_{k,t \in (i,j)}} + \text{LDF}_k \gamma_{T \text{ LDF}_{k,t \in (i,j)}}. \quad (\text{D5}) $$

The LDF$_k$ factor depends on the BVOC compound and can be found in Table 4 of Guenther et al. (2012).

### Light effect $\gamma_{P_{k,t \in (i,j)}}$

PPFD is the flux of photons in the 400–700 nm spectral range of solar radiation that is used for photosynthesis. It is ex-

pressed in µmol m<sup>-2</sup> s<sup>-1</sup> and is calculated from the simulated solar radiation in the grid cell where the tree is located (i, j) as

$$ \text{PPFD}_{t \in (i, j)} = 4.5 \times 0.5 \times \text{SWg}_{t \in (i, j)}, \quad (\text{D6}) $$

where SWg is the global solar radiation (shortwave), 4.5 is a factor to convert the W m<sup>-2</sup> into µmol m<sup>-2</sup> s<sup>-1</sup>, and 0.5 is an approximation of the fraction of the solar radiation energy that is in the 400–700 nm spectral range (Meek et al., 1984).

As no canopy model is used to consider the shadow effects inside the canopy, no distinction between the sunlit and shaded leaves can be done. Therefore, the dependency with the past PPFD that requires this distinction is not included, and the light activity factor is computed as (Guenther et al., 1995)

$$ \gamma_{P \text{ LDF}_{k,t \in (i,j)}} = \frac{C_P \alpha \text{PPFD}_{t \in (i,j)}}{\sqrt{1 + \alpha^2 \text{PPFD}_{t \in (i,j)}^2}}, \quad (\text{D7}) $$

with $\alpha = 0.004$, and $C_P = 1.03$.

### Temperature effect $\gamma_{T_{k,t \in (i,j)}}$

$$ \gamma_{T \text{ LDF}_{k,t \in (i,j)}} = \frac{E_{\text{opt}_{k,t \in (i,j)}} C_{T_2} \exp(C_{T_1 k} x_{t \in (i,j)})}{C_{T_2} - (C_{T_1 k} [1 - \exp(C_{T_2} x_{t \in (i,j)})])}, \quad (\text{D8}) $$

$$ \gamma_{T \text{ LIF}_{k,t \in (i,j)}} = \exp[\beta_k (T_{t \in (i,j)} - T_s)], \quad (\text{D9}) $$

$$ \text{with } x_{t \in (i,j)} = \frac{1}{R} \left( \frac{1}{T_{\text{opt}_{t \in (i,j)}}} - \frac{1}{T_{t \in (i,j)}} \right), \quad (\text{D10}) $$

and $R = 0.00831$. $T_{t \in (i,j)}$ is the leaf surface temperature (K), which is approximated here by the air temperature at 2 m a.g.l. in the horizontal grid cell (i, j) to which the tree t belongs.

$C_{T_2} = 230$, $T_{\text{opt}_{t \in (i,j)}}$, $E_{\text{opt}_{k,t \in (i,j)}}$ are empirical coefficients

$$ T_{\text{opt}_{t \in (i,j)}} = 313 + 0.6 (T_{240_{t \in (i,j)}} - T_s), \quad (\text{D11}) $$

$$ \begin{aligned} E_{\text{opt}_{k,t \in (i,j)}} = C_{\text{eo}_k} \exp[0.05 (T_{24_{t \in (i,j)}} - T_s)] \\ \exp[0.05 (T_{240_{t \in (i,j)}} - T_s)], \end{aligned} \quad (\text{D12}) $$

where $T_{24_{t \in (i,j)}}$ and $T_{240_{t \in (i,j)}}$ are the temperature averages over the past 24 and 240 h, and $T_s = 297$ K is the standard condition for leaf temperature. $C_{T_1 k}$, $C_{\text{eo}_k}$, and $\beta_k$ are BVOC-dependent empirical coefficients that can be found in Table 4 of Guenther et al. (2012).
## D4 Aggregation matrix

<table><caption>Table D3. Aggregation matrix of emitted MEGAN v3.2 species to MELCHIOR2 species.</caption><thead><tr><th>MELCHIOR2/ MEGANv3.2 species</th><th>C<sub>5</sub>H<sub>8</sub></th><th>APINEN</th><th>BPINEN</th><th>LIMONE</th><th>TERPEN</th><th>OCIMEN</th><th>HUMULE</th><th>NO</th><th>CO</th><th>CH<sub>3</sub>OH</th><th>C<sub>2</sub>H<sub>4</sub></th><th>CH<sub>3</sub>CHO</th><th>CH<sub>3</sub>COE</th><th>MEMALD</th></tr></thead><tbody><tr><td>ISOP</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MBO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td></tr><tr><td>MT_PINE</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_ACYC</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_CAMP</td><td>0.0</td><td>0.47</td><td>0.53</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_SABI</td><td>0.0</td><td>0.4</td><td>0.0</td><td>0.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_AROM</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>NO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>SQT_HR</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>SQT_LR</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MEOH</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ACTO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>ETOH</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>ACID</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>LVOC</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.36</td><td>0.0</td><td>0.64</td><td>0.0</td></tr><tr><td>OXPROD</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.9</td><td>0.1</td><td>0.0</td></tr><tr><td>STRESS</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>OTHER</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>CO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr></tbody></table>

## D5 Maps of BVOC emissions

<figure><img src="image_19.png" /><figcaption>Figure D2. Comparison of the 2-month-averaged monoterpane emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.</figcaption></figure>
## D4 Aggregation matrix

<table><caption>Table D3. Aggregation matrix of emitted MEGAN v3.2 species to MELCHIOR2 species.</caption><thead><tr><th>MELCHIOR2/ MEGANv3.2 species</th><th>C<sub>5</sub>H<sub>8</sub></th><th>APINEN</th><th>BPINEN</th><th>LIMONE</th><th>TERPEN</th><th>OCIMEN</th><th>HUMULE</th><th>NO</th><th>CO</th><th>CH<sub>3</sub>OH</th><th>C<sub>2</sub>H<sub>4</sub></th><th>CH<sub>3</sub>CHO</th><th>CH<sub>3</sub>COE</th><th>MEMALD</th></tr></thead><tbody><tr><td>ISOP</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MBO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td></tr><tr><td>MT_PINE</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_ACYC</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_CAMP</td><td>0.0</td><td>0.47</td><td>0.53</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_SABI</td><td>0.0</td><td>0.4</td><td>0.0</td><td>0.6</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MT_AROM</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>NO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>SQT_HR</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>SQT_LR</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>MEOH</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>ACTO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>ETOH</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>ACID</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>LVOC</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.36</td><td>0.0</td><td>0.64</td><td>0.0</td></tr><tr><td>OXPROD</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.9</td><td>0.1</td><td>0.0</td></tr><tr><td>STRESS</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>OTHER</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td></tr><tr><td>CO</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>1.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr></tbody></table>

## D5 Maps of BVOC emissions

<figure><img src="image_21.png" /><figcaption>Figure D2. Comparison of the 2-month-averaged monoterpane emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.</figcaption></figure>
<figure><img src="image_20.png" /></figure>

**Figure D3.** Distribution of monoterpene species emitted and summed over the 2 months (a) for the urban trees in Paris computed with the bottom-up inventory and (b) for the vegetation over Île-de-France region computed with the land use approach.

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

**Figure D4.** Comparison of the 2-month sesquiterpene emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.
<figure><img src="image_22.png" /></figure>

**Figure D3.** Distribution of monoterpene species emitted and summed over the 2 months (a) for the urban trees in Paris computed with the bottom-up inventory and (b) for the vegetation over Île-de-France region computed with the land use approach.

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

**Figure D4.** Comparison of the 2-month sesquiterpene emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.
<figure><img src="image_22.png" /></figure>

**Figure D5.** Comparison of the 2-month nitrite oxide (NO) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.

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

**Figure D6.** Comparison of the 2-month carbon monoxide (CO) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.
<figure><img src="image_24.png" /></figure>

**Figure D5.** Comparison of the 2-month nitrite oxide (NO) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.

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

**Figure D6.** Comparison of the 2-month carbon monoxide (CO) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.
<figure><img src="image_24.png" /></figure>

**Figure D7.** Comparison of the 2-month other VOC (OVOC) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.

## Appendix E: Validation of the reference simulations: meteorology

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

**Figure E1.** Comparison of the temporal variation in the (a) air temperature at 2 m height and (b) photosynthetic photon flux density modeled by WRF (mod) and observed (obs) at the SIRTA observatory site (48.717347° N, 2.208868° E). Heatwave periods are indicated by shaded orange areas.
<figure><img src="image_26.png" /></figure>

**Figure D7.** Comparison of the 2-month other VOC (OVOC) emissions computed with (a) the bottom-up inventory and (b) with the land cover approach in CHIMERE over Île-de-France and greater Paris.

## Appendix E: Validation of the reference simulations: meteorology

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

**Figure E1.** Comparison of the temporal variation in the (a) air temperature at 2 m height and (b) photosynthetic photon flux density modeled by WRF (mod) and observed (obs) at the SIRTA observatory site (48.717347° N, 2.208868° E). Heatwave periods are indicated by shaded orange areas.
This section presents a validation of the surface meteorological fields simulated by WRF-CHIMERE by comparing with measurements performed at SIRTA (Fig. E1 and Table E1) and at seven weather stations operated by Météo-France (Table E2). The meteorological measurements of air temperature (*T*), relative humidity (RH), pressure (*P*), precipitation at 2 m, and wind speed and direction at 10 m above ground level, as well as longwave (LW), global shortwave (SW), and PPFD incident radiations at the surface are compared.

The wind speed and direction observed at 10 m are approximated by the value simulated in the grid cell from 0.15 to 24 m that is supposed to represent the field at the mid-cell altitude (i.e., ≈ 12 m). The meteorological fields are extracted from the horizontal cell of the IDF1 domain which includes the station. Figure E1 shows the comparison of modeled and observed air temperature at 2 m height and PPFD, which are the two meteorological variables used to calculate BVOC emissions, from June to July 2022. They are also compared with statistical indicators (defined in Appendix G) in Table E1, along with the other simulated and observed meteorological variables.

Figure E1 and Table E1 show that the variations in the air temperature at 2 m and PPFD are well modeled, with high correlations and low errors. The temperature is slightly overestimated by the model, especially after 16 July, resulting in an average positive bias of about 1 °C. For PPFD, the daily maximum is overestimated on some days, resulting in a positive bias. As PPFD is computed from the global solar radiation (SW), and the bias on this global solar radiation is lower, this overestimation may also come from the conversion coefficient between PPFD and solar radiation. Some tests have been performed to compare the BVOC emission of the bottom-up inventory calculated with the PPFD / SW ratio measured at SIRTA instead of the ratio used by CHIMERE (2.25) and showed that the impact on BVOC emissions was not significant. Other meteorological variables such as air relative humidity, pressure, wind direction, and incident long-wave radiation are also well modeled (Table E1). The wind speed is slightly overestimated, but this may be due to a difference in the representativity between a punctual wind speed measurement and an average value in the 24 m thick vertical cell. The low rainfall intensity limits the significance of statistical indicator calculations, but the temporal comparison (not shown) demonstrates that the rainy days are well represented by the model, but the intensity of heavy rains is underestimated.

Table E2 shows that the comparison of meteorological variables (*T*, RH, *P*, SW, and wind speed and direction) at seven Météo-France stations leads to very similar conclusions to the comparison at SIRTA.

<table><caption>**Table E1.** Statistical indicators for the comparison of the meteorological variables simulated by WRF and observed at the SIRTA observatory site (48.7°N, 2.2°E) at a 10 min time step. RMSE is for root mean square error, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the meteorological variable.</caption><thead><tr><th>Obs. height</th><th>Variable</th><th>Unit</th><th>Mean obs. su</th><th>Mean mod. su</th><th>RMSE su</th><th>NAD –</th><th>Bias su</th><th>R –</th></tr></thead><tbody><tr><td rowspan="4">2 m a.g.l.</td><td><i>T</i></td><td>°C</td><td>20.5</td><td>21.4</td><td>1.8</td><td>0.03</td><td>0.9</td><td>0.96</td></tr><tr><td>RH</td><td>%</td><td>57.4</td><td>51.1</td><td>11.2</td><td>0.08</td><td>-6.3</td><td>0.88</td></tr><tr><td><i>P</i></td><td>hPa</td><td>998.1</td><td>1000.8</td><td>2.8</td><td>0.00</td><td>2.7</td><td>0.99</td></tr><tr><td>Rain</td><td>mm</td><td>0.008</td><td>0.010</td><td>0.2</td><td>0.83</td><td>0.0</td><td>0.05</td></tr><tr><td rowspan="2">10 m a.g.l.</td><td>Wind speed</td><td>m s<sup>-1</sup></td><td>2.6</td><td>4.5</td><td>2.4</td><td>0.28</td><td>1.9</td><td>0.58</td></tr><tr><td>Wind direction</td><td>°</td><td>179.2</td><td>179.2</td><td>88.3</td><td>0.14</td><td>1.6</td><td>0.68</td></tr><tr><td rowspan="3">Surface</td><td>PPFD</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>524.5</td><td>733.5</td><td>399.7</td><td>0.19</td><td>205.8</td><td>0.91</td></tr><tr><td>SW</td><td>W m<sup>-2</sup></td><td>278.6</td><td>326.0</td><td>153.3</td><td>0.12</td><td>45.3</td><td>0.91</td></tr><tr><td>LW</td><td>W m<sup>-2</sup></td><td>347.5</td><td>338.8</td><td>23.3</td><td>0.02</td><td>-8.2</td><td>0.72</td></tr></tbody></table>
This section presents a validation of the surface meteorological fields simulated by WRF-CHIMERE by comparing with measurements performed at SIRTA (Fig. E1 and Table E1) and at seven weather stations operated by Météo-France (Table E2). The meteorological measurements of air temperature (*T*), relative humidity (RH), pressure (*P*), precipitation at 2 m, and wind speed and direction at 10 m above ground level, as well as longwave (LW), global shortwave (SW), and PPFD incident radiations at the surface are compared.

The wind speed and direction observed at 10 m are approximated by the value simulated in the grid cell from 0.15 to 24 m that is supposed to represent the field at the mid-cell altitude (i.e., ≈ 12 m). The meteorological fields are extracted from the horizontal cell of the IDF1 domain which includes the station. Figure E1 shows the comparison of modeled and observed air temperature at 2 m height and PPFD, which are the two meteorological variables used to calculate BVOC emissions, from June to July 2022. They are also compared with statistical indicators (defined in Appendix G) in Table E1, along with the other simulated and observed meteorological variables.

Figure E1 and Table E1 show that the variations in the air temperature at 2 m and PPFD are well modeled, with high correlations and low errors. The temperature is slightly overestimated by the model, especially after 16 July, resulting in an average positive bias of about 1 °C. For PPFD, the daily maximum is overestimated on some days, resulting in a positive bias. As PPFD is computed from the global solar radiation (SW), and the bias on this global solar radiation is lower, this overestimation may also come from the conversion coefficient between PPFD and solar radiation. Some tests have been performed to compare the BVOC emission of the bottom-up inventory calculated with the PPFD / SW ratio measured at SIRTA instead of the ratio used by CHIMERE (2.25) and showed that the impact on BVOC emissions was not significant. Other meteorological variables such as air relative humidity, pressure, wind direction, and incident long-wave radiation are also well modeled (Table E1). The wind speed is slightly overestimated, but this may be due to a difference in the representativity between a punctual wind speed measurement and an average value in the 24 m thick vertical cell. The low rainfall intensity limits the significance of statistical indicator calculations, but the temporal comparison (not shown) demonstrates that the rainy days are well represented by the model, but the intensity of heavy rains is underestimated.

Table E2 shows that the comparison of meteorological variables (*T*, RH, *P*, SW, and wind speed and direction) at seven Météo-France stations leads to very similar conclusions to the comparison at SIRTA.

<table><caption>**Table E1.** Statistical indicators for the comparison of the meteorological variables simulated by WRF and observed at the SIRTA observatory site (48.7°N, 2.2°E) at a 10 min time step. RMSE is for root mean square error, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the meteorological variable.</caption><thead><tr><th>Obs. height</th><th>Variable</th><th>Unit</th><th>Mean obs. su</th><th>Mean mod. su</th><th>RMSE su</th><th>NAD –</th><th>Bias su</th><th>R –</th></tr></thead><tbody><tr><td rowspan="4">2 m a.g.l.</td><td><i>T</i></td><td>°C</td><td>20.5</td><td>21.4</td><td>1.8</td><td>0.03</td><td>0.9</td><td>0.96</td></tr><tr><td>RH</td><td>%</td><td>57.4</td><td>51.1</td><td>11.2</td><td>0.08</td><td>-6.3</td><td>0.88</td></tr><tr><td><i>P</i></td><td>hPa</td><td>998.1</td><td>1000.8</td><td>2.8</td><td>0.00</td><td>2.7</td><td>0.99</td></tr><tr><td>Rain</td><td>mm</td><td>0.008</td><td>0.010</td><td>0.2</td><td>0.83</td><td>0.0</td><td>0.05</td></tr><tr><td rowspan="2">10 m a.g.l.</td><td>Wind speed</td><td>m s<sup>-1</sup></td><td>2.6</td><td>4.5</td><td>2.4</td><td>0.28</td><td>1.9</td><td>0.58</td></tr><tr><td>Wind direction</td><td>°</td><td>179.2</td><td>179.2</td><td>88.3</td><td>0.14</td><td>1.6</td><td>0.68</td></tr><tr><td rowspan="3">Surface</td><td>PPFD</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>524.5</td><td>733.5</td><td>399.7</td><td>0.19</td><td>205.8</td><td>0.91</td></tr><tr><td>SW</td><td>W m<sup>-2</sup></td><td>278.6</td><td>326.0</td><td>153.3</td><td>0.12</td><td>45.3</td><td>0.91</td></tr><tr><td>LW</td><td>W m<sup>-2</sup></td><td>347.5</td><td>338.8</td><td>23.3</td><td>0.02</td><td>-8.2</td><td>0.72</td></tr></tbody></table>
<table><caption>**Table E2.** Statistical indicators for the comparison of the hourly meteorological variables simulated by WRF and observed at seven Météo-France sites (Montsouris (48.8217°N, 2.3378°E), Longchamp (48.8548°N, 2.2337°E), Melun (48.6103°N, 2.6795°E), Trappes (48.7743°N, 2.0098°E), Versailles (48.8033°N, 2.0900°E), Orly (48.7180°N, 2.3970°E), and Roissy (49.0152°N, 2.5343°E)). RMSE is for root mean square error, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the meteorological variable.</caption><thead><tr><th>Obs. height</th><th>Variable</th><th>Unit</th><th>No. of stations</th><th>Mean obs. su</th><th>Mean mod. su</th><th>RMSE su</th><th>NAD –</th><th>Bias su</th><th>R –</th></tr></thead><tbody><tr><td rowspan="3">2 m a.g.l.</td><td>T</td><td>°C</td><td>7</td><td>20.6</td><td>21.1</td><td>2.2</td><td>0.04</td><td>0.8</td><td>0.91</td></tr><tr><td>RH</td><td>%</td><td>7</td><td>61.9</td><td>52.7</td><td>13.6</td><td>0.10</td><td>-9.4</td><td>0.85</td></tr><tr><td>P</td><td>hPa</td><td>5</td><td>1005.9</td><td>1006.5</td><td>1.2</td><td>0.00</td><td>1.0</td><td>0.99</td></tr><tr><td rowspan="2">10 m a.g.l.</td><td>Wind speed</td><td>m s<sup>-1</sup></td><td>6</td><td>2.7</td><td>4.4</td><td>1.8</td><td>0.20</td><td>1.2</td><td>0.58</td></tr><tr><td>Wind direction</td><td>°</td><td>6</td><td>176.4</td><td>176.8</td><td>101.1</td><td>0.15</td><td>1.4</td><td>0.62</td></tr><tr><td>Surface</td><td>SW</td><td>W m<sup>-2</sup></td><td>7</td><td>273.2</td><td>325.4</td><td>130.8</td><td>0.11</td><td>52.9</td><td>0.94</td></tr></tbody></table>

## Appendix F: Validation of the reference simulations: model to data comparisons of gas and particle concentrations

<table><caption>**Table F1.** Statistical comparison of the observed and simulated concentrations on average over 21 stations in IDF1 (listed in Table A1). Values indicated in bold respect the strictest performance criteria, while those in italics respect the acceptable performance criteria, and those in roman do not respect the performance criteria defined by Hanna and Chang (2012). Correlation coefficients (*R*) are not included in the performance criteria. FB is for fractional bias, MG is for geometric mean bias, NMSE is for normalized mean square error, VG is for geometric variance, FAC2 is for factor of 2, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the species concentration.</caption><thead><tr><th>Species (unit)</th><th>Simulation</th><th>Nb stat.</th><th>Obs. su</th><th>Sim. su</th><th>FB –</th><th>MG –</th><th>NMSE –</th><th>VG –</th><th>FAC2 –</th><th>NAD –</th><th>R –</th></tr></thead><tbody><tr><td rowspan="3">NO<sub>2</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>15.8</td><td><strong>0.05</strong></td><td><strong>1.16</strong></td><td><strong>0.52</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td>REF-TX2</td><td>20</td><td>14.6</td><td>15.7</td><td><strong>0.04</strong></td><td><strong>1.15</strong></td><td><strong>0.52</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>15.6</td><td><strong>0.04</strong></td><td><strong>1.15</strong></td><td><strong>0.53</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td rowspan="3">O<sub>3</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>82.9</td><td><strong>0.19</strong></td><td><strong>1.29</strong></td><td><strong>0.13</strong></td><td><strong>1.31</strong></td><td><strong>0.85</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td>REF-TX2</td><td>12</td><td>68.2</td><td>83.8</td><td><strong>0.21</strong></td><td><strong>1.30</strong></td><td><strong>0.14</strong></td><td><strong>1.32</strong></td><td><strong>0.84</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>84.7</td><td><strong>0.22</strong></td><td><strong>1.32</strong></td><td><strong>0.14</strong></td><td><strong>1.33</strong></td><td><strong>0.84</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td rowspan="3">PM<sub>2.5</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>8.4</td><td><strong>0.17</strong></td><td><strong>1.33</strong></td><td><strong>0.51</strong></td><td><strong>1.68</strong></td><td><strong>0.72</strong></td><td><strong>0.23</strong></td><td>0.41</td></tr><tr><td>REF-TX2</td><td>8</td><td>7.2</td><td>10.4</td><td>0.37</td><td>1.60</td><td><strong>0.81</strong></td><td>2.01</td><td><strong>0.65</strong></td><td><strong>0.26</strong></td><td>0.47</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>12.4</td><td>0.53</td><td>1.88</td><td><strong>1.29</strong></td><td>2.60</td><td><strong>0.57</strong></td><td>0.31</td><td>0.49</td></tr><tr><td rowspan="3">OM (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>1.4</td><td>-0.99</td><td>0.29</td><td><strong>2.24</strong></td><td>7.59</td><td>0.16</td><td>0.50</td><td>0.58</td></tr><tr><td>REF-TX2</td><td>3</td><td>4.3</td><td>2.6</td><td>-0.46</td><td>0.54</td><td><strong>0.78</strong></td><td>2.51</td><td>0.37</td><td>0.35</td><td>0.59</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>3.4</td><td><strong>-0.09</strong></td><td><strong>0.81</strong></td><td><strong>0.58</strong></td><td>1.80</td><td>0.49</td><td><strong>0.26</strong></td><td>0.59</td></tr><tr><td rowspan="3">C<sub>5</sub>H<sub>8</sub> (ppb vol)</td><td>REF</td><td></td><td></td><td>0.09</td><td>-0.89</td><td>0.25</td><td>5.95</td><td>12.27</td><td>0.15</td><td>0.54</td><td>0.54</td></tr><tr><td>REF-TX2</td><td>2</td><td>0.29</td><td>0.09</td><td>-0.89</td><td>0.26</td><td>5.88</td><td>10.31</td><td>0.15</td><td>0.54</td><td>0.54</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>0.09</td><td>-0.89</td><td>0.27</td><td>5.80</td><td>9.23</td><td>0.15</td><td>0.54</td><td>0.55</td></tr><tr><td rowspan="3">MTs (ppb vol)</td><td>REF</td><td></td><td></td><td>0.04</td><td>-0.91</td><td>0.43</td><td>12.53</td><td>2.8 × 10<sup>15</sup></td><td>0.26</td><td>0.60</td><td>0.14</td></tr><tr><td>REF-TX2</td><td>2</td><td>0.09</td><td>0.10</td><td>-0.39</td><td><strong>0.89</strong></td><td>6.11</td><td>4.0 × 10<sup>12</sup></td><td>0.22</td><td>0.58</td><td>0.13</td></tr><tr><td>REF-TX3</td><td></td><td></td><td>0.16</td><td><strong>-0.03</strong></td><td>1.36</td><td>5.49</td><td>1.6 × 10<sup>11</sup></td><td>0.15</td><td>0.59</td><td>0.14</td></tr></tbody></table>
<table><caption>**Table E2.** Statistical indicators for the comparison of the hourly meteorological variables simulated by WRF and observed at seven Météo-France sites (Montsouris (48.8217°N, 2.3378°E), Longchamp (48.8548°N, 2.2337°E), Melun (48.6103°N, 2.6795°E), Trappes (48.7743°N, 2.0098°E), Versailles (48.8033°N, 2.0900°E), Orly (48.7180°N, 2.3970°E), and Roissy (49.0152°N, 2.5343°E)). RMSE is for root mean square error, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the meteorological variable.</caption><thead><tr><th>Obs. height</th><th>Variable</th><th>Unit</th><th>No. of stations</th><th>Mean obs. su</th><th>Mean mod. su</th><th>RMSE su</th><th>NAD –</th><th>Bias su</th><th>R –</th></tr></thead><tbody><tr><td rowspan="3">2 m a.g.l.</td><td>T</td><td>°C</td><td>7</td><td>20.6</td><td>21.1</td><td>2.2</td><td>0.04</td><td>0.8</td><td>0.91</td></tr><tr><td>RH</td><td>%</td><td>7</td><td>61.9</td><td>52.7</td><td>13.6</td><td>0.10</td><td>-9.4</td><td>0.85</td></tr><tr><td>P</td><td>hPa</td><td>5</td><td>1005.9</td><td>1006.5</td><td>1.2</td><td>0.00</td><td>1.0</td><td>0.99</td></tr><tr><td rowspan="2">10 m a.g.l.</td><td>Wind speed</td><td>m s<sup>-1</sup></td><td>6</td><td>2.7</td><td>4.4</td><td>1.8</td><td>0.20</td><td>1.2</td><td>0.58</td></tr><tr><td>Wind direction</td><td>°</td><td>6</td><td>176.4</td><td>176.8</td><td>101.1</td><td>0.15</td><td>1.4</td><td>0.62</td></tr><tr><td>Surface</td><td>SW</td><td>W m<sup>-2</sup></td><td>7</td><td>273.2</td><td>325.4</td><td>130.8</td><td>0.11</td><td>52.9</td><td>0.94</td></tr></tbody></table>

## Appendix F: Validation of the reference simulations: model to data comparisons of gas and particle concentrations

<table><caption>**Table F1.** Statistical comparison of the observed and simulated concentrations on average over 21 stations in IDF1 (listed in Table A1). Values indicated in bold respect the strictest performance criteria, while those in italics respect the acceptable performance criteria, and those in roman do not respect the performance criteria defined by Hanna and Chang (2012). Correlation coefficients (*R*) are not included in the performance criteria. FB is for fractional bias, MG is for geometric mean bias, NMSE is for normalized mean square error, VG is for geometric variance, FAC2 is for factor of 2, NAD is for normalized absolute difference, Bias is for fractional mean bias, *R* is for Pearson correlation coefficient (Appendix G), a.g.l. is for above ground level, and *su* is for the same unit as the species concentration.</caption><thead><tr><th>Species (unit)</th><th>Simulation</th><th>Nb stat.</th><th>Obs. su</th><th>Sim. su</th><th>FB –</th><th>MG –</th><th>NMSE –</th><th>VG –</th><th>FAC2 –</th><th>NAD –</th><th>R –</th></tr></thead><tbody><tr><td rowspan="3">NO<sub>2</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>15.8</td><td><strong>0.05</strong></td><td><strong>1.16</strong></td><td><strong>0.52</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td>REF-TX2</td><td rowspan="2">20</td><td rowspan="2">14.6</td><td>15.7</td><td><strong>0.04</strong></td><td><strong>1.15</strong></td><td><strong>0.52</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td>REF-TX3</td><td>15.6</td><td><strong>0.04</strong></td><td><strong>1.15</strong></td><td><strong>0.53</strong></td><td><strong>1.52</strong></td><td><strong>0.67</strong></td><td><strong>0.24</strong></td><td>0.54</td></tr><tr><td rowspan="3">O<sub>3</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>82.9</td><td><strong>0.19</strong></td><td><strong>1.29</strong></td><td><strong>0.13</strong></td><td><strong>1.31</strong></td><td><strong>0.85</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td>REF-TX2</td><td rowspan="2">12</td><td rowspan="2">68.2</td><td>83.8</td><td><strong>0.21</strong></td><td><strong>1.30</strong></td><td><strong>0.14</strong></td><td><strong>1.32</strong></td><td><strong>0.84</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td>REF-TX3</td><td>84.7</td><td><strong>0.22</strong></td><td><strong>1.32</strong></td><td><strong>0.14</strong></td><td><strong>1.33</strong></td><td><strong>0.84</strong></td><td><strong>0.14</strong></td><td>0.69</td></tr><tr><td rowspan="3">PM<sub>2.5</sub> (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>8.4</td><td><strong>0.17</strong></td><td><strong>1.33</strong></td><td><strong>0.51</strong></td><td><strong>1.68</strong></td><td><strong>0.72</strong></td><td><strong>0.23</strong></td><td>0.41</td></tr><tr><td>REF-TX2</td><td rowspan="2">8</td><td rowspan="2">7.2</td><td>10.4</td><td><strong>0.37</strong></td><td><strong>1.60</strong></td><td><strong>0.81</strong></td><td><strong>2.01</strong></td><td><strong>0.65</strong></td><td><strong>0.26</strong></td><td>0.47</td></tr><tr><td>REF-TX3</td><td>12.4</td><td><strong>0.53</strong></td><td><strong>1.88</strong></td><td><strong>1.29</strong></td><td><strong>2.60</strong></td><td><strong>0.57</strong></td><td><strong>0.31</strong></td><td>0.49</td></tr><tr><td rowspan="3">OM (µg m<sup>-3</sup>)</td><td>REF</td><td></td><td></td><td>1.4</td><td>-0.99</td><td>0.29</td><td><strong>2.24</strong></td><td>7.59</td><td>0.16</td><td>0.50</td><td>0.58</td></tr><tr><td>REF-TX2</td><td rowspan="2">3</td><td rowspan="2">4.3</td><td>2.6</td><td>-0.46</td><td>0.54</td><td><strong>0.78</strong></td><td>2.51</td><td>0.37</td><td>0.35</td><td>0.59</td></tr><tr><td>REF-TX3</td><td>3.4</td><td><strong>-0.09</strong></td><td><strong>0.81</strong></td><td><strong>0.58</strong></td><td>1.80</td><td>0.49</td><td><strong>0.26</strong></td><td>0.59</td></tr><tr><td rowspan="3">C<sub>5</sub>H<sub>8</sub> (ppb vol)</td><td>REF</td><td></td><td></td><td>0.09</td><td>-0.89</td><td>0.25</td><td>5.95</td><td>12.27</td><td>0.15</td><td>0.54</td><td>0.54</td></tr><tr><td>REF-TX2</td><td rowspan="2">2</td><td rowspan="2">0.29</td><td>0.09</td><td>-0.89</td><td>0.26</td><td>5.88</td><td>10.31</td><td>0.15</td><td>0.54</td><td>0.54</td></tr><tr><td>REF-TX3</td><td>0.09</td><td>-0.89</td><td>0.27</td><td>5.80</td><td>9.23</td><td>0.15</td><td>0.54</td><td>0.55</td></tr><tr><td rowspan="3">MTs (ppb vol)</td><td>REF</td><td></td><td></td><td>0.04</td><td>-0.91</td><td>0.43</td><td>12.53</td><td>2.8 × 10<sup>15</sup></td><td>0.26</td><td>0.60</td><td>0.14</td></tr><tr><td>REF-TX2</td><td rowspan="2">2</td><td rowspan="2">0.09</td><td>0.10</td><td>-0.39</td><td><strong>0.89</strong></td><td>6.11</td><td>4.0 × 10<sup>12</sup></td><td>0.22</td><td>0.58</td><td>0.13</td></tr><tr><td>REF-TX3</td><td>0.16</td><td><strong>-0.03</strong></td><td>1.36</td><td>5.49</td><td>1.6 × 10<sup>11</sup></td><td>0.15</td><td>0.59</td><td>0.14</td></tr></tbody></table>