Data screening and quality control
The sampling protocol was considered with the proper care in collecting and preserving the rainwater samples until the complete chemical analysis. In addition, the rainwater samples were excluded that were found any traces, filthy, messy, or contaminated with dust at the sampling site. The data acquired through the chemical analysis were tested through the ion balance method. Data were examined to ensure analytical quality, considering calibration and blank dimensions. The principle of electro-neutrality depicts the eminence of chemical analysis, which can be measured to some extent by the ionic balance in particular occurrences in the samples. The cation/anion ratio in rainwater samples was observed 0.74 and 1.56, with a mean value of 1.10±0.24. There was a significant linear regression (r = 0.60) between the sum of anions and cations (Figure 1). The outcome designates that the quality of data is high and precise enough for further chemical analysis. In the current analysis, the deviation of ≤ 40.4 for the total ion concentration of 169.14 µeq/l is within the acceptable range of 113.45 – 263.33 µeq/l. However, the sample location is located at a high altitude, and such ion-balance deviations may be caused by fluctuations in lower and higher ion concentrations and the presence of weak organic acids. A prevalent organic acid in the atmosphere, coupled with the H+ ion, is likely to evaporate fast in non-preserved water samples and their quantities were lower in rainwater samples [18]. The volume-weighted mean (VWM) concentrations of the observed chemical species in the Shaune Garang catchment rainwater were calculated as follows [24]:


Ci denotes the concentration of a specific chemical species (µeq/l), while Pi and N denote the quantity of rainfall for each event (in mm) and the total number of rainfall events, respectively.
pH and EC variability in rainwater
The pH value of rainwater indicates the presence of major ionic constituents, either from atmospheric gas assimilation or anthropogenic sources. Significant ionic species in rainwater are caused by in-cloud and sub-cloud scavenging processes in the atmosphere [25]. Even though sulfuric and nitric acid is formed from nitrogen and Sulphur oxides, they freely dissolve in cloud water and produce extra hydrogen ions. Furthermore, higher levels of carbon dioxide (CO2) in the atmosphere dissolve in cloud water and produce weak carbonic acid [23]. pH (< 5.6) (i.e., atmospheric CO2 balance) indicates the presence of acidic species in the ionic composition of rainwater. In contrast, a pH value of more than 5.6 specifies the existence of basic species, particularly crustal species, or mineral dust [26]. The pH of rainwater fluctuated from 4.59 to 6.73 and the mean value with the standard deviation of 5.47 ± 0.69 during the sampling period, indicating a mixture of anthropogenic and natural chemical constituents in the catchment (figure S1a). The lowest pH value observed was 4.59 on 3 September 2017, having 7.28 mm rainfall, showing episodic high acidity of rainwater in the catchment. The high pH value was measured on 13 September 2017, when 36.25 mm rainfall was noted with the high alkaline aerosols (Ca2+) presence. Several studies at high altitudes in India and around the world report the average pH value of rainwater 5.10- 6.40. The variation in the pH and EC values of rain events during the study period over the sampling site is presented in figure S1b. In the Himalayas, at Kullu, the pH value of rainwater was reported ranging from 5.16 – 3.36 [21], while at Darjeeling, it varied from 5.0 ± 0.8 [22]. The VWM pH of the rainwater was observed as 5.56 ± 0.29, indicating the mixture of anthropogenic and natural chemical constituents in the rainwater of Shaune Garang catchment. In addition, the other fundamental constituent of rainwater is the specific conductivity, which was considered to check the quality of rainwater in the Shaune Garang catchment. The average specific conductivity observed in this valley was 13 to 36 μS/cm, with an average value of 22.31 ± 7.3. The specific conductivity in this catchment is lower than the reported value of 7 to 57 μS/cm in the Himalayan region [27 28]. To understand the correlation between specific conductivity and the sum of cations and anions in the Western Himalayan region, specific conductivity was plotted against the sum of cations and anions. The findings show a stronger correlation with anions (R2 = 0.89) than with cations (R2 = 0.60), as presented in figure 2.
For determining the rain acidity, linear regression analysis between the ionic concentration of hydrogen ion (H+) in rainwater and rainfall amount is schemed (Figure S2). In general, the analysis results suggested that the concentration of H+ increases with rainfall. The lowest pH value of 4.59 corresponds to a high hydrogen ion concentration (H+= 3.62 µeq/l) with 1.3 mm of rainfall. Furthermore, the highest pH of 6.73 (H+= 1.48 µeq/l) correlates with a little concentration of hydrogen ions (H+= 0.4) in 0.04 mm of rainfall in the Shaune Garang catchment (Figure 3). The significant rain in the middle of the sampling period was a major factor in the dilution of hydrogen ions. Seasonal variations in rainfall amount and air quality have an impact on these.
Chemical constituents of rainwater
A rigorous statistical analysis was done to determine the mean, volume-weighted mean (VWM), minimum, maximum, standard deviation, and standard error for each major ion and pH provided in Table 1. The charge balance error (CBE) analysis of the chemical composition of dissolved ions in rainwater demonstrated the dataset's accuracy. As a result, the following empirical formula was used to calculate it:

