DOI QR코드

DOI QR Code

Chemical characteristics of wet precipitation in urban and mountainous sites of Jeju Island

  • Bu, Jun-Oh (Researcher, Department of Chemistry and Cosmetics, Jeju National University) ;
  • Song, Jung-Min (Research Associate Professor , Department of Chemistry and Cosmetics, Jeju National University) ;
  • Park, Sook-Young (Graduate Student, Department of Chemistry and Cosmetics, Jeju National University) ;
  • Kang, Hee-Ju (Graduate Student, Department of Chemistry and Cosmetics, Jeju National University) ;
  • Kang, Chang-Hee (Professor, Department of Chemistry and Cosmetics, Jeju National University)
  • 투고 : 2019.12.30
  • 심사 : 2020.02.03
  • 발행 : 2020.02.25

초록

Wet precipitation samples were collected in Jeju City and Mt. Halla-1100 site (a site at an altitude of 1100 m on Mt. Halla) during 2011-2013, and their major ionic species were analyzed to examine the chemical composition and characteristics. A comparison of ion balance, electric conductivity, and acid fraction of precipitation revealed correlation coefficients in the range of r = 0.950~0.991, thereby implying the high quality of analytical data. Volume-weighted mean pH and electric conductivity corresponded to 4.86 and 25.5 µS/cm for Jeju City, and 4.98 and 15.1 µS/cm for Mt. Halla-1100 site, respectively. Ionic strengths of the wet precipitation in Jeju City and Mt. Halla-1100 site corresponded 0.3 ± 0.5 and 0.2 ± 0.2 mM, respectively, thereby indicating that more than 30 % of total precipitation was within a pure precipitation criteria. The precipitation with a pH range of 4.5 - 5.0 corresponded to 40.8 % in Jeju City, while the precipitation with a pH range of 5.0 - 5.5 corresponded to 56.9 % in Mt. Halla-1100 site, thereby indicating slightly more weak acidity than that in Jeju city. The volume-weighted mean concentration (µeq/L) of ionic species was in the order of Na+ > Cl- > nss-SO42- > NO3- > Mg2+ > NH4+ > H+ > nss-Ca2+ > PO43- > K+ > CH3COO- > HCOO- > NO2- > F- > HCO3- > CH3SO3- at Jeju City area, while it corresponded to Na+ > Cl- > nss-SO42- > NO3- > NH4+ > H+ > Mg2+ > nss-Ca2+ > PO43- > CH3COO- > K+ > HCOO- > NO2- > F- > HCO3- > CH3SO3- at Mt. Halla-1100 site. The compositions of sea salts (Na+, Cl-, Mg2+) and secondary pollutants (NH4+, nss-SO42-, NO3-) corresponded to 66.1 % and 21.8 %, respectively, in Jeju City and, 49.9 % and 31.5 %, respectively, in Mt. Halla-1100 site. The acidity contributions in Jeju City and Mt. Halla-1100 site by inorganic acids, i.e., sulfuric acid and nitric acid, corresponded to 93.9 % and 91.4 %, respectively, and the acidity contributions by organic acids corresponded to 6.1 % and 8.6 %, respectively. The neutralization factors in Jeju City and Mt. Halla1100 site by ammonia corresponded to 29.8 % and 30.1 %, respectively, whereas the neutralization factors by calcium carbonate corresponded to 20.5 % and 25.2 %, respectively. From the clustered back trajectory analysis, the concentrations of most ionic components were higher when the airflow pathways were moved from the continent to Jeju area.

키워드

1. Introduction

Wet precipitation plays an important role in the washing and elimination of air pollutants from the atmosphere.1-2 Therefore, the chemical compositions of precipitation mainly depends on the major pollutants in ambient air.3 Removal of air pollutants by wet precipitation mainly occurs in two scavenging mechanisms. The first mechanism is termed as in-cloud scavenging (rainout), which mostly absorbs or condenses gaseous pollutants and fine particles around the cloud droplet. The other mechanism is termed as below-cloud scavenging (washout), which is the most efficient pathway to remove fine or coarse particles.1,3

The atmospheric pollutants removed by the processes affect the acidification of precipitation. Specifically, the acidification of precipitation is closely related to the increase in the use of fossil fuels, and it is recognized as a result of rapid industrialization and industrial economy.4 Increases in the use of fossil fuels leads to various atmospheric pollutants including sulfur oxides (SOX) and nitrogen oxides (NOX), which is reported to intensify the acidification of precipitation.4 Specifically, SOX and NOX, correspond to typical pollutants that cause the acidification of precipitation and, are converted into strong acidic substances by ozone, hydrocarbons, hydrogen peroxide, and peroxides (e.g. HCOO•, RCOO•).5

