Spatial Interpolation of the Concentrations of Particulate Matter and Carbon Dioxide of Some Selected Tourist Sites in Srinagar City, Jammu and Kashmir, India Current World Environment

This investigation was carried out to monitor and prepare thematic maps of the monthly spatial variability of the ambient concentrations of particulate matter (PM 1 , PM 2.5 , PM 4 , PM 10 , and TSP) and carbon dioxide (CO 2 ) of some selected tourist sites (viz. Shalimar Garden, Chesmashahi Botanical Garden, Harwan Garden, Nishat Garden, Naseem Bagh, Lal and Sher-e-Kashmir Agricultural and of Kashmir (SKUAST-K) Shalimar campus of Srinagar city, Kashmir valley from November 2019 to April 2020 using the Inverse Distance Weighting (IDW) interpolation technique in Quantum Geographical Information System (QGIS). Considering Srinagar city as a growing city in terms of population, construction, vehicles, etc., the Jammu and Kashmir State Pollution Control Board in a report on managing air quality in Srinagar city identified the pollution sources in the city with their estimated source proportion as follows: vehicular %), and demolition emissions %), minor industrial activities (7-8 %) and other sources (3 %). Thus, Srinagar city’s air quality is deteriorated by these sources. Therefore,this


Introduction
Particulate matter pollution and increase in carbon dioxide concentration are seen as the major troubles man is facing today with respect to health and the environment. In 2016, the health effect institute (IHE) noted that the chronic exposure to ambient PM 2.5 killed about 4.1 million people and caused the disability-adjusted life years' (DALYs') loss to 106 million; the highest total deaths were obtained in India (25%) and China (26%). 1 These data were obtained by the extensive air quality monitoring network stations in urban areas of developed countries and the satellite observation of air quality combined with obtained information from models of global chemical transport to estimate the global PM 2.5 concentration and its consequent effects on humans. 1 PM 2.5 originates from fossil fuel combustion in power plants, congested traffic flow, heating from residential and industrial areas using oil, coal, or wood. 2 PM 10 has adetrimental effect on human health by its blocking and destructive effect on the nasal and bronchial passages, igniting different respiratory-related effects that end up in sickness or death. 3 PM 1 consists of anthropogenically derived particles having dangerous impacts on the human respiratory and cardiovascular systems by directly entering into the blood system. 4 Very small particles (<1μm) pass into the lungs and then the blood stream via the blood barriers and cause some severe health complications. 4 One of the major concerns related to air in today's scenario is the increasing concentration of carbon dioxide and its role in the environment. Atmospheric CO 2 has been increasing by ̴ 3 ppm/ year since 2001 and can be attributed to the increase in the global use of fossil fuels, particularly in China and India. 5 The situation of Kashmir is no different from that of other parts of India and the world. Even Srinagar city is exhibiting increased concentrations of particulate matter and CO 2 due to industrialization and increased vehicular population. 6 The ambient concentrations of particulate matter and carbon dioxide caused by vehicular pollution was monitored from 2019 to 2020 in Srinagar city and the findings showed that during the whole period, average PM 1 concentrations ranged from 15.10 -108.9 μg/m 3 , PM 2.5 (28.70 -577.50 μg/m 3 ), PM 4 (44.50 -780.87 μg/m 3 ), PM 10 (57. .53 μg/m 3 ), TSP (77.77-1410.27 μg/m 3 ) and CO 2 (332.4 -655.0 ppm). 7 The Inverse Distance Weighting (IDW) technique of spatial interpolation that is available in geospatial tools like the QGIS, ArcGIS, etc. estimates cell values by weighting values of geometric data in the neighborhood of each processed cell. In the weighting process, the points situated nearer to the cell centre have more influence or weight on the processed cell. This technique assumes that as distance increases from the sampling site, the influence of the entered variables on the maps decreases with respect to the sampling site. 8 The inverse distance weighting has been proven as the finest technique for the interpolation of air quality assessment data in urban delicate zones. 9 It has also been proven to produce a better comparison between interpolated and measured values of suspended particulate matter, SO 2 , and NO 2 than Kriging during an appraisal of different techniques of interpolation for the parameters of the ambient air quality of Port Blair. 10 The IDW was a very necessary tool to figure out the relations between health effects and air pollution in the assessment of the daily trend of PM 2.5 concentration for the contiguousness of the United States from 2009, at county and census block group level respectively. 11 It was also used to measure and interpolate the concentration of the aerosoloptical thickness (AOT) at 40 locations in the areas of urban Ranchi and to quantify the conditions of the atmosphere at 42 locations in the coalfield zones of Southern Karanpura. 12 The primary use of interpolation techniques has been to map out bed rocks, soils, air quality assessment, and ground and surface water studies. 13,14 Since the inverse distance weighting have such important usefulness in interpolating and mapping air quality parameters, and since no such work has been published for Srinagar city nor given much attention to the in-depth/systematic study of its air quality and the responsible factors for its deterioration, 15 this paper presents six months study IDW maps of the spatial variation and interpolation of particulate matter and carbon dioxide for an idea of the pollutant blanket over the city using some tourists' locations as sampling/monitoring sites.

