Next Article in Journal
GIS-Based Integrated Multi-Hazard Vulnerability Assessment in Makedonska Kamenica Municipality, North Macedonia
Next Article in Special Issue
Spatial and Temporal Patterns of Trace Element Deposition in Urban Thessaloniki: A Syntrichia Moss Biomonitoring Study
Previous Article in Journal
Optimization and Application of Analytical Assays for the Determination of Oxidative Potential of Outdoor and Indoor Particulate Matter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Emission Factors for Particulate Matter and NO2 from Road Transport in Sofia, Bulgaria

1
Faculty of Physics, Department of Meteorology and Geophysics, Sofia University “St. K. Ohridski”, 1164 Sofia, Bulgaria
2
National Institute of Geophysics, Geodesy and Geography-Bulgarian Academy of Sciences (NIGGG-BAS), 1113 Sofia, Bulgaria
3
GATE Institute, Sofia University “St. Kliment Ohridski”, 1113 Sofia, Bulgaria
4
Department of Urban Planning, Faculty of Architecture, University of Architecture Civil Engineering and Geodesy, 1164 Sofia, Bulgaria
5
Environmental Health Division, Research Institute at Medical University of Plovdiv, Medical University of Plovdiv, 4002 Plovdiv, Bulgaria
6
Department Mathematics, Faculty of Transportation Engineering, University of Architecture Civil Engineering and Geodesy, 1164 Sofia, Bulgaria
7
Department of Meteorology, National Institute of Meteorology and Hydrology, 1784 Sofia, Bulgaria
8
Department of Measurements, Metrology, and IT, National Institute of Meteorology and Hydrology, 1784 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
Atmosphere 2024, 15(7), 773; https://doi.org/10.3390/atmos15070773
Submission received: 21 May 2024 / Revised: 21 June 2024 / Accepted: 25 June 2024 / Published: 28 June 2024

Abstract

:
Traffic-related air pollution has a significant impact on the concentration of particulate matter (PM) and nitrogen oxides (NOx) in urban areas, but there are many uncertainties associated with the modeling of PM concentration due to non-exhaust emissions. Bulgarian weather, road surfaces and traffic conditions differ significantly from the UK’s and other EU countries’ averages, which underpin many assumptions in established models. The hypothesis is that the emission factors differ from those used to calculate traffic emissions using the EMIT model. The objective of this work is to adjust the emissions for PM and the relationship between the fractions of NOx and PM using the hourly mean concentrations from road transport and urban background automatic air quality stations in Sofia, Bulgaria. Various already-published and newly developed methods are applied to local observations to derive functions and relations that better represent Bulgarian road and traffic conditions. The ADMS-Urban model is validated and evaluated by comparing pollutant concentrations from simulations using original and adjusted emissions, showing an improvement in results after applying functions and relationships derived from local observations. This work is part of our efforts to improve air quality modeling in urban areas in Bulgaria.

1. Introduction

Traffic-related air pollution (TRAP) is formed by a complex mixture of gases and particulate matter (PM), and it has a significant contribution to urban smog. Vehicle emissions can be categorized into three groups: exhaust, non-exhaust and evaporative. Exhaust emissions from road transport arise from the combustion of fuels such as gasoline, diesel, liquefied petroleum gas (LPG), and compressed natural gas (CNG) in internal combustion engines. Motor vehicles emit a variety of pollutants, including nitrogen oxides (NOx), elemental carbon (EC), PM not exceeding 2.5 μm in aerodynamic diameter (PM2.5), ultrafine particles (UFPs) not exceeding 1 μm in aerodynamic diameter, heavy metals, polycyclic aromatic hydrocarbons (PAHs), and volatile organic compounds (VOCs). Vehicles also produce non-exhaust emissions, with most important sources of airborne PM not exceeding 10 μm in aerodynamic diameter (PM10) likely to be due to the re-suspension of road dust, the abrasion of the road surface and the wear of brakes and tires. Evaporative emissions are a result of vapors escaping from the vehicle’s fuel system. Evaporative emissions are important only for VOCs. Pollutant generation from vehicle emissions depends on the type of vehicle (e.g., light- or heavy-duty vehicles) and age, operating and maintenance conditions, exhaust treatment, type and quality of fuel, wear of parts (e.g., tires and brakes), and engine lubricants used. A number of harmful air pollutants are emitted directly from vehicles, such as “primary” PM and NOx. Others, such as ozone (O3) and “secondary” PM, are formed in the atmosphere from precursors such as NOx and VOCs.
In view of the increasing motorization in the European Union [1,2], which results in growing emissions of exhaust pollutants, coordinated action is being taken across the member states to reduce the harmful effects of vehicle exhausts. One of the key steps in reducing exhaust emissions was the introduction of the EURO standards [3]. Currently, these standards cover NOx, total hydrocarbons (THCs), carbon monoxide (CO), and PM emissions [4]. As a result of the introduction of such air quality regulations and improvements in vehicular emission control technologies, exhaust emissions from motor vehicles, along with the ambient concentrations of most traffic-related pollutants, have decreased steadily over the last several decades in the EU in general and in many of the high-income countries [5]. For instance, NOx and PM2.5 emissions for countries with similar population sizes, like Austria, Denmark, Finland or Portugal, have diminished [6]. Medium-income countries of similar population size such as Bulgaria and other Eastern European countries (e.g., Greece, Hungary, Czechia or Latvia) have had a modest to poor decrease in the generation of exhaust and wear-related non-exhaust emissions. The most probable reason for this is the motorization profile. This profile is characterized by three distinctive features: a higher share of older, predominantly second-hand transport fleets [7], with estimates for Bulgaria of a ratio of around 1 to 10 between new and newly imported second-hand passenger cars for the decade 2010–2020 [8]; increasing motorization rates; and a shift from smaller and older towards newer but bigger vehicles, most often Sports Utility Vehicles (SUVs). Additional national and local fleet data and tendencies, providing background to the current experiment, and data comparing vehicle types, are described in inventory-oriented studies [9] and annexes to a report [10].
The quantification of motor-vehicle emissions is critical in estimating their impact on local air quality and traffic-related exposures and requires the collection of travel-activity data over space and time and the development of emission inventories. Emission inventories are developed based on complex emission models, recently described in [11]. The review paper presents the current state-of-the-art in the field of exhaust emission models and traffic simulation. Emission models could be divided according to their precision scale into macro (regional, national), meso (local), and micro (areas of a dedicated part of a city, road intersections) models [12].
Macroscopic models are mainly based on the number of vehicles; their category, load, average speed, and distance traveled; and weather conditions, among other factors [13]. The most widely used macro-emission models are the European Computer Program to Calculate Emissions from Road Traffic (COPERT) [14] developed by the Joint Research Centre; the Mobile Source Emissions Factor (MOBILE) model, which has been superseded by the Motor Vehicle Emission Simulator (MOVES) [15], both maintained by the U.S. Environmental Protection Agency (EPA); the Handbook of Emission Factors for Road Transport (HBEFA) [16], the product of a common effort by funding agencies and development partners in six countries; and the Comprehensive Emissions Inventory Toolkit (EMIT) [17], supported by the Cambridge Environmental Research Consultants Ltd. (CERC), Cambridge, United Kingdom.
Microscale models need a large amount of data related to the measurements of vehicle parameters, such as acceleration and speed, as well as road parameters such as terrain gradients and position coordinates [18]. The most popular micro-emission models are the Comprehensive Modal Emission Model (CMEM) [19], maintained by the Centre for Environmental Research and Technology; services such as the Environmentally Sensitive Traffic Management (ESTM) [20] solution, developed by Bosch Mobility; and fleet clustering and real consumption models, supported by the Swiss Federal Laboratories for Materials Testing and Research EMPA [21].
The quality of the travel-activity data (such as vehicle miles traveled, number of trips, and types of vehicles) and the complex algorithms used to derive the emission factors suggest the presence of substantial uncertainties and limitations in the resulting emission estimates [22]. In addition to this, it should be noted that estimates of PM emissions have had very limited field valuation and verification [23].
A lot of uncertainty exists in TRAP modeling, especially for PM concentration, due to non-exhaust emissions [24]. An evaluation of existing models for PM emissions from abrasion sources, such as tires, brakes, and road surface wear and resuspension was conducted by Boulter et al. [25] to develop improved prediction methods. The study shows that there are only a few detailed methodologies for predicting emissions of PM from non-exhaust sources and that there is a need for more extensive empirical data. It is also clear from the study that significant quantitative insights can be gained by analyzing measured data from heavily instrumented traffic monitoring sites, and a number of general recommendations for methods of obtaining such data were given.
There are many gaps in the official ground air quality data in Bulgaria. The National Automated System for Environmental Monitoring [26] is operated by the Executive Environment Agency, Ministry of Environment and Water. Only 34 automatic measuring stations (including 4 in forest ecosystems), 5 differential optical absorption spectroscopy systems and 9 points with manual sampling and subsequent laboratory analysis were operated during the first 6 months of 2023 (the last available report) [27]. Only 10 of these stations measure PM2.5 and they are mostly located in cities. Sofia has five automatic and one manual sampling measuring station (only one measures PM2.5), and there are fewer (1–2) in other cities. Another 22 measuring sensor stations are operated by Sofia Municipality [28], and a network of low-cost citizen-owned sensor suites [29] are distributed around the city—both sources of data need additional verification and regular calibration.
Concentrations of PM10 and PM2.5 (and polycyclic aromatic hydrocarbons, PAHs) are elevated in most municipalities in Bulgaria, where air quality stations are deployed. For instance, high hourly average values are sporadically registered in big cities such as Sofia, Plovdiv, Pleven and Varna, and for the period 2015–2019, both the annual average concentration and the allowed number of annual exceedances (18 times in a calendar year) of the hourly limit values of NO2 in the city of Plovdiv were exceeded a number of times [27,30]. In spite of official monitoring not showing a significant problem with NO2 concentration, independent diffusion tube measurements conducted by non-profit civil society organizations show high monthly averaged values [31].
Similar problems with high PM concentrations are registered in other EU countries such as Poland [32] and southeastern European countries that exclusively make up the list of the top five countries with the highest ambient PM2.5 levels in Europe: North Macedonia has the highest exposure in all of Europe, followed by Bosnia and Herzegovina, Serbia, and Montenegro [33]. The seasonal pattern of daily PM10 concentration, with a maximum in the cold period of the year (October to March), is typical for Sofia [34]. This is due to variations in meteorological conditions and the increase in the number and strength of sources. Minimum values are observed during the summer months when intense turbulent mixing and a high atmospheric boundary layer (ABL) results in lower measured concentrations. It is important to note that on a diurnal scale, the pollutant, in particular PM10, mixes well at the urban ABL.
In 2017, the European Court of Justice ruled against Bulgaria for exceeding PM10 (daily and annual) limit values. Bulgaria is facing financial penalties as of December 2020, when the European Commission referred the country to the European Court of Justice for not implementing the 2017 court decision [35]. Additionally, in 2021, two long standing legal cases brought forward by citizens and nongovernmental organizations came to resolution in Plovdiv and Sofia. The court rulings set a 1-year deadline on both cities to achieve air quality norms and implement concrete measures to protect the health of citizens.
In search of a resolution, two atmospheric air quality management programs for Sofia were developed: one for the period 2015–2020, for PM [36,37], and another one for the period 2021–2026, which covers all harmful pollutants [38]. The preliminary findings of the two programs show that there is a significant difference in PM emissions for both periods. From the currently existing emission inventories for Sofia, it is apparent that there is a lack of informative and accessible data on the dynamics and structure of the traffic flow on the major road arteries, and this results in biassed modeling, which in turn may mislead policy- and decision-makers.
Air quality modeling can alleviate some of the deficiencies arising from the lack of official measurement and communication campaigns, but it requires a careful consideration of underlying assumptions and pre-set parameters in computational models. In a previous study [39], it was shown that there is a clear contrast in pollution mechanisms in different parts of Sofia. Transport has a substantial impact on PM concentration in the central parts of the city, where the pollution is mostly local, in proximity to the source, while the contribution of domestic heating from wood and coal burning is highest on the outskirts and in the small towns in Sofia Municipality. These local phenomena must be accounted for. In the same study, some of the deficiencies occurring when conducting an evaluation of the complex pollution processes in urban areas were defined, and concrete steps towards improving the modeling at the urban scale in Bulgaria were formulated. An example of such a step is the development of a high-quality traffic emissions inventory for Sofia [9,40]. Using a tailor-made inventory reduces the uncertainties in existing inventories and helps improve the applicability of established models to the local setting.
The atmospheric dispersion modeling system (ADMS-Urban), which is a comprehensive software suite for modeling air quality in urban areas and close to motorways, roads, and large industrial areas [41,42], was used in this study. The same model was employed in the work of Dimitrova and Velizarova [39], but the present study uses a new high-quality traffic emissions inventory, advanced street canyon modeling, and background concentration, and includes chemical transformations and high-resolution topography description. Despite some deficits of the local Gaussian-type models for air quality studies, the potential of resolving emissions and the physical environment at a fine scale makes them a very useful tool. However, the depiction of the spatial concentration pattern of regional air quality models underrepresents pollution (especially for PM) [43] due to the restriction of their spatial resolution, as they are unable to capture the high concentration gradients nearby roads with heavy traffic and from other sources near the surface. In high-fidelity modeling, accounting for the effect of line sources (e.g., major road arteries with heavy traffic) requires the usage of small spatial steps in the order of tens of meters, enabling us to register the large gradients in proximity of the segment. Pollution concentration decreases rapidly with an increase in the distance from the road, even if the effect of the street canyon is not taken into account for the simulations. It is possible to add further receptor points within the canyon in order to more accurately detect the effect of pollutant retention within the canyon under certain conditions, which is not achievable with the regional models.
Weather, road surfaces and traffic conditions in Bulgaria can also differ significantly from the UK’s and other EU countries’ averages, which underpins many assumptions in established models.
Normal vehicle ageing and wear is likely to be an important factor in transport-related emissions. In Bulgaria, the average age of cars is 16 years old, with almost 70% of all vehicles having a standard of performance below Euro 4, and more than half (56.5%) with no Euro category at all, according to recent investigations [38]. In comparison, the average age of passenger cars is 11.5 years old in the EU, 8 years old in the UK, and 8.6 years old in Switzerland [44]. Fuel type is also important. In [45], the authors show that there is little evidence of NOx emissions reduction from all types of diesel vehicles over the past 15–20 years, and only for petrol passenger cars (including hybrids) was there any evidence for effective NOx control. In this setting, 43.8% of cars in Bulgaria are powered by diesel engines, compared to 50.1% for petrol powertrains [38] (a rate that is slightly higher in comparison to the average split in EU countries) 42.3% for diesel and 52.9% for petrol [44]. There has been recent growth in the demand for electric vehicles in the country, but considering the relatively low purchasing power of the population, there are certain limits to replacing the car fleet with zero-exhaust-emission vehicles.
Road surface wear and dust deposits also vary by type of road and vehicle speed. The typical speed limit in Bulgarian cities is 50 km/h, but on higher-class transit roads, the limit reaches 70 km/h or 80 km/h, which increases dust emissions from the road. In addition, the intensive traffic flow through the central and densely occupied built and intersected parts of the city increases the occurrence and severity of congestions, which further deteriorates air quality. After a significant loss of passengers in the previous two decades, public transportation is slowly recovering part of the lost ridership, with a significant constraint experienced during the pandemic of COVID-19. Coverage by bicycle lanes is much lower than in major European cities, with the existing network of cycle routes being insufficiently integrated, uncomfortable, and sometimes unsafe.
Considering the societal and technological challenges outlined so far, this work has two main contributions:
  • To provide a means of adjusting the emissions of PM and the rates between NO2 and NOx gaseous pollutants, using hourly mean concentration measurements from road transport and urban background automatic air quality stations (AQSs) in Sofia, Bulgaria. Different already-published and new methods are explored and evaluated.
  • To estimate the contribution of PM10 traffic emissions. To identify the contribution of the main groups of PM10 sources in Sofia, the results of a receptor-oriented positive matrix factorization (PMF) analysis for a wide variety of chemical elements [34] are presented. This study covers one year (2019–2020), but selected periods in April and September 2019 are used to validate the newly developed PM10 traffic emissions. These periods were chosen considering several important factors—the absence of domestic heating and additional background sources such as dust transport from the Sahara Desert, lack of precipitation and low wind speed.