The acidification of precipitation leads to destruction of aquatic ecosystems, degradation of crop productivity, destruction of forests, and erosion of buildings. It is typically affected by local and long-range transported pollutants. Therefore, the transboundary of acidic air pollutants including sulfur oxides and nitrogen oxides constitutes a significant environmental controversy between neighboring countries. North America and Northwest Europe have already suffered from serious acid rain and are faced with environmental controversies between neighboring countries. Furthermore, the countries are continuing to examine policy formulation and research to reduce acid rain.6-7

In East Asia, since the mid-1990s, the need for international cooperative research on monitoring of acid rain increased via ongoing consultation with the experts of International Conference on Acid Rain. Subsequently, Acid Deposition Monitoring Network in East Asia (EANET) is in operation since 2001. In Korea, since 1983, the ministry of environment has operated an acid rain measuring network in major cities including Seoul and industrial complexes.7

Jeju Island is a typical background site in Korea with few industrial facilities, low population density, and low pollution sources. Thus, it constitutes an advantageous location to monitor long-range transported air pollutants from the continent of Northeast Asia.4,7 In the study, the ionic components of precipitation samples collected at Jeju City area and Mt. Halla1100 site which are typical urban and mountainous sites in Jeju Island, were analyzed to investigate chemical and pollution characteristics with respect to altitude and airflow transport pathways.

2. Experimental

2.1. Collection of precipitation samples

A total of 147 precipitation samples were collected in the Jeju City area (Jeju National University; 33.26°N, 126.33°E) and 58 precipitation samples were collected in Mt. Halla-1100 site (33.21°N, 126.27°E) during 2011-2013. Automated rain samplers (model SL-4- 001, Shinil Science Co., Gyeonggi, Korea) equipped with a rain sensor and stainless steel bucket (253 mm diameter) were used to collect precipitation samples. The samplers were installed on the roof of a building located in Jeju National University campus and the top of a trailer in Mt. Halla-1100 site. The HDPE bottles and sampling bucket were washed frequently before and after precipitation events with distilled water.

The precipitation samples were transferred to the laboratory, and a few of the samples were used to measure pH and electric conductivity. The remaining samples were partitioned into separate aliquots. An aliquot was preserved in a freezer at -20 ℃ without pretreatment for the analysis of inorganic ions and an other aliquot was stored at the same temperature for the analysis of organic acids after adding a drop of chloroform.4,7-8

2.2. Chemical analysis

The pH value of each precipitation sample was measured via a pH meter (Orion 720A, 81-02 electrode; Thermo Fisher Scientific, Waltham, MA, USA). Two low ionic strength buffer solutions of pH 4.10 and 6.97 at 25 ℃ were used for calibration prior to the measurement. Electric conductivity was measured via a conductivity meter (Orion 3-Star, 013005MD electrode, Thermo Electron, Waltham, MA, USA) at the same temperature.

Ion concentrations of the precipitation samples were analyzed for cations (Na+ , NH4+ , K+ , Mg2+, Ca2+) and anions (SO42− , NO3 , Cl , PO43− , NO2 , F , CH3SO3 , HCOO , CH3COO ) via ion chromatography (Metrohm Modula IC). With respect to the analysis of ionic species, the conditions of ion chromatography (IC) are shown in Table 1. The instrument detection limit (IDL) and coefficient of variation (CV) of the ion chromatography analyses are listed in Table 2.

Table 1. Instrumental conditions for ion chromatography analysis

BGHHBN_2020_v33n1_33_t0001.png 이미지

Table 2. Instrument detection limit (IDL) and coefficient of variation (CV) for the analysis of ionic components (n=7)

BGHHBN_2020_v33n1_33_t0002.png 이미지

2.3. Back trajectory analysis

