Statistical Analysis of Total Column Ozone over Uttarakhand: Environment of Himalaya

Total Column Ozone (TCO) is a critical factor affecting the earth’s atmosphere, especially in the Himalayan region. A comprehensive study of TCO trend analysis and corresponding consequences in the Himalayan atmosphere needs to be analyzed. We statistically examine TCO variability by analyzing the daily TCO dataset of the last 15 years (2005–2019) over the crucial region of the Himalayan environment i.e. Uttara-khand, India. Obtained results indicate that TCO values are at peak during the spring season whereas it shows the least value during the winter season. The highest and lowest value of Coefficient of Relative Variance (CRV) is estimated as 3.14 and 1.09 during winter and monsoon season, respectively. Air mass trajectories have been estimated using Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT), which shows the existence of strong seasonal variability of Ozone corresponding to continental and maritime transportation towards Uttarakhand. Moreover, Least Square Method (LSM) and the Mann-Kendall test estimate a high correlation (86%) for the seasonal and annual trend of TCO with a negative rate. The obtained decreasing rate is very low which indicates recovery of TCO during the study period. Further results imply that the inter-annual oscillation pattern of TCO is similar to Quasi-Biennial Oscillation (QBO) significantly. In addition, a comparative study has been performed for the data measured by two TCO measuring instruments i.e. Ozone Monitoring Instrument (OMI) and Ozone Mapping Profiler Suite (OMPS). TCO values measured from both instruments are highly correlated (96%) with an average relative difference of around 3%. The outcomes of this study are expected to be beneficial for future study of TCO over other crucial regions of Himalayan territory.