Materials and Methods Study Area
Srinagar city located at N 34.08 o , N 74.80 o in the Kashmir valley is internationally known for its remarkable gardens and lakes. The stunning looks of the gardens in Srinagar have attractive hill slides, flowering shrubs and trees, and enthralling water bodies. Notwithstanding the witness of centuries of change, the gardens still attract tourists from all over the world. The study area shown in Figure 1  It weighs about 0.79kg and can be used for up to 24 hours for intermittent sampling. This instrument has a sensitivity that ranges from a high of 0.3μm to a low of 0.5μm and thus monitors the levels of PM 1 , PM 2 . 5 , PM 4 , PM 10 , and TSP simultaneously. 16 The provision of simple and efficient operation is due to the many functions of the rotary dial. It has the capacity of differentiating seven ranges of particles by counting and sizing them and then converting the counted data to mass measurement (µg/m 3 ) by using a proprietary algorithm. It computes for each particle detected a volume then allocates for the transformation a standard density value which is improved by the setting of a K-Factor to make a better measurement accuracy of ±10% to calibration aerosol. 16 A Comet software which is a program that extracts information (alarms, data, settings, etc.) from Met One Instruments' products 17 modifies these K-Factors which, with a reference unit can be analytically obtained. Or a recommended K-Factor setting of 3.0 is used in the case of an unavailable reference unit.This instrument uses the operating principle of particle count to mass conversion using scattered laser light. 18 As the particles enter via the detection chamber which contains a photodetector and samples one particle at a time, the laser light is scattered and detected by the photodetector. The instrument can analyse the scattered light's intensity and deduce the particle's size. In addition, the number of received lights on the photodetector can be used by the instrument to account for the number of particles in the detection chamber.This approach is advantageous because it can use a single detector/analyser to simultaneously sense different diameter particles.
Whereas, the CDM 901-CO 2 monitor works with the principle of NDIR (Non-Dispersive Infrared Radiation),and the NDIR sensors functions on the principle of the absorption of infrared radiation (IR) by a particular gas.The sensor becomes sensitive to a certain gas because of the usage of different bands of absorption in the infrared spectra. At the wavelength of 4.3μm, carbon dioxide strongly absorbs IR radiation, thereby avoiding the bands 2.5 to 2.9 and 5.2 to 7.5μm at which water vapour is absorbed. 19 The path length of the IR that is between the detector and the source orders the gas level the sensor can detect and the Beer's Law, I = I o exp (-αℓ)describes the relationship. From the above formula, the transmitted light via the gas cell is I and the incident light on the gas cell is I o , the sample's absorption coefficient (cm-1 units)is given as α and the optical pathlength of the cell is given as ℓ. 20 As pathlength increases, the signal that is received and the infrared radiation's intensity exponentially decreases with it. The infrared detector produces a higher signal as path length decreases. But, there's a decrease in sensitivity, since a shorter optical path reduces the gas' absorption distance for the radiation. 21 The range of CDM 901 is from 0-2000 ppm. The resolution stands at 1% of the full scale i.e., 20 ppm. The sampling is done by diffusion technique without the need for a pump. The response time for the reading is less than 15 seconds. The inbuilt LCD shows the humidity and temperature range of the ambient air along with the CO2 concentrations. 21