The remainder of the paper is arranged as follows: the materials and methods are described in Section 2. Section 3 presents results from experimental adjustment and modeling, and an overall discussion and conclusions follow in Section 4 and Section 5.

2. Materials and Methods

2.1. Research Question and Hypothesis

The working hypothesis is that modeling at a local scale can be improved by implementing a more reliable inventory of traffic emissions if the relationships used to develop the emission factors are derived using local data that is more representative of Bulgarian conditions. To help test the hypothesis, three main factors are considered: data availability and gaps; the application of different methods to process sparse data; and testing necessary assumptions, conducting experiments, and model validation.
Data at all three traffic AQSs in Sofia for a 12-year period (from 2009 to 2020) were used for local NO2 concentration adjustment, and data for the same period at AQS Pavlovo (traffic station) and AQS Hipodruma (urban background station) were used to estimate the proportional factors between PM10 and the coarse fraction of PM (PM10-PM2.5) ( P M 10 P M 2.5 ). For dispersion modeling and model validation against PM10 obtained through PMF analysis, selected periods in April and September 2019 were used. A conceptual workflow diagram of the study, including the input data, models used and experimental adjustment of the TRAP emissions in Sofia, Bulgaria, is shown in Figure 1.

2.2. Data

2.2.1. Air Quality Data

Hourly measurements for meteorological factors and pollution concentrations are available from the National Environmental Monitoring System of the Executive Environment Agency (EEA) at the Ministry of Environment and Water at six AQSs for the period from 2009 to 2020. Three of the sites located in Druzhba, Nadezhda and Hipodruma are classified as urban background stations, and the other two, in Pavlovo and Orlov Most, are classified as transport stations. The AQS at Orlov Most was moved to a new location in Mladost in November 2015, and as a result, the data sets from both sites are shorter.
Quartile charts for NO2 concentration from the traffic stations are shown in Figure 2. A direct comparison with the established norms for the levels of concentration of NO2 related to the protection of human health, set by European directives and implemented into national legislation [46], can be made on the basis of the average annual limit of 40 µg/m3 and the threshold of an hourly average rate of 200 μg/m3, which may not be exceeded more than 18 times in a calendar year. Concentrations above the hourly average rate are registered at AQS Pavlovo during the first years from 2009 to 2011, and there were eight exceedances at AQS Orlov Most in 2009 (the number of exceedances is below the threshold of 18 incidents within one year). The annual average concentration is close to or above the limit for all sites. These data are used for the calculations described in Section 3.1.1.
Quartile charts for the PM10 concentration from the traffic stations are shown in Figure 3. The annual limit for PM10 is 40 µg/m3, and the threshold of a 24 h average rate is 50 μg/m3, which may not be exceeded more than 35 times in a calendar year. The annual concentration limit is exceeded until 2014 and slowly decreases after. The number of days with PM10 concentrations above the daily limit exceeded 35 for all years. PM2.5 concentration measurements are available only at AQS Hipodruma (an urban background station). These data are used for the calculations described in Section 3.1.2.
Road concentrations strongly depend on transport emissions, which vary with the time of day, the day of the week and the season. Average NO2 and PM10 concentration profiles for weekdays, Saturday, and Sunday (based on the entire 12-year period) for different traffic stations are presented in Figure 4. A diurnal variation shows a well-established binary profile with two peaks in the morning and evening hours. The dynamic is similar for all traffic stations, with the highest concentrations at AQS Orlov Most. During the weekdays (Figure 4a), the first peak concentration is at 9–11, with the second at 19–21. There is a 2 h shift in the evening peak at AQS Mladost, probably due to the location of the station in comparison with the other two. AQS Pavlovo and AQS Orlov Most are located near main transport roads, which connect different parts of the city with the city center. AQS Mladost is close to a big road with somewhat intense traffic that mainly takes local flow, joining one of the major transit traffic axes in the city “Tsarigradsko shoes”, which is also relatively close but at a lower elevation. The morning peak is smoothed on Saturdays (Figure 4b) and almost disappears on Sundays (Figure 4c).

2.2.2. Meteorological Data

The meteorological data for the studied period were provided by the National Institute of Meteorology and Hydrology (NIMH). Both soundings and surface measurements are from the CMO—an experimental base of the NIMH in Sofia and part of the official network of the World Meteorological Organization. The surface observations collect air temperature (T) and relative humidity (f) at 2 m and wind speed (WS) and wind direction (WD) at 10 m, measured with an Vaisala AWS310 [47] meteorological station. The meteorological observations are checked, cleaned, and averaged at 10 min intervals for the entire study period.
Hourly averaged values were used as an input to the ADMS-Urban model, shown in Figure 5. The amplitudes in the fields of T and f are more significant in April compared to September. The temperature is near zero in the early morning hours, increasing to above 20 °C in the afternoon, with a night-time rise at the end of the period (23–26 September). The f decreases during the day and increases above 90% at night for some days, but it never reaches saturation levels (an indication of rain or fog). The wind rose for selected days in April and September, as shown in Figure 6. The prevailing WD is from the southeast, with a low speed of less than 3 m/s.
Aerological soundings were conducted with the Vaisala MW41 [48] sounding system and RS-41SG [49] aerosondes and include data for T, f, WS and WD. The data set consists of midday soundings for periods in April and morning and midday sounding for periods in September. Vertical profiles for days with extremal discrepancies between modeled and observed PM10 concentrations for the two study periods are shown in Figure 6. The largest discrepancy between modeled and observed concentration values is characterized by zero wind velocity (Figure 7a–c). A strong inversion layer is observed in all potential temperature profiles at 06:00 UTC (Figure 7b,e), and a mixed layer up to 1000 m above sea level (ASL) is observed at 12:00 UTC (Figure 7a,c,d,f). Relative humidity is below 80% on all occasions except for the morning hours on the 24, 25 and 26 September. All available profiles are presented in Supplementary Materials S1.

2.3. Modeling

2.3.1. Transport Emission Inventory Modeling

The emissions data for this study were provided by a new traffic emission inventory for the city of Sofia, developed by the authors. Emissions are calculated with EMIT using the activity data described below. To calculate exhaust emissions, data for 2010 were used; from the 2014 urban factor data set, part of the UK National Atmospheric Emissions Inventory (NAEI) was used [50]. This particular year was chosen because it best represents the fleet composition in Sofia, in terms of its average age, for the study period. The fleet approximation is based on comparison between the average age of the fleet, which in Bulgaria is around 8 to 10 years older than that in the UK. In the EMIT input data, most of the vehicles manufactured prior to the introduction of the Euro 4 standard are assumed to lack properly working catalyzers and filters, as is the case observed in practice, although there is no conclusive data available about this well-known and widely discussed phenomenon in the country. It must be noted that this approach only offers a proxy, because vehicles with first registration circa 2010 would have experienced some additional age-related degradation by the time of this study. Thus, the expected emission factors would be higher for all categories, except for newer vehicles, which hold a relatively small share of the fleet.
One of the biggest challenges is the quality of basic, publicly available traffic and fleet composition data. The data come from a variety of sources with poor accessibility, and in some cases, there is questionable integrity and many methodological inconsistencies in the data gathering approaches [9]. This is a typical administrative practice leading to fragmentation, stemming from sectorized institutional silos. In an attempt to overcome this fragmentation, the approach proposed in [40] is applied in this work. In summary, this approach attempts to fill the gaps in the data for traffic along primary and secondary streets by first collecting data for a large number of related features such as road and street categories, population density, functional analysis, spatial syntax, previous traffic, noise and NO2 and then applying ensemble learning to generate insight.
After a series of consecutive data imputations, the annual average daily traffic distribution in the street and road network of the city of Sofia and the metropolitan municipality for 2018 and 2022 (the years with more comprehensive data availability) was estimated. The traffic data modeled for 2018 were used to calculate the exhaust and non-exhaust traffic emissions using EMIT. The chosen year is the previous closest to and most in tune with the previously described meteorological and pollution data from 2019. This period of two years is also just before the big shifts in mobility behavior around the period marked by the COVID-19 pandemic lockdowns and recoveries.
The composition of the vehicle fleet was divided into three groups—heavy, light and motorcycles—according to the available activity data. The non-exhaust emission factor data sets were taken from the European Monitoring and Evaluation Programme (EMEP) [51] for tire, road, and brake wear, and from the Department for Environment, Food and Rural Affairs/Transport Research Laboratory (DEFRA/TRL) for resuspension [25]. Detailed information about the emission factor data sets can be found in Appendix A of the EMIT User Guide [17].
The used approach also implements new route-type data following the NAEI2014 Urban 2014 version 2 detailed taxonomy available in EMIT. It is based on multi-dimensional stratification including fuel types, vehicle types, engine volume and vehicle weight categories, types of EURO standards, and types of cleaning technologies, among others [33]. These data provide the link between the traffic count data/fleet components and vehicle sub-categories with the emission factor. Active traffic data are shown in Figure 8.