INTRODUCTION
Ozone is a robust component of the earth's atmosphere and can be found primarily in two layers of earth atmosphere.A major part i.e. 90% of its mass is found in the stratosphere lying within 16-35 kilometers above the surface of the earth whereas the remaining 10% ozone is found in the troposphere of earth.Stratospheric ozone is considered as shielding layer (good ozone) as it hinders the high energy radiations (2000 0 A to 3000 0 A) from penetrating the earth's main environment.On the other hand, Ultra Violet (UV)-B radiations are biologically harmful solar radiation and the ozone layer prevents these hazardous radiations to reach on the surface of the earth (Dobson et al., 1973;Dobson et al., 1968;Madronich et al., 1998).Most ozone is formed in the tropic region where it moves from lower to higher latitude by the Brewer-Dobson circulation (Weber et al., 2011).TCO is the amount of total ozone (troposphere and stratospheric ozone) containing in a vertical column from the earth surface to the top of the atmosphere at standard temperature and pressure (STP) (David W et al., 2014;Brasseur et al., 1988Brasseur et al., 1997).It is measured using data obtained from ground-based stations and satellites in terms of Dobson unit (DU), which describes the thickness of a layer of pure ozone at STP.One DU is 2.69x10 16 ozone molecules/cm 2 in a vertical column extended till the atmospheric limit or 0.01mm thick layer of pure ozone to the surface of the earth at standard condition (STP).The value of TCO significantly varies from 200 to 300 DU all over the globe.Meanwhile, if its value reduces to 220 DU it is considered as a condition for occurring Ozone hole at any location.A series of reactions of ozone formation and distortion was proposed by Sydney Chapman in 1930.It depicts that the photolysis of Oxygen molecules produces nascent radical in presence of UV radiation which further reacts with oxygen and results in ozone molecule.Moreover, ozone depletion involves the reverse process of conversion of ozone molecules into atomic oxygen due to catalyst reaction [Groves et al., 1980].Apart from this ozone abundance are the key parameters for estimating the formation and depletion of ozone at any location of the globe.The processes of formation and deterioration of ozone are found to be naturogenic as well as anthropogenic.However, both can occur simultaneously (Rubin et al., 2001).The anthropogenic chemicals such as chlorofluorocarbons (CFCs), halocarbons, and other volatile organic compounds can destroy stratospheric ozone significantly by reacting with ozone.(Solomon, S. et al., 1999;Cicerone et al., 1974;Crutzen et al., 1974).Apparently, few studies suggest long term variability of ozone with geophysical parameters, including the Quasi-Biennial Oscillation, solar activity, the stratosphere-troposphere exchange, and volcano eruption (Willett, H.C. et al., 1962;Garcia, R. et al., 1987;Ningombam S S, et al., 2011;Hema Bisht et al., 2014).Despite a small portion of ozone contributes to earth' atmospheric composition, TCO plays a vital role in climate change and ecology (Wayne et al., 1987;Wayne et al., 2000;Ogunjobi et al., 2007).The depletion of the ozone layer results in severe serious health issues i.e. skin cancer, cataract, deficient immunity, respiratory disorders for species (UNEP 1998;UNEP 2003;UNEP 2016).It can also contaminate the marine ecosystem as well (cell-damaging life of phytoplankton and zooplankton).Besides, ozone depletion permits UV-B radiation which may retard the physiological and developmental processes of plants (Schmalwieser et al., 2003;Sivasakthivel et al., 2011).Therefore, ozone layer conservation is mandatory to sustain an ideal ecological surrounding on earth (WMO 1989;WMO 2011, J.I Freijer et al., 2002;R. Venkanna et al., 2015;).
Numerous studies on atmospheric ozone have been carried out all across the world.The Antarctic ozone hole was first observed by the British Antarctic Survey in the 1980s (Farman et al., 1985).Several studies emphasize seasonal depletion in the stratospheric ozone layer in Antarctica and have observed depletion during the spring season (Iwasaka et al., 1987;Rowland et al., 1989;Stolarski et al., 1990).Over the last two decades, minute depletion in ozone was also observed over the Arctic during the late winter and early spring (Brune et al., 1990Goutail, F., et al., 1999).But these studies are limited to the polar region only (Newman et al., 2006).Such studies have been extended for other locations of the globe and established statistical trends of TCO over the northern and southern mid-latitudes (Rowland et al., 1991;Chakrabarty et al., 1998;By Michael L. et al., 2007;P. J. Nair et al., 2013;Ayodeji Oluleye et al., 2013;Kok Chooi Tan et al., 2014;A. Badawy et al., 2017).These studies suggest that the depletion of ozone is not only confined in the Polar Regions but vary significantly in another latitudinal region as well.The study of ozone layer depletion and variation has emerged as a major global scientific and environmental issue in the recent past.In this framework, various international protocols have been designed (UNEP 1987).Montreal Protocol is an international treaty to protect the ozone layer by phasing out the production and limited use of o zone-depleting substances (ODS) (WMO 1989).This agreement was executed on 1Jan 1989 and has been signed by most of the countries including India.This protocol has reduced the risk of further ozone depletion to a large extent (Zubov et al., 2001;Anderson et al., 2000;Guus J. M. Velders et al.,2007;Pinedo Vega et al., 2017).From the perspective of India, a few studies on ozone variation and depletion have been reported to date (Mani et al., 1973;Kundu et al., 1993;Pulikesi et al., 2006;Vazhathottathil Madhu et al., 2014).These studies are limited to industrial areas with high population density.However, there is a paucity of such studies in Himalayan regions i.e. areas with lower population density and enhanced ecological versatility.The ozone climatology is supposed to be crucial in these regions due to climatic disparity.Therefore, a comprehensive study of ozone (TCO variation) over the Himalayan region is needed to monitor the atmospheric ozone concentration in a regular period.In this paper, we have statistically analyzed the spatial and temporal tend of TCO variability using the TCO dataset of the last 15 years i.e. 2005-2019 over the crucial region of the Himalayan environment i.e.Uttarakhand, India.Air mass trajectories have been estimated using Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) to understand the source of air pollutants and corresponding continental and maritime transportation towards Uttarakhand under various climatic conditions has been investigated.We have compared the data meas ured by two TCO measuring instruments i.e.Ozone Monitoring Instrument (OMI) and Ozone Mapping Profiler Suite (OMPS) and corresponding results have been analyzed.

