Skoči na glavni sadržaj

Izvorni znanstveni članak

Analyses of urban and rural particulate matter mass concentrations in Croatia in the period 2006–2014

Amela Jeričević ; Croatian Civil Aviation Agency, Zagreb, Croatia
Vedrana Džaja Grgičin ; Meteorological and Hydrological Service of Croatia, Zagreb, Croatia
Maja Telišman Prtenjak ; University of Zagreb, Faculty of Science, Department of Geophysics, Zagreb, Croatia
Sonja Vidič ; Meteorological and Hydrological Service of Croatia, Zagreb, Croatia
Henk Bloemen ; National Institute for Public Health and the Environment, Utrecht, The Netherlands

Puni tekst: engleski pdf 3.294 Kb

str. 157-181

preuzimanja: 348


Preuzmi JATS datoteku


Air pollution due to particulate matter with an aerodynamic diameter ≤ 10 μm (PM10) and diameter ≤ 2.5 μm (PM2.5) was investigated using available measurements in Croatia. Hourly and daily PM10 and PM2.5 mass concentrations were analysed for urban and rural background stations. Urban and industrial
PM10 measurements from 2006 to 2014 and rural background PM10 and PM2.5 measurements from 2011 to 2014 were used to determine the main characteristics of regional and local air pollution related to particles and to assess the level of air pollution according to the defined limit values (CAFÉ Directive 2008/50/EC). Significant difference between coastal and continental urban and rural PM10 levels was found with high PM10 concentrations in continental areas and considerably lower PM10 values at the coast. Substantially different spatial distribution of rural background PM concentrations with relatively high annual average concentration ~ 30 μm/m3 at e.g. continental station Kopački rit, to only ~ 5 μm/m3 at Hum located on the island of Vis was found. Strong spatial PM gradients from the continent towards the coast are resulting from the distinct differences in climatological characteristics between Croatia’s inlands and the coastline and different physical processes affecting the PM concentrations in the two regions (e.g. sea breeze, deposition, wet removal, resuspension). Furthermore, distinct complex chemical mechanisms that depend on emission rate and PM composition and on meteorological parameters (sun radiation, humidity etc.) have different influence on the formation of secondary organic aerosol in marine and continental boundary layer. Ratio between rural PM2.5 and PM10 mass concentrations (PM2.5/PM10) for all sites was investigated and compared to the results obtained in other European studies. Rural site-specific PM2.5/PM10 ratios ranged from 0.6 to 0.9 during the warmer part of the year while during the colder part of the year ratios were higher and ranged from 0.85 to 0.98. The spatial gradient of rural PM2.5/PM10 values is orientated from the coast to the continent implying that at the coast the PM concentrations are mainly composed of fine particles. Observed hourly temperature, relative humidity and wind values in the period from 2006 to 2014 are analysed together with the corresponding PM concentrations and an initial source identification was conducted based on bivariate polar plots.

Ključne riječi

urban and rural air quality, spatial distribution of PM mass concentration, PM2.5/PM10 ratio, source identification, Croatia

Hrčak ID:



Podaci na drugim jezicima: hrvatski

Posjeta: 1.337 *

1. Introduction

Elevated atmospheric particulate matter (PM) concentrations are associated with significant adverse health effects (Samet et al., 2000; Peters et al., 2001; Pope et al., 2002; Samoli et al., 2005; Anderson, 2009), affect ecosystems, influ- ence visibility and cloud formation in the atmosphere and play an important role in climate change (e.g. Andrea et al., 2005; Jiang et al., 2013). In most European countries including Croatia air quality standards for particulate matter have been introduced in order to protect human health and environment. Current standards are set for PM10 (total mass concentration of particles smaller than 10 μm) and PM2.5 (total mass concentration of particles smaller than 2.5 μm) (European Air Quality Directive 2008/50/EC, further used as Directive).

PM comes from a variety of natural and anthropogenic sources; they can be directly emitted to the atmosphere or formed as secondary pollutants in atmospheric chemical reactions. As a consequence, PM comprises complex mixture of solid and liquid parcels of organic matter, Earth’s crust elements iron (Fe), calcium (Ca), aluminium (Al), silicon (Si), potassium (K), and chlorine (Cl), secondary inorganic aerosols and trace metals. In addition, wind-blown soil and re-suspend- ed dust contribute largely to the coarse particle fraction (Harrison et al., 1999; Putaud et al., 2004; Luhana et al., 2004; Forsberg et al., 2005; Amato et al., 2009; Jeričević et al., 2012). It has been shown that contribution of PM emissions and their origins can be relevant at spatial scales ranging from local to regional and long-range, transboundary transport scales and can be analysed accordingly (e.g. Querol et al., 2004; Juda-Rezler et al., 2011).

Comprehensive analyses on PM data obtained over the past decade at more than 60 regional backgrounds, rural, suburban, urban and kerbside stations across Europe were conducted. Putaud et al. (2010) found that there is no single ratio between PM2.5 and PM10 mass concentrations for all sites although fairly constant ratios ranging from 0.5 to 0.9 were observed at most individual sites. The main constituents of both PM10 and PM2.5 were generally organic matter, sul- phate and nitrate while mineral dust could be a major constituent of PM10 at kerbside sites in Southern Europe. Source apportionment methods and results of particulate matter in Western Europe were reviewed by Viana et al. (2008). Previous European studies mainly did not include data from Eastern and Southeastern (e.g. Croatia, Serbia, Romania, Bulgaria) European countries and generally there is a gap in knowledge on the PM levels and compositions in those areas.

Up to now research on PM in Croatia was mainly focused on urban concen- trations particularly on the mass and chemical analyses of heavy metals and polycyclic aromatic hydrocarbons (PAH) at residential urban background sta- tions in Zagreb area (e.g. Bešlić et al., 2007; Bešlić et al., 2008; Šišović et al., 2008) as well as research on PM and ozone in Rijeka area (Alebić-Juretić et al., 2007). Regional observations obtained with passive samplers during short-term campaigns in Croatia, Serbia and Bosnia and Herzegovina among other coun- tries in Central and Southeastern Europe have been used to analyse PAH levels and conduct the source apportionment analyses based on mass balance tech- nique (Klánová et al., 2007; Dvorská et al., 2011).

