Rainfall variation
Rainfall is the quantity of rain falling to a particular area within a given time. Rainfall varies with the time and area. It is an important factor for the determination of particulates and dissolved pollutants in the atmosphere in a particular area (Migliavacca et al., 2005; Liyandeniya et al., 2020a). Rainfall was calculated using Equation (1) in three sites in mm.
where V is the rainwater collected (mm3), A is the surface area of funnel (mm2), and D is the diameter of the funnel (mm).
Out of a total of 40 samples collected from the University of Peradeniya (UOP) site within the sampling period including both North-east monsoon (NEM) and South-west monsoon (SWM) seasons, five samples were dry only deposition, while the remaining samples were wet deposition (Figure 6). Further, among seven samples collected from the Kandy sampling site, three samples were due to dry deposition while the remaining samples were wet only deposition. In the Polgolla site, two samples were dry only deposition, while five samples were wet deposition. The mean weekly rainfall of UOP, Kandy and Polgolla sampling sites were 40.12 mm, 16.52 mm and 33.15 mm, respectively. The volume weighted average values of water quality parameters and anions concentrations are summarized in Table 1.
Table 1: VWA values of water quality parameters, major anion concentration of three sampling sites.
|
Parameters
|
UOP
|
Kandy
|
Polgolla
|
Water quality parameters/ VWA
|
Mean rainfall/ mm
|
40.12
|
16.52
|
33.15
|
pH
|
7.34
|
7.36
|
7.13
|
Conductivity / µS cm-1
|
40.45
|
53.59
|
37.22
|
Salinity / ppt
|
0.22
|
0.33
|
0.22
|
TDS/ ppm
|
20.05
|
19.64
|
18.22
|
Total hardness / ppm
|
9.31
|
17.74
|
10.73
|
Anion concentration VWA / µeq L-1
|
Cl-
|
70.62
|
82.51
|
141.24
|
NO3-
|
3.47
|
<MDL
|
22.22
|
SO42-
|
34.63
|
48.65
|
38.14
|
3.2 pH Variation
The pH parameter represents the acidity or basicity of aqueous solution (Liljestrand, 1985) . The UOP site recorded an average pH value of 7.34, and further, out of 40 samples, 13 samples showed pH < 7.00, which is usual, due to the dissolution of atmospheric CO2 (Figure 7). The sample collecting period was divided into three categories namely, after fully lockdown (AFLD), normal situation (NS), and partial lockdown (PLD). After the country became normal and within the initial part of PLD, the pH of precipitation was about 6.25 which could be attributed to traffic activities and re-starting regular chemical analysis in laboratories in the proximity to the sampling site (Figure 7). The pH being greater than 7.00 in many samples is indicative of atmospheric contamination due to basic pollutants. Also, Ca is the most abundant element which is generated as soil-borne element in the Kandy city area (Gunawardana et al., 2012). Another indication of atmospheric pollution is the growth of algae in samples having pH close to 9.00.
The average pH values of Kandy and Polgolla sites were 7.36 and 7.13, respectively. Due to the activation of SWM season, the trans-boundary air pollution affects the air quality (Abeyratne & Ileperuma, 2006). Hence, the pH values of depositions at all sites varied between 6.50 and 7.50 during this period (Figure 8). Observation of high pH at the Kandy site can be attributed to high traffic activities with incomplete fossil fuel combustion; such activities are less in the other two sites (Gunawardana et al., 2012; Yatigammana et al., 2022)
Conductivity variation
Dissolved salts and other inorganic substances contribute to increased conductivity in water although organic substances do not much contribute to (Kramer et al., 1996). The average values of conductivity of the UOP, Kandy and Polgolla sampling sites were 40.45 µS cm-1, 53.59 µS cm-1 and 37.22 µS cm-1, respectively. The maximum average conductivity reported at the Kandy site can be attributed to high traffic and construction activities, producing high content of pollutants, some of which would form ions in rainwater thereby increasing the conductivity in the area (Gunawardana et al., 2012; Ileperuma, 2020; Yatigammana et al., 2022) . As UOP and Polgolla sites are in sub-urban areas with low human activities, conductivity of deposition is expected to be low (Liyandeniya et al., 2020a; Madhushani et al., 2023). Nevertheless, prevailing wind due to the large reservoir in the vicinity of the Polgolla site promoting particulate matter to travel through atmosphere complicates the situation having varying conductivities (Figure 9) (Dharmapriya et al., 2023). Moreover, as sample collection periods of all the sites included SWM and NWM seasons, the conductivity of deposition samples were decreased due to the dilution effect because of rainfall (Chathuranga et al., 2020).
Total dissolved solids (TDS) variation
Total dissolved solid is a measure of the dissolved content of all inorganic and organic substances in a liquid (Kramer et al., 1996). The total dissolved solid, usually expressed in ppm, is directly linked with conductivity and salinity ((Rusydi, 2018)). The average values of TDS in UOP, Kandy and Polgolla sites were recorded as 20.05 ppm, 19.64 ppm and 18.22 ppm, respectively. At both Kandy and Polgolla sites, the maximum and minimum TDS values were reported in dry or significantly low wet deposition, and wet deposition periods, respectively. As observed in conductivity measurements, dilution of pollutants in the atmosphere due to SWM activated period leads to decreased TDS values. Regarding the UOP site, the highest TDS of rainwater sample (Figure 10) was recorded as 74.00 ppm during the NEM activation period. During this period, the country had become normal without any lockdown, after the first wave of COVID-19 pandemic, and the number of vehicles traveling in the Colombo-Kandy (A1) Road, which is closer to the UOP site, was normal. Then, the minimum TDS value was reported to be 3.62 ppm in the SWM activation period. During this period, the second wave of COVID–19 pandemic affected the country, and the Western Province was locked down to mitigate the spread of the disease. Consequently, the number of vehicles traveled in the A1 Road became significantly less. This fact supports decreased TDS values of deposition samples.
Salinity variation
Salinity is the total concentration of all dissolved salts, and it mainly contributes to the conductivity of a solution (Lucas Rego Barros Rebello et al., 2020). The average values of salinity in Kandy, Polgolla and UOP sites were recorded as 0.0361 ppt, 0.0310 ppt and 0.0299 ppt, respectively. The highest salinity of 0.0740 ppt was recorded in Polgolla sampling site (Figure 11). The prevailing wind due to the reservoir would be the main cause for varying the salinity at the Polgolla site as compared to the other two sites, and hence, no significant variation in the salinity with rainfall was observed at the UOP and Kandy sites (Dharmapriya et al., 2023).
Hardness variation
Hardness in water is due to the presence of dissolved salts of calcium and magnesium and is expressed as ppm CaCO3. The total hardness occurs due to the presence of bicarbonate, chloride and sulphates of calcium and magnesium ions. The degree of hardness standard as established by the World Health Organization (WHO) is shown in Table 2.
Table 2: Degree of hardness standard established by WHO.
Degree of hardness
|
Range (ppm)
|
Soft
|
< 60
|
Moderately hard
|
60 - 120
|
Hard
|
120 - 180
|
Very hard
|
>180
|
The average values of hardness at UOP, Kandy and Pollgolla sites were determined to be 9.31 ppm, 17.74 ppm and 10.73 ppm, respectively. Hardness values of depositions of all samples at the UOP and Kandy sites were less than 60 ppm, indicating that depositions were soft as per WHO standards (Figure 12). Regarding the Polgolla site, one sample, out of seven, was recorded to have hardness of greater than 60 ppm, while the other samples had hardness < 60 ppm. The maximum hardness value of 83.2 ppm was recorded on the 30th of September 2020 at the Polgolla site. During this period, the SWM was activated, and the anthropogenic activities were normal after the lockdown period of COVID-19 pandemic. Also, limestone quarries located in Digana area of about 13.4 km away from the sampling site, construction activities, and Gneiss quarries located in Kandy area emit Ca2+ and Mg2+ into the atmosphere, contributing to hardness of depositions (Gunawardena et al., 2012; Samaradiwakara et al., 2021; Dharmapriya et al., 2023) . When compared to the volume weighted average hardness values at UOP, Kandy and Polgolla sites reported as 17.32 ppm, 19.69 ppm, and 14.96 ppm, respectively, in a previous publication (Madhushani et al., 2023) it is clear that the hardness values are significantly less in the duration of the present study due to less traffic activities as a result of lockdown periods. However, in Kandy and Polgolla sites hardness values variation is significantly lower, due to the continuous function of limestone quarries located in Digana area after the first wave of COVID 19 pandemic. The area close to Kandy Railway Station has a high level of particulate matter with a high concentration of Ca2+ (Samaradiwakara et al., 2021). The situation in the country after Week 17 became completely normal and the transport activities became regular. Therefore, the sample from the Kandy site in week 17 showed high values for many parameters.
Anion variation
Chemical analysis (Cl-, SO42-, and NO3-) performed for obtaining volume weighted mean anion concentrations on samples of the three sampling sites (UOP site: from 1st week of July 2020 to 4th week of March 2021; Kandy and Polgolla sampling sites: from 3rd week of August 2020 to 5th week of September 2020) indicates that Cl- was the most prominent anion, while NO3- showed the lowest concentration at all three sampling sites (Figures 13 - 16). More specifically, the volume weighted mean (VWM) anion concentration of all the sites follows the order, Cl- > SO42- > NO3-. As samples of UOP included the SWM and NEM, and those of the other two sites included SWM, it is argued that Cl- ions would have come from sea salt through monsoon rains (Ileperuma, 2020; Liyandeniya et al. 2020a). Another fact is that the UOP site, being in the vicinity of chemical laboratories of the Faculty of Science, which emit SO2, NO2, and H2S to the atmosphere, contributes to changing the natural composition of the atmosphere although this site experienced less traffic activities (Dharmapriya et al., 2023; Madhushani et al, 2023). Also, construction activities in the vicinity would also have contributed to enhance the pollution levels in depositions.
Although Kandy site experiences more traffic congestion, more urbanization and less vegetation cover as compared to the Polgolla site, Cl- and SO42- levels, in general, are higher at the Kandy site (Abeyratne et al. 2006; Gunawardena et al.2012). The presence of NO3- at the Polgolla site can be attributed to the Ceylon Tobacco Company located near site. This industry would act as a point source of NO3- (Dharmapriya et al., 2023). Moreover, the digestion of organic waste accumulated in the reservoir in the vicinity of the site would result in emission of SO2, NO2 and H2S (Dharmapriya et al., 2023; Madhushani et al, 2023). Overall, the levels of nitrate anions were much less as compared to the observations reported by Liyandeniya et.al 2020a and Madushani et.al 2023 due to the pandemic.
Statistical analysis
Pearson correlation analysis
The Pearson correlation (r) can be categorized by considering a linear relationship, as identified in Equation 3, between two variables (Alastuey et al., 1999).
where N is the number of pairs of scores, ∑xy is sum of the products of paired scores, ∑x is the sum of x scores, ∑y is the sum of y scores, ∑x2 is the sum of squad x scores and ∑y2 is the sum of squared y scores. The correlation coefficient of greater than zero represents a positive correlation between two variables, while correlation coefficient less than zero represents the negative correlation between the two variables. Then, correlation coefficient near zero is very weak linear relationship and the value of exactly -1 and +1 represent perfect linear relationship. Further, the strength of the relationship is interpreted based on the magnitude of the correlation as follows: 0.00 – 0.19 very weak; 0.20 – 0.39 weak; 0.40 – 0.59 moderate; 0.60 – 0.79 strong; and
0.80 – 1.00 very strong (Chathuranga et al. 2020). Correlation coefficient is associated with a probability (p) value, and the p < 0.05 concludes that there is a significant association between the respective pairs of variables.
Pearson correlation of water quality parameters
Regarding the UOP site, the highest value of 0.906 was reported between salinity and TDS, a very strong positive correlation. The other strong positive corrections observed are conductivity - TDS (0.788), conductivity - salinity (0.756), conductivity - hardness (0.625), and hardness – TDS (0.631), while salinity - hardness showed a moderate positive correlation (0.557). Other pairs of water quality parameters showed weak or very weak correlations at the UoP site. The Pearson correlation coefficients of water quality parameters in bulk deposition at the UOP site are represented in Table 3.
Table 3: Pearson correlation of water quality parameters at UOP site.
|
Rainfall
|
pH
|
Conductivity
|
TDS
|
Salinity
|
pH
p value
|
0.173
0.287
|
|
|
|
|
Conductivity
p value
|
-0.286
0.073
|
0.135
0.406
|
|
|
|
TDS
p value
|
-0.274
0.052
|
0.256
0.110
|
0.788
0.000
|
|
|
Salinity
p value
|
-0.274
0.087
|
0.177
0.276
|
0.756
0.000
|
0.906
0.000
|
|
Hardness
p value
|
-0.448
0.005
|
-0.013
0.938
|
0.625
0.000
|
0.631
0.000
|
0.557
0.000
|
A very strong positive correlation of 0.968 at the Kandy sampling site was observed between salinity and TDS. Rainfall - TDS (-0.743) and rainfall - salinity (-0.700) showed a strong negative correlation. Conductivity - salinity (0.412) and pH - rainfall (0.570) showed moderate positive correlation without a significant relationship. Also, rainfall-conductivity (-0.572) and pH - conductivity (-0.486) showed moderate negative correlation. Other, pairs of water quality parameters showed week or very week correlations at the Kandy site (Table 4).
Table 4: Pearson correlation of water quality parameters at Kandy site.
|
Rainfall
|
pH
|
Conductivity
|
TDS
|
Salinity
|
pH
p value
|
0.570
0.182
|
|
|
|
|
Conductivity
p value
|
-0.572
0.180
|
-0.486
0.268
|
|
|
|
TDS
p value
|
-0.743
0.056
|
-0.194
0.679
|
0.321
0.483
|
|
|
Salinity
p value
|
- 0.700
0.080
|
-0.084
0.858
|
0.412
0.358
|
0.968
0.000
|
|
Hardness
p value
|
-0.104
0.825
|
0.321
0.483
|
-0.015
0.974
|
0.199
0.669
|
0.189
0.684
|
A very strong positive correlation at the Polgolla site was observed between conductivity -TDS (0.999), conductivity – salinity (0.999), TDS - salinity (0.998), conductivity - hardness (0.956), salinity - hardness (0.956), TDS - hardness (0.946) with significant relationship. The other pairs of water quality parameters showed week or very week correlations (Table 5).
Table 5: Pearson correlation of water quality parameters at Polgolla site.
|
Rainfall
|
pH
|
Conductivity
|
TDS
|
Salinity
|
pH
p value
|
0.149
0.750
|
|
|
|
|
Coductivity
p value
|
-0.337
0.460
|
0.102
0.828
|
|
|
|
TDS
p value
|
-0.331
0.469
|
0.076
0.872
|
0.999
0.000
|
|
|
Salinity
p value
|
-0.368
0.417
|
0.098
0.835
|
0.999
0.000
|
0.998
0.000
|
|
Hardness
p value
|
-0.302
0.511
|
0.346
0.447
|
0.956
0.001
|
0.946
0.001
|
0.956
0.001
|
Pearson correlation coefficient of anions
At the UOP site, a moderate positive correlation was observed between NO3- and SO42- with significant relationship. Hence, it is possible that these two anions were generated by the same sources which could be natural or anthropogenic sources such as traffic activities, construction activities, and fossil fuel combustion. Other pairs showed weak or very week correlations. In contrast, the Kandy site showed a very strong positive correlation between Cl- and SO42- with significant relationship. Hence, these two anions would have been released from the same sources. Furthermore, at the Polgolla site, a very strong positive correlation was observed between Cl- - NO3- and SO42-- Cl- with significant relationship, while a moderate positive significant correlation was observed between SO42-- NO3-. Hence, these two anions would have been released by the same source. The Pearson correlation coefficients of anions at the three sampling sites are represented in Table 6.
Table 6: Pearson correlation coefficient of anions at three sampling sites.
|
Cl-
|
NO3-
|
UOP site
|
|
|
NO3-
p value
|
-0.119
0.587
|
|
SO42-
p value
|
0.251
0.249
|
0.433
0.039
|
Kandy site
|
|
|
NO3-
p value
|
|
|
SO42-
p value
|
0.800
0.056
|
|
Polgolla site
|
|
|
NO3-
p value
|
0.918
0.010
|
|
SO42-
p value
|
0.830
0.041
|
0.589
0.218
|