The Eulerian urban dispersion model EPISODE – Part 2: Extensions to the source dispersion and photochemistry for EPISODE–CityChem v1.2 and its application to the city of Hamburg

This paper describes the CityChem extension of the Eulerian urban dispersion model EPISODE. The development of the CityChem extension was driven by the need to apply the model in largely populated urban areas with highly complex pollution sources of particulate matter and various gaseous pollutants. The CityChem extension offers a more advanced treatment of the photochemistry in urban areas and entails specific developments within the sub-grid components for a more accurate representation of dispersion in proximity to urban emission sources. Photochemistry on the Eulerian grid is computed using a numerical chemistry solver. Photochemistry in the sub-grid components is solved with a compact reaction scheme, replacing the photostationary-state assumption. The simplified street canyon model (SSCM) is used in the line source sub-grid model to calculate pollutant dispersion in street canyons. The WMPP (WORM Meteorological Pre-Processor) is used in the point source sub-grid model to calculate the wind speed at plume height. The EPISODE–CityChem model integrates the CityChem extension in EPISODE, with the capability of simulating the photochemistry and dispersion of multiple reactive pollutants within urban areas. The main focus of the model is the simulation of the complex atmospheric chemistry involved in the photochemical production of ozone in urban areas. The ability of EPISODE–CityChem to reproduce the temporal variation of major regulated pollutants at air quality monitoring stations in Hamburg, Germany, was compared to that of the standard EPISODE model and the TAPM (The Air Pollution Model) air quality model using identical meteorological fields and emissions. EPISODE–CityChem performs better than EPISODE and TAPM for the prediction of hourly NO2 concentrations at the traffic stations, which is attributable to the street canyon model. Observed levels of annual mean ozone at the five urban background stations in Hamburg are captured by the model within ±15 %. A performance analysis with the FAIRMODE DELTA tool for air quality in Hamburg showed that EPISODE–CityChem fulfils the model performance objectives for NO2 (hourly), O3 (daily max. of the 8 h running mean) and PM10 (daily mean) set forth in the Air Quality Directive, qualifying the model for use in policy applications. Envisaged applications of the EPISODE–CityChem model are urban air quality studies, emission control scenarios in relation to traffic restrictions and the source attribution of sector-specific emissions to observed levels of air pollutants at urban monitoring stations.

air pollution due to sulfur dioxide (SO 2 ), carbon monoxide (CO), lead (Pb) and benzene have been reduced during the last 2 decades due to emission abatement measures. Tropospheric (ground-level) ozone is a secondary pollutant generated in in photochemical reaction cycles involving two classes of precursor compounds, i.e. nitrogen oxides and volatile organic compounds (VOCs), initiated by the reaction of the hydroxyl (OH) radical with organic molecules. For health protection, a maximum daily 8 h mean threshold for ozone (120 µg m −3 ) is specified as a target value in the European Union, which should not be exceeded at any AQ monitoring station on more than 25 d yr −1 . However, about 15 % of the population living in urban areas is exposed to ozone concentrations above the European Union (EU) target value (EEA, 2015). Traffic is a major source of nitrogen oxides (NO x = NO 2 + NO) and highly contributes to the population exposure to ambient NO 2 concentrations in urban areas because these emissions occur close to the ground and are distributed across densely populated areas. Urban emissions of ozone precursors are transported by local and/or regional air mass flows towards suburban and rural areas, which can be impacted by O 3 pollution episodes (Querol et al., 2016).
Eulerian chemistry-transport model (CTM) systems using numerical methods for solving photochemistry (including chemical reaction schemes with varying degrees of detail) have mainly been used for regional-scale air quality studies. Recent nested model approaches using regional CTM systems have been applied to capture pollution processes from the continental scale to the local scale using between 1 and 5 km resolution and a temporal resolution of 1 h for the innermost domain (e.g. Borge et al., 2014;Karl et al., 2015;Petetin et al., 2015;Valverde et al., 2016). Regional AQ models can give a reliable representation of O 3 concentrations in the urban background, but due to their limitation in resolving the near-field dispersion of emission sources and photochemistry at the sub-kilometre scale, i.e. in street canyons, around industrial stacks and on the neighbourhood level, they cannot provide the information needed by urban policymakers for population exposure mapping, city planning and the assessment of abatement measures.
Urban-scale AQ models overcome the limitation inherent in regional-scale models by taking into account details of the urban topography, wind flow field characteristics, land use information and the geometry of local pollution sources. The urban AQ model EPISODE developed at the Norwegian Institute for Air Research (NILU) is a 3-D Eulerian grid model that operates as a CTM, offline coupled with a numerical weather prediction (NWP) model. EPISODE is typically applied with a horizontal resolution of 1 × 1 km 2 over an entire city with domains of up to 2500 km 2 in size. The Eulerian grid component of EPISODE simulates advection, vertical and/or horizontal diffusion, background transport across the model domain boundaries, and photochemistry. Several sub-grid-scale modules are embedded in EPISODE to represent emissions (line source and point sources), Gaussian dispersion and local photochemistry. In particular, the model allows the user to retrieve concentrations at the sub-grid scale in specified locations of the urban area. Moreover, the EPISODE model is an integral part of the operational Air Quality Information System AirQUIS 2006 (Slørdal et al., 2008).
Part one (Hamer et al., 2019) of this two-part article series provides a detailed description of the EPISODE model system, including the physical processes for atmospheric pollutant transport, the photo-stationary-state (PSS) approximation, the involvement of nitric oxide (NO), NO 2 and O 3 , sub-grid components, and the interaction between the Eulerian grid and the sub-grid processing of pollutant concentrations. Part one examines the application of EPISODE to air quality scenarios in the Nordic winter setting. During wintertime in northern Europe, the PSS assumption is a rather good approximation of the photochemical conversion occurring close to the emission sources. However, when the solar ultraviolet (UV) radiation is stronger, in particular during summer months or at more southerly locations, net ozone formation may take place in urban areas at a certain distance from the main local emission sources (Baklanov et al., 2007). EPISODE in its routine application does not allow for the treatment of photochemistry involving VOCs and other reactive gases leading to the photochemical formation of ozone.
In this part, the features of the CityChem extension for treating the complex atmospheric chemistry in urban areas and specific developments within the sub-grid components for a more accurate representation of near-field dispersion in proximity to urban emission sources are described. Atmospheric chemistry on an urban scale is complex due to the large spatial variations of input from anthropogenic emissions. VOCs related to emissions from traffic are involved in chemical conversion in urban areas. Therefore, it has become necessary to simulate a large number of chemical interactions involving NO x , O 3 , VOCs, SO 2 and secondary pollutants. In order to use comprehensive photochemical schemes in urban AQ models involving VOC interactions, the highest priority for the initial development was to reduce the number of compounds and reactions to a minimum, while maintaining the essential and most important aspects of chemical reactions taking place in the urban atmosphere on the relevant space scales and timescales.
CityChem offers a more advanced treatment for the photochemistry of multiple gaseous pollutants on the Eulerian grid, as well as for dispersion close to point emission sources (e.g. industrial stacks) and line emission sources (open roads and streets).
1. Photochemistry on the Eulerian grid uses a numerical chemistry solver. The available chemistry schemes include (1) EMEP45 , which resulted from an appropriate reduction of the former EMEP (European Monitoring and Evaluation Programme) chemistry scheme (Simpson, 1995); (2) EmChem03-mod, with updated reaction equations and coefficients compared to EMEP45; (3) and EmChem09-mod, which is similar to the current EMEP chemistry mechanism (Em-Chem09; Simpson et al., 2012). EmChem09-mod enables the simulation of biogenic VOCs, such as isoprene and monoterpenes, emitted from urban vegetation.
2. Modifications of the photochemistry in the sub-grid components have replaced the PSS assumption with the EP10-Plume scheme, a compact scheme including inorganic reactions and the photochemical degradation of formaldehyde, using a numerical solver.
3. Modifications of the line source emission model have been made to compute receptor point concentrations in street canyons. A simplified street canyon model (SSCM) is implemented to account for pollutant transfer along streets, including a parameterization of mass transfer within a simplified building geometry at street level.
4. Modifications to the plume rise from elevated point sources allow for a more accurate computation of the plume trajectories. The Meteorological Pre-Processor (WMPP) of the Weak-wind Open Road Model (WORM) is utilized in the CityChem extension to calculate the wind speed at plume height.
Although computational fluid dynamics (CFD) models can be used to solve for local-scale phenomena along point and line emission sources, they are limited to localized applications and are not appropriate for the simulation of dispersion across complex urban areas. In addition, the simulation of the chemical conversions of reactive pollutants using CFD models requires a large amount of computational time (Sanchez et al., 2016).
The EPISODE-CityChem model, which is based on the core of the EPISODE model, integrates the CityChem extension into an urban CTM system. This paper gives a model description of EPISODE-CityChem version 1.2. In the typical setup, EPISODE-CityChem uses downscaled meteorological fields generated by the meteorological component of the coupled meteorology-chemistry model TAPM (The Air Pollution Model; Hurley, 2008;Hurley et al., 2005). TAPM is a prognostic model which uses the complete equations governing the behaviour of the atmosphere and the dispersion of air pollutants. EPISODE-CityChem is coupled offline with the regional-scale air quality model CMAQ (Community Multiscale Air Quality; Byun et al., 1999;Byun and Schere, 2006;Appel et al., 2013) using hourly varying pollutant concentrations at the lateral and vertical boundaries from CMAQ as initial and boundary concentrations.
EPISODE-CityChem has the capability to simulate the photochemical transformation of multiple reactive pollutants along with atmospheric diffusion to produce concentration fields for the entire city on a horizontal resolution of 100 m or even finer and a vertical resolution of 24 layers up to 4000 m of height. The possibility to get a complete picture of the urban area with respect to reactive pollutant concentrations, but also information enabling exposure calculations in highly populated areas close to road traffic line sources and industrial point sources with high spatial resolution, turns EPISODE-CityChem into a valuable tool for urban air quality studies, health risk assessment, sensitivity analysis of sector-specific emissions, and the assessment of local and regional emission abatement policy options.
The paper is organized as follows: Sect. 2 gives an overview of EPISODE-CityChem and a detailed description of the photochemical reaction schemes and modifications of near-source dispersion in the sub-grid components. Section 3 presents tests of the various modules in the CityChem extension. Section 4 describes the application of EPISODE-CityChem within a nested model chain for simulating the air quality and atmospheric chemistry in the city of Hamburg. We assess the performance of EPISODE-CityChem in reproducing the temporal and spatial variation of air pollutant concentrations against data from urban monitoring stations. Model results from EPISODE-CityChem are compared (1) to results from the standard EPISODE model to quantify the total effect of the new implementations and (2) to results from TAPM, acting as reference model for air pollution modelling on the urban scale. Section 5 outlines plans for the future development of the EPISODE-CityChem model, addressing the need for more sophisticated photochemistry, treatment of aerosol formation on an urban scale and further improvements of the source dispersion. A list of acronyms and abbreviations used in this work is given in Appendix A.
2 Development and description of EPISODE-CityChem model extensions EPISODE consists of a 3-D Eulerian grid CTM that interacts with a sub-grid Gaussian dispersion model for the dispersion of pollutants emitted from both line and point sources. We refer to part one (Hamer et al., 2019) for a technical description of the model. The standard EPISODE model simulates the emission and transport of NO x , as well as fine particulate matter with PM 2.5 (particles with diameter less than 2.5 µm) and PM 10 (particles with diameter less than 10 µm) in urban areas, with the specific aim of predicting concentrations of NO 2 , which is the major pollutant in many cities of northern Europe. EPISODE-CityChem solves the photochemistry of multiple reactive pollutants on the Eulerian grid by using one of the following chemical schemes: (1) EMEP45 chemistry, (2) EmChem03-mod or (3) the EmChem09-mod. In the sub-grid components, the PSS assumption involving O 3 /NO/NO 2 is replaced by the EP10-Plume scheme. Dispersion close to point and line sources is modified in the sub-grid component. In the line source sub-grid model, the simplified street canyon model (SSCM) is integrated to calculate pollutant dispersion in street canyons. In the point source sub-grid model the WMPP (WORM Meteorological Pre-Processor) is integrated to calculate the wind speed at plume height. Figure 1 illustrates the modules and processes of the EPISODE model with the CityChem extension. Modules that belong to the CityChem extension are shown in boxes with a magenta frame.
The recommended configuration of model processes in EPISODE-CityChem is given in Table 1. EPISODE-CityChem has been used with this configuration to simulate the air quality and atmospheric chemistry in the city of Hamburg (Sect. 4).
In EPISODE-CityChem, a regular receptor grid is defined, for which time-dependent surface concentrations of the pollutants at receptor points are calculated by summation of the Eulerian grid concentration of the corresponding grid cell (i.e. the background concentration) and the concentration contributions from the Gaussian sub-grid models due to line source and point source emissions. This way, surface concentration fields of pollutants for the entire city at a horizontal resolution of (currently) 100 m are obtained. The modules of the CityChem extension for photochemistry and source dispersion are described in detail in the remainder of this section.