This research is focused on the detailed spatial and temporal analyses of ur- ban and rural PM10 and PM2.5 levels as well as on analyses of the corresponding meteorological parameters in Croatia in the period from 2006 to 2014 in order to investigate local air quality characteristics, compare results with studies con- ducted in other areas in Europe and provide basis for future investigations of PM pollution in Croatia. Data from the rural background measurement in the period from 2011–2014 are investigated for the first time in Croatia that is particularly important for the determination of relevant scientific knowledge of background air quality in Croatia regarding the PM and the contribution to the urban PM levels.

2. Data and methods

2.1. Air quality observations and PM emission sources

The PM measurements were acquired with three different observing ap- proaches (Tab. 1): gravimetric method (BS EN12341:2014, Ambient air. Standard gravimetric measurement method for the determination of the PM10 or PM2.5 mass concentration of suspended particulate matter), β-attenuation monitoring (e.g. Willeke and Baron, 1993) and by the optical particle counters of Grimm 180 instrument. The comparison of PM concentration data obtained by different measurement methods is still considered as a difficult problem. Standard, gravimetric method (BS EN1234:2014) is a classic method of weigh- ing of the deposited mass on a filter that is accepted as a standard reference method against which all other measurement methods are validated (Noble et al., 2001; EC, 2010). Although gravimetric method is a standard used for com- pliance reasons in EU there are numerous papers showing reactions going on between air and deposited particles as well as within the aerosol mass compro- mising these measurements too. Temperature and humidity history of the dust loaded onto filter greatly influences the actual mass on the filter (Allen et al., 1997; Eisner and Wiener, 2002; Pang et al., 2002). It is stated that at 50% relative humidity, aerosol particles can contain up to 30% of water (e.g. Putaud et al., 2004). On the other hand, calibration, temperature and humidity issues can create artefacts for β-attenuation monitoring that has to be taken into ac- count (e.g. Allen et al., 1997; Hauck et al., 2004). Demonstration that other methods than gravimetric yield results equivalent to the reference (EC, 2010) is needed. In Croatia the equivalence study was conducted for one station (Zagreb-1, see Tab. 1) between the gravimetric and β-attenuation monitors for three-year period from 2007 to 2010 and a satisfactory agreement based on cor- relation coefficient was found (Davila et al., 2011). However, the highest un- derestimation by the automatic (β-attenuation) measurements was found dur- ing the winter time when PM10 concentrations generally achieve maximum values. Since rural background measurements obtained with optical counters were not compared against gravimetric measurements the possible underestimation of observed concentrations was taken into account during interpretation of data.

The list of stations and methods for the particulate matter analyses.
Station nameLatitude, LongitudePollutantType of station*Methods of PM analyses
Zagreb-145.8, 16.0PM10, NO2UT(1), (2)
Rijeka-1 45.3, 14.4PM10UT(1), (2)
Rijeka-245.3, 14.5PM10UB(2)
Sisak-145.5, 16.4PM10UI(1),(2)
Kutina-145.5, 16.8PM10UBI(2)
Osijek-145.6, 18.7PM10UB(2)
Kopački rit45.7, 18.8PM10, PM2.5RB(3)
Parg45.6, 14.6PM10, PM2.5RB(3)
Višnjan45.3, 13.7PM10, PM2.5RB(3)
Plitvice44.9, 15.6PM10, PM2.5RB(3)
Polača44.0, 15.5PM10, PM2.5RB(3)
Žarkovica42.6, 18.0PM10, PM2.5RB(3)
Hum43.0, 16.0PM10, PM2.5RB(3)

*Classification of the stations is UT-urban traffic, UB-urban background, UI-urban industrial, RB-rural background, UBI-urban background industrial. Depending on the monitoring sites, concentrations of PM10 were determined by four methods: (1) a reference gravimetric sampler with quartz filter (Standard EN 12341), (2) β -attenuation method (MP101M), and (3) Grimm 180.

Figure 1. Stations in Croatia with the measurements of PM concentrations: a) name, location and type and b) gridded PM10 emissions EMEP4HR emissions on 10 km × 10 km horizontal resolution. The urban traffic stations are marked with blue circle and urban industrial with blue empty circle, urban background stations are red squares, urban background under the influence of industry are red empty squares and rural background stations are green triangles.

Air quality state network (further called state network) operated and maintained by the Meteorological and Hydrological Service of Croatia (MHSC). Locations of all stations are shown in Fig. 1a and details about the stations, classification and measurement programme is given in Tab. 1. In Fig. 1b spatially distributed PM10 gridded emissions (EMEP4HR, Jeričević et al., 2007) were shown to inspect the location of sites in relation to the emission sources. The largest PM10 source is located in the south of Macedonia with emissions ~ 17 × 105 kg/year, following by the sources in Serbia, near Belgrade where emissions ~ 14 × 105 kg/year are recorded, then in Bosnia and Herzegovina with emissions ~ 9 × 105 kg/year and in Budapest, Hungary with emissions ~ 4 × 105 kg/year. The largest PM source with emissions ~ 300 tons/year in Croatia is located in the centre continental part of Croatia and comprises of the cement industry in the city of Našice.

In Croatia the air quality measurements from 13 (urban and rural) stations from the state network were used. Stations: Rijeka-2 and Osijek-1 are urban background stations, Sisak-1 is urban industrial, Kutina-1 is urban background industrial while Zagreb-1 and Rijeka-1 are urban traffic stations. Kutina-1 is under the influence of petrochemical industry while Sisak-1 is mainly impacted by the emissions of oil refineries. City of Rijeka is situated at the Adriatic coast while all other urban stations are continental. In Croatia 12 rural stations as part of the EU-funded PHARE 2006 project “Establishment of Air Quality Monitoring and Management System” have started with measurements in 2011. For this study observed hourly PM10 and PM2.5 concentrations were available from 7 rural background stations in Croatia distributed over the Croatian territory to capture the spatial and temporal concentration variations of different pollutants (Fig. 1a).