Airflow pathway patterns were derived from a five day back trajectory analyses during the study periods. They were calculated via the PC version of HYSPLIT v. 4.0 (HYbrid Single-Particle Lagrangian Integrated Trajectory) developed by the National Oceanic and Atmospheric Administration/Air Resources Laboratory (NOAA/ARL).9 Specifically, HYSPLIT was used to calculate back trajectories using meteorological data with a resolution of 1 × 1° that was generated by via the Global Data Assimilation System (GDAS) model from the National Weathers Service’s (NWS) National Centers for Environmental Prediction (NCEP) (ftp:// arlftp.arlhq.noaa.gov/archives/gdas1).7 The total run time corresponded to 120 h, and the starting heights of the trajectories corresponded to 320 m in Jeju City and 1100 m in the Mt. Halla-1100 site.

3. Results and Discussion

3.1. Quality control of analytical data

It is necessary for the precipitation analysis to measure trace level concentrations of soluble ionic species. Therefore, the analytical data can vary based on the sample treatment and analytical process.7 In precipitation analysis, the comparison methods of ion balance, electric conductivity and acid fraction are effectively applicable to improve the accuracy of the analyzed data. If the analysis of the precipitation components is performed quantitatively, then the ion balance between the sums of cations and anions should be correlated well. Furthermore, the correlation between the electric conductivity calculated by the ion concentrations and the electric conductivity measured by the conductivity meter should be high to significantly improve the accuracy of analytical data.7,10 Additionally, the accuracy of analytical data can be inferred via calculating the acid fractions and their correlations.11

In the study, the ion balances, i.e. the correlation coefficient between the sums of the cation and anion equivalent concentrations corresponded to 0.980 in Jeju City, whereas it was 0.991 in Mt. Halla-1100 site. Additionally, the comparison between the calculated and measured electric conductivities exhibited high correlation coefficients of 0.986 and 0.971 in Jeju City and Mt. Halla-1100 site, respectively. Furthermore, the correlation coefficients between the acid fraction calculated from ionic concentrations and that from pH and conductivity corresponded to 0.962 and 0.950 in Jeju City and Mt. Halla-1100 site, respectively (Fig. 1). The results from three comparison methods confirmed that the quality of analytical data is good and available.

BGHHBN_2020_v33n1_33_f0001.png 이미지

Fig. 1. Evaluation of analytical data by the comparison of ion balance (a), electric conductivity (b), and acid fraction (c).

The percent imbalance (%) between the sums of cationic and anionic equivalent concentrations (Σ[Cation] and Σ[Anion]) was calculated based on the equation (Σ[Cation] − Σ[Anion])/(Σ[Cation] + Σ[Anion]) × 100 to examine the variance of the analytical data (Fig. 2). The results indicated that five (3.4 %) out of the total 147 precipitation samples exhibited imbalances exceeding 25 % in Jeju City. Conversely, the results revealed that six (10.3 %) out of the total 58 precipitation samples exhibited imbalances exceeding 25 % in Mt. Halla-1100 site. The results satisfy a common criterion of relative error, which is typically required for the trace level analysis in experiments.7

BGHHBN_2020_v33n1_33_f0002.png 이미지

Fig. 2. Scattergram of (Σ[Cation] − Σ[Anion])/(Σ[Cation] + Σ[Anion])×100 vs. (Σ[Cation] + Σ[Anion]).

3.2. Concentrations of precipitation components

The volume-weighted mean pH and electric conductivity of the precipitation measured at Jeju City corresponded to 4.86 and 25.5 μS/cm, respectively. Conversely, the values at Mt. Halla-1100 site corresponded to 4.98 and 15.1 μS/cm, respectively. The ionic strength of precipitation corresponded to 0.3 ± 0.5 mM, thereby indicating that 31.3 % of the total precipitation was within the pure precipitation criteria in Jeju City, whereas the ionic strength corresponded to 0.2 ± 0.2 mM in Mt. Halla-1100 site and 46.6 % of the total precipitation samples corresponded to the pure precipitation.10

Frequency of the precipitation with a pH range of 4.5 − 5.0 corresponded to 40.8 % in Jeju City, and this is the most proportional among the total precipitation samples. Conversely, in Mt. Halla-1100 site, the frequency with a pH range of 5.0 − 5.5 corresponded to 56.9 %, thereby indicating a little more basic than in Jeju city. The results suggested that the precipitation acidity can be more affected by altitude and air pollution in Jeju city. Based on the report of the Global Atmosphere Watch (GAW, 2013) on the acidity of precipitation, the pH range of precipitation between 1997 and 2013 corresponded to 4.68 − 4.88, thereby indicating the absence of regional differences or temporal variations.12 The results of the study at two sites were similar to those in previous studies in Jeju Island.

The volume-weighted mean concentration (μeq/L) of the ionic species was in the order of Na+ > Cl > nss-SO42− > NO3 > Mg2+ > NH4+ > H+ > nss-Ca2+ > PO43− > K+ > CH3COO > HCOO > NO2 > F > HCO3 > CH3SO3 at Jeju City area, while it was in the order of Na+ > Cl > nss-SO42− > NO3 > NH4+ > H+ > Mg2+ > nss-Ca2+ > PO43− > CH3COO > K+ > HCOO > NO2 > F > HCO3 > CH3SO3 at Mt.Halla-1100 site (Table 3). The concentration of non-sea salt sulfate (nss-SO42− ) and non-sea salt calcium (nss-Ca2+) are obtained by the equation corresponding to ‘[nss-SO42−] = [SO42−] − 0.251×[Na+]’ and ‘[nssCa2+] = [Ca2+] − 0.044×[Na+]’, respectively.

Table 3. Volume-weighted mean ionic concentrations (μeq/L) of precipitation in Jeju Island and other sites

BGHHBN_2020_v33n1_33_t0003.png 이미지

Overall, the concentrations of sea salt species exceeded those of other species, and next by the secondary pollutant species in two sites in Jeju area. Based on the concentrations of precipitation components, the sea salts (Cl , Na+ , Mg2+) and secondary pollutants (NH4+, nss-SO42− , NO3 ) contributed to the total composition of ionic precipitation species by 66.1 % and 21.8 %, respectively, in Jeju City. The composition of sea salt and secondary pollutants corresponded to 49.9 % and 31.5 %, respectively, in Mt. Halla-1100 site, thereby suggesting that the effect of sea salts is relatively high in low altitude areas (Fig. 3).

BGHHBN_2020_v33n1_33_f0003.png 이미지

Fig. 3. Chemical composition ratios of precipitation components.​​​​​​​

Table 3 shows the ionic concentrations of precipitation obtained from Jeju and other domestic/foreign sites. The results indicate that, the concentrations of sea salt components are higher in island areas, such as Jeju City, Mt. Halla-1100 site, and Okinawa, than in urban areas such as Shenzhen and Shanghai. Conversely, secondary pollutants, such as ammonium, nitrate, and sulfate, exhibit higher concentrations in urban areas than in Jeju City and Mt. Halla-1100 site.12-15

3.3. Acidification and neutralization characteristics

The acidification of precipitation is generally affected by the organic acids (RCOOH) and inorganic acids (H2SO4, HNO3, HCl). Therefore, their acidification contribution rate can be deduced via a comparison of acidic anions, such as SO42− , NO3, Cl , and RCOO- concentrations.7-8 However, Cl is considered to be mainly from the sea salt in Jeju Island, and thus the effect of Cl on the pH of precipitation is negligible.16-17 In the study, the acidification contribution (AC) was calculated via Eq. (1) to confirm the effects of acidic substances on the acidification of precipitation.8

\(A C(\%)=\frac{[X]}{\left[n s s-S O_{4}^{2-}\right]+\left[N O_{3}^{-}\right]+\left[H C O O^{-}\right]+\left[C H_{3} C O O^{-}\right]} \times 100\)       (1)

X = concentration of acidic substance

The acidification contributions of inorganic and organic acids are calculated in the equation, and the results are listed in Table 4. During the study period, the acidification contribution by inorganic acids in Jeju City and Mt. Halla-1100 site corresponded to 93.9 % and 91.4 %, respectively. Conversely, the acidification contribution by inorganic acids at two sites corresponded to 6.1 % and 8.6 %, respectively. The results confirmed that the acidification of precipitation mainly occurred via inorganic acids in Jeju Island.

Table 4. Acidification contributions (%) of acidic anions in precipitation​​​​​​​

BGHHBN_2020_v33n1_33_t0004.png 이미지

The neutralization contribution is calculated by the neutralization factors (NFs) from Eqs. (2) and (3).7,18-19 The results are listed in Table 5.

Table 5. Neutralization factors (NF) for ammonia and calcium carbonate​​​​​​​

BGHHBN_2020_v33n1_33_t0005.png 이미지

\(N F_{N H_{4}^{+}}=\frac{\left[N H_{4}^{+}\right]}{\left[n s s-S O_{4}^{2-}\right]+\left[N O_{3}^{-}\right]+\left[H C O O^{-}\right]+\left[C H_{3} C O O^{-}\right]}\)       (2)

\(N F_{n s s-C a^{2+}}=\frac{\left[n s s-C a^{2+}\right]}{\left[n s s-S O_{4}^{2-}\right]+\left[N O_{3}^{-}\right]+\left[H C O O^{-}\right]+\left[C H_{3} C O O^{-}\right]}\)       (3)

Acidic substances are mainly neutralized by basic substances, such as NH3, MgCO3, and CaCO3, and neutralization is estimated by measuring NH4+ , Mg2+, and Ca2+ concentrations, which correspond to basic cations of precipitation. Specifically, NH3 and CaCO3 are known to mainly contribute to neutralization reaction of acidis substances.16-17 The neutralization factors of the precipitation in Jeju City by ammonia and calcium carbonate corresponded to 0.30 and 0.21, respectively. Similarly, the values for the Mt. Halla-1100 site corresponded to 0.30 and 0.25, respectively. With respect to Jeju Island, the sums of the neutralization factors by ammonia and calcium carbonate approximately corresponded to 51~55 %, and two components were observed to mainly affect the neutralization reaction.

3.4. Effect of airflow transport pathways

Northeast Asia is characterized by rapid economic growth, and the long-range transport of air pollutants is consistently examined since the mid-1980s by organizations such as EANET and LTP.8,20 The long-range transport of air pollutants is very sensitive and constitutes a significant issue between neighboring countries in this area, and thus it is necessary to understand and clarify their transboundary transport mechanism.7 Based on the sampling days at two sampling sites (147 days and 58 days), clustered back trajectory analysis was performed using NOAA’s HYSPLIT 4 model and GDAS meteorological data to confirm transport pathways of air masses (Fig. 4).21

BGHHBN_2020_v33n1_33_f0004.png 이미지

Fig. 4. Clustered 5-day back trajectories of air masses corresponding to precipitation sampling date.​​​​​​​

In the cluster analysis at Jeju City, 147 trajectories were classified into three clusters. Similarly, three clusters were identified at the Mt. Halla-1100 site. Thus, the clustered back trajectories were grouped into three pathways in Jeju City and Mt.Halla-1100 site as follows: In Jeju City, Cluster 1 (Eastern Coast of China & Korean Peninsula) corresponded to 60 %, Cluster 2 (Siberia, Mongolia & China Continent) corresponded to 24 %, and Cluster 3 (East sea & The North Pacific Ocean) corresponded to 16 %. In Mt. Halla-1100 site, Cluster 1 corresponded to 52 %, Cluster 2 corresponded to 21 %, and Cluster 3 corresponded to 27 %. We compare the concentration of ionic components, such as nss-SO42− . NO3 , NH4+ , and nss-Ca2+ by airflow transport pathway, and their components exhibited high concentrations when the airflow moved to Jeju from Siberia, Mongolia, and China continent (Cluster 2) (Table 6). Additionally, the concentrations of ionic components in Cluster 1 and Cluster 2 were high in Jeju City, which exhibited a relatively low altitude, and the concentrations of ionic components were low in Cluster 3. However, in high altitudes, such as the Mt. Halla-1100 site, the concentrations of ionic components in Cluster 3 exhibited high concentrations than that in the case of Cluster 1. The results indicated that the concentrations of components were different based on the airflow transport pathway and altitude.

Table 6. Ionic concentrations (μeq/L) of precipitation in accordance with airflow transport pathways by cluster back trajectory​​​​​​​

BGHHBN_2020_v33n1_33_t0006.png 이미지

4. Conclusions

In the study, the chemical compositions of precipitation collected in urban and mountainous sites of Jeju Island were compared to understand the differences each order. Additionally, their pollution characteristics were investigated relative to the altitude and airflow transport pathways.

The volume-weighted mean pH of precipitation in Jeju City and Mt. Halla-1100 site was weak acidic and corresponded to 4.86 and 4.98, respectively, and the electric conductivity corresponded to 25.5 μS/cm and 15.1 μS/cm, respectively. The precipitation with a pH range of 4.5−5.0 corresponded to 40.8 % in Jeju City, and the precipitation with a pH range of pH 5.0−5.5 corresponded to 56.9 % in Mt. Halla-1100 site, thereby indicating a little slightly weaker acidity than in Jeju city. Furthermore, the ionic strengths of precipitation at two sites corresponded to 0.3 ± 0.5 mM and 0.2 ± 0.2 mM, respectively, thereby indicating the pure precipitation criteria corresponded to 31.3 % and 46.6 %.

The acidification contributions by inorganic acid in Jeju City and Mt. Halla-1100 site corresponded to 93.9 % and 91.4 %, respectively, and the neutralization factors by ammonia and calcium carbonate approximately corresponded to 51~55 %.

The results of clustered back trajectory analysis indicated, that the concentrations of most ionic components were higher when the airflow pathways moved from the Asian continent to Jeju area.

Acknowledgements

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2016R1D1A3B03935687), and supported by the 2019 scientific promotion program funded by Jeju National University

