JpGU-AGU Joint Meeting 2017

講演情報

[EJ]Eveningポスター発表

セッション記号 M (領域外・複数領域) » M-GI 地球科学一般・情報地球科学

[M-GI29] [EJ] データ駆動地球惑星科学

2017年5月20日(土) 17:15 〜 18:30 ポスター会場 (国際展示場 7ホール)

[MGI29-P21] 交換モンテカルロ法を用いた角閃石固溶体の非理想パラメータの最適化

*岡本 敦1桑谷 立2上木 賢太3大森 敏明4福島 孝治5 (1.東北大学大学院環境科学研究科、2.国立研究開発法人 海洋研究開発機構、3.東京大学地震研究所附属高エネルギー素粒子地球物理学研究センター 、4.神戸大学、5.東京大学)

キーワード:amphibole solid solution, activity model, exchange Monte Carlo method

Minerals in rocks possess huge information on dynamics of the Earth’s interiors within their chemical compositions and textures. Recent progress of analytical instruments and computation abilities enable us to open new approaches to extract information from rocks by machine-learning techniques. In this study, we show the results of petrological optimizing problem by using exchange Monte Carlo method (Hukushima and Nemoto, 1996; Okamoto et al., 2015).
Estimates of P-T paths of metamorphic rocks are of special importance for understanding dynamic behaviors of subduction zones and crusts. Based on Gibbs’ phase rule, when the number of compositional variables of a mineral is larger than the degree of thermodynamic freedom of the system, we can obtain pressure and temperature solely from chemical composition of the mineral; this approach is called as Gibbs’ method (Spear 1993). Amphibole is one of the most suitable minerals for this analyses, because (1) amphibole is common in various metamorphic rocks, and (2) it has the complex compositions which compositional variables can compensate the degree of freedom, and (3) it commonly show a growth zoning, which is useful for drawing continuous P-T path during its growth. Okamoto and Toriumi (2001; 2004) applied this method to mineral assemblage of amphibole – plagioclase – epidote – chlorite - water in eight component system, in which thermodynamic freedom, f, is four. By using amphibole compositions, which is written as solid solution of six or seven components, they succeeded in decompression P-T paths from zoning of amphiboles in the Sanbagawa belt.
One of large problems on this thermodynamic analysis using amphibole is the lack of accurate activity model of amphibole due to its complex compositions. Even when we assume simple regular solution model (symmetrical), 21 non-ideal parameters (Margules parameters) are required for mixing of seven endmembers. In binary solution, the Margules parameter is the critical temperature of the solvus, but it is difficult to obtain solvus temperature accurately from complex solid solution like amphibole. Alternatively, Okamoto and Toriumi (2004) proposed the other approach for determining the Margules parameters by using the framework of Gibbs’ method as follows. Since the system described above (thermodynamic freedom = 4, independent amphibole component = 5; one amphibole composition is excess), we can estimate one compositional variable of amphibole itself as well as P, T and the composition of other minerals. Accordingly, we can compare the calculated and observed compositions. Using various amphibole compositions from natural metamorphic rocks, we can optimize the Margules parameters from numerous natural data. However, this optimization problem on such many parameters (21 W’s) would have several local minimums or broad minimum, and the results would be sensitive to the analytical errors of amphibole compositions. To overcome this hardly-relaxing problem, in this study, we construct the optimization scheme by using the exchange Monte Carlo Method (Hukushima and Nemoto, 1996). We succeeded the performance of the algorithm by using the synthetic data. We will discuss the applications to natural system and seek appropriate ways to obtain the realistic P-T paths from zoned minerals.
Hukushima K and Nemoto, K. (1996) J. Phys. Soc. Jpn, 65, 1604
Okamoto, A, Toriumi, M. (2001) Contrib Mineral Petrol 141, 268-286.
Okamoto, A., Toriumi, M., (2004) Contrib Mineral Petrol 146, 529-545.
Okamoto, A., Kuwatani, T., Omori, T., Hukushima, K., (2015) Phys Rev E 92, 042130