2.2. Limit values and data capture

According to the Directive, daily average PM10 concentrations (c_ (PM10 )day) should not exceed 50 μg/m3 more than 35 times in a calendar year while annual average concentrations (c_ (PM10 )year) should be lower than 40 μg/m . For PM2.5 concentrations an average annual limit value amounts 25 μg/m3. In order to determine the baseline of air pollution with PM in Croatia available PM10 and PM2.5 measurements were analysed in the period from 2006 to 2011 in the view of regulated limit values (Directive). Values of c (PM10 )year were calculated only for the urban stations with data capture > 75%. In Tab. 2 the data capture of PM10 measurements at urban stations in Croatia in the period from 2006 to 2014 is shown. Also data capture for rural background stations in the period from 2011 to 2014 is given. Adequate data capture was available at almost all stations. The exception was identified for Zagreb-2 during 2009, 2010 and 2011, Rijeka-2 during 2008, 2011 and 2012, Polača 2011, 2014 and Žarkovica 2014. The calculation of the annual average stations with the number of exceedance higher than 35 were taken into the consideration although data capture was lower than 75%.

Data capture of PM10 measurements* at urban stations in the period from 2006 to 2014 and rural background stations in Croatia in the period from 2011 to 2014.
Station \ Year200620072008200920102011201220132014
Zagreb-1_GRAV 90%90%90%82%90%90%90%90% 90%
Zagreb-399%99%94%98%88% 68%99%99%100%
Kutina-197%97%95%92%93%84%93%81% 90%
Rijeka-1_GRAV46% 83%91%100%
Osijek-189%99%90%100%91%87%99%87% 88%
Kopački rit96%83%78%86%
Plitvice76%81%69% 94%
Višnjan97% 97%83% 79%
Hum96%83% 78%86%

* Measurements at Zagreb-a_GRAV and Rijeka-1_GRAV are conducted with reference gravimetric methods (Standard EN 12341) while all other measurements are performed with automatic, β-attenuation method (MP101M).

2.3. Measurement and analysis of meteorological parameters

Available observed hourly meteorological parameters: temperature, relative humidity and wind speed and direction in the period from 2006 to 2014 for urban and in the period from 2011 to 2014 for rural stations are analysed using com- puter programming language R. The R is a statistical system, a system for statis- tical computation and graphics, and an environment for data analysis and statis- tics. The Openair software with different functions based on R was applied (Carslaw et al., 2012; Carslaw, 2012). The normalised hourly and monthly temperature, relative humidity and PM values are plotted as well as wind roses. Furthermore, bivariate polar plots (e.g. Carslaw et al., 2006; Westmoreland et al., 2007) of wind speed and direction together with the corresponding concentra- tions are analysed. Plots are shown as a continuous surface calculated through modelling using smoothing techniques. These plots have proved to be useful for quickly gaining a graphical impression of potential sources influences at a location. Bivariate polar plots are constructed by portioning wind speed, wind direction and concentration data into wind speed-direction bins and the mean concentration calculated for each bin. Testing on a wide range of data suggests that wind direction intervals at 10 degrees and 30 m/s wind speed intervals capture sufficient detail of the concentration distribution. The wind direction data typically available are generally rounded to 10 degrees and for typical surface measurements of wind speed in the range 0 to 20 to 30 m/s, intervals greater than 30 m/s would be difficult to justify based on a consideration of the accuracy of the instruments. The bins with fewer than 3 measurements were removed,

3. Results

3.1. The analyses of PM10 and PM2.5 concentrations in Croatia in the period from 2006 to 2014

The measurements of urban PM pollutant concentrations in the state network started in 2006 when the annual limit value was exceeded at almost all continental urban sites. In the period from 2007 to 2014 a continuous decreasing trend is found and calculated c_ (PM10 )year were considerably below the annual limit value 40 μg/m3 at all available stations in Croatia (Fig. 2 upper panel). Data are valuable for the analysis although it should be mentioned that non-reference methods were used and the equivalence studies are applied for the reporting pur- 3 poses. Generally c_ (PM10 )year varied between 30 and 35 μg/m at continental urban stations while at coastal urban stations i.e., in the city of Rijeka c (PM10 )year varied between 15 and 28 μg/m3. The number of days (N) with c_ (PM10)day higher than the limit value (50 μg/m3) is shown in Fig. 2 (lower panel). It was found that N was well below the allowed number of days per year (NLV = 35) for the coastal stations in urban Rijeka area while at all continental stations N > NLV for almost all years. From 2012 a decrease in N is found and N < NLV for the majority of the continental stations except Kutina-1 in 2013. During 2006, N was ~ 70 at all continental stations and somewhat higher ~ 90 at Zagreb-3 and Kutina-1, while N ~ 30 was found at most continental stations in 2012.

PM10 and PM2.5 measurements from 7 rural background stations in the period from 2011–2014 were used to assess average regional PM levels in Croatia and conduct the comparison with the average urban values. Calculated c_ (PM10 )year for rural background stations are shown in Fig. 3. For the background stations in the central part of Croatia comprising stations: Plitvice, Parg, Višnjan and Polača, values c_ (PM ) ~ 15 μg/m3 were found for all analysed years

The annual limit value for PM2.5 was exceeded only at continental rural background station Kopački rit (not shown) that is located at the northeastern part of Croatia on the border with Hungary and Serbia indicating higher regional concentrations which is consistent with the gridded emissions (Fig. 4). The lowest c_ (PM10 )year and c_ (PM 2.5 )year among all analysed PM concentrations were found at coastal rural background stations in the southern part of Croatia with yearly average values between 5 μg/m3 at Hum and 10 μg/m3 at Žarkovica station

Figure 2. Annual average concentrations (upper panel) and number of exceedance of the daily limit value 50 μg/m3 (lower panel) at different urban stations in Croatia in the period from 2006 to 2014. Thick lines denote defined limit values according to the Directive, i.e., average annual value of 40 μg/m3 (upper panel) and allowed number of exceedance of the daily concentrations with the defined limit value 35 (lower panel).

