MRO CRISM 다중스펙트럼 데이터의 램버트 반사율 추출 DISORT 기반 대기·열 보정 파이프라인

초록

본 논문은 MRO/CRISM 다중스펙트럼 매핑 모드 데이터를 Lambertian 알베도 스펙트럼으로 변환하기 위한 DISORT 기반 방사전달 파이프라인을 소개한다. 현재 단계에서는 MGS‑TES와 MOLA에서 얻은 기후·지형 정보를 이용해 에어로졸 광학두께, 표면·대기 온도, 압력을 보정 입력값으로 사용한다. 파이프라인은 측정된 I/F와 재구성된 광학각을 입력으로 하여 대기·열 효과를 제거하고, 지질학적 해석에 유용한 알베도 스펙트럼을 산출한다.

상세 분석

이 연구는 화성 궤도선인 MRO에 탑재된 CRISM의 다중스펙트럼 매핑 데이터를 정량적인 지표면 반사율, 즉 Lambertian 알베도로 변환하기 위한 전반적인 방법론을 제시한다. 핵심은 DISORT(Discrete Ordinate Radiative Transfer) 모델을 활용한 전방향 방사전달 계산이며, 이를 위해 입력 파라미터로 에어로졸(먼지·얼음) 광학두께, 표면 압력, 표면 및 하층 대기 온도, 그리고 관측 기하학(태양·관측자 입사각·방출각·방위각)을 필요로 한다. 현재 파이프라인은 이러한 파라미터를 직접 CRISM 데이터에서 추출하기보다는 기존 궤도선(MGS‑TES, MOLA, Viking)에서 구축된 기후·지형 클리마톨로지를 활용한다.

먼저, 에어로졸 광학두께는 TES에서 파생된 먼지·얼음 입자 농도를 CRISM 관측 파장대(400–2500 nm)로 스케일링한다. 이때 파장 의존성을 고려한 Mie 계산 혹은 경험적 스펙트럼 보정이 적용되며, 이는 CO₂ 흡수대(특히 2 µm 근처 H₂O 얼음 밴드)에서 대기 효과를 정확히 제거하는 데 필수적이다. 표면 압력은 MOLA 고도와 TES 하층 대기 온도, 그리고 Viking 착륙지점 기반의 압력-고도 관계식을 결합해 계산한다. 이는 대기 광학두께와 흡수 계수를 압력에 따라 조정하는 데 사용된다.

표면 온도와 하층 대기 온도는 TES의 열복사 스펙트럼을 이용해 월·계절별 평균값을 추출한다. 이러한 평균값을 사용함으로써 단일 관측 시점의 온도 변동성을 완전히 반영하지 못한다는 한계가 있지만, 현재는 전역적인 파이프라인 자동화를 위해 충분히 실용적이다.

DISORT 모델은 입력된 광학두께·압력·온도와 관측 기하학을 바탕으로 다중 산란 및 흡수를 포함한 복합 방사전달을 계산한다. 모델 출력은 대기·열 보정된 지표면 방사율이며, 이를 π로 나누어 Lambertian 알베도로 변환한다. 여기서 I/F는 측정된 복사량을 태양 복사량으로 정규화한 값이며, Lambertian 가정은 실제 표면의 비등방성 반사 특성을 무시하지만, 다중스펙트럼 매핑 해상도(≈200 m)에서는 평균적인 반사 특성을 충분히 대변한다.

파이프라인 검증을 위해 북극 얼음캡, Tyrrhena Terra, Phoenix 착륙지점 등 세 지역의 CRISM 데이터에 적용하였다. 특히 북극 얼음캡에서는 CO₂ 흡수와 H₂O 얼음 밴드가 겹치는 2 µm 근처에서 에어로졸 보정이 없을 경우 알베도 스펙트럼이 크게 왜곡되는 것을 확인했다. 반면, 에어로졸 보정을 포함하면 물 얼음의 고유 스펙트럼이 명확히 드러나며, 이는 광물·얼음 함량 추정에 직접적인 이점을 제공한다.

향후 단계에서는 TES 기반 클리마톨로지를 대체하여 CRISM 자체 관측을 이용해 시공간적으로 변동하는 에어로졸 광학두께·표면 압력·온도를 직접 추정할 계획이다. 이를 위해 다중스펙트럼 데이터의 전파역학적 특성을 활용한 역문제 해결 기법(예: 최적화 기반 대기 역전산)과 머신러닝 기반 보정 모델을 도입할 예정이다. 이러한 발전은 전 지구적 고해상도 알베도 지도 제작과, 화성 표면·대기의 상호작용을 정량적으로 이해하는 데 핵심적인 역할을 할 것으로 기대된다.