Air Pollution is the most detrimental form of pollution in our environment. A significant increase has been seen in atmospheric pollutants sources in Asia while it is considered to be primarily an urban problem in Pakistan. Sheikhupura was selected as study area and data was collected randomly from 15 different points of Sheikhupura city. Data used in this research paper is of secondary type collected from EPD, Sheikhupura. One-way-ANOVA and Interpolation (IDW and Kriging) has been applied for the interpretation and analysis of NOx and SO2 concentration at sampling sites. Results showed that concentration of NOx was very high at areas with high road traffic while that of SOx varied differently. The highest value of NOx found in the study area is 5 ppm while the lowest value is 0.06 ppm, on the other hand the highest value of SOx found in the study area is 1.9 ppm and the lowest value is 0.1 ppm. The most prone areas to NOx in Sheikhupura are Jhang Road, Jinnah Park, Rasul Pura and Landa Bazar. And the areas with high concentration of SOx in Sheikhupura are Ashraf General and Batti Hospital. Proper urban management can help to reduce air pollution.
Keywords: Air pollution; ANOVA; GIS; IDW; Kriging; NOx; SO2; Sheikhupura; Pakistan
From last some years, there is a substantial increase in development, urbanization, motorization associated with increased energy use in the Asian countries. A significant increase has been seen in atmospheric pollutants sources in Asia . Industrial activities, increased population and substantial rise in vehicular use are potential cause of environmental degradation in South Asia . As a result of these factors, air pollution has caused major threats to the environment, life quality, and health effects to the South Asian population. Poor air quality is deteriorating not only environment but also badly affecting health of individuals . A recent report of the OECD says that by the year 2050, air pollution is expected to be the world’s top environmental cause of death, prior to dirty water and poor sanitation . Asian countries are suffering from elevated levels of atmospheric pollutants: Bangladesh , Malaysia , Indonesia , Pakistan , Thailand [7,9], India [7,10], Vietnam  and China . According to , urban air quality badly affects health of individuals in metropolitan areas, and reports that there is an average annual increase of 14.1 % of vehicular emissions in Pakistan from year 1985-2005. Air pollution is considered to be primarily an urban problem in Pakistan .
Concentration of pollution has increased in major cities of Asia like in Lahore, Karachi, Kolkata, Mumbai, Chinghai, Beijing and many other cities. According to WHO report, air pollution is the sole reason for death of more than 2 million people all over the world. 90% of air pollution in urban areas is due to vehicular emissions . Pakistan has the highest concentration of particulate matter and is placed by the WHO in a category called “unhealthy for the sensitive people”. Pakistan’s urban air pollution is among the most severe in the world and it engenders significant damage to human health and the economy . Sheikhupura has developed an inspiring industrial base from the last some years(Javed et al. 2010). Depending on the Census of Manufacturing Industries, large scale industries are higher in Sheikhupura than in Lahore . Total population of Sheikhupura city is approximately 3, 63,000 residents. Population is growing in a faster rate of 1.88% per year. This increased population is associated with the increased traffic resulted in a haphazard situation in Sheikhupura city. According to a survey, 55% of the public use motorbikes and 15% cars as their means of transportation .
Geo-information Sciences is being widely used in decision making regarding transport, land use pattern and environmental quality [14,15]. GIS is a key technology that can be easily used in health impact assessment . Mapping and modeling of pollution and its analysis gives us a better idea of pollution patterns than ordinary techniques . No doubt, geo-informatics has an important role in decision making but still it is important to calculate uncertainties linked with these models . Study of air quality is a complex method because it is quite difficult to collect air samples. Geo-statistics is used for complete analysis and observation of air quality of the study area . Techniques like GIS and RS can be used in order to sort out problems like atmospheric pollution . The designed tool can be used largely by policy makers in residential, industrial and commercial zones .
Sheikhupura in Punjab Province, Pakistan is selected as study area. Sheikhupura is a famous district of Punjab, situated at a distance of 36 km northwest from Lahore. It is well-known for its historical buildings, agricultural surroundings and the best quality rice. It is provincially known as “Qila Sheikhupura”, because of the fort that was constructed by Mughal Emperor Jahangir . Sheikhupura is located on following coordinates: Latitude=31° 42’ 47" N and Longitude=73° 58’ 41" E. According to report of Punjab Development Statistics 2012, estimated population of Sheikhupura district is 2,980,000 covering an area of 3242 Km2. It comes on 3rd ranking for population and 7th for the area as compared to 36 other districts of Punjab province  (Figure 1).
Selection of sampling points
Data on air quality of Sheikhupura city from January to March for 2015 was collected from EPD Sheikhupura. 15 popular and most busy sites of Sheikhupura were observed for concentration of NOx and SOx emissions. Sampling was performed for random sites. Values of data were observed by placing instruments on 2-3 meters away from road sides (Figure 2).
SO2 monitoring data
Sulphur dioxide concentration was collected in glass impingers using sodium tetrachloro-mercurate as absorption solution using Air Sampling Pump (BD Pump, RMN26551 Model, and Manufactured by LaMotte, Maryland, USA). Three replicas for each sampling sites were collected for 24 hours observation and were transferred in to the laboratory in a cold box at 5oC. Samples were analyzed using Colorimetric method by estimating absorbance of SO2 from absorbing reagent using Spectrophotometer at 540 nm.
NOx monitoring data
Nitrogenous Oxides were collected in glass impingers by using sodium hydroxide and sodium arsenate as absorption solutions using Air Sampling Pump (BD Pump, RMN26551 Model, and Manufactured by LaMotte, Maryland, USA.). Transportation and sampling methods are same as SO2.
Data analysis and interpretation
ArcGIS 10, Google Earth and Statistix 10 were used for data analysis and interpretation. Microsoft Excel was used to arrange data according to geographic locations in Geographic Coordinate System (GCS) projection, WGS 1984 datum. ArcGIS interpret data which is geographically referenced to a projection system. Study area was Georef and digitized by using ArcGIS. Then excel sheet was converted into shape file (.shp) in order to perform analysis in GIS. Arranged data for 15 locations is given in table. It shows concentration of NOx and SOx with reference to longitudes and latitudes of sampling sites. Three replicates values are being taken for each of 15 sampling sites (Table 1).
Concentration of SOx and NOx is given in parts per million (ppm) while latitudes and longitudes are given in decimal degrees. One way ANOVA was applied between these parameters to determine whether there is a significant relation between means of three or more unrelated groups or not.
Data was organized in Excel for statistical analysis of two parameters. Each site contains three replica values for each parameter. Mean, Standard Deviation and coefficient of variance were found for each site by using replica values. Then graphs were designed with mean value and difference was finding from the standard deviation value using error bars. Difference of mean value and standard deviation showed that the mean value can deviate positively or negatively from the standard value. Statistix 10 software was used for application of one-way ANOVA on data set. The result of ANOVA shows a significant level.Values are mean of 3 different replicas bearing different letters that are significantly different from each other (at p<0.05) by using the Tukey HSD test. Mean of measured values are further compared with the NEQs.
Spatial analysis of NOx and SOx data for study area was performed by using ArcGIS 10. Interpolation was performed for estimation of values of variables at those sites which are un-sampled. Both IDW and Kriging techniques are used for spatial representation of pollutants at and between sampling sites.
The present research has inquired different pollutants and their sources, spatial patterns and analysis and exposure models in Sheikhupura. Geo-statistical analysis showed that increased concentration of atmospheric pollutants is closely related with the increased traffic but other sources are also of significant consideration for air pollutants like industries and especially brick kilns which are found abundantly in the surroundings of Sheikhupura city. Unmanaged urban transport and urban development, especially an increasing number of industries in Sheikhupura surroundings have increased environmental threats. Geo-statistical results showed that vehicular emission, increased urban settlement and industrial development is a contributing factor in degradation of ambient air quality.
NEQ’s value of NOx is 0.42 ppm (on average hourly value) while NEQs value for SOx is 0.14 ppm (based on average hourly value). The foremost reason for increased concentration of NOx in the air could be the increased traffic from the past few years. According to Transport Planning Unit Road traffic has increased in Sheikhupura up to 50% from the past few years, especially amongst the people having high and medium income. While the cause for increased concentration of SOx emissions could be use of coal and furnace oil in brick kilns, also the use of diesel oil and petrol for transport purposes. The most prone areas of NOx emissions in Sheikhupura are Jhang Road, Jinnah Park, Rasul Pura and Landa Bazar. Jinnah Park and Jhang road are highly populated roads in terms of traffic. While studies at Sarwar Shaheed Road show the lowest concentration of NOx emissions reason for this could be the decreased traffic density at this point. Results demonstrate that the areas with high concentration of SOx in Sheikhupura are Ashraf General Hospital and Batti Chowk. These two areas also have high traffic density but there are many small factories and brick kilns working in this area which also serves to the increased concentration of SOx in these areas. While Sarwar Shaheed Road has the lowest value for SOx emissions reason is again the same it is not commercial area and so contributes less towards SOx or NOx emission.
Concentration of NOx was high in study area in comparison to the NEQs (0.42 ppm based on average hourly value). Nitrogenous oxides have a visible relation with the road traffic emissions as shown in Graph 1. NO is released from the road traffic and industrial source. Trend shows that its concentration is high near roads and decreases near settlements or away from the road. It can be implicit that concentration of NOx decreases by increasing distance from roads and vice versa. Concentration of NOx at Jinnah Park is the highest while at Sarwar Shaheed Road is lowest and also values for both locations are significantly different from other groups. While sampling sites like Rasul Pura, Ayesha Park, Noor Stadium, Ghang Road, Malik Anwar Road showed almost same levels for NOx emissions (Graph 1) (Figures 3 and 4).
The concentration of Sulphur Oxides in different locations of Sheikhupura city varies from Ashraf General Hospital to Rasul pura and it exceeds from the limits of National Environmental Quality standards. Ashraf General Hospital and Batti Chowk shows the highest concentration of SOx emissions and Sarwar Shaheed Road and Damoana Road showed the lowest concentration (Graph 2).
Concentration of SO2 was not as higher in the study area in comparison to the NEQs (0.14 ppm based on average hourly value). SO2 is generated as a result of influx of vehicular emission and industrial pollution. Geo-statistical analysis has proved that SO2 produced in Sheikhupura city is also due to vehicular emission along with sources. It follows the same trends like NOx. With increase in distance from extreme density roads concentration of SO2 decreases. Results showed that traffic could be a contributing factor towards SO2 emissions but is not the only one. Main factor in increase of SOx concentration, other than road traffic could be the burning of coal in brick kilns. Use of coal and furnace oil in brick kilns increases the levels of SOx in air. Concentration of SOx at Batti Chowk and Jinnah Park is the highest while at Sarwar Shaheed Road is lowest and also values for both locations are significantly different from other groups (Figures 5 and 6).
GIS based assessments of an environmental factor have always been an interesting study. GIS along with statistical techniques is best method to find trend in pollution level. GIS is a new and modern tool for assessment and management of environmental problems and is also useful in decision making. Jinnah Park and Rasul Pura is more concentrated with respect to air pollution because of dense settlement. The highest value of NOx found in the study area is 5 ppm while the lowest value is 0.06 ppm, on the other hand the highest value of SOx found in the study area is 1.9 ppm and the lowest value is 0.1 ppm. The most prone areas to NOx in Sheikhupura are Ghang Road, Jinnah Park, Rasul Pura and Landa Bazar. And the areas with high concentration of SOx in Sheikhupura are Ashraf General and Batti Hospital. Proper urban management can help to reduce air pollution. Concentration of NOx is more than SOx because of multiple factors towards its contribution in air like traffic (major source) and industries. Unmanaged urban development and unplanned road infrastructure is major cause of pollution in Sheikhupura. Data was collected within 3 months and there are greater chances that meteorological conditions and data collection and preservation techniques can affects the data collected.