Figure 3. Yearly PM10 averages at rural background stations during period from 2011 to 2014.

Figure 4. The average monthly PM10 concentrations at selected stations in Croatia during the period from 2006 to 2014 at urban stations and from 2011 to 2014 at rural stations. Smooth trend function and 95% confidence intervals are also shown.

The N > NLV at rural stations was found only at Kopački rit in 2011 with N = 55. Over the Croatia there is a decreasing spatial gradient of rural background PM10 and PM2.5 levels from the continent toward the coast which is a consequence of different climatic and dynamical atmospheric processes that influence PM deposition on the coast. The significant decreasing trend in c_(PM10 )year found at continental urban stations that is most likely a consequence of the reduction in contribution of local emission sources was not evidenced at rural background stations where contribution in regional PM levels was almost at the constant level during the analysed period. The contribution of the background PM concentrations to the urban levels was determined through the ratio of rural and urban average annual concentrations. The results are shown in Tab. 3 for the period from 2011 to 2014. Continental rural background station Kopački rit situated near urban station Osijek-1 has the highest average concentrations and indicates that the contribution to urban PM values at the eastern part of Croatia is extremely high with contributions from 80% to over 100% for all years except for 2013. Other continental rural stations situated in the northern and central part of Croatia contributed to the urban PM concentrations ~ 50%. At the coast rural background stations contributed significantly to the urban values with the share of ~ 80%.

Ratio of rural and urban average annual concentrations indicating the contribution of rural background concentrations in the rural PM levels in the period from 2011 to 2014.
Kopački rit/Osijek-196.4179.2041.3399.33
Kopački rit/Zagreb-1100.2788.03 41.98118.11
Kopački rit/Sisak-1 103.9688.06 42.85127.62
Kopački rit/Kutina-1100.7682.6438.34105.93
Plitvice/Kutina-149.51 46.8740.0952.17
Plitvice/Zagreb-149.27 49.9343.9058.18
Plitvice/Kutina-149.51 46.8740.0952.17
Višnjan/Zagreb-151.5654.8146.92 64.35
Višnjan/Rijeka-188.70 74.7081.90 96.05
Plitvice/Rijeka-184.7668.04 76.6286.84

The annual, inter-annual and daily variations of PM concentrations as well as their spatial distribution were further investigated. The yearly course of c_ (PM10 )month exhibited seasonal variations with the highest c_ (PM10 )month values during the winter and late autumn at all continental urban stations (selected stations shown in Fig. 3). A decreasing trend in PM concentrations is found in Zagreb-1 and Kutina-1 while at coastal urban station Rijeka-2 trend was not identified. Significant deviations in c_ (PM10 )day yearly course at continental ur- ban stations within analysed years were not identified, although c_ (PM10 )month levels differed substantially among stations depending on the station location and the influence of local emission sources. The c_ (PM ) ~ 60 μg/m3 was present during the colder part of the year while significantly lower 3 c_ (PM10 )month ~ 25 μg/m were present during the warmer part of the year. At 3 Rijeka-2 almost invariable c_ (PM10 )month ~ 20 μg/m was found. Seasonal and inter annual trend in c_ (PM10 )month at both continental and coastal rural stations that are representative of regional and long-range transport was not observed.

In order to further investigate the portion of PM2.5 in PM10 the monthly PM2.5 / PM10 average ratios and PMcoarse (i.e. PM10 minus PM2.5) concentrations were calculated from the observed hourly concentrations at rural background stations during 2011 (Fig. 5). Calculated yearly PM2.5/PM10 ratio values were ~ 0.7 and 0.8 at all rural background stations indicating that rural PM is mainly consisted of small particles as the contribution of PM2.5 concentrations to PM10 is ~ 80% at all locations. The highest, almost constant, PM2.5 / PM10 ratio ~ 0.9 (Fig. 5 upper panel) was found at Hum where PM concentrations were the lowest indicating that almost all PM is in the form of PM2.5. This is in agreement with the results of other studies that showed generally higher coarse fraction at the urbanized sites compared to the rural background locations (e.g. Querol et al., 2008). Since station Hum is far away from anthropogenic pollution sources and PM levels are mainly influenced by the natural sources such as sea-salt, it is likely that fine particles are formed through secondary aerosol production in the chemical reactions with primary aerosols such as DMS (dimethyl sulphide) (e.g. O’Dowd and De Leeuw, 2007). Seasonal trend in PM2.5 / PM10 ratio was also ascer- tained at other rural background stations in Croatia. During the colder part of the year PM2.5 / PM10 was higher than 0.9 at almost all stations and lower (in the range between ~ 0.8 at Žarkovica and Parg stations to ~ 0.6 at Kopački rit, Polača, Plitvice and Višnjan) during the warmer part of the year with pronounced mini- mum in August 2011. Seasonal trend in PM2.5 / PM10 ratio indicates that during winter the contribution of fine particles in rural background levels is significant while during summer the contribution from coarse particle intensifies. Monthly PMcoarse concentrations (Fig. 5 lower panel) are complementary to the ratios with the lowest PMcoarse values at Hum. Seasonal PMcoarse variations were opposite to PM2.5 / PM10 course with lower during winter and higher PMcoarse during summer. An increase in PM10 coarse concentrations was evident at Kopački rit and Polača during November, 2011 that may be related to the strong northeasterly wind episode when resuspension of crustal material from the ground contributed to the coarse component of the PM mass.

Figure 5. Average monthly PM2.5 / PM10 ratios (upper panel) and PMcoarse mass concentrations (lower panel) calculated from the observed hourly concentrations at rural background stations in Croatia during 2011.

Figure 6. Relation between annual PM2.5 / PM10 ratio and annual average PM concentrations at ru- ral background stations, Croatia, 2011.

Figure 7. Scatter plots between averaged daily PM2.5 mass concentrations and PMcoarse (upper panels) and daily PM10 and PM2.5 mass concentrations (lower panels) at rural background stations Kopački rit (left-red) and Hum (blue-right), Croatia, 2011.