Preparation of IDW Maps of Particulate Matter and CO 2 Concentrations
The inverse distance weighting (IDW) obtainable in QGIS software version 3.16.3 was used to map out the concentration of particulate matter (PM) and carbon dioxide (CO 2 ) by the interpolation of data generated at the seven (7) monitoring sites.
To predict the concentration of pollutants for locations in which sampling was not done, the IDW interpolates the data obtained for the sampled location assuming that areas closer to each other are similar than those farther apart. 22 On each map, a pollutant was interpolated on monthly basis for the six months of monitoring to show the pollutant's spatial variability in the nearby areas of the sampling locations where sampling was not undertaken. Unknown values/concentrations of geographic parameters such as elevation, chemical concentrations, and noise levels of places that are hard, costly, or even impossible to access have been predicted by interpolation techniques with few data for sampled areas/points. 23,14 Flexibility for carrying out interpolations that are efficient and optimal for the spatial variation of any given number of samples is provided by the IDW interpolation technique 4,24,14 The prediction of the concentration of the pollutants in the nearby locations was done by observing the ten (10) partitions by which the maps were scaled.

Statistical Analysis
On the fortnight of each month, three replications for each sampling/monitoring site were taken and the average triplicate data obtained for each month and the average six months of monitoring pollutants load at the seven locations were analysed using the R software version 4.0.4. A one-way analysis of variance (ANOVA) of the randomized complete block design (RCBD)is an experimental design for comparing treatments in blocks with homogenous experimental units was used for the analysis in the R software to get the critical difference (CD) of significance at p ≤ 0.05 to determine the monthly and average six months variation (significant or insignificant) of the pollutants among the different monitoring locations respectively. Also, the average data of the monitored pollutants for the six months of sampling were subjected to a Pearson correlation matrix analysis to know their pattern in a relationship (positive or negative and significant or insignificant).

Results and Discussion
The data in Table 1 shows the monthly mean variation of each monitored parameter at the study sites that were interpolated on monthly basis using the Inverse Distance Weighting (IDW) technique in the QGIS software version 3.16.3. The Inverse Distance Weighting (IDW) maps (Figure 2 -7)show the interpolation of the variation of the concentration of the monitored pollutants at the sampling sites to the unsampled locations in each month to predict concentration values of the monitored pollutants for the unsampled locations. These maps significantly show that the concentrations of pollutants decrease with distance from their respective sampling sites. Thus, the assumption is that each monitored sampling site has an influence locally that diminishes with distance.
On each map, the thick red domains represent the areas with the highest concentration of the pollutants which fades out with distance showing a decrease in the concentration of pollutants with an increase in distance from the study site. The thick blue domains characterize the zones with the lowest concentration of pollutants. The legend on each figure in a descending order i.e., from dense blue to dense red, shows an increase in the concentration of pollutants; and in ascending order i.e., from thick red to thick blue, shows a decrease.
The influence of pollutants on the different nearby locations can be predicted based on how far the scaling area covers the coordinates on the map. The major objective of the investigation is to come up with simple but inclusive maps that show and predict an estimate of the variation and blanket of the monitored pollutants respectively over the city to get an immediate requirement of estimation of pollutants at different sites in Srinagar due to the homogeneity in topography. Areas that are very close to busy roads with tall buildings which create an enclosure that hinders the dispersion of roadside emissions tend to be worse in the concentration of larger diameter particles like PM 10 and TSP. 26 Also, since these larger particles are also often formed by the condensation, coagulation, nucleation, fog evaporation, and cloud droplets, in which dissolution and reaction of gases also occur 27 their concentration will certainly be higher in areas like Naseem Bagh and Lal Chowk. Comparing our data to the data on PM 10 concentration recorded in the revised action plan for the management of air quality in Srinagar city from 2014-15 to 2017-18 shows that November to January i.e., during the winter, are the months with the highest air pollution levels. 6 The monthly spatial variation of CO 2 concentrations shown in Table 1 (f) and Fig. 7 in each month was found to be statistically significant (p ≤ 0.05) between most of the study sites. This might be due to the variation of vehicular movements at the different times of the day at and around the various locations, the burning of biomass, plants' use of CO 2 for photosynthesis, and the weather effect on CO 2 concentration. CO 2 was higher at Shalimar Garden in most of the months (664.27, 662.53, and 612.60 ppm in December 2019, January and March 2020 respectively) due to its location nearby a roadside that experiences high traffic especially in the morning and evening and surrounded by human habitation that do lots of biomass burning for heating purposes especially in the winter months (Nov. 2019 to Jan. 2020). Vehicular pollution and agricultural residue burning are regarded as the main source of emission of CO 2 in north India including Kashmir. 28 Since the study was carried out in winter (Nov. 2019 to Jan. 2020) and spring (Feb. to Apr. 2020) it is certain that both seasons especially the winter had some influence on the increasing concentration of particulate matter and CO 2 . The winter months recorded the highest concentrations of particulate matter and CO 2 at most of the monitoring sites. This is because, during winter, temperatures are low, relative humidity is high, and therefore, the stagnation and less dispersion of particulate matter. 29 Also, vehicles burn more fuel in winter due to cold temperatures and therefore take much time for engines to reach their maximum operating temperature and so enhances the atmospheric CO 2 . Among all the months, particulate matter and CO 2 were lower in April 2020 during the spring and COVID-19 lockdown in all Kashmir. This was due to less/no vehicular movement during the lockdown and the less amount of biomass burning for heating purposes.   Critical Difference (CD) is significant at p ≤ 0.05; SE(m) is the standard error of the mean