2.3.2. Atmospheric Dispersion Modeling in Urban Street Canyons

ADMS-Urban is a second-generation Gaussian dispersion model with an integrated street canyon model that considers the impacts of street canyons on dispersion, turbulence and mixing induced by traffic. It incorporates the latest understanding of the boundary layer structure, using advanced algorithms for the height-dependence of wind speed (WS), turbulence, and stability. ADMS-Urban can model point, area, line, road, volume and grid sources. The model integrates with geographical information system (GIS) tools and emissions databases and supports calculation of concentrations based on a range of averaging periods, including percentiles. The GIS interface is an important visualization tool and provides a means of manipulating the spatial element of the model input data. ADMS-Urban is developed and maintained by the Cambridge Environmental Research Consultants (CERC) Ltd. More details regarding the model are available in [52].
Concentrations of NO2, with a surface lifetime of around 1 day, are primarily related to the dispersion and the chemical transformation of local emissions. The NOx chemical reactions take place over a relatively short time, and in order to calculate NO2 concentrations, NOx chemistry needs to be taken into account. ADMS-Urban has an integral chemistry model, using what is known as the “generic reaction set” [53], which includes the reactions between ozone (O3), VOCs and the oxidation reactions of sulphur dioxide (SO2) that lead to the formation of PM10 and PM2.5 of ammonium sulphate. The background concentrations in this study were taken from the Copernicus Atmosphere Monitoring Service (CAMS) reanalysis [54], which are required by the model to enable the use of the chemical reaction module. Only the dry deposition option was used due to a lack of precipitation during the selected simulated time periods.
The advanced street canyon module was employed in this study. This module provides a more detailed description of the dispersion within the canyon, due to the consideration of a wider range of canyon geometries, including the effects of tall canyons and asymmetries [55]. The model was also used in “network” mode, which allows for the passing of pollutants from the end of one street canyon into another canyon connected to it. The advanced street canyon module requires the specification of additional input parameters, which were obtained via the street canyon tool [55]—a Python plug-in that is designed to export street canyon geometries from within ArcGIS.
By using the time-varying emission factors option, different diurnal profiles for each day of the week and monthly profiles for each source were applied. The temporal profiles were taken from the Dutch Organization for Applied Scientific Research (TNO) [56].
The interested reader may refer to the ADMS—Urban v5.0 user guide [52] and the technical specifications [57] for more information regarding the model settings used in this study.
The modeling domain over Sofia and the locations of the AQSs (where the measurements used in this study were taken from) is shown in Figure 9a. The smaller area in the green rectangle magnified in Figure 9b shows the modeling domain near AQS Mladost employed for emission factors adjustment. Two different sites are used for model validation: AQS Mladost for NO2 and the central meteorological observatory (CMO) in Sofia for PM10.

2.4. Experimental Design of Emissions Adjustment

2.4.1. Study Approach for NO2 Concentration Adjustment

By default, if NO2 emissions are not included explicitly, ADMS-Urban assumes that a certain percentage of the NOx emissions are in fact due to NO2. This percentage is modeled as an empirical relationship between hourly NO2 and NOx based on a third-degree polynomial produced by Derwent and Middleton [58]. To develop the polynomial, the authors of that paper analyzed hourly concentrations of NOx, sorting the data into 10 bins of increasing values in parts per billion (ppb). The corresponding NO and NO2 concentration values were indexed according to the sorted NOx data. The NO2 concentration values corresponding to each bin were then averaged, and a curve was interpolated between NO2 mean concentrations and upper NOx bin limits. The relationship was derived as
[NO2] = 2.166 − [NOx] (1.236 − 3.348A10 + 1.933A102 − 0.326A103),
where the square brackets indicate the hourly average concentration in ppb, and A10 = log([NOx]). This function has limits in the range 9.0 ppb < [NOx] < 1141.5 ppb. Below 9.0 ppb of NOx, the efficiency of [NO2]/[NOx] ratio was limited to 0.723. Above 1141.5 ppb of NOx, the [NO2]/[NOx] ratio was limited to 0.25. The derived polynomial function is based on one year (May 1991 to June 1992) of data collected at a fixed monitoring location in London.
A follow-on study [59], based on a larger amount of data (1991–1997), investigated the validity of new empirical relationships at several sites. The point of significant difference from previous studies was to consider the ratio between NO2 and NOx concentrations as a dimensionless yield quantity: Y = [NO2]/[NOx], where [NOx] = [NO] + [NO2] in ppb, and 0 ≤ Y ≤ 1 is dimensionless by construction. The data were sorted in 10 ppb bins by NOx concentration, and a fourth-order polynomial was fit to the log-transformed upper NOx bin limits (independent variable) and average NO2 values corresponding to the NOx bins. The upper bin limits were used in defining the yield to stabilize the result by avoiding division by small values. The least-squares polynomial fit to the processed data, recommended for urban areas outside London, was reported as
Y2 = −3.083 + 7.472 A10 − 5.116A102 + 1.382A103 − 0.129A104,
The polynomials used in both cases are highly contextual, as they have been fitted on data reflecting the pollution mechanisms in the UK. When modeling pollution dispersion in Sofia, the relationship between NOx and NO2 needs to be adjusted based on local dynamics. To do this, a similar method is used to obtain an empirical relationship between NO2 and NOx using the hourly mean concentration measured at the traffic stations in Sofia. The NOx data are split into 10 ppb bins and log-transformed. The corresponding NO2 concentrations are averaged for each bin, and the bin yield is calculated as the ratio of the average NO2 to the upper NOx bin limit. It is important to note that, due to their fixed end points, the bins contain a variable number of values, so the average NO2 concentration will be a variable-quality estimate of the true mean. Even though, in principle, the standard error of the mean can be computed for each bin and an interval regression can be used to learn the relationship between NO2 and NOx (see, e.g., [60]), a regular point-value regression is used here, with some extensions considered as part of future efforts. Hourly concentrations for NO2 do not order in the same way as NOx values due to a number of variables that influence the formation of NO2, such as turbulent mixing, temperature, light intensity, and the presence of ozone, among others.

2.4.2. Study Approach for PM Emissions Adjustment

The primary emitters of particulate pollution—industry and traffic—generate PM predominantly in the fine particle size range (PM2.5), while emissions from activities such as construction work and those from natural sources, such as pollen, contribute predominantly to coarse particle concentrations (particles greater than 2.5 µm in aerodynamic diameter). As discussed in Section 3.1.2. a major contribution arises from surface soil and dust on paved areas becoming resuspended due to the wind flow or due to turbulence induced by passing road traffic. This makes the estimation of resuspension contributions particularly important for completing a traffic emissions inventory. However, estimating the non-exhaust sources, including the contribution of resuspension, is very challenging due to the lack of direct measurements and various impracticalities related to the necessary road measurements. To address this issue, a methodology to calculate the total traffic emissions of PM10 was developed by [25].
The total emissions attributable to the non-exhaust sources can be expressed as the sum of the abrasion emissions (brake, tire, and road wear) and a resuspended component:
E total = E tire + E brake + E road + E resusp ,
Emissions for tire, brake and road surface wear can be estimated from EMEP emission factors and traffic fleet data. Coupled with an estimate of total PM, Equation (3) can be rearranged to obtain an estimate for resuspension. The approach taken to estimate the contribution of resuspension relies upon several fundamental assumptions. Firstly, it was assumed that PM10 and NOx behave similarly in the atmosphere and therefore the roadside incremental values of the two, which are calculated as the difference in concentration measurements between a traffic station and an urban background station [61,62], are closely related. This enables the estimation of the total PM10 emissions from traffic ( E P M 10 ) based on roadside increments of PM10 (∆PM10) and NOx (∆NOx) and calculated NOx emissions ( E N O x ), as in the following equation:
E PM 10 = E NO x Δ PM 10 Δ NO x
A second assumption in the initial calculations is that the roadside increment of PM2.5 was solely attributable to vehicular exhaust sources and that non-exhaust emissions were largely confined to the coarse particles fraction ( PM coarse = PM 10 PM 2.5 ). Apportioning total PM10 emissions into PM2.5 and PMcoarse fractions according to the ratios between their observed concentrations allows a value for total PMcoarse (i.e., non-exhaust PM) to be calculated as
E PM coarse = β E PM 10
where β can be calculated as the slope of the line of best fit between concentration measurements for PMcoarse and PM10.
Finally, the following relationships for the PM10 and PM2.5 components of the road traffic emissions are used in ADMS-Urban:
E PM coarse = E tire , coarse + E brake , coarse + E road , coarse + E resusp
E PM 2.5 = E exhaust + E tyre , 2.5 + E brake , 2.5 + E road , 2.5
where E x , coarse are the coarse parts of the respective emissions, E x , 2.5 are their fine fraction parts and E PM 2.5 is the total PM2.5 emissions.

2.4.3. Source Apportionment Technique

The majority of non-exhaust PM derives from the resuspension of material already deposited on the road surface (between wheel tracks, on the kerbside or in the pores of the asphalt) due to tire shear, vehicle-generated turbulence and the action of the wind [63]. There are several sources that contribute to road dust, such as particles from brakes, tires, and road wear; particles deposited from exhaust emissions; and particles from nearby environments, including fugitive loading from construction sites, pavements, unpaved roads, dry and wet deposition from the atmosphere, salt residue during freezing periods, traction sand, and the deposition of pollen and plant matter [64]. Because of the heterogeneity of dust sources, it is difficult to distinguish between “direct” wear emissions (tire, brake or road wear) and “resuspended” wear emissions, and it is therefore hard to separate their relative contribution to atmospheric PM levels [65].
Source apportionment (SA) is the technique that relates a source emission (an activity sector or an area) to the ambient air concentration of a pollutant. Positive matrix factorization (PMF, also known as non-negative matrix factorization) is the most frequently used factor analysis technique. It solves a weighted factorization problem with non-negativity constraints using known experimental uncertainties as input data, thereby allowing for the individual treatment of matrices. For SA, measured concentrations and their uncertainties are used to solve the mass balance equation X = G × F + E, where X is the chemical composition matrix, G is the source contribution matrix, F is the factor profile matrix, and E is the matrix of residuals. The most recent versions of the EPA-PMF [66] software solve the factorization problem using the conjugate gradient method and contain routines to estimate the optimal number of factors, to test the rotational ambiguity and introduce constraints. The uncertainty and stability of the solution are estimated by bootstrapping and displacement methods [67].
The results from a SA study in Sofia covering one year (January 2019–January 2020) are presented in [34]. The data were analyzed using the EPA-PMF 5.0 model. More than 200 samples were studied, and PM chemical composition and the application of the PMF method led to the identification of eight factors, six of which have the most significant contribution to the PM10 mass concentration in Sofia. Sources of ambient PM10 have been grouped into eight categories: Resuspension (RES), Secondary (SEC), Biomass burning (BB), Traffic (TR), Industry (IND), Nitrate-rich (N), Fuel-oil burning (FUEL), and Mixed SO42− (SO4). The analysis showed that the resuspension factor is the main contributor to the total PM10 mass (25%), followed by BB (23%), SO4 (19%), SEC (16%), TR (9%), IND (4%), N (4%), and FUEL (0.4%) in Sofia.
The “resuspension” category includes contributions from soil and road dust resuspended into the surface air layer from the wind and the mechanically induced resuspension of air particulate from the vehicles. This factor has seasonal variation, with a minimum at the beginning of 2019 (January–February) and a maximum in the middle of December 2019, when stagnant weather conditions with prolonged inversion and the presence of fog led to high PM10 concentrations. The relatively high impact of resuspension (33%) is connected to the dryer periods. In 2019, significantly lower monthly precipitation amounts were measured. After June and in December 2019 and January 2020, the total monthly precipitation was less than 60% of the climatic norm. The traffic factor accounts for 9% of the measured PM10 and has a weak seasonal pattern, with maximums in the cold periods in Sofia (13 μgm−3). More details related to the method, data sampling, analyses of chemical composition in PM10 and selection of different categories based on specific elements can be found in [34].

3. Results

3.1. Results from Experimental Adjustment

3.1.1. Local NO2 Concentration Adjustment

