Pyrite oxidation flux underestimated due to landscape-mediated microbial sulfate reduction

기술 요약: 경관 매개 미생물 황산염 환원에 의한 황철석 산화 플럭스의 과소평가

문제 제기
황철석의 산화적 풍화는 전 지구적 탄소, 산소 및 황 순환의 핵심적인 구성 요소이다. 이는 대기 중 $O_2$를 소비하고, 용존 황산염($SO_4^{2-}$)을 생성하며, 탄산염 암석을 산성화하여 해양-대기 시스템으로 $CO_2$를 방출한다. 산맥은 이러한 과정의 핫스팟으로 알려져 있으나, 황철석 산화의 순 플럭스(net flux)를 정량화하는 것은 여전히 어려운 과제이다. 전통적으로 하천의 황산염 기원은 안정 동위원소 조성($\delta^{34}S_{SO4^{2-}}$ 및 $\delta^{18}O_{SO4^{2-}}$)을 사용하여 구분된다. 전형적인 해석은 높은 동위원소 수치가 증발암 용해(이는 $O_2$를 소비하거나 $CO_2$를 방출하지 않음)를 나타낸다고 가정하는 반면, 낮은 수치는 황철석 산화를 나타낸다고 본다. 그러나 육상 환경에서의 미생물 황산염 환원(MSR)이 이러한 동위원소 신호를 덮어쓴다면(overprint), 이 접근 방식은 결함이 있을 수 있다. MSR은 황산염을 황화물로 환원시킨 후 다시 산화시키는 과정을 포함하며, 이 과정은 $\delta^{34}S$를 크게 변화시키지 않으면서 $SO_4^{2-}$ 풀(pool) 내의 $^{18}O$를 농축시킬 수 있다. 만약 MSR이 광범위하게 존재함에도 고려되지 않는다면, 현재의 황철석 산화 플럭스 추정치—결과적으로 전 지구적 $CO_2$ 배출량 및 $O_2$ 소비량—는 체계적으로 과소평가되어 있을 수 있다.

연구 방법론
육상 환경에서 MSR의 통제 기제에 대한 메커니즘적 이해의 부족을 해결하기 위해, 저자들은 북서 히말라야의 강가 강(알라크난다 및 바기리티티 하천) 유역에서 통합 연구를 수행하였다. 이 지역은 침식률, 강수량, 온도 및 암석학적 특성이 뚜렷한 구배를 보이는 자연 실험실 역할을 한다.

본 연구는 다중 프록시(multi-proxy) 접근법을 채택하였다:

1. 지구화학적 샘사플링: 고도 구배(1,250–5,300 m)에 따라 하천수 샘플($n=26$)을 수집하였다. 주요 이온 농도는 이온 크로마토그래피를 통해 측정되었다.
2. 동위원소 분석: $\delta^{34}S_{SO4^{2-}}$, $\delta^{18}O_{SO4^{2-}}$, 그리고 $\delta^{18}O_{H2O}$를 측정하였다. 1차 산화 신호와 2차 오버프린팅을 구분하기 위해 산소 편차($\Delta^{18}O = \delta^{18}O_{SO4^{2-}} - \delta^{18}O_{H2O}$)를 계산하였다.
3. 미생물학적 분석: 황산염 환원 박테리아의 프록시인 dissimilatory sulfite reductase의 $\beta$-subunit을 인코딩하는 dsrB 마커 유전자의 풍부도와 총 원핵생물 풍부도(16S rRNA)를 qPCR을 통해 정량화하였다.
4. 지형학적 제약 조건: $^{10}Be$ 농도와 정규화된 채널 경사 지수($k_{sn}$)를 사용하여 유역 평균 침식률을 결정하였다.
5. 암석학적 분할 모델: 주요 이온 기원을 분리하기 위해 몬테카를로 역산 기반 모델을 개발하였다. 결정적으로, 이 모델은 황철석 산화를 모암 암석학과 독립적인 별개의 엔드멤버(end-member)로 취급하였으며, 순환 논리를 피하기 위해 혼합 계산에서 동위원소 데이터를 제외하였다. 이를 통해 황철석 유래 황산염의 분율($f_{py}$)을 독립적으로 정량화할 수 있었다.
6. 전 지구적 종합: 이 프레임워크를 전 지구적 하천 화학 데이터(GLORICH 데이터베이스 포함)에 적용하여 전 지구적 황철석 산화 플럭스를 재평가하였다.

주요 결과

