日本地球惑星科学連合2022年大会

講演情報

[J] 口頭発表

セッション記号 S (固体地球科学) » S-CG 固体地球科学複合領域・一般

[S-CG48] 海洋底地球科学

2022年5月27日(金) 15:30 〜 17:00 105 (幕張メッセ国際会議場)

コンビーナ:沖野 郷子(東京大学大気海洋研究所)、コンビーナ:田所 敬一(名古屋大学地震火山研究センター)、座長:安川 和孝(東京大学大学院工学系研究科エネルギー・資源フロンティアセンター)、田中 えりか(海洋研究開発機構)

16:15 〜 16:30

[SCG48-16] 鉱石生成プロセスの解明に向けた沖縄トラフ伊是名海穴熱水域での反応輸送解析

*冨田 昇平1小池 克明1後藤 忠徳2鈴木 勝彦3笠谷 貴史3 (1.京都大学大学院工学研究科、2.兵庫県立大学大学院生命理学研究科、3.海洋研究開発機構)


キーワード:数値シミュレーション、TOUGHREACT、海底熱水鉱床、鉱物沈殿、化学反応

Seafloor massive sulfide deposits have been attracting attention as new metal resources. Their exploration needs an understanding of the genesis of the deposits. Previously, we conducted a hydrothermal flow simulation of the Izena Hole hydrothermal field, Okinawa Trough and clarified the essential factors of ore generation: temperature distribution, phase separation zone, and fluid flow pattern. However, chemical factors that have a significant impact on the ore generation process, such as fluid chemistry, chemical compositions of fluids and rocks, and chemical reaction were out of our consideration. Aiming at clarifying more detailed ore generation processes by incorporating these chemical factors into the simulation, this study adopts a reactive transport simulation for the Hakurei site in the same hydrothermal field using TOUGHREACT software, which can handle coupled modeling of subsurface multiphase fluid and heat flow, solute transport, and chemical reactions.

A 2-D numerical model covering 2.0 km depth range by 4.5 km width below the seafloor was set for the simulation. Almost all model setting conditions were chosen based on the field observations (i.e., seafloor drilling, heat flux, temperature, and fluid geochemistry measurements). The left side of the model represented a fluid conduit (50 m width) where the hydrothermal fluid ascends from the bottom to the seafloor. The shallow subseafloor consists of alternating layers of pumiceous sediment, hemipelagic sediment, and low-permeable caprock. These layers were underlain by a volcanic basement. The initial and boundary conditions were set as shown in the figure.

To clarify a general pattern of chemical processes as a first step, simple geochemical settings were given by typical rhyolite composition of quartz (39%), oligoclase (26%), K-feldspar (22%), and annite (13%) for the initial rock. Secondary minerals by the water-rock interaction were anhydrite, pyrite, and sphalerite. Because of the limitation of the temperature range of the thermodynamic database (up to 300°C), the hydrothermal fluid source temperature was set as 300°C. Under those conditions, the numerical simulation was conducted for 5,000 years to understand mineral precipitation and accumulation over an order of thousands of years.

As a result of the simulation, the following features were derived: first, hydrothermal fluid ascended through the conduit and discharged from the conduit top, along with the lateral flow within the pumiceous sediments. Second, within the conduit and pumiceous sediments, low pH (approximately 4-5) and high concentrations of K+, H2S(aq), SiO2(aq), Zn+2, and Fe+2 were observed, due to the hydrothermal fluid flow, whereas SO4-2 was recharged to the shallow subseafloor by the seawater percolation. Third, when hydrothermal fluid ascended through the conduit, oligoclase and annite dissolved into the hydrothermal fluid, whereas K-feldspar, anhydrite, pyrite, and sphalerite precipitated from the hydrothermal fluid. Fourth, in the shallow subseafloor, laterally flowing hydrothermal fluid and percolating seawater were mixed and resultantly, pyrite and sphalerite precipitated primarily beneath the mound (the conduit top), with a lateral extension. In addition, anhydrite, quartz, and K-feldspar precipitated in the shallow subseafloor, with much precipitation near the conduit top. Those features generally correspond with the field observations (i.e., drilling and fluid geochemistry).

Because of the above-mentioned thermodynamic database limitation, our simulation temperature is slightly lower than the maximum observed temperature (326°C). Thus, our next step is to develop the thermodynamic database for a wider temperature range, e.g., 0-370°C, and incorporate the phase separation and its effect on fluid chemistry into the simulation.