Study Area
Uttarakhand, a state in northern India crossed by the Himalaya extended from (28.71 0 N to 31.45 0 N and 77.56 0 E to 81.03 0 E) is the study area shown in Fig. 1.It is situated in southern Asia (North of Indo-Gangetic Plain (IGP), East of China, and foothill of Himalayas).It has a total geographic area of 53,483 km 2 of which 86% is mountainous and 65% is covered by forest with an average population density of 189 people/km 2 (Census 2011).It is also surrounded by major Indian industrialized cities i.e.Delhi, Lucknow, and Chandigarh.Along with more industrial development the agricultural activities, burning of crop residues, and forest fire are the primary source of pollutants over Uttarakhand.The mountainous topography and complex land-sea interactions are responsible for phenomenal changes in weather over Uttarakhand.The average temperature varies within a range of 21.8 0 C to 5 0 C with the highest in June while the lowest in December and for few stations the temperature falls below 0 0 C. Maximum rainfall occurs from July to September i.e. early NEM with average precipitation 160-200 cm.Uttarakhand holds a rich heritage of India and is known as Devbhumi i.e. "Land Of God".Its snow-clad peaks, beautiful hill stations, and charming weather make it a favorite destination for tourism.Nanda Devi, Valley of flowers and National Park (1988,2005) are world heritage sites declared by UNESCO are the most fascinating places of Uttarakhand.Uttarakhand is the most popular Himalayan state of India, situated in the natural environment of Himalaya and it plays a vital role in hosting many animals, plants, and rare herbs.Most importantly, Uttarakhand covers a major portion of the Himalayan range and therefore represents a crucial study station (area) to analyze atmospheric phenomena i.e. ozone variation over the Himalayan region.

Data Sources
The daily data of TCO has been acquired from the open resource database of NASA and NOAA.The used data is recorded with two TCO measuring instruments i.e.Ozone Monitoring Instrument (OMI) and Ozone Mapping and Profiler Suite (OMPS).The TCO value is measured in Dobson Unit (DU).The details about the ozone measurement are summarized in Table 1.
The monthly and annual mean for all the stations are derived by taking the cumulative mean of daily data.The principle of measurement of TCO is based on the reflection of solar radiation on certain spectrum bands (Heath et al., 1975;P.K.Bhartia et al., 2013), some unresolved data corresponding to polar night and low light area cannot be measured and recorded as zero (Pinedo-Vega et al., 2017).These unresolved data are few and removed in to account for monthly and annual mean.To examine the seasonal variability of ozone we divide TCO data into five seasons' i.e. winter (Dec-Feb), spring (March-April), summer (May-June), monsoon (July-Sep), and autumn (Oct-Nov).Also, to investigate the variation in ozone profile i.e. the latitudinal variation of TCO, we studied four stations of India i.e.Leh Laddak (34.10N-77.4E),Dehradun (28.37N-77.13E),Banglore (12.58N-77.34E),and Kanyakumari (8.4N-77.32E)from northward to southward for the year 2014.The coefficient of Relative Variance (CRV), correlation coefficient, probability distribution, least square method i.e. regression analysis method, and Mann-Kendall test (MK) are used for statistical analysis and trend analysis of TCO variation over Uttarakhand during a period of 15 years i.e. from 2005 to 2019.

METHODS
a) Percentage relative difference is calculated as; (1) Here 'X' is the value of ozone concentration and N is the number of years.b) A coefficient of relative variation (CRV) for each site has been determined as; (2) Here 'SD' is the standard deviation and 'μ' is the mean for N years.c) Trend analysis and fluctuation of ozone concentration ha ve been estimated by Linear Regression (Least Square method) and Mann-Kendall test for all the selective sites.d) Statistical models are used for the probability distribution (Normal, Lognormal, and Gamma) of TCO.We check their deviation from normality using KST (Kolmogorov-Smirnov D-test) probability distribution over Uttarakhand.