Data at the three traffic AQSs in Sofia for a 12-year period (from 2009 to 2020) were used (see Figure 2) to fit the polynomial regressors. To identify the polynomial that most adequately describes the data, a model selection study was set up with polynomials from the first to the fourth order with all possible combinations of terms. The models are denoted as [i1, i2, …, in], where ij are the powers of monomials having non-zero terms. Thus, for example, [0, 2, 4] denotes y = β 0 + β 2 x 2 + β 4 x 4 . Each polynomial was fit on data from individual stations, as well as data from all stations combined. For the purposes of this study, adequacy is quantified by the coefficient of determination, R2 and the Akaike information criterion (AIC).
Goodness of fit is measured through the adjusted R2, which penalizes the performance score of the model based on its complexity. The AIC is also computed to provide a measure of the complexity of each polynomial with respect to the available training data. In this model selection study, only polynomials with a non-zero intercept were tested due to the nature of the data established in an exploratory data analysis. Note that some polynomials, such as [1, 2, 3] fitted to the data from AQS Orlov Most, also perform well (adjusted R2 = 0.985, AIC = −356.2) without an intercept term, although such cases form a minority. Table 1 shows the results from the model selection investigation.
To verify the value added by increasing the order of the polynomial for each data set, a complexity-oriented F-test was also performed. The F-statistic represents the result of a pairwise hypothesis testing procedure, where the null hypothesis is that the more complex model of the pair does not provide a significantly better fit than the simpler one. The p-value for the test is also reported to help gauge the statistical significance of the test result. The hypothesis testing was limited to full q-degree polynomials, that is, polynomials of degree q, which have non-zero coefficients for all terms of degree q-r with r = {0, 1, 2, …, q}. The results of the F-tests are summarized in Figure 10. Each tile in the figure shows the performance of the polynomials on the data set indicated in the respective title. The main diagonal in each tile contains the adjusted R2 for reference to the overall performance of the model. The upper triangle shows the values of the F-statistic, which are colored in shades of blue, with darker colors signifying more desirable results. The lower triangle contains the p-values for the corresponding F-statistic as the cumulative probability of this statistic under an appropriate F-distribution. The definition of all statistics and metrics are given in Supplementary Materials S2.
The critical significance level was set to 5%, with test results indicated as significant (p-value < 5%) highlighted in green and the rest left in white. It can be quickly seen that despite the increases in R2 for more complex polynomials, the improvements may not be statistically significant to the chosen level and may often be far from it. This is always the case for comparisons between third- and fourth-order polynomials. In one case that did not feature in the comparisons, the [0, 1, 3, 4] model performed better than the [0, 1, 2, 3] one, but these could not be compared in an F-test because both have the same number of coefficients. The minor improvement in R2 (in the order of 10−4) was neglected. In view of these results, a full third-order polynomial was selected for the rest of the work in this paper. To provide a visual assessment of the quality of fit, Figure 11 shows each data set with the corresponding third-degree polynomial superimposed as a black line.
As an additional comparison, the functional maps developed by Derwent and Middleton [58] in Equation (1) and by Dixon et al. [59] in Equation (2) for use in London produce results that severely misrepresent the behavior of NO2 concentration in Sofia. That is, if these expressions are evaluated on the available NOx data, their results are, perhaps unsurprisingly, around three orders of magnitude worse than those from the customized polynomials. Table 2 summarizes the performance, in terms of RMSE, of the two models and compares it to the performance of the models developed specifically for this study. As expected, RMSE scales strongly with sample size, but despite this, the custom polynomials achieve much better results in both the NO2 and yield representations. This serves as a reminder that modeling assumptions must be challenged and relaxed whenever possible to obtain meaningful results. Using ADMS-Urban for other cities will likely require a similar approach.

3.1.2. Local PM Concentration Adjustment

As already mentioned in the introduction, only one AQS in Sofia measures PM2.5, and it is an urban background station. The two traffic AQSs measure only PM10 concentration. The air quality stations at Pavlovo (traffic station) and Hipodruma (background station) were used in this study for two reasons. Firstly, the sites are located in relative proximity to each other (~2.6 km; see Figure 9a), which makes their observation more relevant to the local pollution dynamics. Secondly, all pollutants (except PM2.5 at AQS Pavlovo) and meteorological variables were available at both stations, which becomes useful in estimating PM2.5 concentrations at the traffic station, AQS Pavlovo. Roadside increments of PM10 (∆PM10) and NOx (∆NOx) were calculated using data for 2019. Entries with missing or erroneous values for ∆PM10 and ∆NOx, as well as those with an ∆PM10/∆NOx value larger than 0.1, were removed, corresponding to the data-cleaning steps recommended in [25]. A representative incremental ratio used in Equation (4) was computed as the slope of the zero-intercept line of best fit to the data. Figure 12 depicts the cleaned data and line of best fit. The final ingredient for estimating the total adjusted PM10 emissions in Equation (4) are the total NOX emission (ENOx) for each road source, which were calculated using EMIT.
To estimate the adjusted PMcoarse emissions, using Equation (5), β must first be calculated at the selected traffic station (AQS Pavlovo). In the case of Sofia, this necessitates an estimation of the concentration of PM2.5 at the traffic station, which can be carried out by mapping PM2.5 concentrations from the only AQS that measures them (AQS Hipodruma; background) onto other measured quantities at that station and then applying the map to the traffic station. This is where the second reason for the selection of AQS Pavlovo and AQS Hipodruma—the two measure the same quantities, except PM2.5—plays an important role.
Several linear regression models were developed to serve as a map. Models of the first and second order with all interactions were fit to data for 12 measurands at AQS Hipodruma for the period from 2015 to 2020, due to the lack of data for the majority of measured factors at AQS Pavlovo for the 2009–2014 period, which precludes the modeling efforts for this period. For the first-order model, only temperature was found to be a statistically insignificant factor, whereas a total 25 terms from the second-order model, including the first-order contributions of temperature and relative humidity and second-order contributions of carbon monoxide (CO) and ozone (O3), indicated no statistical significance. However, because all these quantities participate in significant interactions, and the amount of data is sufficient to avoid overfitting (32,859 points factoring missing data), no simplifications to the model were carried out. The two models were validated on one-year worth of data from 2023. The performance of the two models is summarized in Table 3. Due to its better performance statistics, the second-order model was chosen to predict the concentration of PM2.5 at AQS Pavlovo. The results of the correlation tests between predicted and measured values for the test data from 2020 are shown in Figure 13. The two points around (150, 10) in Figure 13b correspond to values for CO, NO2, SO2, PM10, and wind speed and do not seem to correspond to such PM2.5 concentrations in the training set.
Several points need to be mentioned about the model development. First, the approach assumes that the relationship between PM2.5 and other measured quantities will remain the same at traffic stations. There is no empirical evidence to support or refute such an assumption. Second, in contrast to the NO2 work described in Section 3.2, no thorough model selection investigation was conducted. Third, similar to the NO2 study, no error estimation has been included in the modeling process. All of these points are addressed in the outlook for future work.
Once the model was tested, it was used in the calculation of β in Equation (5). The predicted hourly values for PM2.5 were split into years and subtracted from measured PM10 concentrations to determine the concentration of PMcoarse at the traffic station. A zero-intercept line was fit to the PM10 versus PMcoarse data, whose slope is the β coefficient. The values for each year are shown in Table 4.

3.2. Results from Modeling

3.2.1. Transport Emission Inventory Modeling

EMIT [17] is a database tool for storing, manipulating and assessing emissions data from a variety of sources. It can hold data from explicit sources such as major roads, rail and industrial plants. In addition, EMIT can hold data from sources that may be too small to be considered explicitly, and instead, they are treated as average emissions on a regular grid. Source data held in this way are from minor road and commercial and domestic sources. EMIT calculates emissions from source activity data (such as traffic flow, speed and source length for roads, fuel consumption for industrial sources) using up-to-date emission factors. Calculated with EMIT emissions of PM10, PM2.5, NOx, and NO2 using described in Section 2.3.1, data are presented in Figure 14.

3.2.2. ADMS-Urban Model Validation

The model validation is based on comparison of 24 h averaged modeled concentrations of NO2 and PM10 with observations. Measurements for the NO2 concentration at AQS Mladost (42.656′ N, 23.383′ E at 582 m ASL) are used. A comparison of modeled PM10 concentration with observations is performed at the CMO site (42.655′ N, 23.384′ E at 586 m ASL)—the same location where the PM samples were collected and processed using the PMF method (described in Section 2.4.3)—to assess the contribution of traffic emissions. The CMO site is located approximately 200 m southwest of AQS Mladost and 170 m away from the road, while the AQS location is approximately 70 m away from the road, which leads to a small difference between PM10 concentrations measured at both sites.
The output from the EMIT model was used to calculate the concentration fields without applying the methods described in Section 2.4.1 and Section 2.4.2, which were subsequently used to adjust the emissions of NO2, PM10 and PMcoarse. To investigate the contribution of categories related to traffic emissions concentration and to validate the model, two periods without residential heating and Saharan dust intrusion (18–23 April and 20–26 September 2019) were selected. The meteorological conditions for these periods are described in Section 2.2.2.

Model Validation for NO2

The results from two model runs are presented and compared. The NO2 emissions from EMIT are calculated based on the third-degree polynomial developed by Derwent and Middleton [58], referred to as the “original” in this study. As described in Section 2.4.1, a new polynomial was obtained (referred to as “adjusted”), derived from measured concentrations of NO2 and NOx at AQS Mladost, to achieve a more representative relationship between the two fractions. Commonly used statistic measures—mean bias (MB), mean error (ME), normalized (with respect to the measured value) mean bias (NMB), normalized (with respect to the measured value) mean error (NME) and RMSE—are used to estimate the model performance (for definitions of the measures, see Supplementary Materials S2). A comparison of both methods with results from the model output, applying original and adjusted polynomials, is shown in Table 5.
It is difficult to estimate the model performance on a case-by-case basis as the results strongly depend on meteorological conditions. In general, smaller biases and errors are achieved in April in comparison with September. For most days, biases and errors using adjusted polynomial are smaller, except for two days: 21 and 26 September. When using the custom polynomial, RMSE is reduced from 21.78 µg/m3 to 16.67 µg/m3, ME from 20.23 µg/m3 to 13.74 µg/m3 and NME from 80% to 54%.
Although transport is the main contributor to NOx concentration, other sources such as power plants, industrial processes, and the combustion of fossil fuels, wood, and other organic matter can affect the measured concentration at AQS Mladost, depending on the direction of the mean wind flow and the stability of the atmosphere. It must also be noted that the statistical measures were calculated based on a small sample size, which introduces considerable uncertainty.

Model Validation for PM10

The same period was used for validating the model predictions for PM10 concentrations. Observed concentrations of PM10 and those calculated with the PMF model (described in Section 2.4.3) are shown in Table 6, together with the contribution of different categories to the total PM10 concentration (based on data obtained with the PMF model).
All categories were selected based on the chemical composition of PM10 samples, and only RES and TR relate to transport. As the main purpose of this study is to investigate the contribution of transport to air pollution and to provide adjusted emission factors, representative of the conditions in Bulgaria, sources not related to traffic (SEC, BB, IND, N, FUEL and SO4) were not considered.
Numerical results from ADMS-Urban were compared against the sum of the non-exhaust and exhaust traffic sources (RES and TR) for PM10. Model results were extracted at the exact location of CMO, where measurements were collected. Two different runs were performed, one with PM emissions calculated with the original EMIT emission factors and another one using E PM 10 and E PM coarse calculated based on the method described in Section 2.4.2. The 24 h average concentrations, as well as an average for the entire period calculated with both models, are shown in Table 7. For the models, the same statistical measures and nomenclature as those employed for the validation of NO2 are used. There is no significant difference between the results obtained with both runs. The models slightly underpredict PM10 concentration, with MB −1.05 µg/m3 and a positive bias for some days. The ME is 2.27 µg/m3 for the calculation with the original emissions map and 2.25 µg/m3 with the adjusted one. The largest error is for 23 September, following 22 April and 20 and 21 September. The other measures also indicate a similar performance, with an RMSE of 2.92 µg/m3, NMB of −13.30% and NME of 28.55% for the simulation with adjusted emissions, and an RMSE of 2.90 µg/m3, NMB of −13.28% and NME of 28.86% for the simulation with original emissions.
The coefficient calculated for 2019 (see Section 3.1.2) was used to estimate E PM coarse emissions for numerical calculations. The coefficient varies significantly from year to year as many factors (meteorological conditions, traffic flow, and anthropogenic activity, among others) affect air pollution. The emission inventory used in the study was developed for 2018, while the meteorological conditions correspond to 2019 as measurements for the model validation were available for that did not change significantly from 2018 to 2019.

3.2.3. Spatial Distribution of NO2 Concentrations

