How Uncertain Are Estimates of Marine Carbon Dioxide Removal? Insights from a Simplified Model

기술 요약: 해양 이산화탄소 제거 추정치는 얼마나 불확실한가? 단순화된 모델을 통한 통찰

문제 제식 (Problem Statement)
해양 이산화탄소 제거(mCDR) 기술, 특히 해양 알칼리도 증진(OAE) 및 직접 해양 제거(DOR)는 대기 중 이산화탄소를 흡수하기 위해 표층 해수의 이산화탄소 분압을 낮추는 방식에 의존한다. 이러한 개입의 이론적 최대 효율(예: OAE의 경우 알칼리도 1몰당 약 0.8몰의 $CO_2$)은 알려져 있으나, 탄소 제거의 시간적 진화는 물리적 및 생지화학적 과정에 매우 민감하게 반응한다. 현재 mCDR 효능에 대한 추정치는 해양 모델에 크게 의존하고 있지만, 이러한 모델들에 내재된 불확실성—특히 수직 혼합, 해양-대기 가스 교환, 그리고 수괴(watermass) 특성—은 제대로 이해되지 않고 있다. 기존의 글로벌 모델들은 종종 자연적 변동성을 과소평가할 수 있는 고정된 매개변수화를 채택하고 있으며, 그 복잡성으로 인해 CDR 계산을 주도하는 특정 메커니즘을 분리해내는 것이 어렵다. 이러한 정량화된 불확실성의 결여는 mCDR 프로젝트를 위한 모니터링, 보고 및 검증(MRV) 프레임워크와 순 탄소 중립 달성 시점 결정에 근본적인 과제를 제기한다.

방법론 (Methodology)
이러한 불확실성을 체계적으로 조사하기 위해, 본 연구는 1차원(1D) 비생물적(abiotic) 해양 모델을 사용한다. 이 축소된 복잡성의 접근 방식은 수평 방향의 이류(advection), 침강(subduction) 및 생물학적 상호작용은 무시하고 수직 과정(난류 확산 및 해양-대기 가스 교환)만을 분리하여 다룬다. 모델은 총 알칼리도(TA)와 용존 무기 탄소(DIC)를 시간에 따라 진화시키며, 해양-대기 가스 교환을 기반으로 표층 $CO_2$ 플럭스를 계산한다.

본 연구는 두 가지 주요 분석 프레임워크를 활용한다:

1. 민감도 시뮬레이션 (Sensitivity Simulations): 개별 모델 매개변수(예: 혼합층 깊이, 풍속, 가스 전달 계수, 피크노클라인 확산계수, 배경 TA/DIC, 투입 깊이)를 관측된 자연적 변동성과 일치하는 값으로 섭동(perturbation)시킨다. 이를 기준 시뮬레이션과 비교하여, CDR의 기준 시뮬레이션 대비 비율로 정의되는 "CDR 민감도 비율"($r_{CDR}$)을 결정한다.
2. 몬테카를로 앙상블 (Monte Carlo Ensembles): 관측된 해양 변동성을 나타내는 통계적 분포에서 추출한 매개변수를 사용하여 1,000쌍 이상의 시뮬레이션(기준 및 개입 모델)을 동시에 실행한다. 이 접근 방식은 물리적 과정과 수괴 특성의 결합된 불확실성이 CDR 추정치에 어떻게 나타나는지 평가한다.

시뮬레이션은 깊이의 영향을 평가하기 위해 200m와 2,000m의 해양 깊이에서 50년 동안 수행된다.

주요 결과 (Key Results)
연구는 CDR 불확실성을 지배하는 세 가지 뚜렷한 시간적 영역을 식별하였다:

