IS428 AY2018-19T1 Zhuo Yunying

Problem & Motivation

Air pollution is an important risk factor for health in Europe and worldwide. A recent review of the global burden of disease showed that it is one of the top ten risk factors for health globally. Worldwide an estimated 7 million people died prematurely because of pollution; in the European Union (EU) 400,000 people suffer a premature death. The Organisation for Economic Cooperation and Development (OECD) predicts that in 2050 outdoor air pollution will be the top cause of environmentally related deaths worldwide. In addition, air pollution has also been classified as the leading environmental cause of cancer. Air quality in Bulgaria is a big concern: measurements show that citizens all over the country breathe in air that is considered harmful to health. For example, concentrations of PM2.5 and PM10 are much higher than what the EU and the World Health Organization (WHO) have set to protect health. Bulgaria had the highest PM2.5 concentrations of all EU-28 member states in urban areas over a three-year average. For PM10, Bulgaria is also leading on the top polluted countries with 77 μg/m3on the daily mean concentration (EU limit value is 50 μg/m3). According to the WHO, 60 percent of the urban population in Bulgaria is exposed to dangerous (unhealthy) levels of particulate matter (PM10).

This assignment aims to study the following:

• Task 1: Spatio-temporal Analysis of Official Air Quality
• Task 2: Spatio-temporal Analysis of Citizen Science Air Quality Measurements
• Task 3: Find out the relationship of the above analysis with other factors (Local energy sources, Local meteorology, Local topography, Complex interactions between local topography and meteorological characteristics and Transboundary pollution)

Dataset Analysis & Transformation Process

1. Combine all time-series data (e.g. BG_5_9572_2017_timeseries) into one single excel spreadsheet