During the melting period, the error in charge balance between total cation (TZ +) and total anion (TZ -) was found to be (6%), indicating that the dataset is more accurate. Table S1 summarizes the chemical composition of rainwater samples collected from the Shaune Garang glacier during the 2017 study period. The rainwater ionic composition's volume-weighted mean (VWM) was 167.54 µeq/l, indicating low atmospheric pollution concentrations in the Shaune Garang catchment. Mean cation values were 86.13 µeq/l, and anions were 83.01 µeq/l. The outcome demonstrates the supremacy of cations over anions in the catchment. The mean values of ionic species concentrations (169.14 µeq/l) are higher than the median values (157.32 µeq/l), indicating an irregular dispersal of ionic species with skewness to the left. The median and average concentrations show an analogous configuration. Therefore, VWM is used to recognize the greater concentrations during brief periods of rain and to evade clean rain diluting and influencing the rainwater concentration [29]. Among all ionic components, calcium contributed maximum (22.35 %), followed by chloride (17.72 %), sodium (16.39 %), sulphate (13.62 %), nitrate (11.81 %), magnesium (8.97 %), bicarbonate (5.06 %) and potassium (4.08 %). Although, the acidic nature of water generally initiates due to the influence of sulfuric and nitric acid and the neutralization process by cations such as Ca2+and Mg2+ [30]. The input of cations and anions in rainwater during the study period was 48.22 and 51.78 %, correspondingly. The outcome shows that the higher contributions of anions are due to transference from distant sources instead of localized sources [22]. Calcium (Ca2+) was observed as the dominant ionic species with a much lower influence of SO42− and NO3− (13.62 and 11.81 %). Similarly, the sea salts such as Na+ are contributing 16.39 %, and Cl− contributing 17.72 % indicated the transportation from distant sources and influenced by the sea. Notably, SO42− was the second greatest abundant species, and K+ was the least abundant species among the anions during the study period. Figure S4 shows the percentage contribution of measured ionic species in the rainwater of Shaune Garang catchment. The result shows that the maximum contribution in the ionic concentration in the rainwater is of Ca2+ (22.35 %) and minimum from K+ (4.08 %). This might be due to considerable contribution from marine sources in NaCl sea salt [9] and the crustal source in the form of Ca2+and Mg2+. In addition, the contribution from anthropogenic sources SO42- (13.62 %), NO3- (11.81 %), and K+ (4.08 %) is also observed in rainwater composition, which might be due to the resident wood-burning activities [31]. Though, a significant alternative source of the perceived ionic composition in the rainwater might be dust and sea salt transported from other areas. However, the highest VWA ionic species concentration in the rainwater was Ca2+, Cl+, Na+, SO42-, and Mg2+. Volume weighted mean pH was detected 4.59 and reached as high as 6.73 with an increase in Ca2+ (56.23 μeq/l) in rainwater in the catchment. This suggests that Ca2+ be the primary neutralizing agent in rainwater. Though, Mg2+ and K+ can defuse acidity produced due to SO42- and NO3- to regulate the pH of rainwater in the alkaline range.
Ionic ratio of rainwater
Rainwater quality measures a characteristic function such as acidic and alkaline ingredients. The present research found that the ionic strength and concentration in the rainwater during the study period was calculated as 169.14 μeq/l. Although, the acidic nature of water generally initiates due to the influence of nitric and sulfuric acid and the neutralization process by cations containing Ca2++ Mg2+ [21]. Apart from this, the ionic ratio has been calculated to apprehend the comparative involvement of Sulphuric and nitric acid rain formation in the study area.
In addition, fractional acidity (FA) was calculated to understand the acid neutralization capacity of rainwater in the catchment. The correlation between acidic and neutralizing species was evaluated by the subsequent equation [32].