• 황철석 산화의 지배성: 역산 모델은 황철석 산화가 강가 상류의 하천 황산염의 주요 기원임을 나타내며, $f_{py}$는 0.73에서 1.0 사이(평균 $\sim$0.96)이다. 이는 해당 지역의 높은 $\delta^{34}S$ 값을 증발암 용해의 결과로 보았던 기존의 해석과 대조적이다(해당 지역에는 지도상에 표시된 증발암 형성층이 존재하지 않음에도 불구하고).
• MSR 오버프린팅 메커니즘: 강가 상류의 $\Delta^{18}O$ 값(4.6 ~ 14.3‰ 범위)은 1차 황철석 산화의 이론적 최대치($\sim$3.5–6.5‰)를 크게 초과한다. 상승된 $\Delta^{18}O$는 황산염 농도로 정규화된 dsrB 유전자 풍부도와 강한 상관관계($R^2 = 0.71$)를 보이며, 이는 MSR 및 그에 따른 재산화가 $^{18}O$ 농축의 원인임을 확인해 준다.
• 경관 매개 제어: MSR 오버프린팅의 정도는 침식과 강수량에 의해 조절되는 수문학적 통과 시간(transit times)에 의해 제어된다:
  • 침식: $\Delta^{18}O$는 침식률과 음의 상관관계를 갖는다. 낮은 침식률(두꺼운 레골리스, 긴 물의 통과 시간)은 MSR에 필요한 무산소 미세 환경을 촉진한다.
  • 강수량: $\Delta^{18}O$와 강수량의 관계는 비단조적(non-monotonic)이다. 저침식 유역에서 높은 강수량은 불포화대(vadose zone)를 포화시켜 얕은 흐름 경로를 활성화하고 통과 시간을 단축함으로써 MSR을 감소시킨다. 반대로 고침식 유역에서는 낮은 강수량이 얇은 레골리스로 인해 MSR을 제한한다.
  • 최적 조건: MSR 오버프린팅은 낮은 물리적 침식률이 저~중등도의 강수량과 결합될 때 극대화된다.
• 전 지구적 함의: 이러한 통찰을 전 지구적 하천 화학 데이터에 적용한 결과, 황철석 풍화가 지배적인 강들은 이전에 증발암 혼합으로 간주되었던 $\delta^{34}S$ 대 $\Delta^{18}O$ 공간 전체에 걸쳐 분포함을 보여준다. MSR 오버프린팅을 보정했을 때, 황철석 유래 황산염의 전 지구적 분율($f_{py}$)은 약 0.90(더 넓은 산악 데이터셋의 경우 0.72)으로 추정되며, 이는 기존의 동위원소 기반 추정치인 약 0.58보다 현저히 높다.

의의 및 주장
본 논문은 하천 황산염 동위원소를 조절하는 데 있어 육상 MSR의 역할이 체계적으로 간과되어 왔다고 주장한다. 저자들의 주장은 다음과 같다:

1. 동위원소 오해: 황철석과 증발암을 구분하기 위해 $\delta^{34}S$와 $\Delta^{18}O$에 의존하는 기존의 전 지구적 황철석 산화 플럭스 추정치는, MSR에 의한 동위원소 변화를 증발암 용해로 치부함으로써 황철석의 기여도를 과소평가했을 가능성이 높다.
2. 수정된 플럭스 추정치: MSR을 고려할 때, 황철석 풍화 유래의 전 지구적 황산염 플럭스는 2.0–2.5 Tmol yr$^{-1}$로 추정되며, 이는 기존의 동위원소 기반 추정치($\sim$1.3–1.6 Tmol yr$^{-1}$)보다 약 1.5~1.9배 더 높다.
3. 탄소 순환 영향: 황철석 산화가 탄산염으로부터 $CO_2$를 방출하므로, 이러한 과소평가는 산맥 형성 과정이 현재 믿어지는 것보다 대기 중 $CO_2$의 더 큰 원천임을 시사한다. 비록 전 지구적 황 순환이 여전히 순(net) 대기 $O_2$ 축적을 유도할 것이나(해양 퇴적물 내 황철석 매몰이 지배적인 싱크 역할을 하므로), 그 불균형의 규모는 기존 추정치보다 작다.
4. 미래 재구성: 저자들은 향end의 지질학적 탄소 예산 재구성이 황철석 산화 플럭스를 정량화할 때, 육상 MSR의 역할을 반드시 고려해야 한다고 결론짓는다. 이는 결정론적 암석학적 가변성만큼이나 육상 황산염 동위원소 조성 결정에 중요하다.