Back Trajectory Calculation
The origin of air quality and transportation patterns of air mass arriving towards Uttarakhand has been analyzed using Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) backward trajectories.The HYSPLIT Model (Draxler and Rolph, 2011) developed by Air Resource Laboratory (NOAA) and these trajectories have been simulated from the Web site (http://ready.arl.noaa.gov/HYSPLIT_traj.php).Based on the highest concentrations of ozone in the study period, multiple air mass back trajectories for 14th May 2014, 13th June 2014, and 22nd Dec 2014 are stimulated at different elevations: 500 and 1,000 m a.g.l.using HYSPLIT models.The analysis has been performed with the GDAS meteorological dataset at a starting time of 00:00 UTC with a total run time of 120 hours.

RESULTS & DISCUSSION
The daily data of TCO (DU) from January 2012 until December 2019 has been acquired from OMI and OMPS.The monthly mean of TCO from both the datasets has been compared by estimating relative differences based on Equation (1) have been tabulated in Table 3.A range of relative difference 0.51% to 5.29% from 2012 to 2019 has been found which is within the acceptable interval.Moreover, the average relative difference of monthly TCO data via both the datasets is merely 3.5%.Also, we found a high positive correlation between the values of TCO retrieved from both the instruments having a different resolution for the period of 2012-2019.

Monthly Variation
The monthly mean of TCO derived from the cumulative mean of daily data of TCO is found in the range from (297 13.78) DU to (260 8.73) DU for OMI and (300 16.08) DU to (258 6.07) DU from OMPS datasets.The average maximum and minimum value of ozone concentration are measured in May and December respectively from both the instruments over Uttarakhand.Further, a bell-shaped normal distributed curve for monthly ozone mean of 15 years (2005-2019) has been observed as shown in Figure 4. We have observed similar results as the average TCO value from OMI increased, the ozone value measured from OMPS was also increased.Therefore, the quality of OMPS datasets is analogous to OMI data of ozone.
We have statistically analyzed the data and estimated final distributions comply by the mean monthly series of TCO by applying Kolmogorov-Smirnov goodness of fit test of distribution.It is noticed that monthly TCO distribution exhibits lognormal probability distribution in Uttarakhand.The basic statistics (maximum, minimum, average, standard deviation, and distributions) for the Interannual monthly average TCO (2005TCO ( -2019) ) from both data sets are mentioned in Tables 3 and 4. On observing table 3 and 4, we have found that the set of monthly series of TCO follows different distributions and changes from normal, lognormal to gamma distribution which indicates that the distribution is progressively approaching heavier tailed from January to December.Interannual variability of TCO retrieved from the OMI dataset (2005 -2019) showed merely the same oscillating behavior for all the stations along the latitudinal and longitudinal belt of Uttarakhand illustrated in Table 4 depicts the trend of TCO annual and seasonal (retrieved from OMI from 2005 to 2019) for all the stations by the least square method.We obtained a significant negative trend of TCO with rate 0.03 DU/year, 0.162 DU/year, -0.045DU/year, -0.057DU/year, -0.001 DU/year, and 0.033DU/year for Dehradun, Almora, Haridwar, Pithoragarh, U.S.Nagar and Chamoli respectively.However, the annual decreasing rate is very low over Uttarakhand which indicates a significant recovery of TCO during the study period.Also, we observed a negative trend for winters and spring seasons for most of the stations while the increasing trend in summers, monsoon, and autumn seasons for all the stations.Similarly, results obtained from the Mann-Kendall test are illustrated in table 5.The trend by LSM and MK test is highly correlated in all the season as well as annually except monsoon.The correlation coefficient between trend obtained from both the methods has been calculated 0.86 annually, 0.72 in winters, 0.97 in spring, 0.53 in summer, 0.24 in monsoon, and 0.75 in autumn seasons respectively.Moreover, for most of the year, the value of TCO has been observed to fluctuate about its mean value excluding in 2008 and 2015.

Spatial Variation (Coefficient of relative variation (CRV))
The variability of ozone concentration over Uttarakhand has been studied.Annual and seasonal CRV of TCO for six stations i.e.Dehradun, Almora, Haridwar, Pithoragarh, U S Nagar, and Chamoli is shown in Table 5.A significant variation in TCO has been observed among these study areas.Further, it is clear from the observations that the highest value of CRV is obtained in winter while the lowest at monsoon for all the selective stations.However, a gradual ascending trend has been observed for the annual variation of CRV from southern stations to northern stations.These outcomes suggest that annual TCO varies as a function of latitude.On observing these trends, it is a noticeable fact that CRV of ozone concentration decreases gradually from north to south in winter, whereas it increases from west to east in summer and autumn.Also, we found that in monsoon season it decreases from west to east.To have more clear insight, the spatial distribution of CRV is shown in fig. 5.

Latitudinal Variation
To characterize the ozone profile i.e. latitudinal distribution of TCO, we have estimated the monthly mean of ozone concentration for four Indian subcontinent Leh Laddak (34.10N-77.4E),Dehradun (28.37N-77.13E),Banglore (12.58N-77.34E),and Kanyakumari (8.4N-77.32E)for the year 2014.Fig. 6 indicates the legging of two months in maximum ozone value for Leh (high latitudinal area) and Kanyakumari (low latitudinal area).We have statistically analyzed the data and estimated skewness for distribution for northward region i.e.Leh Laddak and Dehradun the value of skewness is found positive (0.515 and 0.183) respectively.Similarly, the skewness is negative i.e. -0.539 and -.600 for Banglore and Kanyakumari respectively.However, a quantitative variation in TCO approx 16% has been calculated from a latitudinal belt of India i.e. from Kanyakumari to Leh Ladakh.The probable reason for this variation might be the processes of formation and transportation ozone from the tropic to the pole region.(Chakrabarty et al., 1979;Nandita D. Ganguly et al., 2005).Also, CRV has been estimated at 6. 58, 4.8, 4.48, and 4.17 for Leh Laddak, Dehradun, Banglore, and Kanyakumari respectively.

Back Trajectory Analysis
Back trajectories are simulated for correlating it with seasonal change of TCO and to understand the air mass transport and source of origin (Stohl, A., 1998) Since marine air mass is relatively clean compared to continental air mass.this oceanic influence, the mass increases the composition of hydroxyl radicals (OH).As ozone is a highly reactive and unstable gas so it reacts with OH radicals and converted into atomic oxygen or molecular oxygen.This causes depletion of O3 to a larger extent and maybe one of the reasons for the reduction of ozone concentration observed during the summer and winter season at this location.This complement our result illustrated in figure 4 i.e. maximum TCO value in spring then further decreases till winters.It concludes that maximum ozone concentrations in different seasons were a consequence of the transport effect of these trajectories.

CONCLUSION
We have statistically analyzed spatial and temporal variability of TCO values over Uttarakhand (Environment of Himalaya) for 15 years i.e. 2005-2019.The daily data of TCO was retrieved from OMI and OMPS datasets having special coverage of 0.25 0.25 and 1 1 respectively.A significant seasonal and annual cycle has been established with the minimum value of TCO in December and maximum in May from both the datasets using daily TCO variability.The statistical outcomes for TCO from both the datasets are summarized.During the study period, we have obtained a monthly mean of TCO data retrieved from OMI was within (297 13.78) DU to (260 8.73) DU from OMI and (300 16.08) DU to (258 6.07) DU from OMPS respectively.Kolmogorov-Smirnov goodness of fit test was applied to find final distributions comply by the mean monthly series of TCO from both the instruments.We found that most of the months follow lognormal probability distribution at Uttarakhand for both the datasets of TCO.Interannual variability of TCO has been performed for OMI datasets for 15 years and it showed merely the similar oscillating behavior for all six stations along the latitudinal and longitudinal belt of Uttarakhand.This fluctuation of TCO from 2005-2019 was an identical quasi-biennial oscillation (QBO) pattern.Further, we have observed an abrupt change in ozone variation in the year 2007 and 2016.This abnormal change might be due to interruption in quasi-biennial oscillation for the respective years.
Seasonal and annual TCO trend has been analyzed using the least square method and Mann-Kendall test for the seasonal and annual values of TCO from OMI datasets for 15 years over six selected sites of Uttarakhand.Corresponding results yield negative trends for annual, winters, and spring season with the maximum in winter while the positive trend for the rest of the seasons.Similar annual and seasonal trends have been estimated from statistical tests i.e.LSM and MK test.These results were highly correlated (0.86%) with each other.Although a negative annual trend has been observed in Haridwar, Pithoragarh, and US Nagar with a decreasing rate of -0.045DU/year, -0.057 DU/year, and -0.001 DU/year respectively but corresponding values are low which suggests a recovered TCO level in Uttarakhand from 2005 to 2019.
Furthermore, annual and seasonal COV pattern has been estimated in Uttarakhand.Annual COV of TCO showed a significant variation as a function of latitude which was increased from low to high latitude stations.The higher value of COV of ozone concentration occurred in winter whereas the lowest values occur at monsoon for all selected sites.Besides, seasonal variation was investigated through backward trajectories in association with continental and maritime during the study period.It showed that maximum ozone concentrations in different seasons are a consequence of the transportation of air mass trajectories.These results are footprints to correlate the TCO variation and climatic factors such as temperature, relative humidity, rainfall, and zonal wind, etc.Furthermore, we observed a positive skewed asymmetric curve of TCO for the higher latitude station of India Leh Laddak (34.10N-77.4E),Dehradun (28.37N-77.13E)and a negatively skewed curve for lower latitude station Banglore (12.58N-77.34E)and Kanyakumari (8.4N-77.32E).This was statistically confirming an earlier peak of maximum ozone concentration from northward to southward stations.
Apart from this, a comparative study for two data sets i.e.OMI and OMPS has been performed in terms of relative difference made between average monthly TCO data retrieved from both the datasets during the period of 8 years (January 2012 to December 2019) at Uttarakhand.The estimated average relative difference was merely 3% and in the range of 0.51% to 5.29%.Indeed, the daily data TCO from 2012-2019 was highly correlated for every year in the range of 68% to 99%.The outcomes of TCO distribution and trends were found similar to both the datasets.These results suggest that although both data sets used in this study i.e.OMI and OMPS have different resolutions but their measurements complement each other.Highly correlated outcomes from both these datasets enhance the validity of the results obtained in this study as well.
In summary, this study provides a detailed statistical interpretation and estimation of TCO variation over the Himalayan region i.e.Uttarakhand, India.Ozone layer depletion, mass trajectories, etc. have been analyzed and discussed using the data set of 15 years (2005-2019) retrieved from OMI and OMPS.The estimated results suggest a crucial variation in TCO has been observed in Uttarakhand over the study period.It is expected that the obtained outcomes would be helpful for the development of further studies of ozone variation (TCO variation) over other geographical regions as well.

Figure 1 :
Figure 1: Map of Uttarakhand with latitude and longitude of the study area.

Figure 2 :
Figure 2: Monthly variation of TCO over Uttarakhand via OMI and OMPS datasets.

Figure 4 .
From 2005 to 2019 ozone variability repeats its ascending and descending pattern at a period of two years.This fluctuation of TCO follows the analogous pattern of the quasi-biennial oscillation (QBO) significantly.Further, we observed an abrupt change in TCO in 2007 and 2016.It suggests this fluctuation might be due to interruption in QBO for the respective years.

Figure 5 :
Figure 5: The Interannual variation of TCO across (a) longitudinal and (b) latitudinal belt of Uttarakhand

Figure 7 :
Figure 7: Back trajectories arriving in Uttarakhand (a) summer season (b) rainy season (c) winter season

Table 1 :
Locations and description of the selected station of Uttarakhand

Table 2 :
Characteristics of the total ozone content datasets.

Table 3 :
We have not observed an average monthly mean of TCO value less than 220 DU (Condition for occurring Ozone hole at any location on the globe).It implies that However, we noted TCO value less then 220DU for one or two days in December 2008 and 2016 and January 2019 from the OMI data set while in December 2008 and January 2019 from the OMPS data set, which can be ignored for accounting it as the ozone hole.For the total number of years considered average ozone maximum and minimum value have observed 395 DU and 212 DU from OMI and 382DU and 208 DU from OMOPS respectively.Comparison of TCO values retrieved from OMI and OMPS datasets.

Table 5 : Probability distribution with statistical parameters of TCO via OMPS datasets (2005-2019) for Uttarakhand
the same years 2015 and 2008 respectively.The behavior is not repeated in any of the seasons (spring, summer, monsoon, and autumn) except winters i.e. the maximum and minimum TCO values from two sensors do not overlap in the same year again.
autumn (264.53 6.87) DU from OMI dataset whereas OMPS shows that the maximum average TOC in the spring season was (297.73 12.16) DU and the minimum was (259.48 6.58) autumn DU.In addition, the sensors recorded annual maximum and minimum TCO measurements in

Table 6 :
Annual and seasonal trend of TCO by least square method from 2005-2019

Table 7 :
Annual and seasonal trend by Mann-Kendall test from 2005-2019

Table 8 : Coefficient of the relative variance of TCO from 2005-2019.
. Further observing the influence of continental and maritime air mass on seasonal variation of surface ozone, backward air trajectories are retrieved from the surface of Uttarakhand fig 7(a), 7(b), and 7(c) during seasons i.e. spring, summer, and winter.These back trajectories imply that the air mass moves from west of Uttarakhand summer and winter seasons.During the spring season (14 May 2014), trajectories are continental and short-range.Air masses at low altitudes (500m, 1000m) appear to be originated from a western direction crossing Afghanistan and Pakistan following through Indian state Punjab and Haryana.Fig. (7a) and influenced by pollutants from nearby industrialized areas.In summer (13 June 2014), trajectories of air mass are marine type and occurred at heights (500m, 1000m) bring ozone and ozone precursors from directions i.e.Sea and crossing some landmass area before reaching Uttarakhand.It is an indication of gasoline combustion originated from the Middle East and transport by the air mass trajectories at different altitudes fig (7b).Trajectories followed by air mass in the last week of June (end of summer) indicated the onset of monsoon from the Arabian Sea.Similarly, during the winter season (22 Dec 2014), trajectories are long-range and marine type transports from, originating from the Caspian Sea (northwestern direction of recep tor site) and crossing Afghanistan Pakistan, and Indian states Punjab and also contributes to the observed law O3 levels due to oceanic influence fig (7c).