The average ratio of H+/ (NO3- + SO42-) measured 0.07 indicates that 93% of the acidity was neutralized in the rainwater during the study period. According to the study conducted in the Kothi, North-Western Himalaya reported that almost 96 % of rainwater was neutralized [21]. The result indicates that rainwater is less acidic in Western Himalayas, particularly in the Shaune Garang catchment than in North-Western Himalaya. Also, a study in Pune and Delhi reported that the average H+/ (NO3- + SO42-) ratio was 0.02 and 0.08, demonstrating that 98% and 92% of the acidity in rainwater was defused by alkaline species [33]. The average ratio of (NO3- + Cl-) / (SO42-) was measured as 2.38, indicating the higher value compared to North-Western Himalaya, mainly due to the insignificant amount of nitric and hydrochloric acid in rainfall. However, in the case of North-Western Himalaya, the ratio is slightly lower because of sulfuric acid [21]. The equal ratio of NO3- / SO42- was measured 0.97 ± 0.60, which suggests the contribution to the acidity in rainwater is dominated by nitric acid and Sulphuric acid. The equal ratio of (Ca2+ + NH4+ / NO3- + SO42-) is mainly used to evaluate the extent of human activity in water chemistry. All the ionic ratios of the observed ionic concentration discussed in this section in the rainwater from Shaune Garang catchment are summarized in Table S2. The plot (Figure 4a) between Ca2+ + NH4+ against NO3- + SO42- shows the positive correlation of all data set throughout the study period and linear spread beyond the 1:1 equiline with a ratio ranging from 0.48 to 2.73 with the average equivalent value of 1.31 ± 0.51. The result indicates that NH4+ and Ca2+ ions are key factors in neutralizing acidity in rainwater through CaSO4 and (NH4)2SO4. The scatter plot between Ca2+ and HCO3− has been shown in Figure 4b. It is imperative from the figures that there is a large variation between Ca2+ and HCO3−. In rainwater, the ionic concentration of Ca2+ is much greater than HCO3− in rainwater during the study. The ratio of NH4+/NO3− extended from 0.42 to 1.00 with the mean value of 0.81 ± 0.21 and the ratio of NH4+/ SO42− ranged from 0.21 to 1.46 with the mean value of 0.73 ± 0.38 signifies the dominancy of compounds NH4NO3 over (NH4)2SO4 in the atmosphere [34]. Furthermore, the result from the ionic composition illustrates that the ammonium nitrate is leading over ammonium sulphate composites in the Shaune Garang catchment during the study period [35].
Neutralization Potential
The factor of acidic and neutralization potential is a significant sign of understanding rainwater's chemical behavior. The difference between rainwater's acid and neutralization potential describes the responsible factors to the entire mechanism in the process. Acid potential (AP) is the summation of Nitrate (NO3− and non-sea salt (NSS) Sulphate (SO42−) concentration and neutralization potential (NP) is the sum of Ca2+, Mg2+, and K+ in the rainwater. The ratio of Acid potential (AP) and neutralization potential (NP) clarify the dominancy of such factors within the system. The ratio of AP/NP measured as (0.74) less than 1 during the study period in the rainwater, indicating that neutralization potential dominates the acidic potential in the Shaune Garang catchment. In addition, the atmospheric neutralization potential of the chemical constituents in rainwater samples can be estimated through the Neutralizing Factors (NF) for the specific parameters. The neutralization factor (NF) measures how well acidic components are neutralized by crustal elements and ammonium ions. Alkaline particles serve a crucial function in modulating the acidity of RW during wet deposition [36]. The two main neutralizing sulfuric and nitric acid agents are calcium and ammonium. The primary supply of calcium is soil dust, whereas the primary source of ammonium is combustion procedures. By computing neutralization factors, the potential of Ca2+ and Mg2+ in neutralization has been established [37]. Neutralization factors can be used to evaluate the ability of key alkaline ions in rainfall to be neutralized (NF). The following equation is used to compute the neutralization factors for ions.

Where: [X] is the concentration of the major ions (Ca2+, Na+, K+, Mg2+, and NH4+) expressed in μeq/L.
The Neutralization factor, which was dominant cations like Ca2+, Na+, K+, Mg2+, and NH4+ in the rainwater of Shaune Garang catchment, was measured as 0.85, 0.63, 0.15 0.34, and 0.39 respectively. The results suggest that the calcium ion (Ca2+) in rainwater is a major controlling neutralizing agent, and potassium is minimum. The particulate matter in the rainwater, which is high in Ca2+ carbonates or bicarbonates, buffers the acidity of cloud-water, which is widespread in India [38-39]. These findings reveal that Ca2+ and Na2+ ions, together with NH4+ ions are the primary neutralization components with the minimal role of K+ in rainfall of Shaune Garang catchment. Table S3 shows similar observations on rainwater neutralizing capacity from several parts of India. The neutralization factor from the previous study demonstrated that the maximum neutralizing capability in rainwater was calcium (Ca2+) and ammonium ions (NH4+) from different parts of India.
Contribution of Sea-salt and non-sea-salt
An effort was made to study the contribution of various ions in rainwater from sea salt (SS) and non-sea salt (NSS). The input of sea salt and non-sea salt to essential ions in rainwater was assessed by relating the Cl−/Na+ ratio in rainwater to saltwater (Table S4). On the other side, the NSS input was calculated by subtracting SS from the total measured ion. The measured ratio of Cl−/Na+ (1.14) was less than the observed seawater ratio (1.16), showing that sources of sea salt considerably impact the Cl in rainwater at the Shaune Garang area. Furthermore, increased K+/Na+, Mg2+/Na+, Ca2+/Na+, and SO42−/Na+ ratios suggest the potential input of additional sources such as soil and nonmarine. The high SO42− values compared to the seawater ratio in the catchment indicate a significant anthropogenic influence. Based on the relative ratios of the major ions, the probable compound formations include NaCl, CaSO4, MgSO4, MgCl2, HNO3, NH4SO4, and (NH4)2SO4 [17 23]. The predominant ionic species in rainwater are thought to indicate the proportional effect of natural and anthropogenic sources. Natural sources include alkaline and insoluble components of mineral aerosols from the earth's crust and seas. The influences from volcanic sources are insignificant at the sample site. The enrichment factor (EF) was calculated using Na+ as the reference element.

where xi is the required major ion and Xi/Na (seawater) is the seawater ratio.
The EF unity infers no enrichment and, as a result, no input from any source other than seawater. EF values greater than one, on the other hand, indicate enrichment of a specific ion relative to non-sea salt sources. It was discovered that the EF of all major ions (Mg2+, K+, Ca2+, and SO42-) was more than one, specifying an important input from sources other than sea salt, i.e., soil and anthropogenic sources in the catchment. The result shows that the only enrichment value of Cl- is less than 1, indicating that around 97.2 % and 1.73 % come from marine and nonmarine sources, respectively. The enrichment factor of Shaune Garang, compared with the other previous studies in India, indicates that the EF value of Kothi, Himachal Pradesh (Table S5) is very similar to the present study. Therefore, it is verified that the predominance of nonmarine contributions influences the Shaune Garang catchment. However, our observation compared with the previous research from Nainital, Uttarakhand indicates that the chemical configuration of rainwater in the Himalayan region is affected by both natural and anthropogenic sources.
Air mass trajectory analysis
Particulate matter concentrations have risen rapidly in developing countries such as India because of high emissions from a variety of human-caused events [40-41]. The researchers studied the influence of aerosols on the earth's radiation budget, weather, rainfall, and cloud formation [42-43]. Aerosol emissions are to blame for rising regional temperatures and have been identified as a key factor in the melting of the Hindu Kush-Himalayan-Tibetan glaciers [44]. The Himalayan region has also been affected by heating caused by another amalgamation of solar energy due to aerosol brown cloud [45]. It can similarly decrease the albedo of snow and glaciers, leading to enhanced melting because of the accumulation of radiation-absorbing aerosols on them. Previous study [46] has estimated that the annual average melting of Himalaya reached from 0.7-0.85 m w.e. per annum at Lahaul/Spiti glaciers during 1999-2004. Moreover, the backward air trajectory modeling significantly showed the impacts up to the Mt. Everest region due to high aerosol emissions from the North-west region. The air mass back trajectories are critical for determining the potential transportation paths of air mass to the station. To understand the other pollution sources, air mass back trajectories were computed for the study area using the NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model [47]. GDAS (Global Data Assimilation System) data was used to generate three-dimensional back trajectories [48] for the monitoring location in Shaune Garang glacier on precipitation days. During July, 80% of air parcels were coming from the northwest direction at relatively higher levels which may carry various air pollutants to the study locations. The speed of the wind and atmospheric boundary layer’s height play an essential part in the diurnal difference of pollutant concentrations and dispersion. The measured ion concentration was approximately two and a half times higher than the air masses passing from the southwest. The bulk of air masses are arriving from the west (the Mediterranean Sea or the mid-west Atlantic Ocean), also known as western disturbances, according to the research. These air masses pass over the Persian Gulf, Iran, Afghanistan, and Pakistan, bringing torrential rainfall to the western Himalayas. Along the dust track, long-distance dust movement interacts with anthropogenic emissions, raising local particulate matter concentrations [49]. The back trajectory analysis (Figure 5) clearly shows that in July, 60% of the air parcels were coming from the Arabian Sea and Bay of Bengal, with the remainder coming from the northwest direction of India. The Arabian Sea disturbance resulted in a significant reduction in ion mass concentrations.
Chemometric analysis
The correlation coefficient matrix was made among the measured chemical components in rainwater to identify the relationship between all chemical species (Table S6). The correlation analysis specifies chemical parameters such as SO42−– NO3–, Mg2+– NO3–, Na+ – SO42−, Ca2+– Cl−, SO42− – Cl− and NO3– –HCO3– have a reasonable correlation. The correlation coefficients between Cl− and SO42− (0.50) indicate that it is influenced by crustal soil. In addition, a minimal correlation (0.22) between bicarbonate and soil-derived Ca2+ was observed among the ions, indicating that the contribution of these sources was influenced by soil dust containing large fractions of CaCO3 originating from the same sources. The correlation among Ca2+ and NO3− was observed to be 0.44 during the study period, suggesting that soil is the most important source of nitrate in rainwater, which is present in the atmosphere in the form of Ca (NO3)2 [21]. Minor contaminants like sulfate (SO42−) and nitrate (NO3−) are significantly correlated (0.67), showing their origin from the same sources. The controlling contribution factor of SO42− and NO3− in acid formation was nearly 54.64 % and 45.38 % respectively in the catchment. Therefore, the findings demonstrate that the high concentrations in the Himalayan area are thought to be caused by airborne primary and secondary particle movement. Apart from this, a good correlation (0.56) observed between Na+ and Cl− suggests that most of Na+ and Cl− components originated from marine sources and were transported with air masses. However, the ionic ratio of Cl−/Na+ in the rainwater compared to the seawater ratio in the previous section (Table S4) indicates a part of chloride from other emission sources.
To identify potential sources of ionic species in rainwater, factor and principal component were used. The correlation matrix was exposed to Bartlett's Sphericity test, which yielded 2 (critical) = 84.93, which was more significant than the critical value 2 (critical = 30.61; df- 45 and p-value 0.0001 with significance level 0.05). Following varimax rotation, factor extraction with an Eigenvalue greater than 1 was used to perform principal component analysis (PCA). The Pearson correlation coefficient between the measured chemical components in rainwater supported the factor components. In principal component analysis, the data were exposed to a Varimax rotation, which optimizes the variance to generate a pattern of loadings on each factor that is as different as possible, allowing for easier interpretation. Factor loadings represent the correlations of each variable with the factor. In this study, four factors were identified, and each variable was assigned a loading with each factor. When the loading of each variable was determined to be greater than 0.50, the significant loading was considered. Table 1 summarizes the loadings with the variance extracted by the factors (Eigenvalue).
Table 1
Principal components (PCs) loading for selected physic-chemical parameters in rainwater.
1
|
F1
|
F2
|
F3
|
F4
|
pH
|
-0.672
|
0.306
|
0.148
|
-0.069
|
EC
|
0.276
|
0.648
|
-0.296
|
-0.458
|
Na+
|
0.299
|
0.763
|
0.101
|
-0.078
|
K+
|
-0.072
|
-0.545
|
0.699
|
0.079
|
Ca2+
|
0.771
|
-0.230
|
0.277
|
-0.072
|
Mg2+
|
-0.541
|
0.452
|
-0.059
|
0.623
|
Cl-
|
0.596
|
0.362
|
0.389
|
0.427
|
SO42-
|
0.423
|
0.646
|
0.480
|
-0.121
|
NO3-
|
-0.685
|
0.534
|
0.250
|
0.007
|
HCO3-
|
-0.643
|
-0.091
|
0.480
|
-0.470
|
Eigenvalue
|
2.924
|
2.492
|
1.363
|
1.037
|
Variability (%)
|
27.894
|
24.981
|
14.640
|
10.857
|
Cumulative (%)
|
27.894
|
52.875
|
67.515
|
78.373
|
The four rotated principal components (PC) were formed during the analysis, including Eigenvalues, percentage of variation, and cumulative percent variability. All the principal components (PC) with its Eigenvalues are considered more than to assess the governing factor in the rainwater. In a scree plot (Figure 6) in the PCA, the cumulative variability was measured as 78.3%, including four-component (PC1 explained 27.89%, PC2 explained 24.98 % PC3 explained 14.64 % PC4 explained 10.85 %). PCA generated four significant factors clarifying 78.37 % of the total data variance. However, Factor 1 explains 27.89 % of the data variance and displays a strong optimistic loading for Ca2+ and Cl-. High loadings of Ca2+ NO3- and moderate loading of SO42- signify that the most important sources of these ions are the burning of fossil fuel and soil dust.
In addition, the first two PC loadings were drawn to understand and interpret the grouping and association amongst the variables. Results clarify that (Factor I and II) explain 52.87 % of the total variance (Figure S3a and S3b). High loadings of a chemical constituent such as Cl–, Ca2+, Mg2+, NO3– and HCO3− demonstrate the natural sources, sea, and soil. Factor II comprises NO3–, K+, and SO42−, indicating that the sources of these chemical constituents from vehicular emissions and biomass burning. Although Na+ and Cl- generally come in the form of sea salt but give the impression of the same factor. Even if SO42− is not originated from the sources of soil or the sea, it seems that it might be due to the marine wind containing SO42− and Ca2+ [19]. Apart from this, as most of the Mg2+ originates from the sea source thus, it appears in this factor. The SO42− over the study region is primarily due to inadequate fuel burning, fertilizer uses, thermal power plants, refineries, and long-range transport.
Comparison of ion concentrations in rainwater at higher altitude regions
The ionic concentration of the rainwater of Shaune Garang catchment was compared with the ionic concentration of rainwater from several study areas around the world, such as China, Nepal, Tibet, and India (Figure 7). However, it is imperative to emphasize that the most relevant parameters which accelerate the formation of acid rain are NO3- and SO42-. In Shaune Garang catchment, high intensities of SO42- (23.83 μeq/l) and NO3- (18.75 μeq/l) were observed during the study period. Therefore, these two parameters are the primary inorganic ions formed from SO42-and NO3- gases during the precipitation event, emitted from different sources such as vehicular emission at high altitudes [50]. The concentration of NO3- in the sampling area (18.75 μeq/l) surpassed the concentrations observed in the study area of Kathmandu (12.75 μeq/l), Dhunche (10.70 μeq/l), Dimsa (8.52 μeq/l), Gosainkunda (4.40 μeq/l), Southern Everest region (0.01 μeq/l), Northern Everest region (1.10 μeq/l), Nam co Tibet (10.30 μeq/l), Lhasa (2.0 μeq/l), Southern Tibet plateau (2.33 μeq/l), Nainital, Uttrakhand (11.9 μeq/l). Therefore, with the known fact that NO3- is a pioneer gas for forming acid rain, it produces nitric acid with soil interaction, causing high nitrification with the metal’s mobilization. Consequently, SO42- originates primarily from vehicle traffic and industrial effluence. In this study, the concentration of Na+ (27.66 μeq/l) and Cl- (31.28μeq/l) were observed low compared to the study area Nainital, Uttrakhand, where Na+ was observed as 49.8 μeq/l and Cl- was 67.3 μeq/l. The ratio (Cl- /Na+) was calculated for the study area and observed the ratio of Shaune Garang catchment was 1.13, indicating the origin in marine sources. According to the European Union, Air Pollution Study Group recommends that the Na+/Cl- ratio between 0.5 and 1.5 indicate a marine source. In the Shaune Garang catchment, values are significantly closer to the range, and it is possible that the Na+ also might have a continental origin. The concentration of HCO3- (8.09 μeq/l) in the Shaune Garang catchment has been lower than the reported values of other locations. This result indicates that the HCO3- ion is also significantly present in rainwater samples of Shaune Garang catchment, and it is largely imitative from soil resuspension and limestone use around the study area [51-52]. Even though the concentration of HCO3- in Lhasa is very high as compared to another study, this might be due to high limestone exploitation at a particular study area. Apart from these parameters, Mg2+ (13.98 μeq/l) was observed high compared to another study, however, K+ (6.56 μeq/l) and Ca2+ (37.13 μeq/l) were almost similar in all studies. These parameters, such as Mg2+, Ca2+, and K+ originate primarily from seawater, soil, and forest fires [53] and heavy agricultural and intense agricultural vehicular traffic [54].