Extensions to the photochemistry
Atmospheric gas-phase chemical reactions are described by ordinary differential equations (ODEs). The ODE set of reactions is considered stiff because the chemical e-folding lifetimes of individual gases vary by many orders of magnitude in the urban atmosphere (from approx. 10 −6 to 10 6 s −1 ; McRae et al., 1982). The non-linear system of the stiff chemical ODEs is solved by the TWOSTEP solver (Verwer and Simpson, 1995;Verwer et al., 1996) using fast Gauss-Seidel iterative techniques, with numerical error control and restart in the case of detected numerical inaccuracies . The solver is applied to chemical reaction mechanisms available in EPISODE-CityChem for photochemical transformation on the Eulerian grid (EMEP45, EmChem03-mod and EmChem09-mod) and in the sub-grid component (EP10-Plume). For solving the EMEP45 scheme, the Gauss-Seidel iterative technique is used for all compounds except for the oxygen atoms and OH, for which reactions are very fast and we use the steady-state approximation instead . The relative error tolerances for the solver are set to 0.1 (10 % relative error) for all chemical compounds, while the absolute error tolerances are set in a range from 2.5 × 10 8 molecule cm −3 to 1.0 × 10 15 molecule cm −3 depending on the compound. Photodissociation rates are specified as a function of the solar zenith angle and cloud cover, as given in Appendix B. The sink terms for the dry deposition and wet removal of gases and particles are presented in Appendix C. The EMEP45 chemistry scheme developed at NILU  contains 45 chemical compounds and about 70 chemical reactions compared to 70 compounds and about 140 reactions in the original EMEP mechanism (Simpson, 1992(Simpson, , 1993Andersson-Sköld and Simpson, 1999). The intention of the development of EMEP45 was to obtain a condensed chemical scheme for urban areas that still captures the key aspects of the photochemistry in the urban atmosphere. The reduction of the EMEP mechanism was guided by the following considerations: first, the new chemistry scheme is applied in rather polluted urban regions. Second, the residence time of the atmospheric compounds in the urban domain is normally limited to less than a day.
The main simplification in EMEP45 compared to the original EMEP mechanism is the neglect of peroxy radical selfreactions. The self-reactions of peroxy radicals, either between the organic peroxy radical (RO 2 ) and hydroperoxyl radical (HO 2 ) or between two organic peroxy radicals, are in competition with the reaction of RO 2 (or HO 2 ) with NO, leading to photochemical ozone formation: At the ambient levels of NO x typical of moderately or more polluted areas, Reactions (R1)-(R3) will be negligible compared with Reaction (R4). Thus, all reactions of organic peroxy radicals of type (R2) and (R3) were omitted in the EMEP45 scheme. However, due to their relevance, the reaction of HO 2 with the methyl peroxy radical (CH 3 O 2 ) and the HO 2 self-reaction (R1) were included. EMEP45 includes a simple four-reaction scheme for the oxidation of isoprene (C 5 H 8 ) with the OH radical. All reaction rates and coefficients in EMEP45 are according to the International Union of Pure and Applied Chemistry (IUPAC) 2001 recommendations .

Development of the EmChem03-mod scheme and
the EmChem09-mod scheme The EMEP45 scheme was updated in recent years at the Helmholtz-Zentrum Geesthacht (HZG). All reaction rate constants were updated in accordance with the default chemistry scheme EmChem09 of the EMEP/MSC-W model . The resulting scheme is called EmChem03-mod and consists of 45 gas-phase species, 51 thermal reactions and 16 photolysis reactions, as listed in Table S1 in the Supplement. The most important technical change compared to EMEP45 is that the new scheme can be dynamically updated and further extended with new chemical reactions and compounds. The chemical preprocessor of the EMEP/MSC-W model, GenChem, developed at the EMEP group , is used to convert lists of input chemical species and reactions to differential equations of the solver in Fortran 90 code. This makes the update and extension of the new scheme entirely flexible.
In the next step, the EmChem09-mod scheme (Table S2) was developed based on the current EMEP chemistry mechanism, EmChem09 , by (1) replacing the detailed isoprene chemistry with the simplified isoprene reaction scheme from EMEP45, (2) adding monoterpene oxidation reactions and (3) including semi-volatile organic compounds (SVOCs) as reaction products which can potentially act as precursors for secondary organic aerosol (SOA) constituents.
EmChem09-mod includes reactions between organic peroxy radicals and HO 2 as well as other organic peroxy radicals; it is therefore appropriate for low NO x conditions in rural and suburban areas of the city domain. With EmChem09mod the chemistry of biogenic volatile organic compounds (BVOCs), emitted from urban vegetation, can be simulated. Two monoterpenes, α-pinene and limonene, are model surrogates to represent slower-and faster-reacting monoterpenes (α-pinene: 5.32 × 10 −11 cm 3 s molecule −1 ; limonene: 1.7 × 10 −10 cm 3 s molecule −1 ; for the OH-reaction, both at 298 K). The scheme considers the OH-initiated oxidation of isoprene, as well as the oxidation of α-pinene and limonene by OH, NO 3 and O 3 . Limonene has two reactive sites (double bonds) allowing for a rapid reaction chain to oxidation products with low vapour pressure. The lumped reaction scheme of α-pinene is adopted from Bergström et al. (2012) and that of limonene is based on Calvert et al. (2000). In total, EmChem09-mod includes 70 compounds, 67 thermal reactions and 25 photolysis reactions.

Development and description of the EP10-Plume chemistry scheme
In the sub-grid components, i.e. the Gaussian models for line and point source dispersion, the PSS assumption involving O 3 /NO/NO 2 was replaced by the EP10-Plume scheme for computation of the chemistry at the local receptor grid points. EP10-Plume includes only the reactions of O 3 , NO, NO 2 , nitric acid (HNO 3 ) and CO, as well as the photochemical oxidation of formaldehyde (HCHO). It contains 10 compounds and 17 reactions; Table S3 provides a list. Only a small portion of NO x from motor vehicles and combustion sources is in the form of NO 2 , the main part being NO. The largest fraction of ambient NO 2 originates from the subsequent chemical oxidation of NO. The only reactions considered to be relevant in the vicinity of NO x emission sources are For conditions in northern Europe, an instantaneous equilibrium between the three reactions relating NO, NO 2 and O 3 is assumed, the so-called PSS, and implemented in the EPISODE model. In EP10-Plume the three reactions are, however, treated explicitly. Reactions occurring with negligible rates at the NO x levels typical of moderately or highly polluted areas were excluded from the scheme. HCHO and acetaldehyde are important constituents of vehicle exhaust gas (e.g. Rodrigues et al., 2012). The photolysis of HCHO is a source of HO 2 radicals.
HCHO also reacts with the OH radical to give two HO 2 radicals. HO 2 competes with ozone for the available NO (Reaction R4), and the reaction between HO 2 and NO results in additional NO-to-NO 2 conversion. Since the generation of HO 2 radicals through HCHO photolysis does not depend on the entrainment of photo-oxidants from the background air, it can trigger the photochemical reaction cycle even in traffic plumes very close to the source. Carbon monoxide (CO) has a lifetime of about 2 months towards OH (at [OH] = 1.2 × 10 6 molecules cm −3 ). Reaction (R10) is therefore not relevant near sources and of very low relevance on the urban scale. For completeness of the OH-to-HO 2 cycling, Reaction (R10) was, however, included in EP10-Plume.

Extensions to the source dispersion
Sub-grid models to resolve dispersion close to point sources and line sources are embedded in the EPISODE model to account for sub-grid variations as a result of emissions along open roads and streets as well as along plume trajectories from elevated point source releases. The sub-grid model for line sources, i.e. open road and urban street traffic, is the Gaussian model HIWAY-2 (Highway Air Pollution Model 2; Petersen, 1980) from the U.S. EPA with modifications. The sub-grid model for point sources, e.g. stacks of industrial plants and power plants, is the Gaussian segmented plume trajectory model SEGPLU (Walker and Grønskei, 1992). SEGPLU computes and keeps a record of subsequent positions of plume segments released from a point source and the corresponding pollutant concentration within each plume segment. The vertical position of the plume segment is calculated from the plume rise of the respective point source. Plume rise for elevated point sources due to momentum or buoyancy is computed based on the plume rise equations originally presented by Briggs (1969Briggs ( , 1971Briggs ( , 1975. A detailed description of the implementation of HIWAY-2 and SEGPLU in the EPISODE model is given in part one (Hamer et al., 2019). In this section, extensions of the sub-grid models for the simulation of dispersion near sources within City-Chem are described.

Implementation of a simplified street canyon model (SSCM) for line source dispersion
In CityChem, a simplified street canyon model (SSCM) to compute concentrations for receptor points that are located in street canyons is introduced. The street canyon model follows in most aspects the Operational Street Pollution Model (OSPM; Berkowicz et al., 1997). A fundamental assumption of this model is that when the wind blows over a rooftop in a street canyon, an hourly averaged recirculation vortex is always formed inside the canyon (Hertel and Berkowicz, 1989). The part of the street canyon covered by the vortex of recirculating air is called the recirculation zone. The concentration at a receptor point located within an urban street canyon is calculated as the sum of the concentration contribution (C line,s ) due to the emissions of the line source s and the urban background concentration, which is taken from the corresponding cell of the Eulerian grid component. The contribution of a line source s is given by the direct contribution (C scdir,s ) from the traffic plume plus a contribution from the recirculation of the traffic plume (C screc,s ) due to the vortex inside the canyon (Berkowicz et al., 1997): The leeward receptor inside a street canyon is exposed to direct contribution from the emissions inside the recirculation zone (unless the wind direction is close to parallel) and a recirculation contribution. For the receptor on the windward side, only emissions outside the recirculation zone are considered for the direct contribution. If the recirculation zone extends through the whole canyon, no direct contribution is given to the windward receptor. The length of the recirculation zone, L rec , is estimated as being twice the average building height of the canyon and limited by the canyon width, W sc . The calculation of the direct and recirculation concentration contributions in this simple approach is adopted from the OSPM following the description in Berkowicz et al. (1997) with certain modifications. Simplifications are made with respect to the street canyon geometry, since only general geometries with average street canyon width and height are used. The rate of release Q s in the street assumes that emissions are distributed homogeneously along the line source segment that is inside the street canyon area, which means emissions are assumed to be distributed homogeneously over the street canyon in the full length and width of the canyon (along the dimension of the respective line source object).
The direct contribution is calculated using a Gaussian plume model. The direct concentration contribution at the receptor point C scdir,s , located at distance x from the line source (i.e. starting from the midline of the street), is obtained by integrating along the wind path at street level. The integration path depends on wind direction, the extension of the recirculation zone and the street canyon length (Hertel and Berkowicz, 1989): where h 0 is a constant that accounts for the height of the initial pollutant dispersion (h 0 = 2 m is used in SSCM), σ w is the vertical velocity fluctuation due to mechanical turbulence generated by wind and vehicle traffic in the street, and u street is the wind speed at street level, calculated assuming a logarithmic reduction of the wind speed at rooftop towards the bottom of the street. Note that the wind direction at street level in the recirculation zone is mirrored compared to the roof-level wind direction. Outside the recirculation zone, the wind direction is the same as at roof level. The vertical velocity fluctuation is calculated as a function of the street-level wind speed and the traffic-produced turbulence by the following relationship (Berkowicz et al., 1997): where α s is a proportionality constant empirically assigned a value of 0.1, and σ w0 is the traffic-induced turbulence, in SSCM assigned a value of 0.25 m s −1 , which is typical for traffic on working days between 08:00 and 19:00 (Central European Time) in situations in which traffic-induced turbulence dominates (Kastner-Klein et al., 2000;Fig. 6 therein). The integration path for Eq. (2) begins from x start , which is defined as the distance from the receptor point at which the plume has the same height as the receptor, which is zero in the case that h 0 is smaller than or equal to the height of the receptor. The upper integration limit is x end , defined by tabular values in Ottosen et al. (2015 , Table 3 therein). The integration is performed along a straight line path against the wind direction. The calculation of the maximum integration path, L max , depends on the wind direction with respect to the street axis, θ street , i.e. the angle between the street and the street-level wind direction (Ottosen et al., 2015).
The recirculation contribution is computed using a simple box model, assuming equality of the inflow and outflow of the pollutant. The cross section of the recirculation zone is modelled as a trapezium with upper length L top and baseline length L base . L top is half of the baseline length, where L base is defined as min(L rec , L max ). The length of the hypotenuse of the trapezium is calculated as L hyp = (L base /2) 2 + H 2 sc , assuming the leeward side edge of the recirculation zone to be the vertical building wall, with the length of the building height. It is further assumed that the slant edge of the recirculation zone towards the opposite street side is not intercepted by buildings.
The recirculation concentration contribution is expressed by the relationship (Berkowicz et al., 1997) where σ wt is the ventilation velocity of the canyon as given by Hertel and Berkowicz (1989), and σ hyp is the average turbulence at the hypotenuse of the trapezium (slant edge towards the opposite street side). For a given receptor point, the concentration contribution from a line source is calculated either by HIWAY-2 or by SSCM. HIWAY-2 does not calculate line source concentration contributions to receptors that are upwind of a line source or receptor points that are very close to the line source. For all windward and leeward receptor points (1) located within a model grid cell defined as a street canyon cell (see below), (2) located close enough to a line source (i.e. within the actual street canyon) and (3) located at a road link with length > 8 m, the concentration contribution from the street is calculated by SSCM. For all windward receptors which do not fulfil these conditions, the concentration contribution is calculated by HIWAY-2.
The complex and diverse geometry of street canyons is approximated by three generic types for which average street canyon geometry properties are applied (Table 2). Street canyons are identified based on the urban land use classes of TAPM. Each line source for which the geometric midpoint is located in a grid cell with urban land use (land use classes 32-35 defined in TAPM) is identified as a potential street canyon. A disadvantage of this method is that some streets and roads, especially in the sparsely built urban areas outside the inner city, will be classified as street canyons despite being open roads with open spaces between buildings.
Furthermore, it is assumed that all buildings at the street canyon line source have the same average building height, H sc , and that there are no gaps between the buildings. The average building heights for the TAPM land use classes were obtained by the intersection of the 3-D city model LoD1-DE Hamburg (LGV, 2014) -which contains individual building heights -with the CORINE (Coordination of Information on the Environment) urban land use information (CLC, 2012). The width of the street canyon, W sc , is defined as twice the width of the (line source) street width W to account for sidewalks and to avoid canyons that are too narrow. The length of the street canyon, L sc , corresponds to the length of the line source within the grid cell.

Implementation of the WMPP for point sources
The wind speed profile function of the meteorological preprocessor WMPP is utilized in the CityChem extension to calculate the wind speed at plume height within the point source sub-grid dispersion model. WMPP replaces the previous routine, which calculated the wind speed at plume height using a logarithmic wind speed profile corrected by the stability function for momentum based on Holtslag and de Bruin (1998). WMPP has been developed as part of NILU's WORM open-road line source model (Walker, 2011(Walker, , 2010 to calculate various meteorological parameters needed by WORM. In the current version of WORM, the profile method is applied using hourly observations of wind speed at one height, e.g. 10 m, and the temperature difference between two heights, e.g. 10 and 2 m, to calculate the other derived meteorological parameters. Given the above input data and an estimate of the momentum surface roughness, WMPP calculates friction velocity (u * ), temperature scale (θ * ) and inverse Obukhov length (L −1 ) according to Monin-Obukhov similarity theory. These quantities are calculated by solving the following three nonlinear equations: where κ is Von Kármán's constant (0.41), g is the acceleration of gravity (9.81 m s −2 ), u is the wind speed difference between heights z u2 and z u1 , where z u2 is e.g. 10 m, and z u1 = z 0m , where the wind speed is zero, so that u = u 10 m − 0 = u 10 m . In the definition of the temperature scale, θ is the difference in potential temperature between heights z t2 and z t1 , which are e.g. 10 and 2 m, respectively, so that we have θ = T 10 m − T 2 m + 0.01, where the +0.01 term is for the conversion from potential temperature to actual temperature. In the definition of the Obukhov length, T ref is a reference temperature, here taken to be the average of T 2 m and T 10 m .
In Eq. (5), the similarity functions ϕ m and ϕ h are defined as follows (Högström, 1996): and Geosci. Model Dev., 12, 3357-3399, 2019 www.geosci-model-dev.net/12/3357/2019/ where Pr 0 = 0.95 is the Prandtl number for neutral conditions and where the empirical coefficients are defined as α m = −19.0, α h = −11.6, β m = 5.3 and β h = 8.2. This set of similarity functions is then used to calculate vertical profiles of temperature and wind speed. The temperature at a height (in metres above the ground) is thus calculated by (8) where z ref = 10 m and c p is the specific heat capacity of air, here set to 1005 J kg −1 K −1 . Similarly, the wind speed at height z (m) above the ground is calculated by In CityChem, WMPP is used in the sub-grid point source model to calculate the wind speed at plume height according to Eq. (9). WMPP can also be used to calculate the convective velocity scale w * and the mixing height h mix , but this is not implemented in CityChem.

Additional modifications
Here we describe the modifications in the CityChem extension to read hourly 3-D boundary concentrations from the output of the CMAQ model and to determine sub-grid concentrations from a regular receptor grid in the surface model layer.

Adapting 3-D boundary conditions from the CMAQ model
CityChem has the option to use the time-varying 3-D concentration field at the lateral and vertical boundaries from the CMAQ model as initial and boundary concentrations for selected chemical species. The adaption of boundary conditions from CMAQ output in the EPISODE model is based on the implementation for boundary conditions from the Copernicus Atmosphere Monitoring Service (CAMS; http://www.regional.atmosphere.copernicus.eu/, last access: 29 July 2019) described in part one (Hamer et al., 2019). The regional background concentrations are adopted for the grid cells (outside the computational domain) directly adjacent to the boundary grid cells of the model domain and for the vertical model layer that is on top of the highest model layer.
The outside grid cell directly adjacent to the boundary grid cell is filled with the CMAQ concentration value for inflow conditions and with the concentration value of the boundary grid cell for outflow conditions, i.e. allowing for a zeroconcentration gradient at the outflow boundary. More details on the treatment of 3-D boundary conditions are given in Appendix D.

Description of the regular receptor grid
In the CityChem extension, a regular receptor grid is defined, for which time-dependent surface concentrations of the pollutants at receptor points are calculated by summation of the Eulerian grid concentration of the corresponding grid cell (i.e. the background concentration) and the concentration contributions from the sub-grid models due to the dispersion of line source and point source emissions. Regular receptor grids with a typical resolution of 100 × 100 m 2 have also been used in earlier versions of EPISODE, but primarily for capturing sub-grid-scale concentration contributions from larger industrial point sources. The establishment of a regular receptor grid is an integral part of CityChem to enable the higher-resolution output required for comparison with monitor data acquired near line sources. Line sources are a major source of pollutant emissions affecting inner-city air quality; thus, the use of the regular receptor grid provides information at much higher spatial resolution than the Eulerian grid output alone. The regular receptor grid in EPISODE-CityChem differs from the downscaling approach by Denby et al. (2014), which allocates sampling points at high density along roads and other line sources but much fewer further away from the line sources. While Denby et al. (2014) interpolate the model-computed high-density set of receptor concentrations to the desired output resolution using ordinary kriging, EPISODE-CityChem gives as output the receptor point concentrations on a regular 2-D grid covering the entire model domain.
The instantaneous concentration C rec in an individual receptor point r * of the receptor grid with coordinates (x r , y r , z r ) is defined as where C m is the main grid concentration of the grid cell (x, y, 1) in which the receptor point is located. The grid (background) concentration C m used in Eq. (10) corresponds to a modified Eulerian 3-D grid concentration, i.e. C(x, y, z), to prevent emissions of point and lines sources from being counted twice. C point,p is the instantaneous concentration contribution of point source p calculated by the point source sub-grid model, and C line,s is the instantaneous concentration contribution of line source s calculated by the line source sub-grid model. Since C rec is not added to the main grid concentration but kept as a separate (diagnostic) variable, the double-counting of emitted pollutant mass is prevented. In the CityChem extension, receptor point concentrations represent the high-resolution ground concentration of a cell with the grid cell area of the receptor grid.
On the 3-D Eulerian grid, time-dependent concentration fields of the pollutants are calculated by solving the advection-diffusion equation with terms for chemical reactions, dry and wet deposition, and area emissions. The hourly 2-D and 3-D fields of meteorological variables and the hourly 2-D fields of area emissions are given as input to the model with the spatial resolution of the Eulerian grid. As the model steps forward in time, an accurate account of the total pollutant mass from the area, point and line sources is kept within the Eulerian grid model component. Emissions from line sources are added to the Eulerian grid concentrations at each model time step.

Test of different model extensions
For the test of the various model extensions, EPISODE was run as a 1-D column model, with vertical exchange as the only transport process. Emissions were injected into the ground cell (grid centre at UTM coordinates: (X) 568500, (Y ) 5935550, 32 N) with an area of 1 × 1 km 2 and flat terrain (15 m a.s.l.). Table 3 shows the general setup for the 1-D column and the specific configuration for the tests. Mixing height, surface roughness and friction velocity were kept constant (h mix = 250m, z 0 = 0.8 m s −1 , u * = 0.12 m s −1 ). Hourly varying meteorological variables included air temperature, temperature gradient, relative humidity, sensible and latent heat fluxes, total solar radiation, and cloud fraction. The test simulations are performed for a period of 5 d, and results were taken as an average of the period.

Tests of the original EMEP45 photochemistry
When the condensed EMEP45 photochemistry was developed, various tests were carried out to compare the condensed mechanism with the standard EMEP chemical mechanism. Results from box model studies with the two chemical mechanisms revealed that there were generally small differences between the full and the condensed chemical mechanisms. Even for conditions more representative of a rural environment, the difference between the standard EMEP and the condensed mechanism was small. For these more rural conditions, the condensed mechanism gave slightly lower levels of NO and NO 2 , while the ozone concentration was almost identical in the two mechanisms. For urban conditions, these differences were expected to be significantly smaller.
The EPISODE model with the condensed EMEP45 mechanism furthermore participated in the CityDelta project (Cu-velier et al., 2007) within which it was applied to the city of Berlin. CityDelta was the first in a series of projects (later named EuroDelta) dedicated to photochemical model intercomparisons. When evaluated against observations of NO 2 and O 3 , the EPISODE model with the EMEP45 chemistry performed favourably when compared to the suite of atmospheric models participating in the CityDelta project .

Test of ozone formation with EmChem03-mod
The ozone-NO x -VOC sensitivity of the EmChem03-mod scheme in the Eulerian model component was analysed by repeated runs with varying emissions of NO x and nonmethane VOCs (NMVOCs) using the daily cycle of mean summer meteorology with clear sky but low wind speed (0.1 m s −1 ). The ozone net production in the runs was taken at the maximum daily O 3 during the simulation.
An area source of traffic emissions of NO x and NMVOCs in the ground cell of the 1-D column was activated in the test. The variation of ozone precursor emissions from the traffic area source was done in a systematic way in order to derive the ozone isopleth diagram (Fig. 2a), which shows the rate of O 3 production (ppb h −1 ) as a function of NO x and NMVOC concentrations. Compound abundances are given in mixing ratios (ppb) for this test to enable comparison with the literature on ozone formation potentials.
The ozone-precursor relationship in urban environments is a consequence of the fundamental division into NO xlimited and VOC-limited chemical regimes. VOC/NO x ratios are an important controlling factor for this division of chemical regimes (Sillman, 1999). VOC-limited chemistry generally occurs in urban centres where NO 2 concentrations are high due to traffic emissions. Rural areas downwind of the city are typically NO x limited (Ehlers et al., 2016).
The "ridgeline" of the ozone isopleth diagram marks the local maxima of O 3 production and differentiates two different photochemical regimes. Below the line is the NO xlimited regime, in which O 3 increases with increasing NO x , while it is hardly affected by increasing VOCs. Above the line is the VOC-limited regime, in which O 3 increases with increasing VOCs and decreases with increasing NO x . The ridgeline in Fig. 2a follows a line of constant VOC/NO x ratio; in the case of EmChem03-mod it is close to the ratio 10 : 1, whereas a slope of 8 : 1 is more typically found (e.g. Dodge, 1977). The traffic NMVOC mixture includes a high share of aromatics (35 %) represented by o-xylene in the model. Due to the high reactivity of the NMVOC mixture, the ridgeline is tilted towards higher VOC/NO x ratios compared to the ozone isopleths for a NMVOC mixture with lower reactivity.
The split into NO x -limited and VOC-limited regimes is closely associated with sources and sinks of odd hydrogen radicals (defined as the sum of OH, HO 2 and RO 2 ). Odd hydrogen radicals are produced in the photolysis of ozone and intermediate organics such as formaldehyde. Odd hydrogen radicals are removed by reactions that produce hydrogen peroxide (Reaction R1) and organic peroxides (Reaction R2). They are also removed by reaction with NO 2 , producing HNO 3 , according to When peroxides represent the dominant sink for odd hydrogen, then the sum of peroxy radicals is insensitive to changes in NO x or VOC. This is the case for the concentrations represented as solid and dash-dotted lines in Fig. 2cd. Doubling NO x emissions from solid lines to dash-dotted lines only marginally changes the peroxy radical sum concentration (Fig. 2d).
When HNO 3 is the dominant sink of odd hydrogen, then the OH concentration is determined by equilibrium between the producing reactions (e.g. photolysis of O 3 ) and the loss reaction (R11); it thus decreases with increasing NO x (Fig. 2c- Plotting the isopleths for the ratio of the production rate of peroxides to the production rate of HNO 3 (Fig. 2b) shows that this ratio is closely related to the split between NO xlimited and VOC-limited regimes. The ratio is typically 0.9 or higher for NO x -limited conditions and 0.1 or less for VOC-limited conditions (Sillman, 1999). The ridgeline that separates the two regimes should be at a ratio of 0.5 (Sillman, 1999), which is the case in Fig. 2b. However, the curves representing the ratio are shifted towards higher NO x mixing ratios compared to the isopleth diagram for the ratio displayed in Sillman (1999, Fig. 8 therein). For instance, for 100 ppbC NMVOCs and 5 ppb NO x , the ratio is below 0.1 (VOC limited) in the isopleth diagram of Sillman (1999), while it is 1.3 (NO x limited) in Fig. 2b. The reason for this discrepancy is the lack of reactions producing organic peroxides (RO 2 H) in EmChem03-mod and thus the reduced removal of odd hydrogen in conditions with a high VOC/NO x ratio. In conditions with NO x below 20 ppbv, EmChem03-mod has an efficiency of NO-to-NO 2 conversion via Reaction (R4) that is too high.

Test of EmChem09-mod photochemistry
The EmChem09-mod scheme was compared to the EmChem03-mod scheme for conditions with relatively low levels of NO x (< 20 µg m −3 ). The configuration of the test was the same as in Sect. 3.1.2, with an area source of traffic emissions of NO x (0.043 g s −1 in the 1 × 1 km 2 ground cell) and varying emissions of NMVOC corresponding to VOC/NO x ratios of 4 : 1, 8 : 1 and 15 : 1. The daily cycle of ozone with EmChem09-mod shows O 3 concentrations which are lower for a VOC/NO x ratio of 4 : 1 (VOC limited) than with EmChem03-mod, similar for a VOC/NO x ratio of 8 : 1 (transition) to EmChem03-mod and higher for a VOC/NO x ratio of 15 : 1 (NO x limited) than with EmChem03-mod showing the ratio of the production rate of peroxides to the production rate of nitric acid, (c) concentration time series of O 3 (black), NO x (red) and OH (grey; second y axis), and (d) concentration time series of H 2 O 2 (blue), HNO 3 (green) and HO 2 +RO 2 concentration (grey, second y axis). Daily concentration cycle as an average from a test run with NMVOC emissions of 695×10 −8 g s −1 m −2 and varying NO x emissions: 1 × 10 −8 g s −1 m −2 (solid lines), 2 × 10 −8 g s −1 m −2 (dash-dotted lines), 38 × 10 −8 g s −1 m −2 (dashed lines) and 55 × 10 −8 g s −1 m −2 (dotted lines). Lines of constant VOC/NO x ratio are annotated with red dashed lines (4 : 1, 8 : 1 and 15 : 1) and the blue dashed line (10 : 1) in panel (a). Note the logarithmic scale of the y axis in panel (b).
( Fig. 3a). Compared to EmChem03-mod, the EmChem09mod scheme includes reactions between organic peroxy radicals and HO 2 , as well as other organic peroxy radicals. In conditions with low levels of NO x , the rates from these reactions will be in competition with the reaction rates of organic peroxy radicals with NO.
The lower O 3 with EmChem09-mod in VOC-limited conditions is related to the competition between organic peroxy radical self-reactions and the reaction with NO, preventing additional NO-to-NO 2 conversion. Compared to EmChem03-mod, the removal of odd hydrogen through Reaction (R11) to form HNO 3 is weakened (Fig. 3b), the formation of H 2 O 2 and organic peroxides is enhanced (Fig. 3c), and the formation of peroxyacetyl nitrate (PAN) is suppressed (Fig. 3d); the latter is due to the competing reaction between the acetyl peroxy radical (CH 3 COO 2 ) and HO 2 , which is not included in EmChem03-mod. As a result, less NO 2 is lost and the NO x concentrations in EmChem09-mod increase compared to EmChem03-mod (Fig. S1), which reduces ozone production in the VOC-limited regime.
The higher O 3 with EmChem09-mod in NO x -limited conditions is related to the much higher production of peroxides and the reduced production of PAN and HNO 3 compared to EmChem03-mod. The NO x concentrations in EmChem09mod are higher, which increases ozone production in the NO x -limited regime.

Test of EP10-Plume sub-grid photochemistry
The photochemistry in the sub-grid component of EPISODE-CityChem was tested for dispersion from a single line source aligned in the SE-NW diagonal of the 1 × 1 km 2 grid cell. The line source was oriented perpendicular to the wind direction, emitting NO x and NMVOCs with a ratio of 1 : 2. The HIWAY-2 line source model was used in the test (SSCM was not activated). Photochemistry tests were made as follows: (1) no chemistry; (2) photochemical steady-state assumption (PSS) for O 3 /NO/NO 2 (default); and (3) with EP10-Plume using the numerical solver. Inside the centre cell, ground air concentrations downwind of the line source were recorded using additional receptor points every 10 m up to a distance of 300 m from the line source.
Comparing O 3 (black lines), NO 2 (red lines) and NO (blue lines) concentrations from the three tests with increasing downwind distance x shows that dilution alone (test with no chemistry; Fig. 4a) leads to a decay of NO, which follows a power function of the form y = ax −b , while O 3 remains constant at the level of the background concentration (30 µg m −3 ).
Applying the PSS reduces O 3 immediately at the line source by reaction (R5) to one-fourth of the concentration without chemistry. At the line source (0 m of distance), PSS converts roughly 15 µg m −3 NO to 21 µg m −3 NO 2 , as deduced from the differences between the no chemistry and the PSS test run. The third option, EP10-Plume, gives very similar results to PSS, with O 3 , NO and NO 2 concentrations deviating by at most 4 % from the solution of the PSS (overlapping lines in Fig. 4b). In EP10-Plume, the line-emitted HCHO during daytime reacts with OH or undergoes photolysis to give HO 2 radicals. However, the odd hydrogen radicals are rapidly removed by Reaction (R11), and the effect of emitted HCHO on O 3 is negligible. It is noted that HCHO accounts for only 2.7 % of traffic NMVOC emissions. Further testing showed that the share of HCHO has to be increased by a factor of 10 or more (for the same VOC emission rate) in order to exceed the PSS concentration of O 3 close to the line source.

Test of SSCM for line source dispersion
Tests with the simplified street canyon model (SSCM; see Sect. 2.2.1) were performed for different roof-level wind speeds (0.5, 1.0, 1.5, 2.0, 4.0 and 6.0 m s −1 ) and compared to results from the HIWAY-2 line source dispersion model. The street canyon was oriented along the SE-NW diagonal of the grid cell, canyon width was 18 m and average building height was 18 m, with no gaps between buildings. Receptor points were placed symmetrically on the northeast side and the southwest side of the canyon 5 m of distance from the street. Time-averaged modelled concentrations of PM 10 , emitted from the line source as a chemically inert tracer, are shown as a function of wind direction and wind speed in Fig. 5 for the northeast side (left) and southwest side (right) receptor. The wind direction dependency at the two receptors is simply shifted by 180 • with respect to the other due to the symmetric arrangement. With SSCM, the leeward concentrations are generally higher than the windward concentrations (grey-shaded areas in the figure). For both models, maximum concentrations are calculated for wind direction close to parallel with the street (135 and 315 • ).
For this specific street canyon, with an aspect ratio (W sc /H sc ) equal to 1, the recirculation zone extends through the whole canyon at high wind speeds and the windward receptor only receives a contribution from the recirculation. At low wind speeds, here at 2 m s −1 or below, the windward side starts to receive a direct contribution because the exten-sion of the vortex decreases at low wind speeds. At wind speeds below 0.5 m s −1 , the vortex disappears and trafficgenerated turbulence determines the concentration levels. Gaussian models are not designed to simulate dispersion in low-wind conditions. Therefore, a lower limit of the rooftop wind speed was placed at 0.5 m s −1 in this test, preventing the test of lower wind speeds. It is, however, obvious from Fig. 5a that the influence of the wind direction on concentrations at 0.5 m s −1 is much reduced compared to higher wind speeds.
Similar to SSCM, the simulation with HIWAY-2 shows a local maximum at the windward side when the wind is perpendicular to the street and a local minimum at the leeward side when the wind is perpendicular. In HIWAY-2, the pollution from traffic is dispersed freely away from the street because it applies to open road without buildings. In SSCM, the leeward side is influenced directly by the traffic emissions in the street and additionally by the recirculated polluted air. HIWAY-2 neglects the contribution of recirculated polluted air. This is also the reason why the baseline contribution (in addition to the urban background) is higher in SSCM.

Test of WMPP-based point source dispersion
The WMPP model code was extensively tested using meteorological observations from a 4-month measurement campaign at Nordbysletta in Lørenskog, Norway, in 2002 (Walker, 2011(Walker, , 2010. WMPP (see Sect. 2.2.2) is used in the plume rise module of SEGPLU (Walker and Grønskei, 1992) for the calculation of the wind speed at (1) stack height, (2) plume heights along the plume trajectory and (3) at final plume height. The modification of the plume rise module is similar to the "NILU plume" parameterization implemented in the WRF-EMEP (Weather Research and Forecasting-European Monitoring and Evaluation) model, as presented in Karl et al. (2014). In comparison with two simple schemes for plume rise calculation, NILU plume gave a lower final plume rise from an elevated point source for all tested atmospheric stability conditions. In neutral conditions, the maximum concentration at the ground (C max ) was found to be roughly twice as high as for the two simple plume rise schemes. In unstable conditions, all plume rise schemes gave similar effective emission heights.
The WMPP integration in the sub-grid point source model for near-field dispersion around a point source was tested in different atmospheric stability conditions and compared to the standard point source parameterization in EPISODE (termed "default" in the following). A single point source was located at the midpoint of the 1 × 1 km 2 grid cell. The dispersion of SO 2 , treated as a non-reactive tracer, released from the point source stack was studied by sampling ground air concentrations from a regular receptor grid with 100 m resolution within a radius of 2 km around the point source. Transport on the Eulerian grid was deactivated so that the test corresponds to a stand-alone test of the Gaussian point source model. Details about the point source and resulting hourly ground concentrations (averaged for 5 d) at the location of maximum impact, C max , for different stability condi- Table 4. Test of the point source dispersion of SO 2 (handled as an inert tracer) for different atmospheric stability conditions in flat terrain at a wind speed of 1 m s −1 . Hourly concentration is given at the location where the maximum is found for the 5 d average within a radius of 2.0 km around the point source. Parameterization of point source: exhaust gas temperature: 20 • C; stack height: 10 m; exit velocity: 5 m s −1 ; stack radius: 0.5 m (circular opening). tions (slightly stable, neutral, slightly unstable, very unstable) are summarized in Table 4.
For WMPP, the maximum impact lies within 250 and 400 m downwind of the point source in neutral and unstable conditions. The effective emission height, H eff , in neutral conditions is about half of that computed by the default parameterization. For WMPP, H eff decreases with enhanced instability (from neutral to very unstable), while C max increases correspondingly. The increase in C max computed by WMPP is 41 %. C max should be roughly inversely proportional to the square of the effective emission height (Hanna et al., 1982); thus, the decrease from 41 m (neutral) to 32 m (very unstable) implies a potential increase in C max by 25 %. The higher increase in C max than expected might be due to continuous wind from one direction (225 • ) and the relatively low wind speed (1 m s −1 ) in the test. For the default, H eff and C max are not affected by changing stability in neutral or unstable conditions; computed C max is a factor of 4-6 smaller than for WMPP. In stable conditions, C max is several orders of magnitude smaller than in neutral and unstable conditions for both plume rise parameterizations. The maximum impact is found at 1700 m of distance from the point source, which is comparable to previous tests with the point source model (Karl et al., 2014).  (Rockel et al., 2008) for the year 2012 using the ERA-Interim reanalysis of the European Centre for Medium-Range Weather Forecasts (ECMWF) as forcing data (Geyer, 2014). Within the northern Europe domain, an inner domain over the Baltic Sea region with 4 km horizontal resolution was nested (Fig. 6a). The 4 km resolved CMAQ simulation of the Baltic Sea region provided the initial and hourly boundary conditions for the chemical concentrations in the Hamburg model domain.
The hourly meteorological fields for the study domain of Hamburg (30 × 30 km 2 ) were obtained from the inner domain in a nested simulation with TAPM (Hurley et al., 2005) with a 1 km horizontal resolution (D4 in Fig. 6b). The meteorological component of TAPM is an incompressible, non-hydrostatic, primitive equation model with a terrainfollowing vertical sigma coordinate for 3-D simulations; it was used for downscaling the synoptic-scale meteorology. The outer domain (D1 in Fig. 6b) was driven by the 3hourly synoptic-scale ERA5 reanalysis ensemble means on a longitude-latitude grid at 0.3 • grid spacing. In addition, wind speed and direction observations at seven measurement stations of the German Weather Service (DWD) are used to nudge the predicted solution towards the observations. TAPM uses a vegetative canopy, soil scheme and an urban scheme with seven urban land use classes (Hurley, 2008) at the surface, while radiative fluxes, both at the surface and at upper levels, are also included. In regions belonging to one of the urban land use classes, specific urban land use characteristics, such as fraction of urban cover, albedo of urban surfaces, thermal conductivity of urban surface materials (e.g. concrete, asphalt, roofs), urban anthropogenic heat flux and urban roughness length, are used to calculate the surface temperature and specific humidity as well as surface fluxes and turbulence. A complete list of the meteorological variables and fields used from TAPM as input to EPISODE-CityChem for the AQ simulations is given in the user's guide for EPISODE-CityChem, which is included in the CityChem distribution.
For a better representation of local features, the coarsely resolved standard land cover classes and elevation heights, which are provided together with TAPM, were updated with 100 m CORINE Land Cover 2012 data (CLC, 2012), and terrain elevation data were adopted from the German Digital Elevation Model (BKG, 2013) at 200 m horizontal resolution.
The procedure to adapt hourly 3-D concentrations of the CMAQ model computed for the 4 km resolution domain as lateral and vertical boundary conditions is described in Sect. 2.3.1 and Appendix D. CMAQ concentrations from the tion (in UTM projection) of the Hamburg study domain, preserving a nesting factor of four (64, 16, 4, 1 km) for the nested model chain. The study domain is located in the southwest part of the 4 km CMAQ domain (CD04); the west border of the study domain is 30 km from the CD04 west border, and the south border is 21 km from the CD04 south border (inset in Fig. 6a). The contribution of the recirculation of NO x from the coarser outer domain to the budget of NO x in the inner domain is very small due to the predominant westerly winds.

Description of the model setup and configuration for Hamburg
The EPISODE-CityChem simulation was performed using the recommended numerical schemes for physics and chemistry, including the new urban parametrization for vertical eddy diffusivity (urban K(z), see part one; Hamer et al., 2019). The segmented plume model SEGPLU with WMPPbased plume rise was used for the point source emissions. The line source model HIWAY-2 with the street canyon option was used for the line source emissions.  (Matthias et al., 2018). Table 6 gives an overview. Spatially gridded annual emission totals for area sources with a grid resolution of 1 × 1 km 2 were provided by the German Federal Environmental Agency (Umweltbundesamt, UBA). The spatial distribution of the reported annual emission totals was calculated at UBA using the ArcGISbased software GRETA (Gridding Emission Tool for Ar-cGIS) (Schneider et al., 2016). Hourly area emissions with a 1 km horizontal resolution for SNAP cat. 03 (commercial combustion), 06 (solvent and other product use), 08 (other mobile sources, not including shipping) and 10 (agriculture and farming) were derived from the UBA area emissions by temporal disaggregation using monthly, weekly and hourly profiles.
For SNAP cat. 02 (domestic heating) the daily average ground air temperature obtained from the TAPM simulation is used to create the annual temporal profiles. The day-today variation of domestic heating emissions is based on the heating degree-day concept implemented in the Urban Emission Conversion Tool (UECT) (Hamer et al., 2019). Domestic heating emissions (SNAP cat. 02) for Hamburg are distributed between 32 % district heating, 40 % natural gas, 14 % fuel oil and 14 % electricity (Schneider et al., 2016).
NMVOC emissions in the UBA dataset were distributed over individual VOCs of the chemical mechanism using the VOC split of the EMEP model  for all SNAP sectors.
A total of 120 point sources were extracted from the PRTR (Pollutant Release and Transfer Register) database (PRTR, 2017) and from the registry of emission data for point sources in Hamburg, representing the largest individual stack emissions.
The line source emission dataset (emissions of NO x , NO 2 and PM 10 ) provided by the city of Hamburg contained 15 816 road links within the study domain. The NO x emission factor from road traffic for the year 2012 was increased by 20 % for all street types because the average NO x emission factor in the new HBEFA (Handbook Emission Factors for Road Traffic) v3.3 for passenger cars is higher by 19.4 % (diesel cars: 21 %) than in HBEFA v3.1 used in the road emission inventory (UBA, 2010). To estimate NMVOC traffic emissions, an average NMVOC/NO x ratio of 0.588, derived from UBA data for SNAP cat. 07, was used.
An NO 2 /NO x ratio of 0.3 was applied to recalculate NO 2 emissions for this study because of the expected higher real-world NO 2 emissions from diesel vehicles. The applied value is higher than suggested by the reported range (3.2-23.5 vol %) of the primary NO 2 emission fraction from vehicular traffic in London (Carslaw and Beevers, 2005) and the NO 2 /NO x ratio of 0.22 for passenger cars in urban areas assumed by Keuken et al. (2012) for the Netherlands. The use of the high NO 2 /NO x ratio for the Hamburg vehicle fleet is consistent with the higher NO 2 /NO x ratio from diesel passenger cars (from 0.12 to > 0.5; Carslaw and Rhys-Tyler, 2012) and the review by Grice et al. (2009), who assumed that Euro 4-6 passenger cars emit 55 % of the total NO x as NO 2 .

Setup of the model evaluation and performance analysis
A 1-year simulation with EPISODE-CityChem was performed for the study domain using the model setup as described in Sect. 4.1.1. Model results were compared to results from the standard EPISODE model to assess the total effect of the new implementations of the CityChem extension.
In the standard EPISODE model, the PSS approximation is used at the receptor points and on the Eulerian grid; the street canyon model and the WMPP module were deactivated. For the simulation with EPISODE-CityChem and standard EPISODE, the boundary conditions from hourly 3-D concentrations of CMAQ were taken as described in Sect. 2.3.1. In addition, EPISODE-CityChem results were compared to results from the air pollution module of TAPM, which is used as a reference model in this study. The TAPM run was performed with the same horizontal resolution (1 km), as well as identical meteorological fields and urban emissions, but with 2-D boundary conditions instead of 3-D boundary conditions. The hourly 2-D boundary concentrations for TAPM were prepared by using the horizontal wind components on each of the four lateral boundaries to weight the CMAQ con- Statistical indicators of the evaluation included the mean (modelled / observed), standard deviation (SD; modelled / observed), correlation coefficient (Corr), root mean square error (RMSE), overall bias (Bias), normalized mean bias (NMB) and index of agreement (IOA). See Appendix E for the definition of the statistical indicators.
The FAIRMODE (Forum for Air Quality Modelling in Europe) DELTA tool version 5.6 (Thunis et al., 2012a(Thunis et al., , b, 2013Pernigotti et al., 2013;FAIRMODE, 2014;Monteiro et al., 2018) was used for the evaluation of model results from air quality simulations for Hamburg. The DELTA tool focuses primarily on the air pollutants regulated in the current AQD (EC, 2008).

Evaluation of downscaled meteorological data
Downscaled meteorological data on temperature, relative humidity, total solar radiation, wind speed and wind direction were examined. Hourly-based data from the meteorological station at the Hamburg airport (Fuhlsbüttel) (operated by DWD) and from the 280 m high Hamburg weather mast at Billwerder (operated by Universität Hamburg) were analysed. Observation data from the DWD station at 10 m of height and from the Hamburg weather mast at 10 and 50 m of height were used in the analysis. TAPM modelled meteorological data from the 1 × 1 km 2 grid cell of the D4 domain, where the stations are located, at the corresponding height were extracted for comparison with observations. Table S4 provides an overview of the statistical analysis of TAPM data.
Hourly temperature predicted by TAPM was in excellent agreement with the observed temperature at both stations and both heights, showing high correlation (Corr ≥ 0.98) and small overall bias (≤ 1.00 • C). Relative humidity also showed good agreement but with lower correlation (Corr = 0.74). Total solar radiation was predicted by TAPM with high correlation (Corr = 0.86) but a high positive overall bias of 27 W m −2 . Situations with reduced solar radiation due to high cloud coverage are often not well captured by TAPM. The IOA for temperature, relative humidity and total solar radiation was 0.99 (average of all observations), 0.86 and 0.92, respectively.
TAPM shows good predictive capabilities for wind speed and direction. Due to the assimilation of wind observations at the DWD Hamburg airport station for nudging wind speed and direction in TAPM meteorological runs, the meteorological performance for wind speed and direction was only compared at the Hamburg weather mast. Modelled hourly data for wind speed at the Hamburg weather mast agreed well with the observations throughout the year at 10 m (Corr = 0.87, Bias = −0.08 m s −1 ) and 50 m of height (Corr = 0.85, Bias = −0.02 m s −1 ); they were within the observed variability. Southwest and west are the most frequent wind directions in Hamburg due to prevailing Atlantic winds, followed by winds from the southeast (Bruemmer et al., 2012). Mean wind direction was computed as a circular average (unit vector mean wind direction) for model and observation data. At the Hamburg weather mast, modelled versus observed wind directions show good agreement at 10 and 50 m of height (IOA ≥ 0.89) with a bias in the mean wind direction of 16.9 and 6.2 • at 10 and 50 m, respectively. The difference is due to a slightly higher frequency of winds from the west predicted by TAPM.

Evaluation of the temporal variation of pollutants
The statistical performance of the models with respect to temporal variation was assessed by comparing modelled concentrations against observed concentrations from the AQ monitoring network of Hamburg. The stations of the monitoring network and the available measurements of pollutant concentrations are listed in Table S5. The monitoring network covers all parts of the city (Fig. 6c). Minimum data availability is required for statistics to be produced at a given AQ monitoring station. In the DELTA tool, the requested percentage of available data over a selected time period (here: 1 year) is 75 %, as defined in the AQD. This has been fulfilled by all stations in Hamburg, except for O 3 and PM 2.5 measurements at two stations. The statistical analysis included all stations for which the data availability criterion was fulfilled. For the comparison, model output at the exact geographic location of the monitoring stations from the model was used. Concentrations of NO 2 and NO were measured at all stations included in this study. The model performance statistics are listed in Table 7 for NO 2 (based on hourly values), in Table 8 for O 3 (based on daily max. of 8 h running mean) and in Table 9 for PM 10 (based on daily mean).
EPISODE-CityChem performs fairly well for NO 2 based on hourly values, with an IOA of 0.70 (average of all stations) and correlation of 0.53 (average of all stations). The average performance at the traffic stations (Corr = 0.57, IOA = 0.73 on average) is better than for the other stations. EPISODE-CityChem performs better than EPISODE and TAPM at the traffic stations. This implies that the use of the street canyon model improves the agreement between the model and observations of NO 2 at traffic stations. In particular, TAPM shows weak correlation at the traffic sites (Corr = 0.29, IOA = 0.52 on average). On the other hand, TAPM shows low bias and relatively good correlation at the industrial sites and stations influenced by emissions from industry and shipping (20VE, 21BI, 61WB, 80KT). For most urban background stations, EPISODE-CityChem tends to underestimate the observation mean. The NMB is slightly negative for stations of the urban background in suburban areas (NMB = −14.2 % on average), which is indicative of NO x emissions from suburban areas that are too low or dispersion of local emissions that is too efficient. Compared to EPISODE, the EPISODE-CityChem model has a lower positive bias at some industrial and background stations (20VE, 21BI, 74BT, 80KT), probably due to the combined effect of the more advanced photochemistry and WMPP used for point source plumes.
EPISODE-CityChem performs better for the O 3 daily maximum of the 8 h running mean than the other two models. The overall performance for ozone is very good, with an IOA of 0.83 (average of all stations), correlation of 0.76 (average of all stations) and small NMB (within ±18 %). The correlation improvement compared to EPISODE is relatively small, and hence the better performance compared to TAPM might be largely from the use of the more comprehensive set of boundary conditions.
The performance of EPISODE-CityChem for PM 10 based on daily means is good, with an IOA of 0.74 (average of all stations). Correlation values are generally satisfactory, with a station average of Corr = 0.60; RMSE values are fairly small, within the range 8.9-12.9 µg m −3 for all stations. For the urban background stations, EPISODE-CityChem and EPISODE perform better than TAPM in terms of correlation and IOA, giving a clear indication of the advantage of using 3-D boundary conditions instead of 2-D boundary conditions. Including the street canyon model leads to an overestimation of the observed daily mean PM 10 concentrations at the traffic stations. However, it cannot be concluded whether the overestimation is due to shortcomings in the street canyon module or due to inaccurate traffic emissions in the respective streets.

Evaluation of the spatial variation of pollutants
Annual mean concentrations of regulatory air pollutants from the EPISODE-CityChem model output were compared to the available observation data (Fig. 7). The model reproduces the spatial variation of NO and NO 2 concentrations (Fig. 7a, b) and the concentration gradients of NO 2 and NO between the urban background (80KT, 51BF, 52NG, 13ST; 61WB, 54BL, 27TA, 74BT), traffic stations (68HB, 64KS, 70MB, 17SM) and industrial stations (21BI, 20VE). For most stations the overall bias of annual mean observed NO 2 is within ±10 µg m −3 . Observed levels of annual mean O 3 at the five urban background stations (13ST, 27TA, 51BF, 54BL, 52NG) are captured by EPISODE-CityChem within ±15 % (Fig. 7c). TAPM overestimates annual mean O 3 by 10 %-25 %, except for the inner-city background station 13ST. Annual mean SO 2 was compared at five stations (Fig. 7d). With the exception of the industrial station 20VE, modelled annual mean SO 2 from EPISODE-CityChem agreed with the observed concentrations within a factor of 2. At station 20VE, modelled SO 2 is about 6 times higher than observed SO 2 . Obviously, the models overestimate the influence of SO 2 emissions from nearby industrial sources. Modelled annual mean SO 2 from EPISODE and TAPM was even higher at this site. The slightly better agreement of EPISODE-CityChem compared to EPISODE might be due to the use of WMPP or to the consideration of the OH reaction of SO 2 . Levels of PM 2.5 and PM 10 in the models are controlled by the primary emission of particulate matter and its atmospheric dispersion, while secondary aerosol formation is not considered in the models. The comparison of the time series of daily means at station 13ST in Fig. S2 indicates that EPISODE-CityChem is able to capture observed peak events with high PM 2.5 and PM 10 concentrations during the winter season (DJF). These events are likely related to the short-or long-range transport of anthropogenic emitted PM or secondary produced inorganic PM. Observed levels of annual means of PM 2.5 are matched by the model within ±43 % (Fig. 7e). TAPM underestimates observed annual mean PM 2.5 and PM 10 at the traffic stations. We note that TAPM has no specific treatment of local dispersion at the street scale. Line source emissions in TAPM are released into the volume of the corresponding grid cell in the surface layer (10 m of height) and therefore immediately diluted.
EPISODE-CityChem overestimated annual mean PM 10 at three traffic stations (17SM, 68HB and 70MB) by 22 %-24 %, while the agreement with observed PM 10 was excellent at the other stations (Fig. 7f). Since EPISODE shows lower modelled PM 10 at the traffic stations than EPISODE-CityChem, we conclude that the overestimation of PM 10 at the traffic stations is due to the street canyon module. How-ever, the estimation of traffic emissions of PM 10 is complicated by coarse particles from non-exhaust emissions, such as tyre abrasion and the resuspension of road dust. In addition, traffic emissions used in the simulations are not based on actual traffic counts, so it is entirely possible that the included PM 10 line source emissions in the respective streets were too high.

Model performance analysis for policy support applications
We assessed whether the model results from EPISODE-CityChem have reached a sufficient level of quality for a given policy support application and compared the outcome to EPISODE and TAPM. The model quality objective in the FAIRMODE DELTA tool has been constructed on the basis of the observation uncertainty and describes the minimum level of quality to be achieved by a model to be fit for policy use (Thunis et al., 2012a(Thunis et al., , b, 2013Pernigotti et al., 2013). The model quality indicator (MQI) provides a general overview of the model performance. The associated model performance criteria (MPC) for correlation, standard deviation and bias can be used to highlight which of the model performance aspects need to be improved. Details on the MQI and MPC are given in Appendix E. Figure 8 shows the model performance evaluation of EPISODE-CityChem in terms of fitness for purpose in the form of scatter diagrams and target diagrams (Thunis et al., 2012a;Monteiro et al., 2018) for NO 2 (hourly), O 3 (daily max. of the 8 h running mean) and PM 10 (daily mean). For the yearly averaged values shown in the scatter diagrams, Table 9. Model performance statistics of EPISODE-CityChem (EPCC), standard EPISODE (EPIS) and TAPM for the temporal variation of PM 10 based on daily mean concentration at all stations with sufficient data available in 2012. Bold numbers represent model results with better performance for Corr, RMSE and IOA per station.  Fig. 8a-c include data from all three models. For EPISODE-CityChem, the data indicate fulfilment (stations lie within green-shaded zone) for all pollutants at all stations. The scatterplot for NO 2 reveals that for TAPM two stations are in the orange-shaded zone (fulfilment but error dominated by bias), and one station does not fulfil the bias MPC. EPISODE-CityChem shows scatter around the 1 : 1 line with more stations lying below than above, indicating a slight tendency to underestimate yearly averaged NO 2 . The scatterplot for O 3 indicates fulfilment for all models, with EPISODE-CityChem and TAPM predictions being closer to the 1 : 1 line than EPISODE. For PM 10 , EPISODE-CityChem and EPISODE gave similar results for most stations (Fig. 7f), and data points therefore overlap in the scatter diagram. EPISODE-CityChem shows values close to the upper measurement uncertainty range at the three traffic stations. TAPM shows an outlier for one traffic station (68HB) with predicted PM 10 that is too low; however, it is still within the green-shaded zone.

Station Model
In the target diagrams ( Fig. 8d-f) the centred RMSE (CRMSE) for each monitoring station is plotted against the normalized bias. The distance from the circle origin gives an estimate for the MQI; MQI values below unity (i.e. within the green-shaded area) are considered to comply with the model quality objective required for policy applications. The EPISODE-CityChem model fulfils the model performance objectives for NO 2 , O 3 and PM 10 in terms of both bias MPC and MQI (max. MQIs for hourly and daily values are NO 2 0.891, O 3 0.506, PM 10 0.820). All stations are located within the green-shaded zone, indicating fulfilment of the RMSE criteria as defined by Eq. (E10), but outside the dashed circle; i.e. the difference between the model and observations is not within the measurement uncertainty range. In the right quadrant of the target diagram (Fig. 8d-f  predicted hourly NO 2 and the predicted daily mean PM 10 is 35 % and 36 %, respectively. Corresponding target diagrams for TAPM and EPISODE are shown in Fig. S3. Notably, for hourly and yearly averaged NO 2 the two other models show weaker performance than EPISODE-CityChem. EPISODE fulfils the MQI for hourly values but exceeds the MQI for yearly averages. TAPM exceeds both MQI criteria, mainly due to the weak performance at the traffic stations (see Sect  Figure 9 depicts spatial maps of the annual mean concentrations of NO 2 , O 3 , total NMVOC, SO 2 , gaseous sulfuric acid (H 2 SO 4 ) and PM 2.5 from the model output of the receptor grid (resolution 100 × 100 m 2 ). With the exception of PM 2.5 , concentrations of the aforementioned compounds are modulated by photochemical reactions in the model simulation. Due to the large temporal, spatial and composi-tional variations of the input from anthropogenic emissions of NO x , NMVOCs and CO within the urban area, the atmospheric chemistry in urban environments is complex. Prevailing winds from the west, on an annual basis, allow for a simplified view of the inflow-outflow pattern for ozone within the study domain of Hamburg. Following the inflow-outflow direction in space from west to east (30 km), the modelled O 3 concentration starts with ca. 50 µg m −3 at the western border, is largely reduced in the inner city to 15-30 µg m −3 and gradually increases again over the eastern part to ca. 40 µg m −3 .
Modelled and observed annual mean NO 2 at the four traffic stations exceeds the annual limit value of 40 µg m −3 in 2012 (see Fig. 7b). The model simulation suggests that there is widespread exceedance of this limit value in the city, mainly at the main streets and road junctions of the inner city and along motorways outside the inner city (Fig. 9a).
Modelled ozone at the outflow border does not reach the level at the inflow border. Within the urban area, trafficrelated emissions of NO destroy much of the O 3 (mainly at night when O 3 is not recycled through the photolysis of NO 2 ), clearly seen as minimum concentrations along the traffic network (Fig. 9b). Thus, the inner urban area pro- vides an efficient sink for ozone, which is qualitatively in accord with the findings of the REPARTEE (Regents Park and Tower Environmental Experiment) measurement campaign carried out in London in the autumn of 2006(Harrison et al., 2012. The photochemical production of O 3 from NO x , NMVOCs and CO, emitted in the urban area, is very limited in the inner city. Main sources of NMVOC in Hamburg are solvent use (SNAP cat. 06) and traffic emissions. NMVOC annual mean concentrations of more than 40 µg m −3 were modelled close to roads in the inner city (Fig. 9c). Ehlers et al. (2016) report the similarity of NMVOC fingerprints in air samples taken in the inner city of Munich and in a road tunnel in Berlin to the fingerprint of petrol cars under cold-start conditions. Our model simulation finds that the loss of NMVOCs by OH-initiated oxidation in the model is inhibited due to the low O 3 concentrations in the inner city. In summer, modelled OH midday maximum concentrations are in the range (0.5-2.0) ×10 6 molecule cm −3 in the inner city. Modelled oxylene, which is the model surrogate compound for the sum of aromatic VOCs, is 5-10 µg m −3 (1.2-2.2 ppbv) at some distance from the roads. This is in the same range as concentrations measured in central London (ca. 2 ppbv for the sum of aromatics; Valach et al., 2015).
SO 2 is an important precursor for secondary aerosol formation. SO 2 emissions in Hamburg are mainly from industrial point sources and ship traffic in the harbour area. The highest yearly averaged modelled SO 2 concentrations are in the range of 20-40 µg m −3 in proximity to the main sources (Fig. 9d).
In the atmosphere, SO 2 reacts with the OH radical and with CH 3 O 2 to give gaseous sulfuric acid (H 2 SO 4 ). The presence of sulfuric acid in gaseous concentrations of 10 6 -10 7 molecule cm −3 is necessary in order to observe new particle formation events in the atmosphere (Zhang et al., 2012). In the model, a constant very low boundary condition (BCON) value (10 −5 µg m −3 ) was chosen for H 2 SO 4 , leading to a reduced sulfuric acid concentration in the boundary cells (Fig. 9e). Towards the inner domain, H 2 SO 4 quickly increases due to the oxidation of SO 2 advected to Hamburg from the regional background. Modelled annual mean H 2 SO 4 peaks in the harbour area with up to 0.018 µg m −3 . On spatial average, the H 2 SO 4 annual mean concentration is 0.013 µg m −3 , corresponding to 8.2 × 10 7 molecule cm −3 , which is higher than typical ambient concentrations in the urban atmosphere. For comparison, reported maximum midday H 2 SO 4 concentrations in Beijing are in the range (0.3-1.1) ×10 7 molecule cm −3 (Wang et al., 2013). Modelled H 2 SO 4 that is too high is explained by the fact that the condensation of sulfuric acid on pre-existing particles is not accounted for in the model. Condensation is the most important atmospheric sink of sulfuric acid and also leads to the formation of particulate sulfate.
Previous studies have shown that levels of PM 2.5 in the urban background are to a large extent controlled by atmo-spheric transport from upwind regions. Modelled PM 2.5 levels in the urban background of Hamburg are on annual average 2.5 µg m −3 higher than the regional background. The urban increment due to road traffic exceeds the urban background by 10-20 µg m −3 in the model simulation (Fig. 9f). Hotspots of PM 2.5 pollution are found within the harbour area to the southwest of the inner city, where several refineries, power plants and industries are located. These industrial emissions are represented as point source plumes in EPISODE-CityChem.

Modelling of ozone formation
To demonstrate the effect of using an advanced photochemical mechanism for modelling ozone concentrations and ozone production, we compared two EPISODE-CityChem summer (JJA) simulations for Hamburg: one using EmChem09-mod chemistry (EmC09-mod) and the other using the PSS assumption (PSSA) on the Eulerian 3-D grid. The implementation of the PSS involves only three reactions for the equilibrium of O 3 /NO/NO 2 (i.e. Reactions R5-R7), while ignoring all other reactions. Spatial maps in Fig. 10a and b show the difference in summer mean concentrations of NO 2 and O 3 between the runs EmC09-mod and PSSA. In EmC09-mod, additional NO-to-NO 2 conversions occur compared to the PSS due to the oxidation of CO and VOCs. The largest impact on NO 2 , an increase by ca. 1.4-2.3 µg m −3 compared to PSS, happens in the western part of the city, along the river Elbe and in the harbour area south of the river. However, the increased NO 2 does not translate into an increase in O 3 in this area because modelled NO 2 is already high (> 35 µg m −3 ), and additional NO 2 decreases the OH concentration via Reaction (R11). Remarkably, photochemical ozone production takes place in the outflow to the southeast of the city and in the urban background southwest of the city. The O 3 concentration difference between EmC09mod and PSSA is in the range 1.6-2.2 µg m −3 in the outflow of polluted air, implying that advanced photochemistry is necessary for a more accurate prediction of O 3 in the urban background. The modelled VOC/NO x ratio in the pollution outflow is ∼ 0.6, indicating VOC-limited conditions. This is in contrast to the general expectation of finding NO x -limited conditions in the outflow of cities; the low VOC/NO x ratio might be associated with an overestimation of NO x by the model, as will be shown below. In the inner city and in the vicinity of highly trafficked streets and motorways, photochemical ozone production is very small (difference < 0.4 µg m −3 ), suggesting that the PSS assumption is valid close to the sources of NO x pollution.
Comparison with measurements was done at two sites: 13ST, an inner-city urban background monitoring site, and 27TA, an urban background station in the SE part of the domain (green stars in Fig. 10a and b). The diurnal cycles of observed median O 3 and NO 2 concentrations at 13ST for summer (JJA) show a minimum of O 3 in the early morning hours between 04:00 and 07:00, as well as a daily maximum of O 3 between noon and 16:00, whereas NO 2 peaks two times during the day (06:00-08:00 and at 18:00-20:00), coinciding with traffic rush hours ( Fig. 10c and d). The diurnal cycles of observations follow a similar pattern at 27TA, but the observed median of NO 2 is roughly half that at 13ST (Fig. 10e  and f).
At 13ST in the afternoons between 14:00 and 19:00, modelled median O 3 concentrations from both EmC09-mod (red line) and PSSA (green line) are below the bandwidth of observed ozone (25th to 75th percentile), and modelled median O 3 is ca. 20 µg m −3 lower than the observed median. At both stations, modelled median NO 2 in the evening (between 16:00 and 20:00) is higher than the bandwidth of observed NO 2 . The prediction of an NO 2 concentration that is too high in the evening at both sites could be due to several reasons: (1) overestimated emissions of NO x in the urban area; (2) inadequate diurnal profile of traffic emissions regarding the afternoon rush hour; and (3) inaccurate mixing height of the nocturnal BL, leading to enhanced accumulation of NO x emissions in the surface layer of the model.

Planned improvements to the EPISODE model
The future development of the EPISODE model with respect to photochemistry and dispersion near sources is outlined in the following. Specifically, the implementation of the photochemistry in relation to sub-grid modelling and the treatment of aerosol formation on the urban scale will be the focus of the planned development for the next versions of EPISODE-CityChem.

Photolysis parameterization development
The procedure for calculating photodissociation rates (j values) has not been changed since the original development of the EMEP45 mechanism and is thus as documented by Walker et al. (2003). This procedure was based on the procedure used in the EMEP oxidant model at that time and needs to be revised and updated. The plan is to update the methodology in accordance with the present version of the EMEP model . In the present EMEP model, the j values are based on precalculated rates from a detailed radiative-transfer model (PHODIS; Kylling et al., 1998) and interpolated between certain fixed cloud fractions.

EP10-Plume development
The photochemical steady state might apply in street canyons because the distance between the source and receptor is short there, and hence only the fastest chemical reactions can have a significant influence on photochemical transformation in the street canyon air. On the timescales applying to a street canyon, CO and hydrocarbons can be treated as inert tracers. However, the photochemical scheme applied in the sub-grid component also needs to consider situations with a larger distance between the source and receptor (within a cell of the Eulerian grid).
The PSS relationship is usually not valid in the urban air because organic peroxy radicals compete with O 3 to convert NO to NO 2 as a result of the oxidation of VOCs in sunlight. EP10-Plume considers only the photochemical degradation of HCHO, but not the oxidation of other reactive VOCs emitted from traffic. It is planned to develop a more detailed chemistry scheme, including VOCs, at the receptor points.
Background ozone concentrations are taken into account in photochemical transformation in the sub-grid component (in both PSS and EP10-Plume). Due to the low ventilation in street canyons, the residence time of pollutants becomes comparable to the timescales of the reactions involved in the PSS assumption, i.e. on the order of tens of seconds. This implies that the exchange rate between the plume from the line source and the background air can become a limiting factor for photochemical transformation in street canyons. The exchange rate is governed by the residence time of the pollutants at the street. However, in the sub-grid line source model, it is assumed that background O 3 is instantaneously mixed with NO from the plume of the line source. The rate of reaction (R5) might therefore be overestimated in the subgrid model, depending on the ambient conditions for plume mixing. Hertel and Berkowicz (1989) take the exchange rate into account for the PSS and suggest that the residence time of pollutants in a street canyon can be approximated by H sc / σ wt , where σ wt is the ventilation velocity. The exchange rate of ozone will be considered in the refinement of the photochemistry of the sub-grid component.

Emissions of biogenic VOC
Biogenic emissions of VOCs from trees might be relevant in the VOC-limited regime. Many urban trees, such as European aspen (Populus tremula) and deciduous oaks, emit large amounts of isoprene (Karl et al., 2009). The ozone formation potential of BVOCs emitted from urban trees is sufficiently high to outperform the ozone uptake capacity of the trees (Grote et al., 2016). Missing emissions of BVOCs from trees in urban green parks and at the roadside might partly explain the underestimation of observed O 3 in summer in the inner city by the model. Isoprene, in the presence of a sufficiently high level of NO x , can contribute substantially to O 3 formation in the urban atmosphere. Further, monoterpenes are considered to be a relevant source of organic peroxy radicals at night (Platt et al., 2002). EmChem09-mod includes chemical reactions for the OH-initiated oxidation of isoprene and for two types of monoterpenes, represented by α-pinene and limonene. It is planned to implement a module for the computation of inline emissions of isoprene and monoterpenes as a function of temperature and solar radiation using data on tree-specific BVOC emission factors by Karl et al. (2009) and high-resolution urban land use information, for example from the Copernicus Urban Atlas.

Secondary formation of particulate matter
Currently, both PM 2.5 and PM 10 are treated as inert tracers with no secondary aerosol formation. The first step towards a better representation of the particulate phase will be the separation of particulate matter into individual particulate chemical compounds. Traffic emissions of fine (PM 2.5 ) and coarse (PM 10 -PM 2.5 ) particulate matter will be separated into (non-volatile) primary organic aerosol (POA) and elemental carbon (EC), thereby assuming that any emission of SVOCs and intermediate-volatility volatile organic compounds (IVOCs) immediately undergoes irreversible condensation to the exhaust particulate matter. A clear advantage of the approach is that it avoids the need to discount SVOCs and IVOCs from the city's particulate matter emission inventory. Mineral dust and sea salt imported to the city from the regional background need to be treated as separate compounds. The main focus is on the formation of secondary inorganic aerosol (SIA). Thus, the Model for an Aerosol Reacting System (MARS; Binkowski and Shankar, 1995) could be used for the calculation of the partitioning between gas and fine-mode particles . However, MARS does not account for mineral dust components and sea salt; the latter could be relevant for the formation of SIA in coastal cities because it increases the water associated with particulate matter.

Ultrafine particles
Among the emerging, yet unregulated, pollutants in cities are ultrafine particles (UFPs; diameter less than 100 nm). Major sources of UFPs in urban environments are motor vehicle exhaust emissions (e.g. Harrison et al., 2011). Emissions from ships contribute to UFP pollution in harbour cities (Pirjola et al., 2014). High concentrations of ultrafine particles are formed during new particle formation events (Salma et al., 2016). Studies have demonstrated the relevance of episodes of new particle formation in urban environments for cities situated in high insolation regions such as southern Europe (e.g. Pey et al., 2008;Dall'Osto et al., 2013;Brines et al., 2015). UFPs are usually evaluated in terms of particle number (PN) concentrations. In urban environments, UFPs dominate the total PN concentrations, but only make a small contribution to particulate matter. A simplified parameterization for the treatment of the dry deposition and coagulation of particles  has already been implemented in EPISODE for modelling PN concentrations in Oslo . The MAFOR aerosol dynamics solver (Karl et al., 2011 will be implemented in EPISODE-CityChem to compute information on the size distribution of UFPs and the total PN. The solver includes nucleation, coagulation, growth due to the condensation of sulfuric acid, and low volatile and/or semi-volatile organic vapours from biogenic and anthropogenic sources.

SSCM development
In addition to the dependence on the wind direction (considered through θ street , the angle between wind direction and road axis), the direct contribution in SSCM is mainly sensitive to the emission intensity of the line source, the streetlevel wind speed and the integration path. The integration path corresponds to the length of the recirculation zone but extends to L max for close-to-parallel wind directions. The length of the recirculation zone (L rec ) depends on the building height along the canyon, while L max is a function of the canyon width for large θ street and a function of the canyon length if θ street is below 45 • . Both dimensions (L rec and L max ) are only roughly estimated in SSCM because it considers only generic street canyon types and not the site-specific street canyon geometry. It is planned to refine SSCM with respect to a better representation of urban street canyon geometry. This could be done, for example, by using spatially resolved information on building height and street canyon width by extracting data from 3-D city building models for each road segment.

WMPP development
Envisaged future improvements to the WMPP are to extend it to use energy budget methods (Thom et al., 1975;van Ulden and Holtslag, 1985;Fritschen and Simpson, 1989;Tunick, 2006) in combination with net surface radiation and heat flux to determine wind and turbulence profiles and mixing height. The method will be based on data such as solar radiation, cloud cover, air temperature, wind speed, relative humidity, precipitation and surface conditions.

Conclusions
The CityChem extension of the urban AQ model EPISODE (Slørdal et al., 2003(Slørdal et al., , 2008Hamer et al., 2019) offers a detailed treatment of the atmospheric chemistry in urban areas and a more advanced treatment for dispersion close to point emission sources, such as industrial stacks, and line emissions sources, such as open roads and street canyons. EPISODE consists of a 3-D Eulerian grid CTM with embedded Gaussian dispersion models that track the sub-scale dispersion of pollutants from line and point emission sources until the fine-scale variability becomes unimportant. The EPISODE-CityChem model, which is based on the core of the EPISODE model, integrates the CityChem extension into an urban CTM system with the capability to simulate the photochemistry and dispersion of multiple pollutants on urban scales. Photochemistry on the Eulerian grid is computed using a numerical chemistry solver. Photochemistry in the sub-grid components is solved with a compact reaction scheme, replacing the photo-stationary-state assumption. The integration of SSCM in the sub-grid line source model results in higher concentrations in street canyons because it considers the reduced ventilation inside the canyon and the recirculation of the traffic plume. The integration of WMPP in the sub-grid point source model for dispersion around an elevated point source increases the maximum ground concentration of an inert tracer by a factor of 4 to 6 in neutral and unstable conditions compared to the standard parameterization in EPISODE.
The EPISODE-CityChem model takes into account the fact that long-range transport contributes to urban pollutant levels by using hourly varying pollutant concentrations at the lateral and vertical boundaries from the CMAQ model (Byun and Schere, 2006) as initial and boundary concentrations. The model reads meteorological fields generated by the prognostic meteorology component of TAPM (Hurley, 2008;Hurley et al., 2005) but can also use meteorological fields constrained by observations.
The performance of EPISODE-CityChem was evaluated with a series of tests to study the basic functionalities of the CityChem extension and with a first application to the air quality situation in the city of Hamburg, Germany. The ability to reproduce the temporal variation of major regulated pollutants at AQ monitoring stations in Hamburg was compared to that of the standard EPISODE model and TAPM (AQ model) using identical meteorological fields and emissions. EPISODE-CityChem performs better than EPISODE and TAPM for the prediction of hourly NO 2 concentrations at the traffic stations, which is attributable to the street canyon model. EPISODE-CityChem was in better agreement with the observed O 3 daily maximum of the 8 h running mean than the other two models. For daily mean PM 10 at urban background stations, EPISODE-CityChem and EPISODE gave better results than TAPM, largely due to the use of hourly 3-D boundary conditions from CMAQ. The performance analysis with the FAIRMODE DELTA tool for air quality in Hamburg showed that EPISODE-CityChem fulfils the model performance objectives for NO 2 (hourly), O 3 (daily max. of the 8 h running mean) and PM 10 (daily mean) set forth in the AQD, qualifying the model for use in policy applications.
The effect of using an advanced photochemical mechanism (EmChem09-mod) compared to the PSS assumption for modelling ozone concentrations and ozone production was investigated in summer (JJA) simulations for Hamburg. Photochemical ozone production was found to take place in the outflow of polluted air from the city, implying that advanced photochemistry is necessary for a more accurate prediction of O 3 in the urban background. However, the modelled daily maximum O 3 in summer afternoons was ca. 25 % lower than observed at an inner-city urban background station. In addition, the model predicted an NO 2 concentration in the summer evenings at two urban background stations that is too high. Further investigation of the high modelled evening NO x in summer will require a sensitivity analysis of the various source categories contributing to the NO x levels in the inner city, which remains as a task for future studies.
BVOC emissions from urban parks and forests might partly explain the underestimation of observed ozone in summer. There is evidence that the contribution of BVOC emissions to ozone formation can be up to 60 % during heat waves in densely populated areas (Churkina et al., 2017), depending on the type and amount of urban vegetation. In the future, BVOC emissions in urban areas are expected to become even more important in ozone formation if anthropogenic NMVOC emissions continue to decline as a result of technological progress (Wagner and Kuttler, 2014). The implementation of BVOC emissions as a function of temperature and solar radiation will be the focus of coming developments of the CityChem extension.
Envisaged applications of the EPISODE-CityChem model are urban air quality studies, emission control scenarios in relation to traffic bans introduced in German cities and the source attribution of sector-specific emissions to observed levels of air pollutants. The model can also be utilized in the evaluation of air pollution exposure and in the assessment of adverse health impacts. Features of the model that facilitate its application to urban AQ in cities worldwide include integrated utilities for input preparation and output processing, moderate computational demand, photochemistry options (ozone formation studies), high spatial and temporal resolution, and demonstrated fitness for policy use.
Code and data availability. The source codes of the EPISODE-CityChem model version 1.2 and the preprocessing utilities are accessible in release under the RPL licence at https://doi.org/10.5281/zenodo.3063356 .
A tar package with example data for a 1-month simulation and the user's guide are included in the release. All preprocessing tools are written in Fortran 90. Software requirements for the utilities and the EPISODE-CityChem model are installation of the GCC-GFortran Fortran 90 compiler (version 4.4. or later) and the NetCDF library (version 3.6.0 or later).
The following datasets are available for download from the HZG FTP server upon request: where H is the influence function for heat transfer, z 0 is the surface roughness (for momentum) and L is the Monin-Obukhov length. The quasi-laminar layer resistance is calculated according to the parameterization given by Simpson et al. (2003). The canopy resistance, i.e. deposition due to the capture of pollutants by the surface, is currently only considered by a minimum value, i.e. R c = R c,min = 1 m s −1 . The parameterization of the canopy resistance is complex, since it depends on both surface characteristics and the chemical characteristics of the depositing gas. The dry deposition velocity of particles, V dry,p , is calculated as (Simpson et al., 2003) where υ s is the gravitational settling velocity and the other terms are as for gases. Equation (C4) involves the assumption that all deposited particles stick to the surface so that the surface resistance becomes zero. The dry deposition velocity of atmospheric aerosols depends on their sizes. The current formulation distinguishes between PM 2.5 and PM 10 , which are presently assigned the particle diameters of 0.3 and 4 µm. All the resistances are integrated over the aerosol sizes, assuming a lognormal particle size distribution with geometric standard deviations of 2.0 and 2.2 µm for PM 2.5 and PM 10 , respectively.

C2 Wet deposition
Wet deposition is described as a sink term within the advection-diffusion equation and can be parameterized by dC i /dt = − · C i , where C i is the grid concentration of a gaseous or particulate species and is the scavenging coefficient (s −1 ). Wet scavenging is different from zero in grid cells in which precipitation (rainfall or snowfall) occurs. The chosen crude approach for representing wet deposition treats in-cloud scavenging in the same way as below-cloud scavenging. Further, the cloud base is assumed to be at the model top, which means that scavenging occurs throughout the entire 1-D model column for which the precipitation rate in the surface grid cell is greater than zero. For the short-term estimation of near-ground concentrations in urban areas, belowcloud scavenging is expected to be the dominant wet removal process. A more accurate treatment of below-cloud scavenging requires knowledge of the cloud base height (which is not standard output of TAPM) in order to limit wet deposition to the model layers that are actually affected by raining clouds and to separate between in-cloud and below-cloud scavenging. The scavenging of gases is calculated as (Simpson et al., 2003) C i,wet = −C i W sub P r H sub + ρ w , where W sub is the sub-cloud scavenging coefficient for gases, supplied as a constant value by the model user, P r (kg m 2 s −1 ) is the precipitation rate, H sub is the scavenging depth (corresponding to the total vertical depth of the model) and ρ w is the water density (1000 kg m −3 ). Precipitation is a 2-D surface field, either from observations of precipitation rate or computed by TAPM. The wet deposition rate of particulate compounds is calculated as (Simpson et al., 2003) where V dr is the raindrop fall speed (V dr =5 m s −1 ), A = 5 m 3 kg −1 s −1 is an empirical coefficient when the Marshall-Palmer size distribution is assumed for raindrops and E is the tabulated size-dependent collection efficiency of aerosols by the raindrops.

Appendix D: Treatment of boundary concentration data
The utility BCONCC v2.1 (included in the CityChem distribution) is used to produce EPISODE-format binary files containing hourly varying 3-D boundary concentrations for the most relevant chemical compounds. BCONCC makes use of the IOAPI version 3.1 library (https://www.cmascenter.org/ ioapi, last access: 29 July 2019) to access CMAQ output files. The IOAPI (Models-3/EDSS Input/Output Application Programming Interface) provides a variety of data structure types for organizing the data and a set of data access routines. The preparation of chemical boundary conditions from the CMAQ model output is done in two steps. First, the city's 3-D domain extent plus one grid cell to each side is cut out from the CMAQ model grid, interpolating the hourly concentrations to the horizontal main grid resolution of EPISODE using bilinear interpolation. Second, EPISODE-format binary files for boundary conditions (BCON files) containing background concentrations of all individual CityChem compounds are created for the defined model domain in the required input format. Linear interpolation is used to convert concentrations from the vertical layers of the CMAQ model to the vertical layers of the EPISODE-CityChem model. Temperature and pressure from the METCRO3D file (meteorological input file of the CMAQ simulation) are used to convert the concentrations of gaseous compounds from mixing ratios to mass-based concentrations.
The background concentrations are adopted for the grid cells directly adjacent to the grid cells of the model domain (with nx × ny cells per model layer) and also for the vertical model layer that is on top of the highest model layer. Boundary conditions from CMAQ concentrations are created for the gas-phase chemical compounds: O 3 , NO, NO 2 , H 2 O 2 , SO 2 , HCHO, CO, N 2 O 5 , HNO 3 , PAN and the individual VOCs. Boundary conditions for PM 2.5 include primary aerosol components: EC, POA, sea salt (NaCl) and mineral dust. They also include secondary inorganic aerosol (SIA) components: sulfate (SO 2− 4 ), ammonium (NH + 4 ), ni-trate (NO − 3 ) and SOA (PM 2.5 was defined including modes I and J of the CMAQ aerosol components). Since the focus of the AQ study is mainly on photochemistry and fine particulate mass, boundary conditions for PM 10 were approximated as [PM 10 ] = [PM 2.5 ] × 1.5.
The entrainment of O 3 and PM 2.5 from the regional background into the model domain and their effect on the concentrations inside the domain was studied with a numerical experiment using the model setup for Hamburg as described in Sect. 4.1.1. A constant concentration offset (BCON offset) was added to the hourly CMAQ concentrations at the lateral and vertical boundaries. In a series of test runs, the BCON offset of O 3 was increased from 0 to 60 µg m −3 in steps of 10 µg m −3 , and the BCON offset of PM 2.5 was increased from 0 to 30 µg m −3 in steps of 5 µg m −3 between the runs. A linear relationship was found between the monthly mean concentration (July 2012) of O 3 and PM 2.5 in the grid cell in which the inner-city urban background station 13ST is located and the BCON offset (Fig. D1). Fitting a linear regression model of the form y = a + bx to the data gave a slope of 0.66 and 0.93 for O 3 and PM 2.5 , respectively. Since PM 2.5 is treated as a chemical non-reactive tracer in the model, the reason for a slope smaller than 1 is removal by dry and wet deposition within the study domain. For ozone, the addition of an offset to the concentrations at the boundaries does not fully propagate into the grid cell concentration at station 13ST due to removal by dry deposition, photolysis by sunlight and the chemical reaction with NO x emitted in the urban area.

Appendix E: Statistical indicators and model performance indicators
In the statistical analysis of the model performance the following statistical indicators are used: overall bias (Bias), normalized mean bias (NMB), standard deviation (SD), root mean square error (RMSE), correlation coefficient (Corr) and index of agreement (IOA).
The overall bias captures the average deviations between the model and observed data and is defined as follows: where M and O stand for the model and observation results, respectively. The overbars indicate the time average over N time intervals (number of observations). The normalized mean bias is given by The root mean square error combines the magnitudes of the errors in predictions for various times into a single mea- sure and is defined as where subscript i indicates the time step (time of observation values). RMSE is a measure of accuracy to compare prediction errors of different models for particular data and not between datasets, as it is scale dependent. The correlation coefficient (Pearson r) for the temporal correlation is defined as SD M and SD O are the standard deviation of the model and observation data, respectively. The standard deviations are The index of agreement is defined as An IOA value close to 1 indicates agreement between modelled and observed data. The denominator in Eq. (E6) is referred to as the potential error.
The model performance criterion (MPC) for dispersion models is the minimum level of quality that has to be achieved for use in policy support related to AQ regulations. The MPC implemented in the FAIRMODE DELTA tool has been constructed on the basis of the observation uncertainty (Thunis et al., 2012a).
The uncertainty of a single observation value U 95 (O i ) is expressed as where u RV r represents the relative measurement uncertainty estimated around a reference value, RV, for a given time averaging, e.g. the hourly or daily limit values of the Air Quality Directive (AQD). The fraction of uncertainty around the RV is given by α 2 . Most commonly, the expanded uncertainty is scaled by using a value of 2 for the coverage factor, k, to achieve a level of confidence of approximately 95 %.
The root mean square of the observation uncertainty (RMS U ) is then A model quality indicator (MQI) is defined as the ratio between the model-observation bias and a quantity proportional to the observation uncertainty as with β = 2 in the DELTA tool. The model quality objective (MQO) is fulfilled when the MQI is less than or to equal 1.
A characteristic of the MQI is that errors in Bias, SD and Corr are condensed into a single indicator value, as follows: The MQI as described by Eq. (E10) is used as the main indicator in the target diagram (Thunis et al., 2012a). In the normalized target diagram, it represents the distance between the origin and a given station point. The normalized bias (first term on the right-hand side of Eq. E11) is used for the y axis, while the centred root mean square error (CRMSE) (sum of the two last terms on the right-hand side of Eq. E11) is used to define the x axis. More details on the normalized target diagram can be found in Thunis et al. (2012a).
Author contributions. MK was mainly responsible for the development of the CityChem extension, development of the research questions, most of the writing, evaluating air concentration data, and creating a framework for data processing, visualization and plotting. SEW drafted the overall structure of the paper, developed EMEP45 and WMPP, which became part of the CityChem extension, and implemented TWOSTEP for solving photochemistry on the Eulerian grid. SS contributed to the development of EMEP45. MOPR prepared input datasets for the air quality study of Hamburg, performed TAPM simulations to produce air concentration data and meteorological input data, evaluated meteorological data from TAPM, tested early versions of the EPISODE-CityChem model, and assisted with the setup and use of the DELTA tool. All co-authors contributed to the writing of the paper and discussion of the model results.
Competing interests. The authors declare that they have no conflict of interest.