서 론
농약류에는 유기인계(organophosphorous), 유기염소계(organochlorine), 카바메이트계(carbamate), 피레스로이드계(pyrathroid), 클로로니코티닐계(clothianidin), 페녹시계(phenoxy), 아닐린계(aniline) 등 다양한 종류가 있으며, 살충제(insecticide), 살균제(germicide), 제초제(herbicide), 생장조정제(growth regulator), 유인제(attracting agent), 기피제(repellent), 전착제(spreader) 등으로 쓰이고 있다.1 본 연구에서는 국내에 유통되었거나 유통되고 있는 농약 중 출하량2 및 독성3−13을 고려하여 살충제인 메소밀, 알디캅, 카바릴, 카보퓨란과 제초제인 2-methyl-2,4-dichlorophenoxyacetic acid(2,4-D), 4-chlorophenoxy-acetic acid(MCPA), molinate를 분석항목으로 선택하였다. 이와 같이 다양한 효능을 보이는 많은 농약들이 동물대상의 연구에서 발암성, 생식독성, 만성신경독성, 면역독성과 같은 만성 건강장해에 대한 충분한 근거를 보고하고 있고, 이러한 만성 건강영향은 사람에서도 기능장애나 암 등의 만성질환이 발생 가능하다는 점에서 중요성이 있다.3 따라서 농약 노출에 대한 전반적인 실태조사와 관리가 지속적으로 이루어져야 하며, 합성 원료인 quinoline과 비스페놀-A의 경우 각종 산업단지 등에서 널리 쓰이고 있으므로 이들에 대한 하천수 내 관리도 필요하다.
현재 2,4-D는 논 제초제로 주로 사용하고 있고, 인간에 대한 독성으로 경구, 경피, 흡입 경로에 따른 신경독성 등이 언급되고 있으며,4 인체 발암가능성 물질 2B로 분류되었다.5 메소밀은 맹독성 살충제로, 쥐의 급성 경구독성 LD50(암컷)가 30 mg/kg이다.6 알디캅은 살충제로, 포유동물에 대한 급성독성은 매우 강하여 쥐에 대한 급성 경구독성 LD50가 0.9 mg/kg이다.7 그리고, 인체에 많은 양에 노출되면 시력 저하, 메스꺼움, 두통 등을 일으키고, 호흡체계를 마비시킬 수 있다.8 MCPA는 호르몬형 선택성 제초제로 쥐에 대한 급성경구독성 LD50가 700~1,330 mg/kg이고, 노출 시 빈혈, 근육약화 위장 장애, 간 손상 등이 보고되었다.9 Molinate는 논 제초제로, 수생생물에 강한 독성을 나타내고, 갑상선 기관의 독성작용이 보고되었으며, 구역질, 설사, 복통, 발열, 결막염 등을 유발할 수 있다.10 카바릴은 접촉살충제로서 인간에게 발암가능성 있는 물질이며, 쥐에 대한 경구독성은 250~850 mg/kg이다. 카보퓨란은 살충제로, 쥐에 대한 급성 경구독성 LD50는 8~14 mg/kg, 개에 대한 급성 경구독성 LD50는 19 mg/kg으로 높은 독성을 보이며, 신경독성이 있다.11 합성원료의 경우, quinoline은 생산 산업간 중간체로 발생되고 급성 흡입노출은 눈, 코, 목을 자극하고, 두통, 현기증 등을 일으킬 수 있으며, 만성섭취시 간 손상이 의심된다. 미국 EPA에서 발암물질 그룹C로 분류하였다.12 또한 비스페놀-A는 플라스틱 제조 및 산화방지제의 원료 등으로 사용되고 내분비계장애물질로 적은 양에서도 신경 발달에 문제를 일으킬 수 있으며, 급성노출시, 피부, 눈, 점막, 상부기도에 자극성이 있다.13
최근 연구에서는 하천수 안의 농약류를 검출하기 위한 다양한 전처리 과정 및 기기분석법이 시도되고 있다. 전처리 과정으로 회전 디스크 흡착 추출법(rotating disk sorptive extraction, RDSE),14 고체상 미량추출법(solid phase microextraction, SPME),15,30 고체상 추출법(solid phase extraction, SPE),15,18,20,21,23,24,26,32 단일 방울 미량추출법(single-drop microextraction, SDME),16 직접 담금 고체상 미량추출법(direct immersion, DI)-SPME,17,19 직접 주입법(direct injection, DI),22 교반 막대 흡착 추출법(stir bar sorptive extraction, SBSE),25 액액 미량추출법(liquid-liquid microextraction, LLME),27 액체상 미량추출법-분산 액액 추출법(liquid phase microextractiondispersive liquid-liquid microextraction, LPME-DLLME),28 초음파 이용 에멀젼화(ultrasound-assisted emulsification, USAE),29 헤드 스페이스 고체상 미량추출법(headspace-SPME, HS-SPME),31 직접 주입법(directly injection, DI)33 등이 있다.
그리고, 고감도 기기분석을 위해서 기체 크로마토그래피 질량분석법(gas chromatography mass spectrometery, GC-MS)14−19 또는 액체 크로마토그래피 질량분석법(liquid chromatography mass spectrometery or tandem mass spectrometery, LC-MS or LC-MS/MS)20−26을 주로 사용하고 있고, 이 외에도 GC-불꽃 이온화 검출기(GC-flame ionisation detector, GC-FID),27 GC-불꽃 광도 검출기(GC-flame photometric detector, GC-FPD),28 GC-전자포획 검출기(GC-electron capture detector, GC-ECD),29−31 LC-다이오드 배열 검출기(LC-diode array detector, LC-DAD),32 이온 크로마토그래피-형광 검출기(ion chromatography-fluorescence detector, IC-FLD)33 등을 활용하여 분석하였다.
하천수 내 극미량의 분석물질을 검출하기 위해서는 높은 농축계수와 효과적인 방해물질 제거과정 등의 고감도 분석기술이 필요하다. 그래서, 미국 환경부(USEPA) method34−36 등 여러 연구문헌 조사를 통한 농약류 분석법을 토대로 합성 원료도 동시에 분석할 수 있는 최적의 SPE-LC-MS/MS법을 개발하였다. 이 분석법으로 농약류 7종 및 합성원료 2종 각각에 대해 극미량까지 분석 가능한 낮은 정량한계뿐 아니라 신뢰도 높은 정도관리 결과를 얻었으며, 하천수 10개 지점에 적용하여 분석을 실시하였다.
재료 및 방법
재료
본 연구에 사용된 용매는 Merck(Darmsadt, Germany)의 잔류농약 분석용 시약을 이용하였고, 내부표준물질(internal standard, IS)을 포함한 각종 표준물질은 Sigma-Aldrich(Taufkirchen, Germany)의 고농도 표준용액을 희석하여 사용하였다. 분석용 메탄올과 아세토나이트릴은 Burdick & Jackson(Muskegon, Michigan, USA)의 시약을 사용하였다. 암모늄 아세테이트는 Sigma-Aldrich, 개미산(formic acid)과 염산은 Junsei(Tokyo, Japan)의 분석용 시약을 사용하였다. 정제수는 초순수 장치 milli-Q(millipore corp., milford, MA, USA)를 사용하여 정제하였으며, HLB(hydrophilic lipophilic balance)와 MCX (mixed-mode cation exchange) cartridge는 Waters(Milford, MA, USA)의 Oasis® 시료 추출용 6 cc 제품을 사용하였다. 거름종이와 실린지 필터는 Advantec(Tokyo, Japan)사 제품을 사용하였고, 모든 초자 기구는 세제, 증류수, 아세톤, 메탄올로 연속하여 세척한 후 350 ℃에서 4시간 이상 가열한 후 오염이 없는 곳에서 방냉한 것을 사용하였다.
시료채취
시료는 인근 하천인 금강의 10개 지점에서 채수하였다. 시료용기는 미리 세척·건조한 Environmental Sampling Supply(San Leandro, CA, USA)사 제품의 테플론 마개가 있는 1 L 시료채취용 갈색유리용기로서 유기물이 용출되지 않음을 증명하는 인증서를 포함하고 있다. 채수한 시료는 4 ℃ 온도를 유지하기 위해 현장에서 얼음을 넣은 아이스박스에 보관하였으며, 채수한 시료는 실험실로 이송하여 분석 전까지 냉장 보관하였다.
농약류 7종 동시분석을 위한 전처리 과정
0.45 μm 거름종이를 사용하여 거른 후 시료 500 mL를 6N 염산으로 pH 2.0으로 조절 후 내부표준물질(카보퓨란-d3)을 첨가하였다. HLB 카트리지를 진공감압장치(vacuum manifold)에 장착한 후 증류수 4 mL와 메탄올 4 mL를 흘려주고 다시 증류수 4 mL, pH 2로 조절된 증류수 5 mL를 차례로 통과시켜 준비하였다. 그리고 HLB 카트리지에 시료를 10 mL/min의 속도로 통과시켰다. 시료가 통과된 HLB 카트리지를 건조시킨 다음 메탄올 6 mL로 용리시켰다. 이 용리액을 질소농축기를 사용하여 100 μL까지 농축 시킨 다음 30% 메탄올 100 μL로 재 용해하여 0.2 μm 실린지 필터로 여과시킨 후 2 mL 바이알에 옮겨 LC-MS/MS로 분석하였다.
합성원료 2종 동시분석을 위한 전처리 과정
0.45 μm 거름종이를 사용하여 거른 후 시료 500 mL를 NaHCO3:K2CO3 (2:1) 완충용액을 사용해 pH 10으로 조절 한 다음 내부표준물질(quinoline-d7, 비스페놀-A-d16)을 첨가하였다. MCX 카트리지를 진공감압장치(vacuum manifold)에 장착한 후 증류수 4 mL와 메탄올 4 mL를 흘려주고 다시 증류수 4 mL, pH 10으로 조절한 증류수 5 mL를 차례로 통과시켜 준비하였다. 그리고 MCX 카트리지에 시료를 10 mL/min의 속도로 통과시켰다. 시료가 통과된 MCX 카트리지를 건조시킨 다음 메탄올 6 mL로 용리시켰다. 이 용리액을 질소농축기를 사용하여 100 μL까지 농축 시킨 다음 30% 메탄올 100 μL로 재 용해하여 0.2 μm 실린지 필터로 여과시킨 후 2 mL 바이알에 옮겨 LC-MS/MS로 분석하였다.
LC-MS/MS 분석
LC-MS/MS는 Agilent(Palo Alto, CA, USA) 1200 series인 LC와 Agilent 6460 series인 triple quadruple instrument를 사용하였고, ESI 모드로 분석하였다. 컬럼은 50 mm × 2.1 mm, 1.8 μm pore size의 Agilent사 Eclipse Plus C18을 사용하였으며, 이동상(mobile phase)의 유속(flow rate)은 0.2 mL/min으로 일정하게 유지하였다. 이동상 A는 농약 분석 시 0.1% 개미산(formic acid, 이동상 A1), 합성원료 분석 시 10 mM 암모늄 아세테이트(이동상 A2)를 선택하였다. 이동상 B로는 아세토나이트릴을 사용하였고, 이동상 구배 조건은 분석시작시 이동상 B를 0%로 하고, 1분 후 5%, 10.5분에는 100%, 11분에는 100%, 12분에는 0%가 되도록 설정하였다(Table 1). 정성 및 정량분석을 위한 다중반응 모니터링(multiple reaction monitoring, MRM)은 Table 2와 같다.
Table 1.LC-MS/MS conditions of pesticides and synthetic materials
Table 2.LC-MS/MS MRM conditions
결과 및 고찰
방법검출한계 및 정량한계
검출한계(the low limit of detection, LOD) 및 정량한계(the limit of quantitation, LOQ)는 증류수가 아닌 매질효과를 고려한 분석물질이 검출되지 않은 하천수를 이용하여 0.005 또는 0.01 µg/L가 되도록 분석물질을 첨가한 시료 7개를 준비하고, 각 시료를 실험절차와 동일하게 추출하여 표준편차를 얻은 후, 표준편차에 3.14를 곱한 값을 방법검출한계, 10을 곱한 값을 정량한계로 계산하였다. 위에서 구한 농도가 실제 시료에서 검출될 수 있는지를 검증하기 위해 분석물질이 검출되지 않은 하천수를 이용하여 위에서 구한 정량한계로 첨가한 후 실험절차와 동일하게 추출하여 검출하였을 때 signal/noise(S/N)비가 3배 이상인지 확인하였다. 농약류의 LOD와 LOQ 분석결과, 메소밀은 각각 0.0008 µg/L, 0.003 µg/L이었다. Molinate는 0.002 µg/L, 0.007 µg/L이었고, MCPA는 0.002 µg/L, 0.006 µg/L이었다. 또한 카바릴과 알디캅은 똑같이 0.001 µg/L, 0.004 µg/L이었고, 2,4-D는 0.003 µg/L, 0.009 µg/L이었다. 가장 낮은 LOD과 LOQ 값을 보인 카보퓨란은 각각 0.0006 µg/L, 0.002 µg/L이었다. 합성 원료류의 LOD와 LOQ 분석결과, quinoline은 0.004 µg/L, 0.011 µg/L이었고, 비스페놀-A는 0.003 µg/L와 0.009 µg/L이었다.
검정곡선과 직선성
검정곡선 작성 및 검증은 정량범위 내의 5개 농도에 대해 검정곡선을 작성하고 얻어진 검정곡선 결정계수(r2)가 0.98 이상이어야 하며 결정 계수가 허용범위를 벗어나면 재 작성 하도록 하였다. 농약류의 검정곡선과 직선성 분석결과, Table 3과 같이 검출되는 농도범위 내에서 모든 검정곡선의 직선성(r2)이 0.997 이상으로 좋은 결과를 나타내었다.
Table 3.The LOD, LOQ and linear equation results for the analysis of seven pesticides and two synthetic materials in surface water
정밀도와 정확도
정밀도(precision)와 정확도(accuracy)는 분석물질이 검출되지 않은 하천수에 방법검출한계의 3~100배 정도의 농도수준이 되도록 표준물질을 동일하게 첨가한 시료(농약류의 경우 0.01과 0.05 µg/L, 합성원료의 경우 0.1과 0.5 µg/L)를 5개씩 준비하고 시험 분석절차와 동일하게 진행하여 구하였다. 농약류의 경우, 2.4~16.4%의 정밀도와 83.5~110.0%의 정확도를 보였다. 합성원료의 경우는 4.7~8.4%의 정밀도와 84.0~107.0%의 정확도를 보였다(Table 4).
Table 4.Intra-day laboratory precision and accuracy results for the analysis of analytes in surface water (n=5)
회수율
회수율(recovery)은 분석물질이 검출되지 않은 하천수에 방법검출한계의 3~100배 정도의 농도수준이 되도록 표준물질을 동일하게 첨가한 시료(농약류의 경우 0.01과 0.05 µg/L, 합성원료의 경우 0.1과 0.5 µg/L)를 5개씩 준비하고 시험 분석절차와 동일하게 진행하였다. 그 다음 시험 분석절차를 거치지 않고 같은 표준물질의 양을 첨가한 다음 그것만 직접 분석한 것과 비교하여 백분율을 구하였다. 농약류의 경우, 57.8~98.6%의 평균 회수율과 3.8~19.6의 상대표준편차(%)를 보였다. 합성원료의 경우는 56.2~76.7%의 평균 회수율과 3.9~15.9의 상대표준편차(%)를 보였다(Table 5).
Table 5.Recovery test results for the analysis of seven pesticides and two synthetic materials in surface water (n=5)
하천수 분석결과
하천수 총 10지점에 대한 시료를 채취하여 SPE-LC-MS/MS로 분석한 결과, 농약류인 알디캅과 카바릴은 전 지점에서 모두 검출되지 않았다. 농약류에서 메소밀의 경우, 0.012 μg/L(평균 0.003 μg/L) 농도로 검출되었고, 검출빈도는 10%였으며, molinate의 경우, 0.008~0.009 μg/L(평균 0.005 μg/L) 농도범위로 검출되었고, 검출빈도는 20%였다. 그리고 MCPA의 경우, 0.007~0.073 μg/L(평균 0.029 μg/L) 농도범위로 검출되었고, 검출빈도는 90%였으며, 2,4-D의 경우, 0.010~0.101 μg/L(평균 0.031 μg/L) 농도범위로 검출되었고, 검출빈도는 70%였다. 또한 카보퓨란의 경우는 0.023~0.465 μg/L(평균 0.191 μg/L) 농도범위로 검출되었고, 검출빈도는 90%였다. 합성원료에서 quinoline의 경우, 0.013 μg/L(평균 0.006 μg/L) 농도범위로 검출되었고, 검출빈도는 10%였으며, 비스페놀-A의 경우, 0.013~0.189 μg/L(평균 0.050 μg/L) 농도범위로 검출되었고, 검출빈도는 90%였다. 표준물질을 첨가한 후와 10개의 채취지점 중 한 지점에서의 검출된 분석항목에 대한 크로마토그램을 Fig. 1과 Fig. 2에 나타내었다.
Figure 1.LC-MS/MS chromatogram of the extract from surface water sample spiked in the concentration of 100 μg/L (aldicarb: 7.782 min, methomyl: 5.636 min, molinate: 10.167 min, MCPA: 8.997 min, carbaryl: 8.821 min, 2,4-D: 8.972 min, carbofuran: 8.613 min, quinoline: 7.735 min, bisphenol A: 8.971 min).
Figure 2.LC-MS/MS chromatograms from Gim-Je (sampling site) quantified in concentration of 0.012 μg/L (methomyl), 0.035 μg/L (MCPA), 0.015 μg/L (2,4-D), 0.100 μg/L (carbofuran) and 0.189 μg/L (bisphenol A).
결 론
본 연구에서는 하천수 내 농약류 및 합성원료를 극미량까지 분석할 수 있는 고감도 SPE-LC-MS/MS 동시 분석법을 개발하였다. 이 분석법으로 농약류인 카보퓨란의 경우 정량한계 0.002 μg/L까지 달성하였으며, 정량한계 수준에서의 각 분석물질들의 피크의 S/N비가 10배 이상으로 명확하였다. 높은 농축계수을 가지며 방해물질의 효과적인 제거가 용이한 SPE 전처리를 통한 정도관리로 신뢰성을 확인하였으며, 정량에 활용되는 검정곡선은 실제 하천수를 이용하였기 때문에 매질의 영향을 반영한 분석결과를 얻을 수 있게 되었다. 하천수 10개 지점의 분석결과, 농약류인 카보퓨란, MCPA, 2,4-D, molinate, 메소밀 등 5개 항목과 합성원료인 quinoline과 비스페놀-A 등 2개 항목을 검출하였다. 본 연구의 분석법은 하천수 모니터링에서 신뢰도 높은 분석방법으로 사용할 수 있을 것으로 판단된다.
References
- Korean Crop Protection Association. Pesticide Guide; Korean Crop Protection Association: Seoul, Korea, 2012.
- Korean Crop Protection Association. Annual Report of the Pesticides; Korean Crop Protection Association: Seoul, Korea, 2012.
- Kwon, Y. J.; Kang, T. S.; Kim, K. R.; Lee, K. S.; Ju, Y. S.; Song, J. C. Korean J. of Rural Med. 2004, 29(2), 265.
-
US EPA, Reregistration Eligiblility Decision for 2,4-D 2005(
http://www.epa.gov/pesticides/reregistration/ REDs/ 24d_red.pdf ). -
IARC, Monographs Supplement 7 1987 (
http://monographs.iarc.fr/ENG/Monographs/suppl7/Suppl7-48.pdf ). -
US EPA, Reregistration Eligiblility Decision for Methomyl 1998 (
http://www.epa.gov/pesticides/reregistration/REDs/0028red.pdf ). -
US EPA, Reregistration Eligiblility Decision for Aldicarb 2007 (
http://www.epa.gov/pesticides/reregistration/REDs/aldicarb_red.pdf ). -
US EPA, Illegal Pesticide Products 2012 (
http://www.epa.gov/pesticides/health/illegalproducts/#pet ). -
USDA, EXTOXNET (Extension Toxicology Network) 1993 (
http://pmep.cce.cornell.edu/profiles/extoxnet/haloxyfopmethylparathion/mcpa-ext.html#6 ). -
USDA, EXTOXNET (Extension Toxicology Network) 1996 (
http://extoxnet.orst.edu/pips/molinate.htm ). -
NPIC, Carbaryl(General Fact Sheet) 2003 (
http://npic.orst.edu/factsheets/carbgen.pdf ). -
US EPA, Quinoline-Hazard summary 2000 (
http://www.epa.gov/airtoxics/hlthef/quinolin.html#ref4 ). - Korea Food and Drug Administration. Hazardous material series; Korea Food and Drug Administration: Osong, Korea, 2010.
- Giordano, A.; Richter, P.; Ahumada, I. Talanta 2011, 85, 2425. https://doi.org/10.1016/j.talanta.2011.07.087
- Bonansea, R. I.; Ame, M. V.; Wunderlin, D. A. Chemosphere 2013, 90, 1860. https://doi.org/10.1016/j.chemosphere.2012.10.007
- Pinheiro, A. S.; Rocha, G. O.; Andrade, J. B. Microchem. J. 2011, 99, 303. https://doi.org/10.1016/j.microc.2011.05.019
- Tankiewicz, M.; Morrison b, C.; Biziuk, M. Talanta 2013, 107, 1. https://doi.org/10.1016/j.talanta.2012.12.052
- Quintana, J.; Martí, I.; Ventura, F. J. Chromatogr. A 2001, 938, 3. https://doi.org/10.1016/S0021-9673(01)01168-2
- Filho, A. M.; Santos, F. N.; Pereira, P. A. P. Microchem. J. 2010, 96, 139. https://doi.org/10.1016/j.microc.2010.02.018
- Shi, Z.; Hu, J.; Li, Q.; Zhang, S.; Liang, Y.; Zhang, H. J. Chromatogr. A 2014, 1355, 219. https://doi.org/10.1016/j.chroma.2014.05.085
- Xu, C.; Armstrong, D. W. Analytica Chimica Acta 2013, 792, 1. https://doi.org/10.1016/j.aca.2013.05.054
- Reemtsma, T.; Alder, L.; Banasiak, U. J. Chromatogr. A 2013, 1271, 95. https://doi.org/10.1016/j.chroma.2012.11.023
- Corcia, A. D.; Nazzari, M., Rao, R.; Samperi, R.; Sebastiani, E. J. Chromatogr. A 2000, 878, 87. https://doi.org/10.1016/S0021-9673(00)00263-6
- Gervais, G.; Brosillon, S.; Laplanche, A.; Helen, C. J. Chromatogr. A 2008, 1202, 163. https://doi.org/10.1016/j.chroma.2008.07.006
- Assoumani, A.; Margoum, C.; Chataing, S.; Guillemain, C.; Coquery, M. J. Chromatogr. A 2014, 1333, 1. https://doi.org/10.1016/j.chroma.2014.01.063
- Dujaković, N.; Grujić, S.; Radišić, M.; Vasiljević, T.; Laušević, M. Analytica Chimica Acta 2010, 678, 63. https://doi.org/10.1016/j.aca.2010.08.016
- Farajzadeh, M. A.; Khoshmaram, L. Food. Chem. 2013, 141, 1881. https://doi.org/10.1016/j.foodchem.2013.05.088
- Xiong, J.; Hu, B. J. Chromatogr. A 2008, 1193, 7. https://doi.org/10.1016/j.chroma.2008.03.072
- Su, Y. S.; Jen, J. F. J. Chromatogr. A 2010, 1217, 5043. https://doi.org/10.1016/j.chroma.2010.06.006
- Boussahel, R.; Bouland, S.; Moussaoui, K. M.; Baudu, M.; Montiel, A. Wat. Res. 2002, 36, 1909. https://doi.org/10.1016/S0043-1354(01)00372-4
- Dong, C.; Zeng, Z.; Yang, M. Wat. Res. 2005, 39, 4204. https://doi.org/10.1016/j.watres.2005.08.004
- Llasera, M. P. G.; Lez, M. B. Wat. Res. 2001, 35(8), 1933. https://doi.org/10.1016/S0043-1354(00)00478-4
- Subhani, Q.; Huang, Z. P.; Zhu, Z. Y.; Liu, L. Y.; Zhu, Y. Chin. Chem. Lett. 2014, 25, 415. https://doi.org/10.1016/j.cclet.2013.11.014
-
US EPA, Method 1699 2007 (
http://water.epa.gov/scitech/methods/cwa/bioindicators/upload/2008_01_03_methods_-method_1699.pdf ). -
US EPA, Method 508.1 1995 (
http://water.epa.gov/scitech/methods/cwa/bioindicators/upload/2007_11_06_methods_method_508_1-2.pdf ). -
US EPA, Method 507 1989 (
http://water.epa.gov/scitech/methods/cwa/bioindicators/upload/2007_11_06_methods_method_507.pdf ).
Cited by
- Monitoring of Pesticides in the Yeongsan and Seomjin River Basin vol.34, pp.4, 2015, https://doi.org/10.5338/KJEA.2015.34.4.37
- Simultaneous Analysis of 13 Pesticides in Groundwater and Evaluation of its Persistent Characteristics vol.50, pp.5, 2015, https://doi.org/10.7745/kjssf.2017.50.5.434