• Create a new "Station" Column (indicating station code i.e. 9421) and "Year" Column based on the year for each of the time-series data
• Use excel to combine the rest of time-series files based on all the common columns (e.g. Countrycode, Namespace, AirQualityNetwork and etc.)
• Based on the analysis of the existing combined time-series data from 2013 to 2018, there is a drastic difference in the level of aggregation across the years. As shown in the table below, Year 2016 has a combination of hourly air quality readings and daily air quality readings while in 2017, there is a mixture of hour and var readings as for certain days, readings are not measured at one-hour interval continuously. As such, the analysis on air quality readings will be based on "Day" averaging time from year 2013 to 2016 as it is impossible to lower the aggregation level of "Day" to "Hour" readings in 2016. On the other hand, the analysis on air quality readings will be based on "Hour" averaging time for 2017 and 2018. For the data in 2017,average readings will be taken if the readings for any specific days are done on "Var" basis. Hourly analysis is also more accurate for both 2017 and 2018 as there are missing data on specific months. (In 2017, only Nov and Dec data are available while in 2018, only Jan to Sep's data are available)

• Due to the standardization in averaging time, the values under "Concentration", "DatetimeBegin" and "DatetimeEnd" have been adjusted accordingly while other column values remain unchanged.
• There are quite a number of duplicated readings in the dataset. These duplicated readings are removed during the transformation process to avoid unequal weightage.
• Due to the high variation in raw data and small dataset, excel is used for the transformation process.

2. Merge metadata file with combined time-series data

• According to the source of scrapped data (http://discomap.eea.europa.eu/map/fme/AirQualityExport.htm),"the join between time-series files and the metadata file should be made using the Countrycode (or Namespace) and SamplingPoint". Thus, metadata file and time-series data time are merged via Vlookup function in excel based on SamplingPoint (since Countrycode are all "BA" for both data file).
• Upon further inspection, these two datasets have a number of common columns with same values. These columns include "Countrycode, Namespace, AirQualityNetwork, AirQualityStation, AirQualityStationEoICode, SamplingProcess, AirPollutantCode, AirPollutant" and thus repeated columns are removed.

Task 2: Spatio-temporal Analysis of Citizen Science Air Quality Measurements

• Geohash represents station locations. However, Tableau is not able to interpret geohash as geographic data. Before Air Tube data is imported to Tableau for analysis, geohash needs to be decoded into geographical coordinates. Due to the sheer size in Air Tube datasets(data_bg_2017.xlsx and data_bg_2018.xlsx), R packages ("devtools", "tidyverse" and "geohas") as indicated by Prof Kam will be used to transform the geohashed raw data to a new csv file containing the coordinates of the locations.
• Next, 2017 and 2018 datasets that have been transformed will be combined based on the R programming code below.

• The combined datasets comprises of readings across the whole Bulgaria area which is beyond the scope of our concern (i.e. Sofia). To maintain the consistency in comparison, Inclusion filter function is used to filter in only longitudes and latitudes that are in Sofia city. Due to the difficulty in comparing with the real geographical area of Sofia, an existing map on Sofia City (from Open Street Map) was used to filter through the Tableau filter function. The folowing shows map of Sofia and the filtering function.

Task 3: Find out the relationship of the air quality analysis with other factors

• The Meteorological data scrapped is currently in crosstab format which is not suitable for analysis on Tableau. Hence, the dataset has to be transformed to columnar format. Hence, pivoting will be done to transform the data.
• As the data includes average, minimum and maximum readings for different types of measurements and are presented on separate columns. The format is not suitable to be processed on Tableau if there is a need to introduce filters in the dashboard. Hence the dataset is transformed to as shown below.

Interactive Visualization

The interactive visualization can be accessed here:

Official Air Quality in Sofia City

This dashboard shows the daily and hourly PM10 concentration in the 4 measured stations. A user can select a day on the heat map and investigate PM10 concentration at different times of the day.

Sensor Coverage, Performance and Operation

This dashboard shows sensors and their readings in Sofia City. When selecting a sensor, the line charts will display the sensor readings at different times of the day. A user thus can tell when the sensors start not working properly.

Air Quality Measures in Sofia City

This dashboard helps user see the flow of air pollution from day to day and analyse the impacts of local topography and meteorology on air quality.

Relationship between Local Meteorology and Air Pollution

The last dashboard helps user find the relationship between local meteorology and air pollution. A user can select a few days and compare the patterns among temperature, pressure, rainfall, humidity, wind and air quality.

Interesting & Anomalous Observations

Task 1: Spatio-temporal Analysis of Official Air Quality

1.1 Characterize Recent and Past Situation of Air Quality in Sofia City The figures above shows the PAST average daily concentration heat map for year 2013 to 2016 as daily readings are not computed for year 2017 and year 2018. The calendar heat map visualization allows users to appreciate the trend in average daily concentration from Jan to Dec at one glance. The calendar map also includes filter by Sampling location, air quality range as well as Year. With reference to EU Air Quality classification based on research, the air quality range includes 0-20:Good, 20-40:Moderate, 40-50:Normal, 50-75:Unhealthy, 75 and above: Very Unhealthy.The heat map’s colour intensity is specified within a range such that dark red represent unhealthy and very unhealthy air quality range and vice versa. Trends Analysis: From year 2013 to year 2016, the air quality exceeds EU limit (50 µg/m3) at the end of each year and at the beginning of each year. Specifically, the air quality reached unhealthy and very unhealthy range for majority of the days In January, February, November and December across all four years. In the middle of year (i.e. Quarter 2 and Quarter3), the air quality is much better and are generally below the EU limit value of 50 µg/m3. There isn’t any clear trend in the change in daily concentration of PM10 across the days in a month throughout 2013 to 2016.

1.2 A Typical Day in Sofia City This calendar heat map shows the hourly concentration of PM 10 throughout a typical day across different months from Nov 2017 to Sep 2018. The rightmost column “Average” shows the overall average readings of hourly concentration. Filter on sampling location and air quality range could be utilized to narrow down the scope of comparison. Trends Analysis: For the most RECENT average hourly concentration of PM 10 calendar map from Nov 2017 to Sep 2018, a typical day in Sofia city does not experience drastic change in air quality on average over the months. However, the average hourly concentration across these 11 months are generally higher in the range of 30+ µg/m3 before 9am and after 6pm which are after working hours. The average hourly concentration is lower in the range of 20+ µg/m3 from 9am to 6pm. The air quality is considered good throughout the day on average. However, the air quality readings are exceptionally high throughout the day (24 hours) in both November 2017 and January 2018. The air quality has reached unhealthy and very unhealthy range. This signifies that there might be external factors that influence the air quality to be extremely poor in those months.

1.3 Further Analysis on Sampling Points TThis dashboard shows the distribution of sampling points in Sofia city as well as the readings concentration, distance to Kerb, altitude and distance to building for these sampling points. Trends Analysis: All the sampling points follow the same cyclical pattern for the time-series analysis of average daily concentration of PM10 pollutants. The cyclical pattern repeated across the years from 2013 to 2018. Across the years, the mean concentration of PM10 pollutants is the highest for Orlov Most followed by Nadezhda. Orlov Most is near to residential or industrial building while Nadezhda is near not as near as Orlov Most to buildings. The variation in concentration is highest for Hipodruma (as seen by the high range and outliers). There might be association with other factors since the distance to kerb, altitude and distance to building is similar to most of other sampling points.

1.4 What anomalies do you find in the official air quality dataset? Inconsistent measurement on air quality data: The consistency occurs in the way they measure the air quality as hourly, var and daily concentration values are all available. Incomplete data for robust time-series analysis: There are a few months (Jan 2017 to Oct 2017) that do not have any air quality data. There is a need to find out why. Inconsistent coverage for Sampling Points: For Mladost (60881), only has one year of air quality reading available while for Orlov Most (9484), only have readings from year 2013 to 2015.

How do these affect your analysis of potential problems to the environment? Need to consider the appropriate measurement time method when comparing the air quality concentration for time-series analysis The data points for sampling points might be too small (only 4 points with consistent readings) and thus affect the accuracy of analysis. The poor air quality might be related to (1) the surrounding regions and facilities around the sampling points, (2) cyclical patterns across the months in different years and (3) Consumption behaviours of residents living in Sofia city

Task 2: Spatio-temporal Analysis of Citizen Science Air Quality Measurements

2.1 Sensor Coverage, Performance and Operation 2.1.1 Sensor Coverage The most recent sensor distribution (on 7 April 2018) shows that the sensors are not evenly spread out across the entire city. They are saturated in the center of the city. Nevertheless, the sensors did cover the central part of the city well. 2.1.2 Sensor Performance and Operation However, the sensors were not always working properly. As shown in the graphs above, fora example, on 3 February 2018, the sensor at (42.625, 23.321) was working properly. Its readings of pressure, humidity and temperature fluctuate within reasonable ranges though the day. However, on 4 February 2018, the sensor temperature reading had a sharp drop from 7 degree to -139 degree at 9AM, which is not reasonable. At the same time, pressure and humidity readings also increased sharply, suggesting that the sensor might have a problem on 9AM, 4 February 2018. The suspect was approved by the static sensor readings on 5 February 2018, suggesting that the sensor was not working. Similar sensor behaviors can also be seen on other sensors such as the sensor located at (42.686, 23.347) after 14 February 2018 and the sensor located at (42.715, 23.161) after 7 February 2018.

2.2 Air Pollution Most of the time, pollution started from the northern part of the city. As time passes, the pollution can spread in different directions. Case 1: Pollution spread to the south Case 2: Pollution spread to the south first and then back to the north

Task 3: Factors Affecting Air Quality in Sofia City

Factor 1: Local Energy Sources The diagrams above shows that in cities like Sofia and Blagoevgrad, P1 concentration is relatively high. This might be because these cities are along the main gas pipeline system. Since the cities are rich in hard coal, the gas and energy generated in these cities might be from burning of hard coals, which damages the air.

Factor 2: Local Meteorology With the given Meteo-data, which describes the temperature, wind speed, humidity, pressure and rainfall in Sofia Airport (42.6537, 23.3829 ), we can investigate the relationships between meteorology and air pollution. Since station STA-BG0052A(42.6665, 23.4002) is closest to Sofia Airport compared with other three stations, I use data in station STA-BG0052A for comparison. The diagrams above compare the patterns of temperature, wind speed, humidity, pressure and rainfall in Sofia Airport (42.6537, 23.3829 ) with the pattern of PM10 concentration in Station STA-BG0052A(42.6665, 23.4002) in 2014. It shows that PM10 concentration tends to be high in a low temperature, low wind speed, high humidity and high pressure environment. However, the relationship between rainfall and PM10 concentration is not very obvious. Data in 2016 also shows the same findings.

Factor 3: Local Topography As mentioned in Task 2, pollution usually starts from the northern part of Sofia city. This may be due to the fact that Sofia's altitude is low in the north and high in the southwest. PM10 particles from the north can flow into Sofia city much more easily than those from the southwest.

Factor 4: Complex interactions between local topography and meteorological characteristics The diagram shows that even though the pressure was high, humidity was high and temperature was low across the city, the pollution level is high in the northern part of the city mainly because of the low altitude in the northern part.

Factor 5: Transboundary Pollution Pollution does spread from one place to another as discussed in Task 2. The graph below also displays that when air quality in Sofia is poor, the cities around it also have places with poor air quality.

References

Understanding the current issues of poor air quality in Bulgaria https://www.eea.europa.eu/publications/air-quality-in-europe-2018/at_download/file

Map on the Bulgarian Coal Resources and Energy Infrastructure hhttps://www.researchgate.net/figure/Map-on-the-Bulgarian-coal-resources-energy-infrastructure-and-largest-power-plants_fig1_257941554

Coordinates of Sofia's Thermal Plants http://www.wikiwand.com/en/List_of_power_stations_in_Bulgaria#/Thermal