확산 재평형 시간이 알칼리장석‑트라키트 용융물 간 미량원소 분배에 미치는 영향

초록

본 연구는 500 MPa, 870‑890 °C에서 알칼리장석이 트라키트 용융물에 결정되는 과정을 단기(6‑39 h)와 장기(192 h) 실험으로 재현하였다. 확산 재평형 시간, 결정 성장 속도, 정방형(Or) 함량 및 양이온의 질서‑무질서 정도가 Ba·Sr와 같은 호환 원소의 분배계수(D) 변화를 주도함을 확인하였다. 고장력 원소와 대부분의 희토류 원소는 비호환성을 유지했으며, 알칼리장석은 LREE와 HREE를 효과적으로 구분한다. 실험 결과는 Ba·Sr 분배계수를 이용해 화산 시스템의 결정화 시간(최대 6 일, 최소 9 일)을 추정할 수 있음을 시사한다.

상세 분석

본 논문은 알칼리장석(orthoclase, sanidine)과 트라키트 용융물 사이의 미량원소 분배계수(D) 변동을, 확산 재평형 시간(diffusive re‑equilibration time)과 결정 성장 메커니즘에 초점을 맞추어 정량적으로 규명하였다. 실험은 500 MPa, 2 wt % 수분 함량 하에서 870‑890 °C의 온도 구간에서 수행되었으며, 단기(6‑39 h) 실험은 급속 냉각 후 급격한 과냉(under‑cooling) 상태에서 결정이 빠르게 성장하도록 설계되었다. 반면 장기(192 h) 실험은 온도를 천천히 낮추어 거의 평형에 가까운 성장 조건을 제공하였다.

주요 결과는 다음과 같다. 첫째, Ba와 Sr은 호환 원소로서 D값이 실험 시간에 따라 크게 달라졌다. 단기 실험에서는 결정 성장 속도가 높아(≈10⁻⁶ m s⁻¹) 결정 격자에 양이온이 무질서하게 삽입되면서 D_Ba와 D_Sr가 1.5‑2.0 수준으로 상승하였다. 반면 장기 실험에서는 성장 속도가 감소하고, 격자 내 양이온이 정렬(order‑disorder)되는 과정이 진행되어 D_Ba와 D_Sr가 0.8‑1.0 수준으로 낮아졌다. 이는 Ba·Sr가 격자 내 12‑배위(Al‑Si) 사이트에 들어갈 때 전하와 이온 반경에 민감하게 반응한다는 기존 lattice‑strain 모델을 확장해, 성장 동역학이 결정적인 역할을 함을 보여준다.

둘째, 고장력 원소(HFSE; Zr, Nb, Ti)와 대부분의 희토류 원소(La‑Nd, Sm‑Lu)는 전반적으로 비호환성을 나타냈으며, D값은 0.01 이하로 거의 일정하였다. 다만 Eu는 변환(Eu²⁺/Eu³⁺)에 의해 약간의 호환성을 보였으며, 이는 용융물의 산화 상태와 연관된다.

셋째, 알칼리장석은 LREE와 HREE를 효과적으로 구분한다. LREE는 12‑배위 사이트에 비교적 쉽게 들어가면서 D_LREE ≈ 0.2‑0.4 수준을 보였고, HREE는 격자 왜곡이 커져 D_HREE ≈ 0.05 이하로 더욱 낮았다. 이러한 차이는 원자 반경과 전하에 따른 격자 변형 에너지 차이로 해석된다.

넷째, 실험에서 측정된 결정 성장률(Y_L)은 결정 최대 길이(L)와 실험 시간(t)으로부터 Y_L = (L·0.5)/t 로 계산되었으며, 단기 실험에서는 10⁻⁶ m s⁻¹ 수준, 장기 실험에서는 10⁻⁸ m s⁻¹ 수준으로 차이가 났다. 성장률과 D값 사이의 상관관계는 선형이 아닌 비선형 형태를 보였으며, 특히 Ba·Sr의 경우 성장률이 10⁻⁷ m s⁻¹를 초과하면 D가 급격히 상승한다는 임계값이 존재한다는 점을 제시한다.

다섯째, 실험 결과를 Campi Flegrei 화산암에 적용했을 때, 현장 알칼리장석의 Ba·Sr D값을 이용해 결정화가 시작된 후 최대 6 일(비평형)에서 최소 9 일(평형에 근접) 사이의 시간을 추정할 수 있었다. 이는 화산 활동 전·후의 마그마 저장소에서 결정화가 얼마나 빠르게 진행되는지를 정량화하는 새로운 타이머 역할을 할 수 있음을 의미한다.

전반적으로 본 연구는 미량원소 분배가 단순히 온도·압력·조성에 의존하는 것이 아니라, 결정 성장 속도와 확산 재평형 시간이라는 동역학적 인자에 의해 크게 변조될 수 있음을 실험적으로 입증하였다. 이는 기존의 평형 기반 분배 모델을 보완하고, 화산·심성 마그마 시스템에서 결정화 속도와 시간 규모를 역추적하는 데 중요한 방법론적 기반을 제공한다.