3.1 General features of SO2 concentration and meteorology: Brief details of background meteorological conditions and daily variation of SO2 concentration in the atmosphere of SGV have represented in Fig. 2. The SO2 concentration exhibit a strong variability within the range from 1.10 μg/m3 (On Jul 15, 2018) to the highest SO2 concentration as 17.39 μg/m3 (On May 29, 2019) and also showing an increasing trend (Upslope = 0.011 ± 7.1E-4) between this period. However, after May 29, 2019, we observed a declining trend of SO2 with the downslope of -0.17 ± 0.02 until Jul 15, 2019 over SGV. The average daily SO2 concentration was 5.09 ± 1.43 μg/m3, which was lower than the nearest previous reported locations, i.e. Dehradun (12.84 ± 1.18 μg/m3) & Haridwar (25 ± 2.51 μg/m3) during the year 2009 and 2011 (Chauhan et al. 2010; Deep et al. 2019).
The temperature (Fig. 2b) shows a declining trend from Jul 15, 2018 (30.32 ℃) to Dec 28, 2018 (8.78 ℃) and increasing onwards up-to Jul 15, 2019 (28.52 ℃). The highest daily temperature was reported on Jul 1, 2019 (34.87 ℃) at SGV. Temperature also shows different correlation nature with SO2 in different seasons. SO2 shows significant (P-value = < 0.05) and negative correlation with the temperature in the M-2018 (r = -0.48 & p- value= 7.93E-4) and PoM-2018 (r = -0.23 & p- value= 0.01). Whereas, W-2019 (r = 0.43 & p- value= < 0.05) and PrM-2019 (r = 0.32 & p- value= < 0.05) have positive and significant correlation with SO2 at SGV (Table 1). In W-2018 low temperature supports the SO2 due to less oxidation and shallow boundary layer height (Gaur et al. 2014).
The relative humidity (Fig. 2c) shows several peaks throughout the observation with a vertical bell shape type pattern. In Jul-Aug- Sep, 2018, the high humidity levels were recorded on Jul 28, 2018 (93.62%), Aug 4, 2018 (84.70%), Sep 23, 2018 (93.58%) and Sep 24, 2018 (93.04%). Whereas, on Jan 22, 2019 (87.88%), 24 Jan (77.56%), 7 Feb (81.91%), 14 Feb (88.69%), 21 Feb (88.98%) and 4 Jul 2019 (84.11%) respectively. Humidity has positive and significant correlation with SO2 in the M-2018 (r = 0.26 & p- value= 0.01) & PoM-2018 (r = 0.24 & p- value= 0.10) season and negatively significant in the PrM-2019 (r = -0.33 & p- value= <0.05) during study (Table 1). (Chaudhuri and Dutta 2014) also observed the variation of pollutants (SO2, NO2, PM10 & CO and O3) over the Kolkata during year 2011 to 2012 and reported r= -0.03 and r = -0.09 with Humidity and temperature. However (Sharma et al. 2010) observed the a seasonal variation of SO2, NO, NO2, & NH3 and found good correlation with temperature & humidity in the W-2008 (r= -0.92 & 0.92), PrM-2008 (r= -0.83 & 0.86) and PoM-2008 (r= -0.63 & 0.52) over Delhi during 2008.
The wind direction and wind speed have represented in the Fig. 2d & 2e. In Jul-Aug 2018 and Jul, 2019 the winds were coming from South East direction with a lower wind speed below 5 m/s. Whereas, In the Sep-Oct-Nov-Dec of 2018, Jan-Feb -early Mar-May of 2019, the air mass was originating from the southwest direction above 5 m/s. Wind speed has significant impact on SO2 concentration during different season. Wind speed has a negative correlation (r= -0.29 & p-value = < 0.5) in the M-2018 and positive correlation (r= 0.28 & p-value = < 0.5) with SO2 during observation. However wind direction shows weak correlation with SO2 during M-2018 (r= 0.03 & p-value = 0.79), PoM-2018 (r= 0.05 & p-value = 0.75), W-2019 (r= 0.05 & p-value = 0.75), PrM-2019 (r= 0.13 & p-value = 0.21), and M-2019 (r= 0.07 & p-value = 0.62) respectively Table 1). (Gupta et al. 2004) also reported a negative correlation (r= -0.88) between SO2 and wind speed during 1997 to 2000 over Mumbai. Whereas (Turalioǧlu et al. 2005) also reported negative correlation between SO2 and wind speed (r= -0.49) along with temperature (r= -0.75) and Humidity (r= 0.028) over Erzurum, Turkey. Furthermore, the simultaneous seasonal variations of wind speed, wind direction and SO2 concentration have also explained in section 3.6.
3.2 Diurnal and monthly variations of SO2 over SGV: In the M -2018 season (Fig 3a), the maximum SO2 concentration was 3.66 ± 2.05 μg/m3 at the 1900 LT and the minimum was reported of 3.01± 1.31 μg/m3 in the morning time at 0500 LT. The SO2 variation is almost constant, i.e. changed by only 0.64 μg/m3 with an average value of 3.39 ± 1.64 μg/m3. The SO2 starts fluctuating from 0000 LT and drops at 0700 LT, after that SO2 was increasing up to 1200 LT may be due to local anthropogenic and photochemical activities as described in Fig. 3a (Husar and Patterson 1980). Furthermore, a dip was observed at 1400 LT due to change in traffic volume and SO2 is increasing up to 6.45 ± 3.49 μg/m3 at 1900 LT due to regional and long-range transport of pollutants (Gautam et al. 2018; Sandeep et al. 2020) over SGV (also explained in section 3.8). In the PoM season, The highest SO2 concentration was 5.75 ± 1.83 μg/m3 reported at 1300 LT which may be attributed to the effects of hefty traffic contribution on the national highway (NH) -58 (Fig. 3b). We also recorded another second peak of SO2 (5.54 ± 2.23 μg/m3) at 2000 LT during the same observation period that is showing the effects of local and regional biomass burning in Punjab, Haryana and transport of pollutants (also explained in section 3.7 & 3.8). In the case of W-2019, the highest SO2 concentration 6.42 ± 1.79 μg/m3 has observed at 1800 LT. Similarly, we have recorded high SO2 concentration at 1900 LT as well as second highest at 1300 LT due to the effect of long-range transport pollutant and photochemical production of SO2 during PrM-2019 season (Fig. 3b). The identical deeps in SO2 concentration were observed at 0700-0800 LT and 0500 LT during M -2018, PoM-2018, W-2019, PrM-2019 and M -2019 probably may be due to low anthropogenic activities in early morning hours at the observation site (Fig. 3a & b). In the M-2018, the SO2 concentration (3.39 ± 1.64 μg/m3) was low as compare to PoM-2018 (5.25 ± 1.74 μg/m3), W-2019 (5.71 ± 1.56 μg/m3), PrM-2019 (6.20 ± 3.09 μg/m3) and M-2019 (5.63 ± 3.09 μg/m3) season. (Gaur et al. 2014) observed high seasonal concentration during W-2009 due to shallow boundary layer height and weak oxidation of sulfate and observed a similar valley shape pattern during afternoon hours in the PrM-2009 season as well as M-2009 has almost constant variation.
The monthly variation of SO2 is represented by using box plots (Fig. 3c), to understand the standardized statics and distribution. The low SO2 concentrations were found in Jul 2018 (1.07 ± 0.82 μg/m3), and Aug 2018 (3.67 ± 0.80 μg/m3) may be due to the process of extensive wet deposition in M-2018 (Pandey et al. 2005; Rai et al. 2010; Gaur et al. 2014). However, in Jul 2019 the SO2 concentration was high (4.19 ± 0.68 μg/m3) as compared to Jul 2018 because of fewer precipitation events and bio mass burning activities (section 3.8). Whereas, the high SO2 were reported in May 2019 (7.58 ± 3.30 μg/m3), followed by Jun 2019 (6.51 ± 2.08 μg/m3), Feb 2019 (6.11 ± 0.65 μg/m3) and Nov 2018 (5.79 ± 0.71 μg/m3) during the study. The mean SO2 concentration during May and Jun 2019 was higher than median levels, i.e. SO2 concentration has clustered at a higher level due to forest fire. The SO2 levels mainly in Oct-Nov-Dec 2018 were due to extensive biomass burning activities in the Punjab Haryana and rest of the IGP region (Mittal et al. 2009; Gadi et al. 2011). Whereas, in the case of Jan and Feb 2019, a higher level of SO2 was observed due to shallow boundary layer dynamic over the monitoring site (Gaur et al. 2014; Deep et al. 2019). However, predominately high SO2 concentrations have observed in May 2019 and Jun 2019 may be attributed due to extensive fire forest activities over the Uttarakhand (Jha et al. 2016; Yarragunta et al. 2020). The Central Pollution Control Board and WHO stands have also given in Table 2 and compared with other locations (Table 3).
3.3Weekdays andWeekend analysis: Many researchers examined the variation of atmospheric pollutants on weekdays and weekends to get an idea of different pollutant sources over another part of the world’s due to change in vehicular emission (Gour et al. 2013; Filonchyk et al. 2016; Radin Mohamed et al. 2016; Özden Üzmez 2018). We have considered days from Monday to Saturday, weekend days have considered only for Sunday, and other holidays for our study for weekdays. One of the primary sources of SO2 in the SGV is due to vehicular transport activities for holy pilgrims (Badrinath, Kedarnath & Hemkund sahib) and hill stations (Auli, Chopta, etc.) every year (Semwal and Upreti 2019). The weekend effect has been defined as the decline of pollutant concentration in the atmosphere due to low anthropogenic activities (Cleveland et al. 1974). The difference between weekdays and weekends clearly shows the different levels of transportation and other anthropogenic activities over SGV. In Aug 2018, Sep 2018, Feb 2019, Apr 2019 and May 2019 the weekend effect can be easily observed with the percentage change between weekdays and weekends as 42.77 %, 38.99%, 30.09%, 26.82%, 16.30%, respectively (Fig . 4). These reductions of SO2 may be due to heavy rainfall, snowfall, summer vacations as well as change in traffic volume. However, some significant difference was also reported in Jul 2018 (-5.49%), Oct 2018 (-3.18%), Nov 2018 (-9.70%), Dec 2018 (-10.16%), Jan 2019 (-1.39%) Mar 2019 (-5.97%) and Jul 2019 (-27.01%) due to the significant vehicle moments for pilgrims/ hill stations (Chopta Auli, Badrinath & Kedarnath) along with firecracker activities in Diwali. In May 2019, the high mean than median indicated the presence of high SO2 concertation during fire forest activity (Fig. 4). The notable weekend effect of SO2 and other pollutants (O3, NO2, etc.) have already observed by many researchers due to change in traffic volume, industrial and transportation activities in the different part of the world (Riga-Karandinos et al. 2006; Gour et al. 2013; Özden Üzmez 2018).
3.4 Impact of Forest fire on SO2: The MODIS and FIRMS data have been used to understand the daily variation of fire activities/fire count (Over Uttarakhand) and ground-based SO2 (Fig. 5a). Furthermore, the significant fire days have been identified and considered as the high fire activity period (HFAP) based on the methodology suggested by Yarragunta et al. (2020). If the three-day rolling mean/ running mean exceeds the overall median values, then it has considered as HFAP. In the first HFAP (May 6 to May 14, 2019), the SO2 and fire count shows a remarkable variation in the range of 4.81 to 7.95 μg/m3 along with daily fire counts in the range of 8 to 98. The second HFAP (May 21 to Jun 4, 2019) shows a sharp increment in the fire count (4 to 150) and SO2 (5.26 to 17.39 μg/m3) with the highest values of 17.39 μg/m3 on May 29, 2019, which is very close to daily WHO standards (20 μg/m3 Daily average). However, a slightly lower SO2 concentration has observed in the third HFAP (Jun 10, 2019, to Jun 18, 2019). After, Jun 20, 2019, the SO2 and fire count deceased may be due to wet deposition activities. The SO2 concentration has reported 3.61 times higher than the lowest SO2 concentration in May 2019 as well it is only 2.63 μg/m3 less than WHO standards.
3.5 Satellite versus ground-based SO2: The satellite-based SO2 data set were extracted from the National Aeronautics and Space Administration (NASA)'s Giovanni platform by using Modern-Era Retrospective analysis for Research and Applications (MEERA) -2 model between the 77 -79 E longitude and 29 -31 N latitude above SGV (Ma et al. 2020). The satellite-based SO2 surface mass concentration (SMC) has compared with ground-based monitored SO2 concentration (C) over SGV (Fig. 6).
A moderate (r = 0.36 to 0.67) and positive correlation (r = 0.46) was also observed between the SO2 SMC and SO2 C during observation (Taylor 1990; Ma et al. 2020). The variation trend matches with the Jul, Aug, Sep, Oct and Nov 2018 along with May, Jun and Jul 2019. The average SO2 SMC & SO2 C were 3.04 ± 0.76 μg/m3 & 5.25 ± 1.84 μg/m3 reported over the observation site. The SO2 SMC was 1.68 times lower than SO2 C during the study. The satellite-based SO2 observations are very rare over Uttarakhand. (Naja et al. 2014) performed the ground-based and satellite comparison as well as model simulation of SO2 in the nearest station over Nainital, Uttarakhand.
3.6 Wind dependency of SO2: The seasonal variation of SO2, wind speed and wind direction can be visualized by using the bivariate Polar plots (Fig. 7) in the R (open-air Package) open platform (Carslaw and Ropkins 2012; Grange et al. 2016). The SO2 dominated in the southeast (SW) to northwest (NW) direction in the range of 5 to 6 μg/m3 with corresponding wind speed ranging from 4 m/s to 8 m/s in PoM-2018 season (Fig. 7b). The minimum SO2 concentration has dispersed from the northeast (NE) direction below 3.5 μg/m3. However, in the case of the W-2019 season, the SO2 was dominantly present in the South (S) to northwest (NW) direction in the range of 9.07 to 10.25 μg/m3 and low values have sprayed in the range of 3.20 to 4.37 μg/m3 at observation site (Fig. 7c).
However, in the PrM-2019 season, the SO2 concentration was dominated in the north-east (NE) to the south-east (SE) in the range of 6.5 to 9.5 μg/m3 along with the 2 to 14 m/s. Whereas, low values of SO2 have been reported in the south (S) to north (N) direction (Fig. 7d). Now M -2019 follows almost similar trends, but the concentration was high. (Fig. 7e). The wind speed less than 5 m/s indicates the local transport of pollutants over Nainital (Dumka et al. 2015) Therefore, SGV has been affected by the local and transportation of pollutants from different directions in different seasons. The high wind speed also supports the high SO2 concentration over Srinagar valley possibly due to long-range transport of air mass to SGV (The origin of possible source and air mass have also been explained in section 3.8).
3.7 Air back mass trajectory and Cluster analysis: The seasonal air back mass trajectory (AMBT) has been extracted from Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model & Global Data Assimilation System (GDAS) data set and further plotted to investigate the pathway of air mass (Stein et al. 2015). Seven days AMBTs have been calculated for local time 1200 LT at 500 m, 1000 m & 1500 m (AMSL) and plotted by using MeteoInfo software with TrajSet (version 1.5) plugin (Rico-Ramirez et al. 2007; Zachary et al. 2018). Fig. 8a suggested mass air arrival from the Bay of Bengal (BOB) and the Arabic ocean via IGP to Srinagar. However, at 1500 m, the air mass also contributed from the Sahara region in the M-2018.
In the case of the PoM-2018 and W-2018 seasons, the predominant air mass is coming from the Gulf region, Western Pakistan, Afghanistan, Eastern Rajasthan, Haryana and Punjab to the monitoring site (Fig. 8b & Fig 8c). But in the case of PrM-2019 season the transportation of air mass is from Iran, Iraq, western Pakistan, Rajasthan, Punjab via Haryana to Srinagar Garhwal in the as well from western Uttar Pradesh (Fig. 8d). A similar trend has been observed for the M -2019 season (Fig. 8e). This result influences the previous high SO2 observation during PoM-2018, W-2019 seasons (Fig. 3 & 5). The SO2 concentration at the monitoring has also been influenced by long-range transportation as well as regional emission of pollutants. Air back mass trajectory was well documented and examined (Gautam et al. 2018; Sandeep et al. 2020) also investigated the long-range and significant contribution of the pollutants from IGP to SGV.
Now to quantify the contribution of air mass we also performed cluster analysis based on the GDAS data set and TrajSet platform over the observation site (Rozwadowska et al. 2010; Sharma et al. 2016; Gautam et al. 2021). The southeasterly air mass from BOB has a remarkable contribution of 59.83% along with a contribution of 10.68 % from southeasterly air mass originating from Pakistan and the Arabian ocean to the receptor site during M-2018. However, local air mass from Punjab and Rajasthan also contributes up to 29.49% of the transportation (Fig. 9a). In the PoM-2018, the air mass has a dominant contribution of 40.43% from Central Asian and Northern African countries along with 59.56% contribution from Afghanistan side (Fig. 9b). Now, W-2019 & PrM-2019 have remarkable contributions in the range of 61.54% to 72.29% from the Central Asian, Northern African countries along with Afghanistan and Pakistan. But it should be also noticed that Punjab and Rajasthan have 20.60 % to 37.85% contribution during transportation of air mass (Fig. 9c-d). The local contribution of 31.16 % and 26.81% from BOB and Central Asia have observed during M-2019 (Fig. 9e). The air mass from the Arabian Ocean region has only a 14.79% contribution. (Gautam et al. 2021) has reported the significant contribution of air mass from the Central Asian, and neighbors countries such as Afghanistan and Pakistan through southwesterly wind over Himalayan cloud Observatory, Badshahithaul (Located 80 km from observation site).
3.8 Fire Spot over Uttarakhand and India: The MODIS and FIRMS data sets have plotted in the origin software (Serial Number GL3S4-6089-7,609,063) to represent the biomass burning activities (BBAs) over Uttarakhand (Demonstrated by green box) and India in Fig. 10 (Kharol et al. 2012; Shaik et al. 2019; Yarragunta et al. 2020). In M-2018, the BBAs are very rare over Uttarakhand and the rest part of India due to rainfall events during M-2018 (Fig. 10a). The Punjab, Haryana, Himachal Pradesh (HP) and Jammu & Kashmir (J& K) have significant BBAs and therefore, air mass helps to transport the pollutant to SGV during PoM-2018 (Previously explained in Fig. 8b & 9b of section 3.7). Dehradun, Haridwar and Udham Singh Nagar of Uttarakhand (Fig. 1) have also contributed to the local BBAs as depicted in Fig. 10b (Shaik et al. 2019). The impact of BBAs also reflects in the relatively high concentration SO2 concentration as compare to M-2018 (Fig. 3b). In W-2019, The BBAs dominantly occur in Uttarakhand due to solid wood & coal burning for warming and cooking purpose. BBAs also contribute to local SO2 emissions along with the long-range transport of pollutants (Fig. 9c & 10c). Whereas, In the PrM-2019 season (Fig. 10d), the BBAs have a predominant impact on SO2 may be due to agricultural waste burning and intense forest fire in Uttarakhand, J& K, and HP (Yarragunta et al. 2020). According to Fig. 10e, BBAs have a dominant influence as compare to the M-2018, which results in the relatively high SO2 concentration previously describe in Fig. 3a during M-2019. The BBAs in the PrM and PoM in the IGP also significantly contributes to the higher level of pollutant at high altitudes of the Himalayas (Kumar et al. 2011; Bhardwaj et al. 2018).