Mor-After-Eve-Mor-After-Eve
The average daytime (morning, afternoon, and evening) concentrations of the monitored pollutants for the monitoring sites are depicted in Table 2.
From the data, it can be observed that for each monitored pollutant at most of the monitoring sites and the average mean daytime concentrations of pollutants, the morning and evening concentrations were highest. For example, the average mean daytime concentration of PM 1 was highest in the morning (67.07 µg/m 3 ) followed by the evening (57.62 µg/m 3 ) and then the afternoon. The same trend was followed by PM 2.5 and PM 4 . But, with the larger particles (PM 10 and TSP) and CO 2 , the average mean daytime concentration was highest in the evening followed by the morning and then the afternoon. This can be attributed to a lot of vehicular movements in the morning as people drive to their different working places and in the evening as they drive from the same. Also, traffic congestion in Srinagar city is observed mainly in the morning and the evening times of the day. During the winter months (November -January), the high biomass burning for heating purposes that are observed around the various monitoring sites plus the lesser atmospheric wind to disperse particles and CO 2 , and the high level of humidity keeps particles particularly the larger ones stagnated in the ambient air for long period especially in the morning and evening. It has been proven by many research data that pollutants transport from a long-range, lower atmospheric mixing, and less dispersion of pollutants, especially during the winter are responsible for higher pollution as pollutants are always stagnant in the lower atmosphere. 15,32,33 30 The average six months spatial variation of the concentration of the monitored pollutants was statistically significant (p ≤ 0.05) between most of the monitoring sites.
The correlation of particulate matter and CO 2 given in Table 4 shows a positive correlation even though non-significant between the monthly concentration of most pollutants. Nevertheless, PM 1 , PM 2.5 , and PM 4 show a significant (p ≤ 0.01) strong positive correlation with each other. Also, PM 4 and PM 10 , and likewise PM 10 and TSP had a significant (p ≤ 0.05) strong positive correlation with each other. The positive correlation shows how the increase of any of these parameters will cause a corresponding increase of the other. Particulate matter concentration has been found to have a mutual correlation with outdoor/ambient environment. 31

Conclusion
From the data and IDW maps obtained, it was clear that there was a significant monthly spatial variation of the concentration of all the monitored pollutants between most of the monitoring sites and their surrounding unsampled areas. Also, the spatial variation of the average six months concentration of monitored pollutants was statistically significant (p ≤ 0.05) between most of the monitoring sites. This might be attributed to the location of the monitoring sites, changes in traffic flow, different seasonal (winter and spring) activities like biomass burning, and weather conditions like atmospheric temperature, wind speed, precipitation, and relative humidity. This informs us that the concentration of particulate matter and carbon dioxide varies on monthly basis with distance from one location to another in Srinagar city. Also, the data proves that the ambient air in Srinagar city experiences increasing pollution in the morning and evening than in the afternoon. And that all the parameters/pollutants correlate positively with each other, that is, they all increase with each other simultaneously even though most are insignificantly positively correlated. The sources of pollution with their estimated source proportion in the city have been noted by the Jammu and Kashmir State Pollution Control Board in a report on managing air quality in Srinagar city as follows: vehicular emission (65-75 %), dust from bad roads (10-15 %), biomass and garbage burning (10-20 %), construction and demolition emissions (5-8 %), minor industrial activities (7-8 %) and other sources (3 %). 6 Thus causing the deteriorating condition of the city's ambient air.It is therefore concluded that the poor air quality of Srinagar city varies with distance as depicted by the data and IDW maps with respect to the monitoring sites and the monitored pollutants. Thus, giving an idea of the pollutants blanket over the city. In the future, a more expanded study with regards to this work is recommended to be carried out covering the whole Srinagar district.