EMIT provides emission estimations for each segment of primary and secondary roads in the study domain. A comparison of the spatial distribution of NO2 using original and adjusted emissions and the difference in concentration patterns for concentration fields averaged over the entire period is shown in Figure 15. The emissions calculated with the original map provide lower concentrations over the primary roads and higher concentrations over the secondary roads. The adjusted emissions lead to an increase in the NO2 fraction in close proximity to the roads and a reduced concentration within one hundred meters along the roads with intensive traffic, which is in agreement with experimental observations. The difference in simulated concentrations using original and adjusted emissions is up to 20 µg/m3 over the road and 10 µg/m3 far away from it.
Maps with a 24 h average difference between model results obtained using original and adjusted emissions are shown in Supplementary Materials S3. The same pattern as described for the mean field is typical for all days, except for the 23 and the 24 September, when there is no difference over the main road in the area: Tsarigradsko Shose. As can be seen from Table 5, the maximum concentrations for the period were measured at AQS Mladost for these days, at 45.33 µg/m3 and 43.82 µg/m3, respectively. Simulations using adjusted emissions also provide the highest concentrations for the same dates, but there is no correspondence in the 24 h average concentration values obtained using the original emissions. These two days are characterized by very high relative humidity (near saturation) during the afternoon hours and a decrease in incoming solar radiation and temperature. Vertical profiles (see Supplementary Materials S1.) show a strong inversion layer and calm conditions (wind speed less than 1 m/s).

3.2.4. Spatial Distribution of Concentrations of Different PM Fractions

Similar average concentration field maps for PM10, using original and adjusted emissions and the difference in concentration patterns, are shown in Figure 16. Adjusted emissions increase the concentration over the road segments with the most intensive traffic. The difference in simulated concentrations using original and adjusted emissions is up to 5 µg/m3 over the road and is insignificant away from it. Maps of a 24 h average concentration difference between model results obtained using original and adjusted emissions are shown in Supplementary Materials S4. Differences over the road segments with the most intensive traffic are higher for the same days, 23 and 24 September, with a concentration approximately 8 µg/m3 higher for the simulation using adjusted emissions.
The same pattern can be found in the spatial distribution of PM2.5 using original and adjusted emissions. The concentration patterns for the two models, averaged over the study period, and their difference are shown in Figure 17. Adjusted emissions increase the concentration over the road segments with the most intensive traffic. The difference in simulated concentrations using original and adjusted emissions is up to 5 µg/m3 over the road and insignificant away from it. The main reason is that small particles are associated with exhaust traffic emissions and coarse particles with non-exhaust emissions. Maps with a 24 h averaged difference between model results obtained using original and adjusted emissions are shown in Supplementary Materials S5. Differences over the road segments with the most intensive traffic are higher for the same days, 23 and 24 September, with a concentration approximately 7 µg/m3 higher for the simulation using adjusted emissions.
The concentration maps for PMcoarse differ from the other fractions (Figure 18). This fraction is more widely dispersed when using the original emissions, and the values are higher far from the road. In the simulations with adjusted emissions, the coarse fraction is mainly dispersed near the roads due to the faster deposition of heavier particles. Maps with a 24 h average difference between model results obtained using original and adjusted emissions are shown in Supplementary Materials S6.

4. Discussion

Due to the scarcity of reliable experimental measurements supporting the evaluation of emission factors that represent traffic, roads and driving conditions in Bulgaria, advanced statistical models were developed, taking into account the influence of various factors: geographic and urban morphology, artificial surface and canyons, road infrastructure, demographics and public life, etc. The available data were processed on several different levels in order to achieve an optimal representation of traffic conditions in Bulgaria. To the authors’ best knowledge, the approach proposed here is being applied for the first time in that particular context.
Sofia was selected for this pilot study for the following reasons:
  • The city of Sofia is the biggest urban area in the country, and despite the efforts made during the last few decades, the citizens in the capital are still exposed to high levels of PM10 as well as other pollutants, especially NO2. The latter is mostly related to traffic, but the evidence has somehow stayed hidden due to gaps in the monitoring of areas with more intensive traffic [26].
  • Sofia as a study object is a challenging and complex urban system because of its geographical setting and due to the fast expansion of the city and the recent tendencies in compact and car-oriented urban development. Unfavorable in terms of air quality are the rising density and share of impervious surfaces, urban street canyon formation, and the interruption of “green wedges”, which have an important role in the ventilation of the city. In the previous two decades, investment priority was given to the construction of new transport infrastructure for untapping some “bottlenecks,” but this has induced demand, encouraging higher motorization and more intensive travel by private cars while increasing competition for the narrow space with the alternative mass or active mobility options [68,69,70].
  • Despite the relative scarcity of reliable observational data, Sofia is well ahead in terms of experimental infrastructure in comparison to other cities in Bulgaria. Matching the different sources of data and information is a challenging endeavor in the pursuit of better knowledge, which can contribute to more appropriate decision-making.
Our previous work shows many uncertainties related to local-scale modeling, and the transport emission inventory has the highest ambiguity in air quality research [39]. The working hypothesis was that modeling at a local scale can be improved by implementing a more reliable inventory of traffic emissions if the relationships used to develop the emission factors are derived using local data more representative of Bulgarian conditions.
Our results support this hypothesis, showing an improvement in simulated concentrations for the main pollutants after model validation. Estimating NO2 concentration with a customized polynomial, based on local measurements, shows a significant improvement with increasing near-road concentration, most notably during stable conditions (temperature inversion). As already mentioned, there are many uncertainties related to TRAP modeling, especially for PM concentration due to non-exhaust emissions. This was the reason for developing PM10 and PMcoarse emissions following an established methodology but using local observations to better represent real Bulgarian road conditions and traffic structures. Emissions for PM2.5 are obtained by subtracting PM10 and PMcoarse emissions. Although model validation did not show a significant difference at the point of comparison, a significant difference was recorded above and near the major roads between the simulations using the original and adjusted emissions. For the average PM10 concentration field, there was an approximately 20% increase in value over the road with the most intensive traffic, and on days with a strong inversion (e.g., 23–25 September 2019), the difference reached 8 µg/m3. Adjusted emissions show a small contribution of PMcoarse concentrations away from the road, mainly due to the higher deposition velocity of this fraction, but there was a significant increase in PM2.5 concentration of approximately 30% for the average field and 7 µg/m3 on days with stable conditions.
These results raise a very important question about the missing observational data for PM2.5 in Sofia and other urban areas. This pollutant is measured only in one urban background site in Sofia. The situation is the same with other Bulgarian cities. The concentration of PM2.5 is closely related to health risk and mortality mainly from cardiovascular diseases (including ischemic heart diseases and strokes) [71]. Bulgaria faces severe air pollution, with particulate matter levels being amongst the highest in the EU [72]. However, the lack of measurements for the concentration of PM2.5, and therefore the inability to validate numerical model predictions for this pollutant, make health risk assessment difficult and uncertain.
Furthermore, the EU Commission proposed to revise the Ambient Air Quality Directives as part of the European Green Deal. The revision aligns the air quality standards more closely with the recommendations of the World Health Organization (WHO). According to the WHO’s Air Quality Guidelines [73], the global burden of disease associated with air pollution exposure takes a massive toll on human health worldwide: exposure to air pollution is estimated to cause millions of deaths and lost years of life annually. The recommended values are for significant reduction in annual and 24 h air quality levels for PM2.5, PM10 and NO2, pollutants closely related to TRAP emissions. A higher share of population in larger cities, which is a stable trend in Bulgaria as well, is closely related to the increase in the number of cars and human activities. Poor road quality and conditions and the large number of ageing vehicles in use contribute to the systematic exceedance of EU standards and the subsequent financial sanctions for the country.
Another very important issue is the location of the so-called traffic air quality stations. On paper, there are only four of these stations in Bulgaria—two in Sofia, one in the city of Plovdiv and one in the city of Stara Zagora. Despite it being designated as one, AQS Mladost does not conform to the EU requirements for positioning of traffic stations (e.g., traffic-oriented sampling probes have to be no more than 10 m from the kerbside) [74], and consequently, it is not representative of such a device. The lack of reliable measurements near the roads leads to severe biases in the data and the impossibility of properly validating any transport model. More precise measurements with high temporal resolution at true traffic sites are needed to assess the real traffic impacts; volume and structure; TRAP; and associated health and comfort effects for the population.

5. Conclusions

The aim of the present study was to develop a high-quality inventory of traffic emissions using an established methodology but applying local observations, in order to adjust emissions for PM as well as the rates between NO2 and NOx and between PM2.5 and PM10, to better represent traffic conditions in Bulgaria.
The main outcomes from the conducted study are as follows:
  • Descriptive statistics based on observations at three transport stations in Sofia over an extended period (2009–2020) show a trend in the average annual concentrations of NO2 and PM10. A slightly negative trend (decrease in averaged annual concentration) was registered at AQS Pavlovo. The highest values were measured at AQS Orlov Most, but these decreased significantly after the station was moved to a new location in Mladost. The lowest recorded concentrations in 2019 and 2020 are most likely due to the imposed lockdown during the COVID-19 pandemic. Hourly concentrations of NO2 (averaged for the entire period) for weekdays, Saturday, and Sunday are also presented as evidence of the diurnal profiles of traffic-related emissions.
  • A new inventory of traffic emissions has been developed for the city of Sofia based on publicly available traffic and fleet data, as well as a dedicated data collection campaign to fill the gaps in secondary street traffic data. Various processing steps were applied to align diverse geometry, attribute data, and acquisition methods, utilizing the advantages of ensemble learning. This activity data were used with the EMIT model to calculate the traffic emissions. The emissions were then exploited to simulate air pollution in a specific area by the ADMS-Urban model.
  • The polynomial relationship between NOx and NO2 was adjusted based on local dynamics at three AQSs in Sofia for a 12-year period (from 2009 to 2020), to provide a more adequate estimation of NO2 from the ADMS-Urban model.
  • Models of the first and second order with all interactions were fit to data for 12 measurands at AQS Hipodruma (where PM2.5 observations are available) for the period from 2015 to 2020. The two models were validated on one-year worth of data from 2023. Due to its better performance statistics, the second-order model was chosen to predict the concentration of PM2.5 at the chosen traffic station, AQS Pavlovo. Data from both transport and urban background stations were used to calculate the roadside increments and the adjusted emissions of PM10 and PMcoarse.
  • The ADMS-Urban model was validated and evaluated by comparing pollutant concentrations from simulations using original and adjusted emissions, showing an improvement in results after applying functions and relationships derived from local observations.
This work is part of our effort to improve modeling at the local scale and part of a framework methodology to support the management and planning of healthy urban environments and lifestyles. The methodology was developed for the city of Sofia as a pilot study but will be applicable to any other urban area in Bulgaria.
There are a number of points in the present study that will be addressed in the future. The presented work can also benefit from several extensions in the statistical and data processing parts. The first and perhaps most challenging deficiency to overcome is the lack of any sort of uncertainty in the data used in setting up the physics-based models and training the various data maps. This includes measurement and replication uncertainty, model form error, and statistical model uncertainty. Close collaboration with AQS manufacturers, operators, and service personnel as well as careful data collection and formal validation procedures can help quantify and potentially reduce all this uncertainty. Another point of improvement is the validation or elimination of some key assumptions, chiefly the one about mapping PM2.5 from background to traffic stations. In the authors’ opinion, this can only be satisfactorily achieved by equipping the AQSs with relevant sensors and eliminating statistical modeling altogether. Sofia Municipality operates a set of equivalent sensors under the AirThings initiative [28], which could be used as an interim measure to at least validate some of these assumptions.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/atmos15070773/s1, Figure S1: Vertical profiles of wind speed and direction, potential temperature and relative humidity (f) for selected dates of the present study; Figure S2: Definitions of statistics; Figure S3: Maps of NO2 showing the difference in concentration patterns derived from the ADMS-Urban model using original and adjusted emissions; Figure S4: Maps of PM10 showing the difference in concentration patterns derived from the ADMS-Urban model using original and adjusted emissions; Figure S5: Maps of PM2.5 showing the difference in concentration patterns derived from the ADMS-Urban model using original and adjusted emissions; Figure S6: Maps of PMcoarse showing the difference in concentration patterns derived from the ADMS-Urban model using original and adjusted emissions.

Author Contributions

Conceptualization and methodology, R.D. and E.H.; software and validation, M.V., A.B., D.B. and P.O.H.; formal analysis and investigation, R.D., A.B., D.B., P.O.H. and M.V.; resources, A.B., D.B., M.V., E.H. and O.G.; data curation, M.V., P.O.H., A.B., D.B. and O.G.; writing—original draft preparation, R.D., A.B., D.B., P.O.H., M.V., E.H. and O.G.; writing—review and editing, R.D. and P.O.H.; visualization, A.B., D.B., P.O.H. and M.V.; supervision, R.D.; project administration, R.D.; funding acquisition, R.D. and P.O.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Bulgarian National Science Found (BNSF), under project “Development of a methodology for assessing air quality and its impact on human health in an urban environment”, grant number: KP-06-H54/2, 15 November 2021. P.O. Hristov is funded by the Bulgarian National Scientific Fund under National Scientific Program “Petar Beron i NIE”, agreement no. KP-06-DB/3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used to support the findings of this study are available from the corresponding author upon request due to privacy.

