Rainfall-Driven Variability of Physicochemical Water Quality Parameters in the Meghna River, Bangladesh: A Preliminary Monthly Assessment Using Correlation, Regression, and Principal Component Analysis

문제 제기

방글라데시와 같은 몬순 중심의 열대 지역에서 하천의 수질 화학적 특성은 강우에 의해 상당한 계절적 변동을 보입니다. 이러한 변동은 두 가지 주요 메커니즘에 의해 특징지어집니다: 희석(총 용존 고형물(TDS)과 같은 보수적 물질의 농도를 감소시킴) 및 유출/이동(지표 유출수를 통해 화학적 산소 요구량(COD) 및 생화학적 산소 요구량(BOD)과 같은 유기물 부하를 증가시킴). 방글라데시에서 이러한 계절적 변동에 대한 질적 증거는 존재하지만, 메그나 강(Meghna River)의 강우와 수질 매개변수 간의 상관관계를 구체적으로 다룬 정량적 연구는 제한적입니다. 또한, 이 특정 유역에 대해 방글라데시 환경부(DoE)의 수질 기록과 위성 유도 강우 데이터(NASA POWER)를 통합하는 연구도 부족한 실정입니다. 기존 모니터링의 결정적인 제약 요인은 불충분한 데이터 커버리지입니다. 2023년의 경우, 완전한 월간 수질 데이터가 6개월치만 확보되었으며, 이는 통계적 한계를 인정하면서도 기초 자료(baseline)를 구축하고자 하는 연구 설계의 필요성을 시사합니다.

연구 방법론

본 연구는 2023년 메그나 가트(Meghna Ghat) 관측소에서 수행된 예비 월간 평가를 실시하였습니다. 분석에는 강우와 수질 데이터가 모두 확보된 6개월(1월, 3월, 7월, 9월, 10월, 11월)의 데이터셋을 활용하였습니다.

• 데이터 출처:
  • 수질: pH, 용존 산소(DO), BOD, 총 용존 고 sólidos(TDS), COD에 대한 월간 데이터를 방글라데시 DoE로부터 확보하였습니다.
  • 강우: 남아시아 유역에서 검증된 GPM IMERG 제품을 활용하여, NASA Prediction Of Worldwide Energy Resources (POWER) 데이터베이스로부터 월간 강우 데이터를 추출하였습니다.
• 분석 기법:
  1. 기술 통계: 모든 매개변수의 평균, 표준 편차, 범위를 계산하였습니다.
  2. 피어슨 상관계수 및 선형 회귀: 강우와 각 수질 매개변수 간의 쌍별 관계를 분석하였습니다. 단순 선형 회귀를 사용하여 결정 계수($R^2$)를 결정하였습니다.
  3. 종합 수질 지수(WQI): 가중치를 두지 않은 등가중치 WQI를 구성하였습니다. 각 매개변수의 하위 지수($Q_i$)는 기준값(예: BOD 6 mg/L, TDS 1000 mg/L, COD 40 mg/L)을 바탕으로 계산되었으며 백분율로 스케일링되었습니다. 최종 WQI는 5개 하위 지수의 평균이며, 값이 높을수록 수질이 나쁨을 의미합니다.
  4. 주성분 분석(PCA): 6개 변수(강우 + 5가지 수질 매개변수) 간의 결합 변동의 지배적인 모드를 식별하기 위해 6개월간의 표준화된 데이터(z-score)를 대상으로 PCA를 수행하였습니다.

주요 결과

분석 결과, 작은 표본 크기($n=6$)로 인해 강력한 통계적 추론은 어려웠으나, 강우가 서로 다른 물리화학적 매개변수에 미치는 영향에 대한 뚜렷한 패턴이 드러났습니다.

• 강우와의 상관관계:
  • COD: 강우와 매우 강한 양의 상관관계($r = 0.99, p < 0.001$)를 보였습니다. 회귀 모델은 강우가 COD 변동의 97.2%를 설명함을 나타냈으며($R^2 = 0.972$), 이는 강우 시 유출수가 유기 오염물질을 이동시키는 유출 효과(wash-off effect)를 시사합니다.
  • BOD 및 TDS: 강우와 중등도의 음의 상관관계(BOD $r = -0.71$, TDS $r = -0.62$)를 보였습니다. 이는 유량이 증가함에 따라 용존 고형물과 유기물 부하의 농도가 낮아지는 희석 효과와 일치합니다. 그러나 표본 크기로 인해 이러한 상관관계는 $p < 0.05$ 수준에서 통계적으로 유의미하지 않았습니다.
  • pH 및 DO: 강우와 약하고 유의미하지 않은 상관관계($r = -0.15$ 및 $r = -0.12$)를 보였습니다.