참고문헌

  1. J. H. Seinfeld and S. N. Pandis, 'Atmospheric Chemistry and Physics: From Air Pollution to Climate Change', 3rd Ed., John Willey & sons, New Jersey, 2016.
  2. L. Wang, Z. Shen, D. Lu, Q. Zhang, T. Zhang, Y. Lei and H. Xu, Atmospheric Pollution Research, 10, 395-403 (2019). https://doi.org/10.1016/j.apr.2018.08.011
  3. A. Roy, A. Chatterjee, S. Tiwari, C. Sarkar, S. K. Das, S. K. Ghosh and S. Raha, Atmos. Res., 181, 44-53 (2016). https://doi.org/10.1016/j.atmosres.2016.06.005
  4. J. O. Bu, J. M. Song, S. H. Shin, W. H. Kim and C. H. Kang, J. Korean Soc. Environ. Anal., 20(1), 1-10 (2017).
  5. C. H. Kang, W. H. Kim, S. B. Hong, K. H. Lee, M. S. Hong and S. G. Shim, J. KOSAE, 15(5), 555-566 (1999).
  6. N. K. Kim, Y. P. Kim, C. H. Kang and K. C. Moon, J. KOSAE, 20(1), 119-128 (2004).
  7. J. O. Bu, J. M. Song, L. K. Kim, S. Y. Kim, W. H. Kim and C. H. Kang, Bull. Korean Chem. Soc., 40(2), 157-162 (2019). https://doi.org/10.1002/bkcs.11661
  8. J. O. Bu, J. M. Song, S. H. Shin, W. H. Kim and C. H. Kang, Anal. Sci. Technol., 29(1), 19-28 (2016). https://doi.org/10.5806/AST.2016.29.1.19
  9. R. R. Draxler and G. D. Rolph, HYSPLIT (HYbrid Single-Particle Lagrangian Integrated Trajectory) Model access via NOAA ARL READY Website, http://ready.arl.noaa.gov/HYSPLIT_traj.php, 2013.
  10. R. Sequeira and F. Lung, Atmos. Environ., 29(18), 2439-2447 (1995). https://doi.org/10.1016/1352-2310(95)00161-Q
  11. M. Kramer, M. Schule and L. Schutz, Atmos. Environ., 30(19), 3291-3300 (1996). https://doi.org/10.1016/1352-2310(96)00070-2
  12. Climate Change Information Center, 'Report of Global Atmospheric Watch 2013', Korea, 2013.
  13. H. Sakihama, M. Ishiki and A. Tokuyama, Atmos. Environ., 42(10), 2320-2335 (2008). https://doi.org/10.1016/j.atmosenv.2007.12.026
  14. D. Y. Huang, Y. G. Xu, B. Zhou, H. H. Zhang and J. B. Lan, Environ. Monit. Assess., 171, 429-439 (2010). https://doi.org/10.1007/s10661-009-1289-7
  15. K. Huang, G. Zhang, C. Xu, Y. Wang and A. Tang, Atmos. Res., 89, 149-160 (2008). https://doi.org/10.1016/j.atmosres.2008.01.006
  16. Bu, J. O. 'Seasonal Characteristics of Rainwater Components Jeju Area during 2009-2010' Master's degree, Jeju National University, Korea, 2013.
  17. Lee, H. J. 'Acidity and Pollution Characteristics of Precipitation in Jeju City during 2015-2016' Master's degree, Jeju national University, Korea, 2017.
  18. N. Zhang, Y. He, J. Cao, K. Ho and Z. Shen, Atmos. Res., 106, 50-60 (2012). https://doi.org/10.1016/j.atmosres.2011.11.006
  19. S. Tiwari, D. M. Chate, D. S. Bisht, M. K. Srirastava and B. Padmanabhamurty, Atmos. Res., 104-105, 128-138 (2012). https://doi.org/10.1016/j.atmosres.2011.09.006
  20. K. J. Kim, J. O. Bu, W. H. Kim, Y. S. Lee, D. R. Hyeon and C. H. Kang, J. KOSAE, 29(6), 818-829 (2013). https://doi.org/10.5572/KOSAE.2013.29.6.818
  21. J. M. Song, J. O. Bu, J. Y. Lee, W. H. Kim and C. H. Kang, Asian J. Atmos. Environ., 11(4), 313-321 (2017). https://doi.org/10.5572/ajae.2017.11.4.313