Acknowledgments

R. Dimitrova acknowledged the provided access to the cluster PHYSON funded by the National Roadmap on Research infrastructures: 2020–2027; National Centre of HPC and Distributed Systems (NIS-3318, SU).

Conflicts of Interest

The authors declare no conflicts of interest.

Abbreviations

The following abbreviations are used in this manuscript:
TRAPTraffic-Related Air Pollution
PMParticulate Matter
LPGLiquefied Petroleum Gas
NOxNitrogen Oxides
ECElemental Carbon
PM2.5Particulate Matter Not Exceeding 2.5 μm in Aerodynamic Diameter
UFPsUltrafine Particles Not Exceeding 1 μm in Aerodynamic Diameter
PAHPolycyclic Aromatic Hydrocarbons
VOCVolatile Organic Compounds
PM10Particulate Matter Not Exceeding 10 μm in Aerodynamic Diameter ()
THCTotal Hydrocarbons
COCarbon Monoxide
COPERTEuropean Computer Model to Calculate Emissions from Road Traffic
MOBILEMobile Source Emissions Factor Model
MOVESMotor Vehicle Emission Simulator
HBEFAHandbook Emission Factors for Road Transport
EMITComprehensive Emissions Inventory Toolkit
CMEMComprehensive Modal Emission Model
ESTMMultimodal Traffic Simulation Software
EMPASwiss Federal Laboratories for Materials Testing and Research
NAEIUK National Atmospheric Emissions Inventory
NO2Nitrogen Dioxide
AQSAutomatic Air Quality Stations
ADMS-UrbanAir Quality Management & Assessment System
CERCCambridge Environmental Research Consultants Ltd.
O3Ozone
SO2Sulphur Dioxide
CAMSCopernicus Atmosphere Monitoring Service
TNONetherlands Organization for Applied Scientific Research
CMOCentral Meteorological Observatory
EEAExecutive Environment Agency
EMEPEuropean Monitoring and Evaluation Programme
Defra/TRLDepartment for Environment, Food and Rural Affairs/Transport Research Laboratory
RMSERoot-Mean-Square Error
IAIndex of Agreement
SASource Apportionment
PMFPositive Matrix Factorization
EPAEnvironmental Protection Agency
RESResuspension
SECSecondary
BBBiomass Burning
TRTraffic
INDIndustry
FUELFuel-Oil Burning
COCarbon Monoxide
NONitrogen Monoxide
NMBNormalized Mean Bias
NMENormalized Mean Error

