마리아나 해구 환경 요인 변동성 분석

초록

본 연구는 R 프로그래밍과 GIS를 활용해 마리아나 해구의 지형 변화를 좌우하는 주요 환경 요인들을 정량적으로 평가한다. 침전물 두께, 경사각, 깊이, 그리고 네 개의 판 구조(필리핀, 태평양, 마리아나, 캐롤라인) 간의 상호작용을 변수로 설정하고, 요인분석(Factor Analysis)을 수행해 핵심 인자를 도출하였다. 결과는 침전물 두께와 경사각이 깊이 변동에 가장 큰 영향을 미치며, 판 경계와의 공간적 연관성이 뚜렷함을 보여준다.

상세 분석

이 논문은 심해 지형학 연구에 통계적 모델링과 공간분석을 결합한 점이 가장 큰 특징이다. 먼저, 연구자는 마리아나 해구를 가로지르는 4개의 판(필리핀, 태평양, 마리아나, 캐롤라인)의 경계와 위치 정보를 GIS로 디지털화하고, 각 판별로 대표적인 관측 지점을 선정하였다. 이 지점들에서 얻은 심도, 좌표, 경사각, 침전물 두께 등 12여 개의 환경 변수를 데이터베이스화한 뒤, R의 ‘psych’, ‘FactoMineR’, ‘ggplot2’ 등 패키지를 이용해 전처리와 정규화를 수행하였다.

전처리 단계에서는 결측값을 다중대체법(MICE)으로 보완하고, 변수 간 스케일 차이를 최소화하기 위해 Z-스코어 변환을 적용하였다. 이후 상관행렬을 계산해 변수 간 다중공선성을 확인했으며, VIF(Variance Inflation Factor) 값이 5 이하인 변수만을 최종 요인분석에 투입하였다. 요인분석은 주성분 분석(PCA) 기반의 탐색적 요인분석(EFA)으로 진행되었으며, 스크리 플롯과 고유값(eigenvalue) >1 기준으로 3개의 주요 요인을 추출하였다.

첫 번째 요인은 ‘심도·경사 복합 요인’으로, 깊이와 경사각이 높은 양의 부하값을 보이며, 이는 해구의 급격한 절벽 구간과 직접 연관된다. 두 번째 요인은 ‘침전물 두께 요인’으로, 침전물 두께와 지질학적 층위가 높은 부하를 나타내어, 퇴적물 압축과 물리적 강도가 해구 형성에 미치는 영향을 강조한다. 세 번째 요인은 ‘판 경계·위치 요인’으로, 관측 지점의 위도·경도와 판 경계와의 거리 변수가 높은 부하를 보여, 판 구조적 스트레스가 지역별 지형 변동을 조절한다는 결론을 뒷받침한다.

요인 점수는 각 관측 지점에 대해 계산되어 군집분석(k-means)과 결합했으며, 결과는 판 경계 근처가 높은 ‘심도·경사 복합 요인’ 점수를, 판 내부는 상대적으로 낮은 점수를 나타냈다. 또한, Euler‑Venn 다이어그램을 활용해 네 개 판이 공유하는 환경 변수와 고유 변수를 시각화함으로써, 판 간 상호작용이 복합적인 환경 요인 네트워크를 형성한다는 점을 명확히 했다.

통계적 검증 측면에서는 요인 적재량의 신뢰성을 확인하기 위해 부트스트랩(1000회) 재표본추출을 수행했으며, 모든 요인에 대해 95% 신뢰구간이 양의 부하값을 포함함을 확인했다. 또한, 요인 회전(Varimax) 후 각 요인의 설명력은 총 변동성의 68%를 차지했으며, 이는 해구 지형 변동을 설명하는 데 충분한 모델 적합도를 의미한다.

결과적으로, 본 연구는 R 기반 요인분석이 심해 지형학 데이터의 복합성을 해소하고, 주요 환경 요인을 정량적으로 도출하는 데 유효함을 입증한다. 특히, 침전물 두께와 경사각이 깊이 변동을 주도한다는 점, 그리고 판 경계와의 공간적 연관성이 지형 변동에 부수적인 영향을 미친다는 점을 과학적으로 뒷받침한다. 이러한 통합적 접근은 향후 심해 지형 변화 예측 모델링과 해양 자원 관리에 중요한 기초 자료가 될 것이다.