In Fig. 6 the relation between annual PM2.5 / PM10 ratio and annual average PM10 concentrations at rural background stations is shown. Ratio of PM2.5 / PM10 is systematically decreasing with the increase in PM10 mass concentrations. The lowest PM10 concentrations with the highest PM2.5 / PM10 ratios were found for Hum. The lowest PM2.5 / PM10 ratio was present at Kopački rit with the highest PM10 concentrations mainly affected by regional transport from anthropogenic sources in the area.

Scatter plots (based on Putaud et al. (2010) but for daily average values) were determined between the daily average PM10 and PM2.5 mass concentrations and PM2.5 and PMcoarse. These are shown for two most distinct stations in rural back- ground network, i.e. Kopački rit and Hum (Fig. 7). Generally, PM10 and PM2.5 concentrations were exceptionally well correlated (R2 > 0.9) at all rural back- ground sites. Correlations between PM2.5 and PM10 concentrations indicate that PM10 and PM2.5 have the significant shared influence of meteorology, dispersion and transport as in Putaud et al. (2010). However, correlation between PM2.5 and PMcoarse was not established showing that there was no relationship between PM2.5 and PMcoarse concentrations which is valid for all background sites. This means that PMcoarse concentrations can generally not be inferred from PM2.5 mea- surements and that PM2.5 and PMcoarse are controlled by different sources and processes.

3.2. The analyses of meteorological parameters in relation to PM10 concentrations in Croatia in the period from 2006 to 2014

3.2.1. Temperature, relative humidity and PM concentrations

In Figs. 8 and 9 normalised average hourly and monthly values of tempera- ture, relative humidity and PM concentrations are shown for selected rural and urban coastal and continental stations. Normalising is achieved by dividing the value by its mean value. Also shown on the plots is the 95% confidence interval in the mean. Daily courses in temperatures indicate higher variability at conti- nental stations while at coastal rural and urban stations average hourly temper- ature values exhibited lower variability. The largest differences between daily temperature maximum and night time temperature minimum values was ob- served at Plitvice, Parg and Kopački rit which are surrounded by forest. The smallest differences in hourly values during a day was observed at Hum and Žarkovica. Maximum average hourly temperatures were present in Polača ~ 20 °C while at Parg and Plitvice the lowest average daily temperatures ~ 12 °C were observed. Annual course was similar at all stations with minimum values during the colder part of the year and maximum values in July and August.

Average hourly relative humidity daily course was similar at rural and ur- ban continental stations with minimum values during early afternoon and maxi- mum during night while at Hum (~ 76%) and Žarkovica (~ 65%) humidity was high and almost constant during the day which is a consequence of the sea influence. Similar was observed at urban Rijeka station but with more variability in hourly relative humidity values due to urban influences. Maximum values in hourly humidity was observed at 7 local standard time (LST) at Kopački rit ~ 9% at Plitvice and ~ 90% and at Parg ~ 85%. The lowest hourly humidity was ob- served at Žarkovica. Annual course at Plitvice, Parg and Kopački rit indicate higher values during the colder part of the year a decrease in March, almost con- stant values during the warmer part of the year and an increase in autumn. At other stations average monthly humidity continuously decreased from maximum in February toward a minimum in August and again a sharp increase in autumn occured. Annual course of average monthly relative humidity at urban continen- tal station was similar to those at Plitvice, Parg and Kopački rit stations.

Figure 8. Normalised average hourly (left) and monthly (right) values of temperature (red line), rela- tive humidity (green line) and PM concentrations (blue line) are shown for selected stations: Žarkovica – coastal rural (top) and Rijeka-2 – coastal urban station (bottom) calculated from the observations in the period from 2011 to 2014.

Figure 9. Same as Fig. 8, but for Kopački rit – continental rural station (top) and Kutina-1 – urban/ industrial station (bottom).

Daily course in hourly PM concentrations indicate two maximums at all urban and at rural Žarkovica and Polača stations one at 7 LST and the second in the afternoon around 6 PM. Hourly maximum values differed significantly between urban and rural stations with the highest values ~ 45 μg/m3 at Kutina-1 while at Polača ~ 16 μg/m3 and at Žarkovica ~ 8 μg/m3 were observed. Different daily course was found at Parg, Plitvice and Višnjan where a continuous increase in hourly PM concentrations was detected staring at 6 LST from minimum val- ues ~ 13 μg/m3 and reaching a maximum value ~ 17 μg/m3 around 18 LST after which a decease started. At Hum and Kopački rit the daily course was opposite with maximum values around 6 LST and minimum in the afternoon around 15 LST. However daily variability at Hum was significantly lower than in Kopački rit. Differences in daily maximum and minimum PM values was ~ 1 μg/m3 at Hum and 16 μg/m3 at Kopački rit. Yearly course was similar at all rural and ur- ban stations with maximum values during the colder part of the year and minimum during the warmer part of the year except for rural coastal station Žarkovica where yearly course was opposite with maximum in July and urban coastal station Rijeka-2 where maximum also during June and July.

3.2.2. Wind roses and bivariate polar plots

In Fig. 10 wind roses for Kutina-1 in the period from 2006–2014 and Hum for the period between 2011-2014 are shown. The city of Kutina-1 is situated in the continental, plain part of Croatia called Slavonija in Pannonian valley that is mainly characterized with lowland and flat horizontal surface. Wind rose indi- cate the prevalence of mainly light north-easterly wind with the wind speeds ranging from 0–2 m/s. Wind rose show that the most frequent wind is from NW and SE directions covering the whole range of wind speed values up to 28 m/s and much stronger than in Kutina. This is expected as coastal station Hum is situated in a complex orography where the impact of strong winds sirocco (SE) and etesian (NW) are present. The wind rose is a very useful way of summarising meteorological data particularly for showing how wind speed and wind direction conditions vary by year.

Figure 10. Wind roses for Kutina-1 (left) in the period from 2006 to 2014 and Hum (right) in the pe- riod from 2011 to 2014.