References

  1. EUROSAT, Database—Transport—Eurostat—European Commission. Available online: https://ec.europa.eu/eurostat/web/transport/data/database (accessed on 13 May 2024).
  2. EUROSAT, Database—Transport—Eurostat—European Commission. Available online: https://ec.europa.eu/eurostat/web/products-eurostat-news/w/ddn-20240117-1 (accessed on 13 May 2024).
  3. Selleri, T.; Melas, A.D.; Joshi, A.; Manara, D.; Perujo, A.; Suarez-Bertoa, R. An Overview of Lean Exhaust deNOx Aftertreatment Technologies and NOx Emission Regulations in the European Union. Catalysts 2021, 11, 404. [Google Scholar] [CrossRef]
  4. Selleri, T.; Gioria, R.; Melas, A.D.; Giechaskiel, B.; Forloni, F.; Villafuerte, P.M.; Demuynck, J.; Bosteels, D.; Wilkes, T.; Simons, O.; et al. Measuring Emissions from a Demonstrator Heavy-Duty Diesel Vehicle under Real-World Conditions—Moving Forward to Euro VII. Catalysts 2022, 12, 184. [Google Scholar] [CrossRef]
  5. Frey, H.C. Trends in Onroad Transportation Energy and Emissions. J. Air Waste Manag. Assoc. 2018, 68, 514–563. [Google Scholar] [CrossRef]
  6. European Environmental Agency. National Emission Reduction Commitments Directive Emissions Inventory Data. 2024. Available online: https://www.eea.europa.eu/data-and-maps/dashboards/necd-directive-data-viewer-7 (accessed on 17 June 2024).
  7. Mehlhart, G.; Merz, C.; Akkermans, L.; Jordal-Jorgensen, J. European Second-Hand Car Market Analysis, Final Report. 2011. Available online: https://www.oeko.de/oekodoc/1114/2011-005-en.pdf (accessed on 13 May 2024).
  8. National Statistical Institute. Theme 7. SUSTAINABLE TRANSPORT: 7.6. Share of New Registered and Registered Brand New Vehicles. 2016. Available online: https://www.nsi.bg/sites/default/files/files/data/timeseries/SDI%207.6_en.xls (accessed on 17 June 2024).
  9. Burov, A.; Brezov, D. Transport Emissions from Sofia’s Streets—Inventory, Scenarios, and Exposure Setting. In Environmental Protection and Disaster Risks. EnviroRISKs 2022. Lecture Notes in Networks and Systems; Dobrinkova, N., Nikolov, O., Eds.; Springer: Cham, Switzerland, 2023; Volume 638, pp. 223–233. [Google Scholar] [CrossRef]
  10. Za Zemiata. Spatial-Based Scenarios for the Introduction of Low-Emission Zones for Transport in the Metropolitan Municipality Guidelines for Related Measures/Appendixes/. 2024. Available online: https://www.zazemiata.org/wp-content/uploads/2023/02/AirNEZ_web_landscape_FINALANNEXES_2024.pdf (accessed on 17 June 2024).
  11. Mądziel, M. Vehicle Emission Models and Traffic Simulators: A Review. Energies 2023, 16, 3941. [Google Scholar] [CrossRef]
  12. Zhang, R.; Wang, Y.; Pang, Y.; Zhang, B.; Wei, Y.; Wang, M.; Zhu, R. A Deep Learning Micro-Scale Model to Estimate the CO2 Emissions from Light-Duty Diesel Trucks Based on Real-World Driving. Atmosphere 2022, 13, 1466. [Google Scholar] [CrossRef]
  13. De Nunzio, G.; Laraki, M.; Thibault, L. Road Traffic Dynamic Pollutant Emissions Estimation: From Macroscopic Road Information to Microscopic Environmental Impact. Atmosphere 2021, 12, 53. [Google Scholar] [CrossRef]
  14. COPERT Computer Model to Calculate Emissions from Road Traffic. Available online: https://web.jrc.ec.europa.eu/policy-model-inventory/explore/models/model-copert (accessed on 8 March 2024).
  15. Description and History of the MOBILE Highway Vehicle Emission Factor Model. Available online: https://www.epa.gov/moves/description-and-history-mobile-highway-vehicle-emission-factor-model (accessed on 8 March 2024).
  16. The Handbook of Emission Factors for Road Transport (HBEFA). Available online: https://www.hbefa.net (accessed on 8 March 2024).
  17. Cambridge Environmental Research Consultants Ltd. Atmospheric Emissions Inventory Toolkit (EMIT). Version 3.9. User Guide. 2022. Available online: https://www.cerc.co.uk/environmental-software/assets/data/doc_userguides/CERC_EMIT3.9_User%20Guide.pdf (accessed on 8 March 2024).
  18. Mao, F.; Li, Z.; Zhang, K. A Comparison of Carbon Dioxide Emissions between Battery Electric Buses and Conventional Diesel Buses. Sustainability 2021, 13, 5170. [Google Scholar] [CrossRef]
  19. Comprehensive Modal Emission Model (CMEM). Available online: https://www.cert.ucr.edu/cmem (accessed on 8 March 2024).
  20. Environmentally Sensitive Traffic Management (ESTM). Available online: https://www.bosch-mobility.com/en/solutions/air-quality-solutions/environmentally-sensitive-traffic-management (accessed on 8 March 2024).
  21. Swiss Federal Laboratories for Materials Testing and Research (EMPA), Future Mobility. Available online: https://www.empa.ch/web/empa/future-indiviual-mobility (accessed on 8 March 2024).
  22. West, J.L.; Mobley, J.D.; Deslauriers, M.; Feldman, H.; Frey, C.; Rojas-Bracho, L.; Wierman, S.S.; Werner, A.S. NARSTO emission inventory assessment. In Proceedings of the American Association of Aerosol Research 2005, Atlanta, GA, USA, 7–11 February 2005. [Google Scholar]
  23. Health Effects Institute. Panel on the Health Effects of Traffic-Related Air Pollution. In Traffic-Related Air Pollution: A Critical Review of the Literature on Emissions, Exposure, and Health Effects; HEI Special Report 17; Health Effects Institute: Boston, MA, USA, 2010. [Google Scholar]
  24. Thorpe, A.; Harrison, R.M. Sources and properties of non-exhaust particulate matter from road traffic: A review. Sci. Total Environ. 2008, 400, 270–282. [Google Scholar] [CrossRef]
  25. Boulter, P.G.; Thorpe, A.J.; Harrison, R.M.; Allen, A.G. Road Vehicle Non-Exhaust Particulate Matter: Final Report on Emission Modelling, CPEA23/SPU82 Project Report. 2007. Available online: https://uk-air.defra.gov.uk/assets/documents/reports/cat15/0706061624_Report2_Emission_modelling.PDF (accessed on 8 March 2024).
  26. Executive Environment Agency, System for Informing the Population about the Quality of Atmospheric Air. Available online: https://eea.government.bg/kav (accessed on 8 March 2024). (In Bulgarian)
  27. Executive Environment Agency. National Report for the Condition and Preservation of the Environment: Pollutant Emissions and Quality of the Atmospheric Air. 2023. Available online: https://eea.government.bg/bg/soer/2023/1Air.pdf (accessed on 8 March 2024). (In Bulgarian)
  28. AIRTHINGS—A Pilot System of Sensors for Monitoring Atmospheric Air. Available online: https://air.sofia.bg (accessed on 8 March 2024). (In Bulgarian).
  29. Network of Low-Cost Citizen-Owned Sensor Suites. Available online: https://airsofia.info (accessed on 8 March 2024). (In Bulgarian).
  30. European Clean Air Centre. State of Air Quality in Bulgaria: Overview. 2020. Available online: https://cleanaircentre.eu/wp-content/uploads/2019/10/State-of-AQ-Bulgaria.pdf (accessed on 8 March 2024).
  31. Popov, I.; Hlebarov, I.; Raeva, D. Is There Nitrogen Dioxide Pollution of the Air in Sofia? 2022. Available online: https://www.zazemiata.org/wp-content/uploads/2022/05/Za-Zemyata-Doklad-NO2-Online.pdf (accessed on 8 March 2024). (In Bulgarian).
  32. Zareba, M.; Weglinska, E.; Danek, T. Air pollution seasons in urban moderate climate areas through big data analytics. Sci. Rep. 2024, 14, 3058. [Google Scholar] [CrossRef] [PubMed]
  33. Health Effects Institute. Trends in Air Quality and Health in Southeast Europe: A State of Global Air Special Report; Health Effects Institute: Boston, MA, USA, 2022; ISSN 2578-6881; Available online: https://www.stateofglobalair.org/sites/default/files/documents/2022-05/soga-southeast-europe-regional-report_1.pdf (accessed on 17 June 2024).
  34. Hristova, E.; Veleva, B.; Georgieva, E.; Branzov, H. Application of Positive Matrix Factorization Receptor Model for Source Identification of PM10 in the City of Sofia, Bulgaria. Atmosphere 2020, 11, 890. [Google Scholar] [CrossRef]
  35. Air Quality: Commission Decides to Refer Bulgaria to the Court of Justice over Its Failure to Comply with Previous Judgement. Available online: https://ec.europa.eu/commission/presscorner/detail/en/ip_20_2150 (accessed on 8 March 2024).
  36. Stefanov, B.; Dombalov, I.; Pelovski, Y.; Kozarev, N.; Ilieva, N.; Sokolovski, E.; Borisov, D.; Nikolova, Y.; Georgieva, A.; Angelova, S.; et al. Atmospheric Air Quality Management Program of Sofia Municipality for the Period 2015–2020—Reduction of Emissions and Reach of the Established Norms for Fine Particles PM10. 2015. Available online: https://www.sofia.bg/documents/20121/129306/Project_Program_KAV_Sofia_2015-2020.pdf (accessed on 8 March 2024). (In Bulgarian).
  37. Kozarev, N.; Dombalov, I.; Sokolovski, E.; Ilieva, N.; Predyov, I. Program for Supplementation of the Program “Atmospheric air Quality Management Program of Capital Municipality for the Period 2015–2020—Reduction of Emissions and Achievement of the Established Norms for Fine Particles PM10”, by Indicators of Fine Particles with Size up to 2.5 Microns and Polycyclic Aromatic Hydrocarbons. 2019. Available online: https://air.sofia.bg/media/Documents/R.523-Pr1.pdf (accessed on 8 March 2024). (In Bulgarian).
  38. Sofia Municipality. Program for Improving the Air Quality on the Territory of Sofia Municipality for the Period 2021–2026. 2019. Available online: https://www.sofia.bg/documents/20121/992246/Пpoeкт+нa+Пpoгpaмa+KAB+2021-2026.pdf/ (accessed on 8 March 2024). (In Bulgarian).
  39. Dimitrova, R.; Velizarova, M. Assessment of the contribution of different particulate matter sources on pollution in Sofia City. Atmosphere 2021, 12, 423. [Google Scholar] [CrossRef]
  40. Brezov, D.; Burov, A. Ensemble Learning Traffic Model for Sofia: A Case Study. Appl. Sci. 2023, 13, 4678. [Google Scholar] [CrossRef]
  41. Carruthers, D.J.; Edmunds, H.A.; McHugh, C.A.; Singles, R.J. Development of ADMS-Urban and Comparison with Data for Urban Areas in the UK. In Air Pollution Modeling and Its Application XII. NATO—Challenges of Modern Society; Gryning, S.E., Chaumerliac, N., Eds.; Springer: Boston, MA, USA, 1998; p. 22. [Google Scholar] [CrossRef]
  42. Stocker, J.; Hood, C.; Carruthers, D.; McHugh, C. ADMS-Urban: Developments in modelling dispersion from the city scale to the local scale. Int. J. Environ. Pollut. 2012, 50, 308–316. [Google Scholar] [CrossRef]
  43. Georgieva, E.; Syrakov, D.; Prodanova, M.; Etropolska, I.; Slavov, K. Evaluating the performance of WRF-CMAQ air quality modelling system in Bulgaria by means of the DELTA tool. Int. J. Environ. Pollut. 2015, 57, 272. [Google Scholar] [CrossRef]
  44. Average Age of Registered Passenger Cars in Selected European Countries as of 2019. Available online: https://www.statista.com/statistics/974713/passenger-car-average-age-europe/ (accessed on 25 July 2023).
  45. Carslaw, D.C.; Rhys-Tyler, G. New insights from comprehensive on-road measurements of NOx, NO2 and NH3from vehicle emission remote sensing in London, UK. Atmos. Environ. 2013, 81, 339–347. [Google Scholar] [CrossRef]
  46. Regulation 12 of 15 July 2010 for the Limit Levels of Sulphur Dioxide, Nitrous Dioxide, Fine Particulate Matter, Lead, Benzene, Carbon Monoxide and Ozone in the Atmosphere. Available online: https://lex.bg/bg/laws/ldoc/2135691821 (accessed on 12 February 2024). (In Bulgarian).
  47. VAISALA Product Documentation Portal. Available online: https://docs.vaisala.com/v/u/B211729EN-E/en-US (accessed on 8 March 2024).
  48. VAISALA Product Documentation Portal. Available online: https://docs.vaisala.com/v/u/B211221EN-L/en-US (accessed on 8 March 2024).
  49. VAISALA Product Documentation Portal. Available online: https://docs.vaisala.com/v/u/B212709EN-A/en-US (accessed on 8 March 2024).
  50. UK National Atmospheric Emissions Inventory (NAEI). Available online: https://naei.beis.gov.uk/about (accessed on 8 March 2024).
  51. The Co-Operative Programme for Monitoring and Evaluation of the Long-Range Transmission of Air Pollutants in Europe (EMEP). Available online: https://www.emep.int (accessed on 8 March 2024).
  52. Cambridge Environmental Research Consultants Ltd. ADMS-Urban. Urban Air Quality Management System. Version 5.0. User Guide. 2020. Available online: https://www.cerc.co.uk/environmental-software/assets/data/doc_userguides/CERCADMS-Urban5.0_User_Guide.pdf (accessed on 25 July 2023).
  53. Venkatram, A.; Karamchandani, P.; Pai, P.; Goldstein, R. The development and application of a simplified ozone modeling system (SOMS). Atmos. Environ. 1994, 28, 3665–3678. [Google Scholar] [CrossRef]
  54. Copernicus Atmosphere Monitoring Service. Available online: https://atmosphere.copernicus.eu (accessed on 8 March 2024).
  55. Cambridge Environmental Research Consultants Ltd. ADMS Street Canyon Tool. Version 2.0. User Guide. 2018. Available online: https://www.cerc.co.uk/environmental-software/assets/data/doc_userguides/CERC_Street_Canyon_Tool_User_Guide.pdf (accessed on 8 March 2024).
  56. Builtjes, P.J.H.; van Loon, M.; Schaap, M.; Teeuwisse, S.; Visschedijk, A.J.H.; Bloos, J.P. Project on the Modelling and Verification of Ozone Reduction Strategies: Contribution of TNO-MEP, TNO-Report, MEP-R2003/166; Netherlands Environmental Assessment Agency: Apeldoorn, The Netherlands, 2003. [Google Scholar]
  57. Technical Specification Documents. Available online: https://www.cerc.co.uk/environmental-software/technical-specifications.html (accessed on 8 March 2024).
  58. Derwent, R.G.; Middleton, D.R. An empirical function for the ratio [NO2]:[NOx]. Clean Air 1996, 26, 57–60. [Google Scholar]
  59. Dixon, J.; Middleton, D.R.; Derwent, R.G. Sensitivity of nitrogen dioxide concentrations to oxides of nitrogen controls in the United Kingdom. Atmos. Environ. 2001, 35, 3715–3728. [Google Scholar] [CrossRef]
  60. Neto, E.d.A.L.; de Carvalho, F.d.A.T. Constrained linear regression models for symbolic interval-valued variables. Comput. Stat. Data Anal. 2010, 54, 333–347. [Google Scholar] [CrossRef]
  61. Harrison, R.M.; Yin, J.; Mark, D.; Stedman, J.; Appleby, R.S.; Booker, J.; Moorcroft, S. Studies of the coarse particle (2.5–10 µm) component in UK urban atmospheres. Atmos. Environ. 2001, 35, 3667–3679. [Google Scholar] [CrossRef]
  62. Harrison, R.M.; Jones, A.M.; Barrowcliffe, R. Field study of the influence of meteorological factors and traffic volumes upon suspended particle mass at urban roadside sites of differing geometries. Atmos. Environ. 2004, 38, 6361–6369. [Google Scholar] [CrossRef]
  63. Grigoratos, T.; Martini, G. Brake wear particle emissions: A review. Environ. Sci. Pollut. Res. 2015, 22, 2491–2504. [Google Scholar] [CrossRef] [PubMed]
  64. Amato, F.; Pandolfi, M.; Viana, M.; Querol, X.; Alastuey, A.; Moreno, T. Spatial and chemical patterns of PM10 in road dust deposited in urban environment. Atmos. Environ. 2009, 43, 1650–1659. [Google Scholar] [CrossRef]
  65. Denby, B.R.; Kupiainen, K.J.; Gustafsson, M. Chapter 9—Review of road dust emissions. In Non-Exhaust Emissions; Amato, F., Ed.; Academic Press: Cambridge, MA, USA, 2018; pp. 183–203. [Google Scholar]
  66. Norris, G.; Duvall, R.; Brown, S.; Bai, S. EPA Positive Matrix Factorization (PMF) 5.0 Fundamentals and User Guide. 2014. Available online: https://www.epa.gov/sites/default/files/2015-2/documents/pmf_5.0_user_guide.pdf (accessed on 26 September 2021).
  67. Mircea, M.; Calori, G.; Pirovano, G.; Belis, C. European Guide on air Pollution Source Apportionment for Particulate Matter with Source-Oriented Models and Their Combined Use with Receptor Models; Publications Office of the European Union: Luxembourg, 2020. [Google Scholar] [CrossRef]
  68. Madzhirski, V.; Dimitrova, E. Challenges and Lessons Learned during the Initial Implementation of EU Urban Mobility Policy in Bulgarian Cities in the period 2007–2017. Corvinus Reg. Stud. 2019, 4, 83–97. [Google Scholar]
  69. Sustainable Urban Mobility Plan of Sofia. 2019. Available online: https://sofiaplan.bg/wp-content/uploads/2022/05/SUMP-21-05-2019.pdf (accessed on 13 May 2024).
  70. Program for Sofia, I_Territorial Scope and Analysis of the State. 2021. Available online: https://sofiaplan.bg/wp-content/uploads/2021/09/I_Tepитopиaлeн-oбхвaт-и-aнaлиз-нa-състoяниeтo-2.pdf (accessed on 13 May 2024).
  71. Mensah, G.A.; Fuster, V.; Murray, C.J.L.; Roth, G.A. Global Burden of Cardiovascular Diseases and Risks, 1990–2022. J. Am. Coll. Cardiol. 2023, 82, 2350–2473. [Google Scholar] [CrossRef]
  72. Health Effects Institute. Trends in Air Quality and Health in Bulgaria: A State of Global Air Special Report. Boston, MA, 2020. Available online: https://www.stateofglobalair.org/sites/default/files/documents/2022-06/soga-southeast-europe-bulgaria-report.pdf (accessed on 5 January 2024).
  73. World Health Organization. WHO global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulphur Sulfur Dioxide and Carbon Monoxide. Executive Summary. Geneva: World Health Organization. Available online: https://iris.who.int/handle/10665/345334 (accessed on 13 May 2024).
  74. Directive 2008/50/EC of the European Parliament and of the Council of 21 May 2008 on Ambient Air Quality and Cleaner Air for Europe; Annex III, Section C: Microscale Siting of Sampling Points. Available online: https://eur-lex.europa.eu/eli/dir/2008/50/oj (accessed on 13 May 2024).