• 수질 지수(WQI):
  • 종합 WQI는 27.58(7월, 몬순)에서 51.11(3월, 건기) 사이의 범위를 나타내어, 우기에 수질이 더 좋음을 보여주었습니다.
  • 강우와 WQI 사이의 상관관계는 중등도의 음의 상관관계($r = -0. -46$)였으나 통계적으로 유의미하지 않았습니다 ($p = 0.36$).
• 주성분 분석(PCA):
  • 처음 두 개의 주성분이 전체 분산의 **85.3%**를 설명하였습니다.
  • PC1 (분산 58.9%): 대조적인 로딩 패턴을 보였습니다. 강우와 COD는 음의 방향으로 로딩된 반면, BOD와 TDS는 양의 방향으로 로딩되었습니다. 이는 상반된 계절적 행동을 확인시켜 줍니다: 즉, 강우는 COD를 높이는 반면(유출), BOD와 TDS는 낮춥니다(희석).
  • PC2 (분산 26.4%): pH에 의해 지배되었으며, 이는 pH가 강우 기반의 경사도와 비교적 독립적으로 변함을 시사합니다.

주요 기여

1. 데이터 통합: 본 연구는 지상 기상 관측망이 희박한 지역을 위한 실용적인 방법론으로서, DoE 수질 기록과 NASA POWER 위성 강우 데이터를 메그나 강에 성공적으로 통합하였습니다.
2. 정량적 기초 자료 제공: 2023년 메그나 가트 관측점에서 월간 강우와 특정 물리화학적 매개변수(pH, DO, BOD, TDS, COD) 간의 관계에 대한 최초의 정량적 평가를 제공하였습니다.
3. 방법론적 투명성: 본 논문은 적은 표본 크기($n=6$)를 사용하는 데 따른 통계적 추론의 한계를 명확히 밝히고 있습니다. 저자들은 PCA와 상관관계 분석을 확정적인 증명이 아닌, 다변량 구조를 시각화하고 향후 연구를 위한 가설을 생성하기 위한 탐색적 도구로 활용하였습니다.
4. 계절적 특성화: 결과는 "혼합 희석-유출(mixed dilution–wash-off)" 가설을 뒷받침하며, 강우가 보수적 물질(TDS)은 희석하는 동시에 유기 오염물질(COD)은 이동시킨다는 점을 보여줍니다.

의의 및 주장

저자들은 본 연구를 단기적인 강우-수질 역학에 대한 결정적인 규명이 아닌, 엄격하게 예비 탐색적 연구이자 기초 자료로 규정합니다.

• 신중한 주장: 저자들은 강한 COD 상관관계($r=0.99$)를 더 광범위하고 지속적인 역사적 데이터에 의한 확인이 필요한 가설로 취급해야 한다고 명시했습니다. 다른 상관관계에서 통계적 유의성이 나타나지 않은 것은 작은 표본 크기 때문이라고 설명했습니다.
• 과학적 가치: 주요 의의는 데이터가 부족한 지역에서 위성 유도 강우 데이터를 사용하여 수질 추세를 분석하는 것의 타당성을 입증하고, 보다 빈번한 이벤트 기반 모니터링의 필요성을 강조하는 데 있습니다.
• 향후 방향: 본 연구는 향-후 연구가 월간 집계에 의존하기보다는, 여러 해에 걸친 월간 모니터링을 완료하고 수문학적 반응을 더 잘 포착할 수 있는 이벤트 기반 샘플링 설계(사전, 도중, 사후 샘플링)를 채택해야 한다고 결론짓습니다.

요약하자면, 본 논문은 데이터 부족에도 불구하고 엄격하게 결과를 해석하며, COD에 대한 잠재적인 유출 메커니즘과 TDS/BOD에 대한 희석 메커니즘을 식별하는 투명하고 데이터 중심적인 예비 평가를 제공합니다.