The bivariate polar plots indicating the variation of observed PM concentra- tions by wind speed and wind direction are plotted for Kopački rit and Hum in Fig. 11. For Kutina-1 the potential sources of PM concentrations are situated SE from the monitoring site where highway and petrochemical industry are located. In Hum bivariate polar plot indicates eastern wind sector as the most fequent for the occurence of the higher concentrations which are ~ 6 μg/m3 and significantly lower than in Kutina-1. The similar analysis of wind roses and bivariate polar plots at all sites (not shown) indicated that the occuence of higher concentrations is not usually related to the most frequent wind sectors. Similarly as for Hum at Žarkovica the highest PM concentrations are occurring during SE winds with wind speeds ~ 5m/s. At urban coastal station Rijeka-2 the highest concentrations are originating from south, while in urban continental station Osijek-1 two situa- tions with high concentrations are identifiend: one taking place during the low wind speed conditions indicating local influences and another occuring with higher wind speds from N and NE sectors. In Kopački rit, rural station relatively near Osijek-1 higher concentrations occurred during the light wind conditions. However it should be pointed that measurements in Kopački rit indicated signifi- cant underestimation of winds from northern sector which might be the conse- quence of inadequate siting conditions and will be further investigated. At Zagreb-1 station that is dominantly under the influence of traffic, higher mean concentrations are present also during the low wind speeds and winds from NW and E, but the highest concetration are clearly under the influence of local sources i.e. traffic. At Parg low wind speeds < 5 m/s from SE sector and in Plitvice NE winds with speeds in the range between 8 and 10 m/s mostly contributed to the highest observed concentrations while at Polača those were pesent primarely during the calm wind conditions.

Figure 11. Same as Fig. 10, but for bivariate polar plots.

4. Discussion and conclusions

The aim of the research was dedicated to the analyses of the urban and rural PM10 and PM2.5 concentrations in Croatia. The analyses of air pollution with PM were conducted due to: high observed PM levels, large regional emission sources and limited regional research studies. Datasets on PM mass concentrations are acquired with different measurement approaches: gravimetric method, β-attenuation and by the optical particle counters of Grimm 180 instrument from 12 sites among which 7 new rural background stations with the measurements of PM10 and PM2.5 in Croatia. Main conclusions are:

  1. Calculated annual average c_ (PM10 )year were considerably below the de- fined limit value of 40 μg/m3 at all available urban and rural stations in Croatia and decreasing trend in average annual values as well as in the number of excess days was observed. High level of PM2.5 up to 80% at con- tinental rural background station comparable to the PM2.5 level at urban background site was found which is in agreement with the research re- sults established for Central Europe (e.g. Juda-Rezler et al., 2011). European aerosol phenomenology studies conducted by several authors (Van Dingenen et al., 2004; Putaud et al., 2004; Putaud et al., 2010) have shown that the annual background PM10 and PM2.5 concentrations for con- tinental Europe are strongly affected by the regional aerosol transport. Significant contribution up to 76% of the long range transport in the urban PM2.5 concentration was found in Finland (Pakkanen et al., 2001; Karppinen et al., 2004).

  2. The pronounced spatial gradient in PM10 and PM2.5 concentrations de- creasing from the north to the south, namely from the continent toward the coast for urban and rural sites in Croatia was revealed both in average concentrations and the number of exceedance. Characteristic geographic and climatological conditions as well as intrinsic atmospheric processes influencing PM concentrations in the air such as: local wind circulation, vertical diffusion and deposition were considerably different at coast than inlands.

  3. As in other parts in Europe, a significant correlation between the PM2.5 and PM10 average mass concentrations calculated at 7 rural background sites in Croatia was found indicating that both concentrations are govern by the same meteorological processes. Also PM2.5 / PM10 ratios at rural sta- tions differ from each other to a larger extent during the warmer part of the year depending on the level of PM concentrations. There is no relation- ship between PM2.5 and PMcoarse at rural background stations meaning that those are controlled by different processes. The content of PM2.5 in PM10 was high at all rural background stations and a gradient of PM2.5 / PM10 oriented from the south to the north, i.e., from the coast to the continent, inverse to the gradient of PM10 and PM2.5 concentration, was identified.

  4. Rural site-specific PM2.5 / PM10 ratios ranged from 0.6 to 0.9 during the warmer part of the year while during the colder part of the year ratios ranged from 0.85 to 0.98 which is in accordance with conclusions of Putaud et al. (2010). They found that PM2.5 / PM10 ratios increase with the increase of PM10 levels at single station which is generally the case in Europe during the colder part of the year. The seasonal variability in ra- tios is less pronounced for the stations with lower PM concentrations in- dicating the same PM2.5 composition and contributing sources during the year. However, our results lead to conclusion that the contribution of PM2.5 in PM10 is higher and less variable during the year for the lower PM10 concentrations at rural background stations and that ratio is not undependable of PM10 levels. This is reasonable since higher contribu- tion of the PMcoarse is found mainly at kerbside urban traffic stations in Europe.

  5. Daily courses in average hourly temperatures indicated higher variability at continental than at coastal stations while annual course of temperature was similar at all stations with minimum values during the colder part of the year and maximum values in July and August opposite to annual courses of relative humidity and PM. The only exception in yearly PM courses was found at Žarkovica where maximum was present during sum- mer. This might be a consequence of the summer forest fires frequently taking place at coast. Average daily course of hourly relative humidity was similar at rural and urban continental stations with minimum values during early afternoon and maximum during night while at Hum (~ 76%) and Žarkovica (~ 65%) humidity was high and almost constant during the day which is a consequence of the sea influence. Bivariate polar plots indi- cated location of sources which will be further investigated in future studies.


This work has been supported by Croatian Science Foundation (grants »CATURBO« No. 09/151)



Alebić-Juretić A.; Cvitaš T.; Kezele N.; Klasinc L.; Pehnec G.; Šorgo G. (2007): Atmospheric particulate matter and ozone under heat-wave conditions: Do they cause an increase of mortality in Croatia?, B. Environ. Contam. Tox., 79, 468-471, DOI:


