Coupled atmospHere Interior modeL Intercomparison (CHILI). I. Evolutionary Modelling -- Primordial Magma Oceans of Earth and Venus

기술 요약: 결합된 대기-내부 모드 상호 비교 (CHILI). I. 진화적 모델링 – 지구와 금성의 원시 마그마 바다

문제 제기
지구와 금성은 초기에 용융된 "마그마 바다" 상태로부터 서로 다른 진화적 결과를 나타내지만, 이들의 초기 열적 및 지질학적 역사는 여전히 불분명하게 남아 있다. 수치 시뮬레이션은 이러한 역사를 재구성하고 외계 행성 관측을 해석하는 주요 도구이지만, 기존의 결합된 내부-대기 코드들은 이질적인 근사법과 단순화를 채택하고 있다. 이러한 불일치는 행성의 냉각 시간 규모, 대기 조성 및 휘발성 성분 예산을 예측하는 데 상당한 불확실성을 초래한다. CHILI 프로젝트는 태양계의 지구형 행성이라는 "정답(ground truth)"을 기준으로 이 모델들을 벤치마킹하여, 모델의 핵심적인 민감도와 한계를 식별하는 것을 목표로 한다.

방법론
본 연구(CHILI의 Part I)는 8개의 서로 다른 진화 코드를 통제된 방식으로 상호 비교한다: GOOEY, NEONGOOEY, PROTEUS, PACMAN, LINCS, MOAI, PlanAtMO, 그리고 MOAI의 초기/후기 변형 모델이다. 이 모델들은 "표준" 지구 및 금성 시나리오에 대해 결합된 내부 및 대기의 진화를 시뮬레이션한다.

• 초기 조건: 모든 시뮬레이션은 항성 연령 50 Myr(백만 년)에서 완전히 용융된 상태로 시작하며, 고정된 핵 반경(행성 반경의 55%)과 동일한 초기 휘발성 재고량($4.7 \times 10^{20}$ kg의 수소(지구 해양 3개 분량에 해당) 및 $2.73 \times 10^{20}$ kg의 탄소)을 가진다.
• 파라미터 그리드: 표준 사례를 넘어, 본 연구는 휘발성 성분 예산에 대한 민감도를 평가하기 위해 수소(1, 5, 10 ocean) 및 탄소 재고량을 변화시키는 9가지 시나리오 그리드를 탐색한다.
• 모델 제약 조건: 시뮬레이션은 구름과 에어로졸을 무시하며, 고정된 10%의 본드 알베도(Bond albedo)를 적용하고, (MOAI를 제외하고) 레일리 산란을 비활성화한다. 종료 기준은 모델마다 다르지만 일반적으로 맨틀 용융 분율이 특정 임계값(1%에서 5% 사이) 아래로 떨어질 때 발생한다.
• 물리적 변수: 코드들은 복사 전달(line-by-line vs. correlated-k vs. gray), 대기 탈출(분별 탈출 vs. 비분별 탈출, 에너지 제한형 vs. 확산 제한형), 맨틀 대류(경계층 이론 vs. 혼합 길이), 유변학(rheology), 그리고 휘발성 성분 분배(평형 vs. 동적 용융 포획) 처리 방식에서 차이를 보인다.

주요 기여 및 결과
상호 비교 결과, 모든 모델이 지구의 마그마 바다의 급격한 고결(solidification)에 대해 질적으로는 일치하지만(표준 사례의 경우 약 4 Myr 이내), 냉각 시간 규모와 대기 진화에 있어서는 상당한 양적 차이를 보였다.

1. 냉각 시간 규모 및 고결:

   • 지구: 표준 지구 모델들은 하데안(Hadean) 지르콘의 경험적 제약 조건과 일치하는 4 Myr 이내의 고결을 예측한다. 그러나 냉각 시간 규모는 초기 휘발성 재고량보다 특정 모델 물리(예: 맨틀 점성도, 용융 곡선)에 더 민зу하게 반응한다.
   • 금성: 금성 시나리오는 더 큰 변동성을 보인다. 몇몇 모델(LINCS, PROTEUS, PlanAtMO)은 행성이 전 지구적 복사 평형(폭주 온실 효과) 상태에 진입하여, 대기 탈출이 온실가스를 제거할 때까지 최대 49 Myr 동안 마그마 바다를 유지하는 "Type II" 거동을 예측한다. 다른 모델들(MOAI, GOOEY)은 직접적인 고결("Type I")을 예측한다.

