비와 염분의 무작위 상호작용 토양 염류 장기 예측 모델

초록

이 논문은 건조·습윤 염분 침적과 강우에 의한 간헐적 세척을 고려한 최소주의 확률 모델을 제시한다. 다중 포아송 잡음으로 구동되는 단일 확률 미분 방정식으로 토양 염질량과 농도의 장기 확률밀도함수를 analytically 도출하고, 강우 빈도에 따라 두 가지 상이한 염류 거동 구역이 존재함을 밝힌다.

상세 분석

본 연구는 토양 수분과 염분 양을 동시에 기술하는 전통적인 2차 방정식 체계를 하나의 확률 미분 방정식(SDE)으로 축소함으로써 수학적·계산적 복잡성을 크게 낮췄다. 핵심은 강우 사건을 ‘점프’ 형태의 포아송 프로세스로 모델링하고, 이 점프가 염분 농도에 곱셈적 영향을 미친다는 가정이다. 따라서 SDE는 dM = (a − b M)dt + c M dN(t) 형태를 띠며, 여기서 M은 토양 내 총 염질량, a는 건조·습윤 염분 침적 속도, b는 자연 증발·식물 흡수에 의한 감소율, c는 강우에 의한 세척 효율, N(t)는 평균 λ인 포아송 과정이다. 이 식을 Itô 해석이 아닌 Stratonovich 해석으로 풀어, 장기 정상 상태 확률밀도함수(PDF)를 정확히 구한다.

분석 결과는 두 개의 파라미터 구역을 제시한다. 첫 번째 구역에서는 λ가 증가함에 따라 평균 염질량 ⟨M⟩이 거의 일정하거나 감소한다. 이는 강우가 충분히 자주 발생해 염분이 효과적으로 씻겨 나가므로 토양 염류가 축적되지 않음을 의미한다. 두 번째 구역에서는 λ가 일정 수준을 초과하면 ⟨M⟩이 급격히 상승한다. 여기서는 강우가 너무 자주 발생해 토양 수분이 지속적으로 포화 상태에 머무르면서 염분 침적이 지속적으로 이루어지고, 세척 효율(c)이 상대적으로 낮아 누적 효과가 지배한다. 이러한 전이점은 건조 기후에서 강우 빈도가 소폭 감소할 경우 토양 염류가 급격히 악화될 수 있음을 시사한다.

또한 파라미터 민감도 분석을 통해 식생(뿌리 흡수율), 토양 투수성, 염분 침적 속도 등과 같은 환경·생물학적 요인이 장기 염류 거동에 미치는 상대적 중요성을 정량화하였다. 특히 λ와 c의 비율이 전이 구역을 결정하는 핵심 지표임이 밝혀졌다.

이 모델은 기존의 결정론적 물리‑화학 모델이 놓치기 쉬운 ‘우연성’과 ‘점프 사건’의 영향을 명시적으로 포함함으로써, 기후 변화 시나리오 하에서 비·염분 상호작용을 보다 현실적으로 예측할 수 있는 기반을 제공한다.