마그마 해양의 결정층 침식·침전 메커니즘 연구

초록

본 논문은 실험적 아날로그 모델을 이용해 대류 마그마 해양 내에서 부유·침강하는 결정(플라스틱 비드)의 안정성, 침식 속도, 침전 속도를 규명한다. 침식 법칙은 레일리‑로버츠 수와 셸즈 수 두 차원수만으로 기술되며, 침전 속도는 스톡스 속도와 선형적으로 스케일링한다. 얻어진 스케일링 관계는 차후 행성 미행성체의 마그마 해양 진화 모델에 적용될 예정이다.

상세 분석

이 연구는 두 가지 핵심 차원수, 즉 내부 가열에 의한 대류 강도를 나타내는 레일리‑로버츠 수(Ra_H)와 입자‑유체 상호작용을 정량화하는 셸즈 수(Θ) 를 중심으로 실험적 스케일링을 구축한다. 실험 장치는 마이크로파 흡수에 의해 부피 가열되는 고점도 유체(글리세롤‑에틸렌글리콜 혼합물)와 플라스틱 비드(밀도 차이로 부양·침강 전이 가능)를 사용한다. 비드의 부양·침강 전이 온도(T_inv≈37.4 °C)를 기준으로 부유층(lid)과 침전층(cumulate)의 형성을 제어한다.

침식 메커니즘은 초기 대류가 급속히 발생한 뒤, 부유층 바닥에 베드로드(bed‑load) 형태의 입자층이 형성되고, 이 층 위에서 사구( dune) 형태의 돌출이 나타나는 과정으로 관찰된다. 사구가 성장하면 전단 응력이 증가해 입자 재현탁이 촉발되고, 부유층 두께가 점진적으로 감소한다. 실험 결과는 침식 속도가 Θ_c · Ra_H^−1/2 형태의 함수로 표현될 수 있음을 보여준다(Θ_c는 임계 셸즈 수).

침전 측면에서는 비드가 액체보다 무거워지는 온도 구간에서 바닥에 침전이 시작된다. 침전 속도는 입자의 스톡스 낙하 속도(v_s)와 거의 일치하며, 대류 흐름이 침전 속도에 미치는 영향은 1차 스케일링에서 무시할 수 있을 정도로 작다. 이는 기존 연구(예: Solomatov & Stevenson 1993)의 결과와 일치하지만, 본 실험은 대류와 침전이 동시에 진행되는 비정상 상태(transient)에서도 동일한 스케일링이 유지됨을 실증한다.

열전달 모델은 내부 가열에 의한 온도 스케일 ΔT_H = H h²/λ_f와 고 Prandtl 수(P_r≫1) 가정 하에, 상부 열경계층(TBL)의 두께와 온도 강하가 Ra_H^−1/4 로 스케일링되는 식을 도입한다. 부유층 두께 δ(t)와 기저 온도 T_lid 사이의 열플럭스 Q_s,lid는 대류층에서 방출되는 열플럭스 Q_s,conv와 에너지 보존식(ρ c_p ∂T_bulk/∂t = H – Q_s,conv)으로 연결된다. 실험 데이터에 기반한 상수 C_T≈3.59를 사용해 Q_s,conv ∝ (T_bulk – T_lid)^{4/3}·Ra_H^{−1/3} 형태의 스케일링을 도출하였다.

결과적으로, 부유층의 침식·소멸은 셸즈 수와 레일리‑로버츠 수의 조합으로 예측 가능하고, 침전은 스톡스 속도에 의해 지배된다. 이러한 두 스케일링 법칙은 행성 규모의 마그마 해양에서 결정층(플라그이오클라스·올리빈 등)의 부양·침강을 예측하는 데 직접 적용될 수 있다.