2. 대기 조성 및 휘발성 성분 분배:

   • 탄소 종 분화: 초기 대기는 CO/CO$_2$가 지배적일 것으로 예측된다. 그 비율은 맨틀의 산화 환원 상태($f_{O2}$)와 표면 온도에 크게 의존한다. PROTEUS(고정된 $f_{O2}$)와 같은 모델은 CO$_2$가 풍부한 대기를 예측하는 반면, MOAI(동적 산화 환원)는 CO가 지배적인 대기를 예측한다.
   • 수소 저장: 모델 간 불일치의 핵심적인 원인은 고체 맨틀 내 수소 저장 처리 방식이다. PACMAN과 MOAI는 "동적 용 melt trapping"을 포함하여 상당량의 수소를 고결된 맨틀에 격리함으로써, 더 얇고 CO$_2$가 지배적인 대기를 형성한다. 고체상 저장을 무시하는 PROTEUS는 거대한 H$_2$O 지배적 대기(>450 bar)를 예측한다.
   • 탈출: 대기 탈출은 주요한 불일치 지점이다. LINCS는 상당한 수소 손실(Type II 거동)을 예측하는 반면, GOOEY는 무시할 만한 수준의 손실을 예측한다. 연구 결과, 탈출이 일부 모델에서 맨틀의 산화 과정을 유도하긴 하지만, 내부 지질 역학 및 복사 전달에 비해 열 진화 시간 규모를 조절하는 역할은 부차적임을 발견했다.

3. 복사 전달 및 지질 역학:

   • OLR 변동: 출구 장파 복사(OLR)는 복사 전달 방식(예: 충돌 유도 흡수, 성층권 온도 프로파일 포함 여부)의 차이로 인해 모델 간 수십 배의 차이를 보인다. OLR-T$_{surf}$ 관계의 이러한 차이는 초기 탄소 예산의 변화보다 냉각률의 발산을 더 강력하게 주도한다.
   • 맨틀 유변학: 점성 파라미터화와 용융 곡선(고상선/액상선)의 차이는 열 수송 속도에 상당한 영향을 미친다. 예를 들어, PROTEUS는 침강하는 쇄설물(cumulates)에 의한 대류 억제와 특정 점성 법칙의 처리로 인해 더 느린 냉각을 예측한다.

의의 및 주장
본 논문은 현재의 마그마 바다 모델 "zoo(집단)"가 단순히 알려지지 않은 초기 조건 때문이 아니라, 주로 서로 다른 물리적 파라미터화 때문에 서로 다른 진화 경로를 나타낸다고 주장한다.

• 모델 민감도: 본 연구는 휘발성 분배(특히 심부 수소 저장), 대기 탈출 효율, 맨틀 점성도, 그리고 용융 곡선이 모델 간의 발산을 일으키는 가장 결정적인 민감도임을 식다.
• 외계 행성에 대한 시사점: 지구와 금성이 벤치마킹을 위한 유일한 "정답"이기 때문에, 이러한 불일치는 외계 행성 관측(예: JWST)을 해석하는 능력을 제한한다. 저자들은 이러한 특정 물리적 처리에 대한 검증이 의미 있는 외계 행성 데이터 해석을 위한 전제 조건이라고 주장한다.
• 향-후 방향: 논문은 미래의 모델 개선이 정확한 복사 전달(연속 흡수 포함), 견고한 melt trapping 및 유변학 파라미터화, 그리고 조석 가열 및 구름 형성과 같은 미모델링된 물리학의 포함에 집중되어야 한다고 결론짓는다. 저자들은 이러한 발견이 CHILI 프로젝트의 Part II(TRAPPIST-1 외계 행성) 및 Part III(정적 대기 상호 비교)를 위한 토대를 마련했다고 언급한다.

저자들은 단일 모델이 "옳다"고 단정하지 않고, 현재의 상호 비교가 근본적인 물리 법칙의 결정론적 해결책보다는 인식론적 불확실성을 강조한다는 점을 인정하며 신중한 입장을 유지한다. 본 연구의 주요 가치는 향후 관측 및 이론적 노력을 안내하기 위해 이러한 불확실성을 정량화하는 데 있다.