Allen G.; Sioutas C.; Koutrakis P.; Reiss R.; Lurmann F. W.; Roberts P. T. (1997): Evaluation of the TEOM® method for measurement of ambient particulate mass in urban areas, J. Air Waste Manag. Assoc., 47, 682-689, DOI:


Amato F.; Pandolfi M.; Escrig A.; Querol A.; Alastuey A.; Pey J.; Perez J.; Hopke P. K. (2009): Quantifying road dust resuspension in urban environment by Multilinear Engine: A comparison with PMF2, Atmos. Environ., 43, 2770-2780, DOI:


Anderson H. R. (2009): Air pollution and mortality: A history, Atmos. Environ., 43, 142-152, DOI:


Andrea M. O.; Jones C. D.; Cox P. M. (2005): Strong present-day aerosol cooling implies a hot future , Nature, 435, 1187-1190, DOI:


Bešlić I.; Šega K.; Čačković M.; Klaić Z. B.; Vučetić V. (2007): Influence od weather types on concentrations of metallic components in airborne PM10 in Zagreb, Croatia, Geofizika, 24, 93-107.


Bešlić I.; Šega K.; Čačković M.; Klaić Z. B.; Bajić A. (2008): Relationship between 4-day air mass back trajectories and metallic components in PM10 and PM2.5 particle fractions in Zagreb air, Croatia, B. Environ. Contam. Tox., 80, 270-273, DOI:


Carslaw D. C.; Beevers S. D.; Ropkins K.; Bell M. C. (2006): Detecting and quantifying aircraft and other on-airport contributions to ambient nitrogen oxides in the vicinity of a large interna- tional airport, Atmos. Environ., 40, 5424-5434, DOI:


Carslaw D. C.; Ropkins K. (2012): OpenAir – An R package for air quality data analysis, Environ. Modell. Softw., 27–28, 52-61, DOI:


Carslaw D. C. (2012): The OpenAir Manual – Open-source tools for analysing air pollution data, Manual for version 0.6-0 , King’s College London available at


Davila S.; Šega K.; Bešlić I. (2011): Comparison of gravimetric and β-attenuation methods for PM10 concentration measurement at the monitoring station Zagreb 1, in: Zbornik sažetaka Sedmog hrvatskog znanstveno-stručnog skupa “Zaštita zraka ‘11”, edited by: Šega, K., ZagrebHrvatsko udruženje za zaštitu zraka, 78-79.


EC (2008): 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, available at


EC (2010): Guide to the demonstration of equivalence of ambient air monitoring methods, Report by an EC Working Group on Guidance for the Demonstration of Equivalence 92 pp, available at


Eisner A. D.; Wiener R. W. (2002): Discussion and evaluation of the volatility test for equivalen- cy of other methods to the federal reference method for fine particulate matter, Aerosol Sci. Technol., 36, 433-440.


Harrison R. M.; Shi J. P.; Jones A. M. (1999): Multisite study of particle number concentrations in urban air, Evniron. Sci. Technol., 39, 6063-6070, DOI:


Hauck H.; Berner A.; Gomiscek B.; Stopper S.; Puxbaum H.; Kundi M. (2004): On the equiva- lence of gravimetric PM data with TEOM and beta-attenuation measurements, J. Aerosol Sci., 35, 1135-1149, DOI:


Forsberg B.; Hansson H. C.; Johansson C.; Areskoug H.; Persson K.; Järvholm B. (2005): Comparative health impact assessment of local and regional particulate air pollutants in Scandinavia, Ambio, 34, 11-19.


Jeričević A.; Kraljević L.; Vidič S.; Tarrason L. (2007): Project description: High resolution envi- ronmental modelling and evaluation programme for Croatia (EMEP4HR), Geofizika, 24, 137-143.


Jeričević A.; Kraljević L.; Grisogono B.; Fagerli H.; Večenaj Ž. (2010): Parameterization of verti- cal diffusion and the atmospheric boundary layer height determination in the EMEP model, Atmos. Chem. Phys., 10, 341-364, DOI:


Jeričević A., Ilyin I., Vidič S. (2012): Modelling of heavy metals: Study of impacts due to climate change, in:National Security and Human Health Implications on Climate Change, edited by: Fernando, H. J., Klaić, Z. B., and McCulley, J. L. Springer, Dordrecht, The Netherlands 175-189.


Jiang H.; Liao H.; Pye H. O. T.; Wu S.; Mickley L. J.; Seinfeld J. H.; Zhang X. (2013): Projected effect of 2000–2050 changes in climate and emissions on aerosol levels in China and associated transboundary transport, Atmos. Chem. Phys. Disc., 13, 6501-6551, DOI:


Juda-Rezler K.; Reizer M.; Oudinet J.-P. (2011): Determination and analysis of PM10 source apportionment during episodes of air pollution in Central Eastern European urban areas: The case of wintertime 2006, Atmos. Environ., 45, 6557-6566, DOI:


Karppinen A.; Härkönen J.; Kukkonen J.; Aarnio P.; Koskentalo T. (2004): Statistical model for assessing the portion of fine particulate matter transported regionally and long range to urban air, Scand. J. Work Environ. Health, 30(2), 47-53.


Luhana L.; Sokhi R.; Lynne W.; Mao H.; Boulter P.; McCrae I.; Wright J.; Osborn D. (2004): Measurement of non-exhaust particulate matter, in: Deliverable 8 from EU Project PARTICULATES – Characterization of Exhaust Particulate Emissions from Road Vehicles, European Commission 5th Framework Programme 96 pp, available at port/roadsafety_library/publications/particulates_d8.pdf


Noble C. A.; Vanderpool R. W.; Peters T. M.; McElroy F. F.; Gemmill D. B.; Wiener R. W. (2001): Federal reference and equivalent methods for measuring fine particulate matter, Aerosol Sci. Technol., 34, 457-464, DOI:


O’Dowd C. A.; De O’Dowd G. (2007): Marine aerosol production: A review of the current knowledge, Philos. T. Roy. Soc. A, 365, 1753-1774, DOI:


