지중해 심층 엽록소 최대층의 시공간 변동과 미소조류 동역학 확률 모델

초록

본 연구는 지중해 남부의 저혼합 해역에서 미소조류(피코식물성 플랑크톤)의 분포를 설명하기 위해, 빛·영양소 경쟁과 환경 변동을 포함한 확률적 반응‑확산‑탐색 모델을 구축하였다. 모델 시뮬레이션 결과는 시실리 해협 두 관측지의 실측 데이터와 높은 일치성을 보였으며, 향후 올리고트로픽 해역의 식물성 플랑크톤 예측에 활용 가능함을 제시한다.

상세 요약

이 논문은 지중해의 특수한 물리‑생물학적 환경, 즉 표층에서 강한 혼합이 일어나지 않아 심층에 영양소와 빛이 제한되는 상황을 배경으로 한다. 저자는 피코식물성 플랑크톤(pico‑phytoplankton)의 성장·사멸을 기술하기 위해, 반응‑확산‑탐색(Reaction‑Diffusion‑Taxis) 방정식에 두 가지 핵심 요소를 도입하였다. 첫째, 빛과 영양소에 대한 동시 경쟁을 나타내는 인트라스페시픽(intraspecific) 경쟁 항을 비선형 성장 함수 형태로 구현했으며, 이는 광합성 효율과 영양소 흡수율을 각각 깊이에 따라 감소하는 함수로 매개한다. 둘째, 환경 변수(예: 일사량, 온도, 영양염 농도)의 무작위 변동을 반영하기 위해 곱셈형 노이즈(multiplicative noise)를 포함시켰다. 이 노이즈는 확률 미분 방정식 형태로 삽입되어, 평균적인 환경 조건 위에 순간적인 변동을 겹쳐 놓는다. 모델은 1‑차원 수직 좌표(z)를 기준으로, 확산 계수 D와 탐색 계수 χ를 통해 물질의 수평 혼합과 빛‑영양소 구배에 따른 이동을 동시에 고려한다. 수치 해석은 안정적인 Crank‑Nicolson 스킴과 Stochastic Euler‑Maruyama 방법을 결합해 수행했으며, 파라미터는 현장 측정값(광 투과도, 영양소 농도, 온도 프로파일)과 문헌값을 기반으로 보정하였다. 시실리 해협의 두 관측지(동부와 서부)에서 얻은 심층 엽록소 최대층(DCM) 깊이와 피코식물성 플랑크톤 농도 프로파일을 모델 출력과 비교했을 때, 평균 오차가 5 % 이하로 매우 낮았다. 특히, 노이즈 강도가 증가할수록 DCM 위치가 변동성을 보이며, 실제 관측에서 나타나는 연간·계절적 변동을 재현하는 데 중요한 역할을 함을 확인했다. 이러한 결과는 전통적인 결정론적 모델이 놓치기 쉬운 환경 변동성의 영향을 정량화하고, 저영양 해역에서 피코식물성 플랑크톤이 어떻게 공간적으로 집중되는지를 설명한다는 점에서 학문적·실용적 의의를 가진다.