1. 초기 급속 평형 단계 ($t \ll \tau_{mix}$): 초기 몇 년 동안, CDR 추정치는 초기 해양-대기 가스 교환 시간 척도($\tau_{air-sea}$)를 제어하는 매개변수(혼합층 깊이, 풍속, 가스 전달 계수, 배경 탄산염 화학 성분(TA, DIC), 온도)에 매우 민감하다. 이 단계에서 CDR 추정치는 매개변수 편향에 따라 최대 2배까지 차이가 날 수 있다.
2. 혼합층으로부터의 알칼리도 손실 ($t \approx \tau_{mix}$): 수직 혼합이 알칼리도를 표층에서 제거함에 따라, 불확실성은 해양-대기 교환 시간 척도와 혼합 시간 척도($\tau_{mix}$) 사이의 상호작용에 의해 지배된다. 피크노클라인 확산계수($\kappa_{pyc}$)에 대한 불확실성이 결정적인 역할을 하며, 이로 인해 CDR 추정치의 변동성은 투입 후 2~3년 사이에 정점에 달한다(최대 30%).
3. 심해 혼합 ($t \gg \tau_{mix}$): 장기적으로, 심해로의 혼합 시간 척도가 해양-대기 평형 시간 척도와 비슷하거나 더 빠를 경우 불확실성은 수십 년간 지속된다. 본 연구는 피크노클라인 확산계수가 장기적인 CDR 불측성을 결정하는 지배적인 요인임을 발견했다. 이는 해당 계수가 평형에 도달한 알칼리도의 분율과 평형에 도달하지 않은 심층수가 표층으로 돌아오는 속도 모두를 제어하기 때문이다.

정량적 결과:

• 불확실성의 크기: 초기 몇 년 동안 CDR 추정치는 기준 시뮬레이션과 최대 100%(2배)까지 차이가 날 수 있다. 2050년이 지나면 이 차이는 1020%로 좁혀지지만, 심해 투입의 경우 여전히 상당한 불확실성이 존재한다.
• 투입 깊이: 혼합층 아래(심지어 20~40m 아래라도)에 알칼리도를 투입하면 효율이 현저히 낮아지고 불확실성이 높아지며, 일부 추정치는 50년이 지난 후에도 기준 시뮬레이션보다 훨씬 낮은 수준에 머문다.
• 순 중립 시점: (전과정 평가 비용을 고려한) 개입이 순 탄소 음성(net carbon-negative)이 되는 데 걸리는 시간은 매우 가변적이다. 기준 사례는 약 1.4년이 소요된다고 제안하지만, 몬테카를로 앙상블은 넓은 분포를 보여준다: 약 30%의 시뮬레이션은 1년 이내에 순 중립에 도달하지만, 50%는 2년 이상 걸리며, 약 4%는 50년 이상 걸린다. 이러한 "긴 꼬리(long tail)" 현상은 주로 심해 투입 또는 조기 평형을 방해하는 빠른 혼합과 관련이 있다.
• 깊이 의존성: 불확실성은 200m 깊이의 해양보다 2,000m 깊이의 해양에서 더 오래 지속된다. 200m의 경우 전체 수층이 시뮬레이션 기간 내에 평형을 이룰 수 있기 때문이다.

의의 및 주장 (Significance and Claims)
본 논문은 mCDR 추정치의 불확실성이 단순히 일시적인 것이 아니라, 해양-대기 가스 교환 대비 내부 혼합의 시간 척도로 인해 수십 년간 지속될 수 있다고 주장한다. 저자들은 이러한 결과가 물리적 과정이 제대로 규명되지 않은 지역 및 글로벌 mCDR 시뮬레이션에 근본적인 과제를 제시한다고 주장한다.

구체적으로, 본 연구는 다음을 강조한다:

• MRV에 대한 시사점: 탄소 제거의 시간적 진화에 대한 불확상성은 순 탄소 중립 시점에 직접적인 영향을 미치며, 이는 탄소 크레딧 체계에서 매우 중요하다. 불확실성이 잘 정량화되지 않는다면 탄소 크레딧 부여가 성급하거나 부정확할 수 있다.
• 모델의 한계: 현재의 해양 모델들은 상당한 공간적, 시간적 편향을 겪을 수 있는 난류 및 가스 교환 매개변수화에 의존하는 경우가 많다. 본 연구는 이러한 특정 물리적 과정(혼합 및 가스 전달)을 검증하는 것이 일반적인 기후 예측보다 mCDR에 있어 더 중요하다고 제안한다.
• 투입 전략: 결과는 효율을 극대화하고 불확실성을 최소 최소화하기 위해 알칼리도를 반드시 표층 혼합층 내에 투입해야 할 필요성을 강조한다.

저자들은 자신의 1D 모델이 수평 이류와 생물학적 피드백을 제외하고 있다는 점을 언급하며 연구 범위를 겸손하게 설정하였으나, 비생물적 시스템에서도 물리적 과정만으로 발생하는 불측성이 충분히 크다는 점을 들어, mCDR 개입이 수십 년의 척도에서 진정으로 탄소 음성이 되기 위해서는 심도 있는 추가 조사가 필요하다고 결론짓는다.