Pakkanen T. A.; Loukkola K.; Korhonen C. H.; Aurela M.; Mäkelä T.; Hillamo R.; Aarnio P.; Koskentalo T.; Kousa A.; Maenhaut W. (2001): Sources and chemical composition of atmo- spheric fine and coarse particles in the Helsinki area, Atmos. Environ., 35, 5381-5391, DOI:


Pang Y., Eatough N. L., Eatough D. J. (2002): PM2.5 semivolatile organic material at Riverside, California: Implications for the PM2.5 federal reference method sample, Aerosol Sci. Technol., 36, 277-288, DOI:


Peters A.; Dockery D. W.; Muller J. E.; Mittleman M. A. (2001): Increased particulate air pollution and the triggering of myocardial infarction, J. Am. Heart Assoc., 103, 2810-2815, DOI:


Pope C. A.; Burnett R. T.; Thun M. J.; Calle E. E.; Krewski D.; Ito K.; Thurston G. D. (2002): Lung cancer, cardiopulmonary mortality, and long-term exposure to fine particulate air pollution, JAMA – J. Am. Med. Assoc., 287, 1132-1141, DOI:


Putaud J.-P.; Van Dingenen R.; Alastuey A.; Bauer H.; Birmili W.; Cyrys J.; Flentje H.; Fuzzi S.; Gehrig R.; Hansson H. C.; Harrison R. M.; Herrmann H.; Hitzenberger R.; Huglin C.; Jones A. M.; Kasper-Giebl A.; Kiss G.; Kousa A.; Kuhlbusch T. A. J.; Loschau G.; Maenhaut W.; Molnar A.; Moreno T.; Pekkanen J.; Perrino C.; Pitz M.; Puxbaum H.; Querol X.; Rodriguez S.; Salma I.; Schwarz J.; Smolik J.; Schneider J.; Spindler G.; ten Brink H.; Tursic J.; Viana M.; Wiedensohler A.; Raes F. (2010): A European aerosol phenomenology—3: Physical and chem- ical characteristics of particulate matter from 60 rural, urban, and kerbside sites across Europ, Atmos. Environ., 44, 1308-1320, DOI:


Putaud J. P.; Raes F.; Van Dingenen R.; Brüggemann E.; Facchini M. C.; Decesari S.; Fuzzi S.; Gehrig R.; Hüglin C.; Laj P.; Lorbeer G.; Maenhaut W.; Mihalopoulos N.; Müller K.; Querol X.; Rodriguez S.; Schneider J.; Spindler G.; ten Brink H.; Tørseth A.; Wiedensohler A. (2004): A European aerosol phenomenology—2: Chemical characteristics of particulate matter at kerbside, urban, rural and background sites in Europe, Atmos. Environ., 38, 2579-2595, DOI:


Querol X., Alastuey A., Viana M. M., Rodriguez S., Artíñano B., Salvador P., Garcia do Santos S., Fernandez Patier R., Ruiz C. R., de la Rosa J., Sanchez de la Campa A., Menéndez M., Gil J. I. (2004): Speciation and origin of PM10 and PM2.5 in Spain, J. Aerosol Sci., 35, 1151-1172, DOI:


Querol X.; Pey J.; Minguillón M. C.; Pérez N.; Alastuey A.; Viana M.; Moreno T.; Bernabé R. M.; Blanco S.; Cárdenas B.; Vega E.; Sosa G.; Escalona S.; Ruiz H.; Artíñano B. (2008): PM speciation and sources in Mexico during the MILAGRO-2006 Campaign, Atmos. Chem. Phys., 8, 111-128, DOI:


Samet J. M.; Dominici F.; Curriero F. C.; Coursac I.; Zeger S. L. (2000): Fine particulate air pollution and mortality in 20 U.S. cities, 1987–1994, N. Eng. J. Med., 343, 1742-1749, DOI:


Samoli E.; Analitis A.; Touloumi G.; Schwartz J.; Anderson H. R.; Sunyer J.; Bisanti L.; Zmirou D.; Vonk J. M.; Pekkanen J.; Goodman P.; Paldy A.; Schindler C.; Katsouyanni K. (2005): Estimating the exposure-response relationships between particulate matter and mortality within the APHEA multicity project, Environ. Health Perspect., 113, 88-95, DOI:


Šišović A., Bešlić I., Šega K., Vađić V. (2008): PAH mass concentrations measured in PM10 particle fraction, Environ. Int., 34, 580-584, DOI:


Van Dingenen R.; Raes F.; Putaud J. P.; Baltensperger U.; Charron A.; Facchini M. C.; Decesari S.; Fuzzi S.; Gehrig R.; Hansson H.-C.; Harrison R. M.; Hüglin C.; Jones A. M.; Laj P.; Lorbeer G.; Maenhaut W.; Palmgren F.; Querol X.; Rodriguez S.; Schneider J.; ten Brink H.; Tunved P.; Tørseth K.; Wehner B.; Weingartner E.; Wiedensohler A.; Wåhlin P. (2004): A European aerosol phenomenology—1: Physical characteristics of particulate matter at kerbside, urban, ru- ral and background sites in Europe, Atmos. Environ., 38, 2561-2577, DOI:


Viana M.; Kuhlbusch T. A. J.; Querol X.; Alastuey A.; Harrison R. M.; Hopke P. K.; Winiwarter W.; Vallius M.; Szidat S.; Prévôt A. S. H.; Hueglin C.; Bloemen H.; Wåhlin P.; Vecchi R.; Miranda A. I.; Kasper-Giebl A.; Maenhaut W.; Hitzenberger R. (2008): Source apportion- ment of particulate matter in Europe: A review of methods and results, J. Aerosol Sci., 39, 827-849, DOI:


Westmoreland E. J.; Carslaw N.; Carslaw D. C.; Gillah A.; Bates E. (2007): Analysis of air qua- lity within a street canyon using statistical and dispersion modelling techniques, Atmos. Environ., 41, 9195-9205, DOI:


Willeke K.; Baron P. A. (1993): Aerosol Measurement: Principles, Techniques, and Applications. First Edition Van Nostrand ReinoldNew York 977 pp.

This display is generated from NISO JATS XML with jats-html.xsl. The XSLT engine is libxslt.