주거지 토양 납 농도에 대한 페인트·휘발유 기여도 비선형 회귀 모델링

초록

본 연구는 1987년 미네소타 토양 납 조사와 1990년 전국 토양 납 조사 데이터를 활용해, 납 기반 페인트와 무연 휘발유가 주거지 토양 납 농도에 미치는 상대적 영향을 시간에 따라 추정하는 새로운 비선형 회귀 모델을 제시한다. 모델은 각 원천별 기여도를 정량화하고, 지역별·시대별 납 오염 패턴을 해석함으로써 환경 정책 및 잠재적 소송에 활용 가능한 과학적 근거를 제공한다.

상세 분석

이 논문은 기존의 선형 회귀 접근법이 납 오염 원천을 시간에 따라 구분하는 데 한계가 있다는 점을 지적하고, 비선형 회귀 모델을 도입함으로써 두 주요 원천(납 기반 페인트와 무연 휘발유)의 기여도를 동시 추정한다. 모델은 토양 납 농도를 (C(t)=\alpha + \beta_1 P(t)^{\gamma_1} + \beta_2 G(t)^{\gamma_2}) 형태로 설정한다. 여기서 (P(t))는 해당 연도까지 사용된 페인트량(연도별 주택 건축·리모델링 데이터 기반), (G(t))는 연도별 휘발유 소비량(국가 통계)이며, (\gamma_1,\gamma_2)는 각각의 비선형 효과를 나타내는 지수이다.

모델 파라미터는 비선형 최소제곱법(Levenberg‑Marquardt 알고리즘)으로 추정했으며, 잔차 분석을 통해 이분산성 및 자기상관을 검증하였다. 결과는 두 데이터셋 모두에서 (\gamma_1)과 (\gamma_2)가 1보다 작아, 초기 단계에서 급격히 증가하던 납 농도가 시간이 지남에 따라 포화 현상을 보인다는 것을 시사한다. 특히 1970년대 이전에 건축된 주택이 밀집된 지역에서는 페인트 기여도가 60~70%에 달했으며, 1980년대 이후 차량 운행이 급증한 교외 지역에서는 휘발유 기여도가 45%까지 상승한다.

민감도 분석에서는 (\beta) 계수의 변동이 전체 모델 예측에 큰 영향을 미치지 않으며, (\gamma) 값이 변할 경우 기여도 비율이 급격히 바뀌는 것을 확인했다. 이는 비선형 지수가 원천별 시간 가중치를 결정하는 핵심 변수임을 의미한다. 또한, 교차 검증을 통해 모델의 일반화 능력이 우수함을 입증했으며, 기존 선형 모델 대비 평균 절대 오차(MAE)가 15% 감소하였다.

법적 측면에서는 모델이 제공하는 “시간‑공간적 기여도” 추정치를 근거로, 특정 토양 오염 사건에 대한 책임소재를 보다 정밀하게 가려낼 수 있다. 예를 들어, 1995년 건축된 주택 인근 토양에서 높은 납 농도가 발견될 경우, 모델은 해당 시점 이전에 축적된 페인트 기여도가 30% 이하임을 보여주어, 휘발유 배출이 주요 원인임을 입증할 수 있다.

이러한 접근법은 정책 입안자에게도 유용하다. 지방 자치단체는 모델 결과를 토대로 오래된 주택이 밀집된 구역에 페인트 제거·재도장 프로그램을 우선 배치하고, 교외 고속도로 주변에서는 차량 배출 저감 대책을 강화하는 등 자원 배분을 최적화할 수 있다.