Figure 1. Conceptual workflow diagram of the input data, models used and experimental adjustment of the TRAP emissions in Sofia, Bulgaria.
Figure 1. Conceptual workflow diagram of the input data, models used and experimental adjustment of the TRAP emissions in Sofia, Bulgaria.
Atmosphere 15 00773 g001
Figure 2. Box plots for individual traffic stations by year (based on hourly concentrations of NO2 for a 12-year period): for AQS Pavlovo (a) and for AQS Orlov Most, moved to a new location in Mladost (b). The box edges represent the first and third quartiles, the line inside the box is the median, and the cross is the mean value. The lower whisker indicates the thresholds for outliers, which are in turn shown as points. Red dots indicate exceedance of hourly threshold.
Figure 2. Box plots for individual traffic stations by year (based on hourly concentrations of NO2 for a 12-year period): for AQS Pavlovo (a) and for AQS Orlov Most, moved to a new location in Mladost (b). The box edges represent the first and third quartiles, the line inside the box is the median, and the cross is the mean value. The lower whisker indicates the thresholds for outliers, which are in turn shown as points. Red dots indicate exceedance of hourly threshold.
Atmosphere 15 00773 g002
Figure 3. Box plots for individual traffic stations by year (based on hourly concentrations of PM10 for a 12-year period): for AQS Pavlovo (a) and for AQS Orlov Most, moved to a new location in Mladost (b). The box edges represent the first and third quartiles, the line inside the box is the median, and the cross is the mean value. The lower whisker indicates the thresholds for outliers, which are in turn shown as points. Red dots indicate exceedance of hourly threshold.
Figure 3. Box plots for individual traffic stations by year (based on hourly concentrations of PM10 for a 12-year period): for AQS Pavlovo (a) and for AQS Orlov Most, moved to a new location in Mladost (b). The box edges represent the first and third quartiles, the line inside the box is the median, and the cross is the mean value. The lower whisker indicates the thresholds for outliers, which are in turn shown as points. Red dots indicate exceedance of hourly threshold.
Atmosphere 15 00773 g003
Figure 4. Diurnal variation of mean concentrations (averaged over the entire period) for weekdays (a), Saturday (b), and Sunday (c) for NO2 and PM10.
Figure 4. Diurnal variation of mean concentrations (averaged over the entire period) for weekdays (a), Saturday (b), and Sunday (c) for NO2 and PM10.
Atmosphere 15 00773 g004
Figure 5. Temperature (a) and relative humidity (b) hourly observations for the selected periods.
Figure 5. Temperature (a) and relative humidity (b) hourly observations for the selected periods.
Atmosphere 15 00773 g005
Figure 6. Wind rose for the study period. The concentric circles show the contour of constant number of measurements.
Figure 6. Wind rose for the study period. The concentric circles show the contour of constant number of measurements.
Atmosphere 15 00773 g006
Figure 7. Vertical profiles of wind speed and direction, potential temperature and relative humidity for selected days and hours. The days with the highest biases and errors comparing modeled and observed PM10 concentrations are selected on 22 April at 12 UTC (a), and 23 September at 06 UTC (b) and 12 UTC (c); and days with the lowest biases and errors comparing modeled and observed PM10 concentration on 23 April at 12 UTC (d), 25 September at 06 UTC (e) and 12 UTC (f).
Figure 7. Vertical profiles of wind speed and direction, potential temperature and relative humidity for selected days and hours. The days with the highest biases and errors comparing modeled and observed PM10 concentrations are selected on 22 April at 12 UTC (a), and 23 September at 06 UTC (b) and 12 UTC (c); and days with the lowest biases and errors comparing modeled and observed PM10 concentration on 23 April at 12 UTC (d), 25 September at 06 UTC (e) and 12 UTC (f).
Atmosphere 15 00773 g007
Figure 8. Active traffic data using the ensemble learning algorithm of [28]. Clockwise from top left: traffic density, vehicle speeds, presence of canyons and longitudinal slopes of sections.
Figure 8. Active traffic data using the ensemble learning algorithm of [28]. Clockwise from top left: traffic density, vehicle speeds, presence of canyons and longitudinal slopes of sections.
Atmosphere 15 00773 g008
Figure 9. Modeling domain and locations of air quality stations (a); modeling domain for emission factors adjustment with both sites used for model validation (distance to the road is approximately 70 m for AQS Mladost and 170 m for CMO) (b).
Figure 9. Modeling domain and locations of air quality stations (a); modeling domain for emission factors adjustment with both sites used for model validation (distance to the road is approximately 70 m for AQS Mladost and 170 m for CMO) (b).
Atmosphere 15 00773 g009
Figure 10. Results from the pairwise F-tests for polynomials for all data sets.
Figure 10. Results from the pairwise F-tests for polynomials for all data sets.
Atmosphere 15 00773 g010
Figure 11. Scatter plot of hourly data at different stations in Sofia, showing yield versus total oxides of nitrogen.
Figure 11. Scatter plot of hourly data at different stations in Sofia, showing yield versus total oxides of nitrogen.
Atmosphere 15 00773 g011
Figure 12. Computing a representative incremental ratio for PM10 emissions estimation. Cleaned data for 2019 (blue dots).
Figure 12. Computing a representative incremental ratio for PM10 emissions estimation. Cleaned data for 2019 (blue dots).
Atmosphere 15 00773 g012
Figure 13. Scatter plot of predictions versus measurements of PM2.5 at AQS Hipodruma (background AQS) for 2020, for (a) first order and (b) second-order linear regression. In the case of a perfectly predicting model, the blue dots should lie on the 45° red, dashed line.
Figure 13. Scatter plot of predictions versus measurements of PM2.5 at AQS Hipodruma (background AQS) for 2020, for (a) first order and (b) second-order linear regression. In the case of a perfectly predicting model, the blue dots should lie on the 45° red, dashed line.
Atmosphere 15 00773 g013
Figure 14. Estimated pollutant emissions. Clockwise from top left: PM10, PM2.5, NO2, NOx with EMIT using the active data shown in Figure 7.
Figure 14. Estimated pollutant emissions. Clockwise from top left: PM10, PM2.5, NO2, NOx with EMIT using the active data shown in Figure 7.
Atmosphere 15 00773 g014
Figure 15. NO2 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Figure 15. NO2 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Atmosphere 15 00773 g015aAtmosphere 15 00773 g015b
Figure 16. PM10 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Figure 16. PM10 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Atmosphere 15 00773 g016
Figure 17. PM2.5 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Figure 17. PM2.5 concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Atmosphere 15 00773 g017
Figure 18. PMcoarse concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Figure 18. PMcoarse concentration map averaged over the entire study period using (a) original and (b) adjusted emissions; (c) difference in concentration patterns for the two types of emissions.
Atmosphere 15 00773 g018
Table 1. Summary of model selection study. Best performing models are bolded for each data set.
Table 1. Summary of model selection study. Best performing models are bolded for each data set.
All StationsPavlovoMladostOrlov Most
Active TermsR2 adj.AICR2 adj.AICR2 adj.AICR2 adj.AIC
[0, 1]0.92−249.40.86−130.60.84−102.80.94−268.2
[0, 1, 2]0.92−247.90.89−136.80.96−142.70.95−275.0
[0, 2]0.91−242.60.89−138.00.90−119.10.90−238.9
[0, 1, 2, 3]0.99−369.10.98−201.00.98−166.00.99−372.4
[0, 2, 3]0.93−257.20.88−136.00.94−135.00.96−296.1
[0, 1, 3]0.92−247.40.88−134.50.95−137.20.95−279.0
[0, 3]0.88−220.10.88−134.90.93−131.10.85−209.9
[0, 1, 2, 3, 4]0.99−371.20.98−200.30.98−164.00.99−370.4
[0, 2, 3, 4]0.98−348.40.98−201.00.98−164.50.99−367.4
[0, 1, 3, 4]0.99−355.90.98−202.20.98−165.40.99−369.5
[0, 1, 2, 4]0.99−363.70.98−202.10.98−165.90.99−371.8
[0, 3, 4]0.95−276.70.89−140.20.94−132.00.98−329.3
[0, 2, 4]0.94−262.20.88−136.20.94−132.00.97−306.6
[0, 1, 4]0.92−247.60.87−132.60.94−132.40.95−283.2
[0, 4]0.83−198.60.85−126.50.94−134.00.78−188.2
Table 2. Comparative performance between two established models and the polynomial developed in this study. RMSE in µg/m3.
Table 2. Comparative performance between two established models and the polynomial developed in this study. RMSE in µg/m3.
Air Quality Monitoring SiteAccording to [59]
(Yield)
This Study
(Yield)
According to [58]
(NO2)
This Study
(NO2)
All stations10.90.0138457.26.0
Pavlovo5.40.0183476.44.5
Orlov Most10.90.0138457.56.0
Table 3. Estimated coefficients for the linear regression model applied for AQS Hipodruma.
Table 3. Estimated coefficients for the linear regression model applied for AQS Hipodruma.
ModelNumber of TermsAdjusted R2RMSE TrainingRMSE Test
First-order130.86210.5 µg/m38.29 µg/m3
Second-order910.9237.84 µg/m36.20 µg/m3
Table 4. Proportional factors β between P M 10 and P M c o a r s e for the period 2015–2020 at a traffic and a background station.
Table 4. Proportional factors β between P M 10 and P M c o a r s e for the period 2015–2020 at a traffic and a background station.
Year β t r a f f i c β b a c k g r o u n d
20150.340.38
20160.350.46
20170.290.33
20180.320.33
20190.360.45
20200.380.47
Table 5. ADMS-Urban model validation for NO2 concentration comparing results obtained with original and adjusted polynomial function. All units are µg/m3.
Table 5. ADMS-Urban model validation for NO2 concentration comparing results obtained with original and adjusted polynomial function. All units are µg/m3.
DateAQS
Mladost (Measured)
ADMS-
Urban
(Original)
ADMS-
Urban
(Adjusted)
Bias
(Original)
Bias
(Adjusted)
Error
(Original)
Error
(Adjusted)
18 April201914.8920.8117.065.912.175.912.17
19 April 201913.8132.2410.4918.43−3.3218.433.32
20 April 201912.4240.085.2627.66−7.1527.667.15
21 April 20198.0126.664.4018.65−3.6118.653.61
22 April 20197.9617.043.799.09−4.179.094.17
23 April 20198.3415.0914.046.755.706.755.70
20 September 201934.949.1411.17−25.79−23.7725.7923.77
21 September 201930.158.463.76−21.69−26.3921.6926.39
22 September 201929.829.1210.60−20.69−19.2220.6919.22
23 September 201945.3319.7024.13−25.63−21.2025.6321.20
24 September 201943.8215.7231.66−28.10−12.1628.1012.16
25 September 201940.898.5118.82−32.38−22.0732.3822.07
26 September 201939.0916.9211.42−22.16−27.6622.1627.66
Average25.3418.4212.81−6.92−12.53 20.2313.74
Table 6. Comparison of obtained with PMF model (PMF) and measured (meas) concentrations, and contribution of different categories * to the total PM10 (PMF) concentration. All units are µg/m3.
Table 6. Comparison of obtained with PMF model (PMF) and measured (meas) concentrations, and contribution of different categories * to the total PM10 (PMF) concentration. All units are µg/m3.
DatePM10 (PMF)PM10
(meas)
RESTRBBINDFUELNSO4SEC
18 April 2019 20.3815.235.143.286.490.050.250.004.930.37
19 April 201921.4117.584.913.874.670.110.240.002.804.84
20 April 201918.4418.973.542.316.960.000.240.000.714.99
21 April 201919.5920.666.262.324.220.680.190.000.007.08
22 April 201927.1324.166.833.575.161.440.310.003.076.79
23 April 201928.6327.553.343.319.402.180.220.176.873.13
20 September 201919.1721.826.942.081.410.080.220.302.795.34
21 September 201919.2024.065.932.684.090.150.160.076.080.04
22 September 201925.7328.004.172.325.342.580.210.029.951.15
23 September 201929.9132.427.844.844.041.090.030.289.572.20
24 September 201927.7333.882.803.847.742.270.020.006.185.12
25 September 201929.7330.792.903.282.802.000.050.0012.086.74
26 September 201919.6722.542.821.260.000.430.140.1210.156.13
* Resuspension (RES), Secondary (SEC), Biomass burning (BB), Traffic (TR), Industry (IND), Nitrate-rich (N), Fuel oil burning (FUEL), and Mixed SO42− (SO4).
Table 7. ADMS-Urban model validation for PM10 related to transport comparing results obtained with original adjusted emission factors and PMF model at CMO site. All units are µg/m3.
Table 7. ADMS-Urban model validation for PM10 related to transport comparing results obtained with original adjusted emission factors and PMF model at CMO site. All units are µg/m3.
DatePMF Model
at CMO
ADMS-
Urban
(Original)
ADMS-
Urban
(Adjusted)
Bias
(Original)
Bias
(Adjusted)
Error
(Original)
Error
(Adjusted)
18 April2019 8.428.869.040.430.620.430.62
19 April 20198.788.138.24−0.65−0.550.650.55
20 April 20195.858.018.062.162.222.162.22
21 April 20198.575.835.87−2.75−2.702.752.70
22 April 201910.405.855.89−4.55−4.514.554.51
23 April 20196.656.566.63−0.09−0.020.090.02
20 September 20199.024.975.20−4.05−3.824.053.82
21 September 20198.615.375.28−3.24−3.333.243.33
22 September 20196.497.927.871.431.381.431.38
23 September 201912.686.446.19−6.24−6.496.246.49
24 September 20196.647.517.350.870.710.870.71
25 September 20196.186.596.420.420.240.420.24
26 September 20194.076.746.722.672.652.672.65
Average7.876.836.83−1.05−1.052.272.25
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Velizarova, M.; Dimitrova, R.; Hristov, P.O.; Burov, A.; Brezov, D.; Hristova, E.; Gueorguiev, O. Evaluation of Emission Factors for Particulate Matter and NO2 from Road Transport in Sofia, Bulgaria. Atmosphere 2024, 15, 773. https://doi.org/10.3390/atmos15070773

AMA Style

Velizarova M, Dimitrova R, Hristov PO, Burov A, Brezov D, Hristova E, Gueorguiev O. Evaluation of Emission Factors for Particulate Matter and NO2 from Road Transport in Sofia, Bulgaria. Atmosphere. 2024; 15(7):773. https://doi.org/10.3390/atmos15070773

Chicago/Turabian Style

Velizarova, Margret, Reneta Dimitrova, Petar O. Hristov, Angel Burov, Danail Brezov, Elena Hristova, and Orlin Gueorguiev. 2024. "Evaluation of Emission Factors for Particulate Matter and NO2 from Road Transport in Sofia, Bulgaria" Atmosphere 15, no. 7: 773. https://doi.org/10.3390/atmos15070773

APA Style

Velizarova, M., Dimitrova, R., Hristov, P. O., Burov, A., Brezov, D., Hristova, E., & Gueorguiev, O. (2024). Evaluation of Emission Factors for Particulate Matter and NO2 from Road Transport in Sofia, Bulgaria. Atmosphere, 15(7), 773. https://doi.org/10.3390/atmos15070773

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop