Geochemical modelling of ore-forming processes in the Börzsöny Mountains, N-Hungary

기술 요약: 뵈르죄니 산맥(Börzsöny Mountains)의 광화 작용 과정에 대한 지화학적 모델링

문제 제기
헝가리 북부의 중앙 뵈르죄니 산맥은 네오기(Neogene) 칼크-알칼리 화산 활동 중에 형성된 금 함유 환원성 포피리 광상(reduced porphyry deposit)을 보유하고 있다. 자철석, 적철석, 안하이드라이트(anhydrite)를 특징으로 하는 일반적인 산화형 포피리 구리 광상(PCD)과 달리, 이 시스템은 황산염과 자철석이 결여된 환원된 광물 조합(휘몰리브데나이트, 휘철석, 황동석)으로 정의된다. 이전의 지질학 및 암석학 연구를 통해 온도, 산화-환원 조건, 황 자유도(sulfur fugacity)의 변화를 포함하는 3단계 광화 모델이 확립되었으나, 복잡한 공생 단계(paragenetic sequence)를 구동하는 구체적인 열역학적 제어 기작은 정량적으로 규명되지 않은 상태였다. 본 연구는 다수의 유체 펄스(fluid pulses)와 물리화학적 변동이 어떻게 광물 침전을 제어했는지 이해하기 위해, 이 환원된 열수 시스템의 진화를 시뮬레이션해야 할 필요성을 다룬다.

방법론
저자들은 열수 시스템의 열역학적 진화를 시뮬레이션하기 위해 PHREEQC 3.0을 이용한 지화학적 모델링을 채택하였다. 단일 모델로는 관찰된 전체 공생 단계를 재현할 수 없다는 점을 인식하여, 두 가지 상호 보완적인 열역학적 배치(batch) 모델을 개발하였다:

• 모델 1 (환원 단계): 주요 황화물 조합을 형성하는 초기 고온(400–150 °C) 환원 단계를 시뮬레이션한다. 이 모델은 bl-1kb.dat 데이터베이스(1 kbar 압력에 대해 보정됨)를 사용하며, Thermoddem 데이터베이스와 WORM(Water-Organic-Rock-Microbe) 포털에서 유래한 커스텀 열역학 데이터를 사용하여 누락된 종(예: 휘 arsenic-pyrite, 방연석, 섬아연석)을 포함한다. 모델은 중등도 환원된 마그마-열수 유체와 안데사이트-데사이트질 모암 사이의 반응을 추적한 후, 서로 다른 유체 펄스(철-풍부, 황-풍부)의 주입을 추적한다.
• 모델 2 (산화적 오버프린트): 후속되는 저온(<300 °C) 산화 사건을 시뮬레이션한다. 이 모델은 LLNL.dat 데이터베이스를 사용하며, 모델 1의 출력을 초기 평형 상(initial equilibrium phases)으로 포함한다. 이 모델은 점진적인 산화, 황 함량 감소, 그리고 다양한 pH가 기존 환원 광물의 안정성에 미치는 영향을 테스트한다.

시뮬레이션은 관찰된 광물 서열을 재현하기 위해 유체의 조성(주요 이온, 기저 금속, 황 종)을 반복적으로 조정하고 pH, pe, 산소 자유도($fO_2$), 황 자유도($fS_2$)와 같은 매개변수에 대한 민감도 분석을 수행하였다.

주요 결과
모델링은 관찰된 공생 단계를 성공적으로 재현하였으며, 이는 광 형성이 단순한 냉각 과정이 아니라 서로 다른 물리화학적 특징을 가진 다수의 유체 펄스를 필요로 했음을 확인시켜 주었다:

1. 주요 환원 조합 (모델 1):

   • 400 °C에서 초기 환원 유체로부터 **휘arsenic-pyrite(arsenopyrite)**가 침전되었다.
   • 시스템이 150 °C로 냉각됨에 따라 **방연석-I(galena-I)**과 **섬아연석-I(sphalerite-I)**이 결정화되었고, 이어 소량의 **황동석-I(chalcopyrite-I)**이 형성되었다.
   • 낮은 황 자유도를 가진 철-풍부 유체의 후속 펄스가 250 °C에서 발생하여 황동석-I을 불안정하게 만들고 **휘철석(pyrrhotite)**의 침전을 촉진했다.
   • 세 번째의 고온(400 °C), 황-풍부 유체 펄스는 방연석-II와 섬아연석-II의 형성을 유발했다.

2. 산화적 오버프린트 (모델 2):

   • 점진적인 산화와 황 결핍은 300~200 °C 사이에서 황동석-II와 휘arsenic-pyrite-II의 침전을 초래했다.
   • 기존의 휘철석은 **황철석(pyrite)**과 **적철석(hematite)**으로 치환되었다.
   • 시스템이 진화함에 따라 **사이드라이트(siderite)**와 **방해석(calcite)**이 더 낮은 온도(150–100 °C)에서 침전되었다.
   • **마르카사이트(marcasite)**는 모델의 평형 조건 하에서는 침전되지 않았으며, 황철석 침전을 인위적으로 금지했을 때만 형성되었다. 이는 평형 모델이 준안정상(metastable phases)을 포착하는 데 갖는 한계를 보여준다.

주요 기여

• 다단계 진화의 검증: 본 연구는 뵈르죄니 산맥의 환원성 포피리 시스템이 온도, 산화-환원 상태, 유체 조성의 동시적인 변동을 포함하는 뚜xt별 유체 펄스를 통해 진화했다는 열역학적 증거를 제공한다.
• 제어 기작: 모델은 광물 안정성에 대한 특정 제어 요인을 식별하였다:
  • 휘arsenic-pyrite의 안정성은 온도와 용존 탄산염 함량에 의해 제어된다.
  • 고온에서의 방연석 및 섬아연석 형성은 용존 기저 금속 비율과 황 함량에 의해 제어된다.
  • 황동석은 환원된 고온 환경에서 매우 불안정하며, 저온(<150 °C) 또는 산화적 오버프린트 중에만 침전되는 것으로 나타났다.
• 저구리 함량에 대한 설명: 주요 환원 단계에서 황동석의 불안정성은 산화형 시스템에 비해 환원성 포피리 광상에서 나타나는 특징적인 낮은 구리 함량을 설명하는 열역학적 근거를 제공한다.

의의 및 주장
저자들은 본 연구 결과가 준안정 침전(예: 마르카사이트)에 관한 내재적 한계에도 불구하고, 환원성 포피리 관련 열수 시스템을 조사하는 데 있어 PHREEQC 기반 지화학 모델링의 적용 가능성을 입증한다고 주장한다. 본 연구는 산화-환원 조건만으로는 광물 서열을 설명하기에 불충분하며, 온도와 유체 화학의 동시적인 변화가 필수적임을 확인한다.

또한, 본 논문은 네오기 화산 과정과 광 형성 사이의 밀접한 유전적 연관성을 확립한다. 연구는 마그마 상승 과정 중 하부 지각 내의 흑연질 변성 퇴적암의 동화(assimilation)에 의해 광화 용융물의 환원이 촉진되었음을 시사하며, 이는 초기 화산 단계의 가넷(garnet) 조성을 통해 입증된다. 이후의 산화적 오버프린트는 화산 활동의 최종 단계에 기인한다. 궁극적으로, 이 연구는 중앙 뵈르죄니 광화를 복잡한 다중 펄스 열수 사건에 의해 형성된 금 함유 환원성 포피리 광상으로 분류하는 것을 뒷받침한다.