第32回日本情報地質学会講演会 講演要旨集 015-016頁 2021年

地球統計学による岡山市北西部日応寺花崗岩体中のα石英の構造シフトパラメター分布のモデリング
山川純次*

Modeling of a distribution of structural shift parameter of alpha-Quartz in
the Nichioji granite body on the Northwest part of Okayama city by theGeostatistics
Junji Yamakawa** 岡山大学 学術研究院 自然科学学域 Academic Field of Natural Science and Technology, Okayama
University, Japan, 3-1-1, Kita-ku ward, Tsushima-Naka, Okayama 700-8530, Japan. E-mail:ufeq0285@cc.okayama-u.ac.jp

キーワード：クリギング法，結晶構造解析，α石英，花崗岩，日応寺花崗岩体
Key words : Kriging method, Crystal structure analysis, alpha-Quartz, Granite, Nichiouji
granite body

１．はじめに岡山市北西部に産出する日応寺花崗岩体から分離したα石英の構造シフトパラメター(DI)の分布を地球統計学を用いてモデリングした.その結果，α石英の DI が岩体の南西部から北東部に向けて漸進的に増加する傾向を示すモデルが得られた.このモデルは岩体周辺の地質構造および地史と調和的であると考えられたので報告する.

２．データと手法
α石英を含む花崗岩のサンプルは日応寺花崗岩体の露頭から採集した．各サンプルから分離したα石英の DI はRIETAN-fp (Izumi and Monma, 2007) を使用して解析した．DI の空間構造はバリオグラムを用いて推定した．花崗岩体における DI の分布はドリフトを持つと仮定し普遍クリギング法 (Mathoron, 1973) を使用してモデリングした．
ドリフトには地質境界からのバッファ距離モデルおよび標高データを同時に用いた．予測グリッドデータとバッファ
距離モデルは 20 万分の 1 日本シームレス地質図 V2 (産総研地質調査総合センター, 2021) による岩体境界データから生成した．標高データは基盤地図情報数値標高モデル
10m メッシュ (国土地理院, 2021) から生成した．バッファ
距離モデルと標高データはα石英の DI データと予測グリッドデータそれぞれに結合した．全てのデータは日本測地系 2000 第 IV 系に投影した．クリギング法は R (Ikaha andGentleman, 1996) と gstat パッケージ (Pebsma, 2004)を用いて実行した．得られた構造シフトパラメターモデル
は Google earth (Google, 2021) で表示した．

第 1 図 DI と共変量の相関．左: DI と標高の相関．右: DI とバッファ距離の相関．

岩体内部の 1422 点について DI を普遍クリギング法で予測した結果, DI = 1.131〜5.723, Var = 0.1946〜0.3601 となった．DI 分布モデルを Google earth で表示した結果を
第 2 図に示す．

３．結果
日応寺花崗岩体のサンプル (19 個) から分離したα石英の DI は 2.06〜4.84 となった．α石英の DI と各ドリフトの相関を第 1 図に示す．

第 2 図 Googe earth で表示した DI 分布モデル．

16

４．議論
日応寺花崗岩体は岡山県中・南部分布する白亜紀花崗岩類である広島型花崗岩類に属する狭義の花崗岩である．北部および東部のルーフペンダントはペルム系と三畳系，夜久野岩類を含む先第三紀堆積岩類である 石英の斑状〜集斑状結晶が丸みを帯びており半深成岩である可能性が指摘されている (濡木・浅見・光野, 1979)．77.1 +/- 2.5 M.y.のK-Ar 黒雲母年代を示す (柴田, 1979)．西部は，より古いと推定されているトーナル岩から花崗岩までの岩相変化を示す塁帯深成岩体である日近花崗閃緑岩体と隣接している(多賀・村田・草地, 2004)．中央部では 34.0 +/- 2.7 Ma のFT 年代が得られている富良層および 27.4 +/- 1.9 Ma と27.1 +/- 1.5 Ma の FT 年代が得られている津高層と不整合する (鈴木, 2003)．日応寺花崗岩体の北端には日応寺衝上断層 (今村, 1966) が存在し岡山空港が位置している．
α石英の DI 標高と正の相関を示し岩体境界からのバッファ距離と負の相関を示した．
普遍クリギング法で得られた DI 分布モデルでは岩体南東部から北西部にかけて DI が漸進的に高くなった．これはこの方向にマグマが陥入し既存の岩体と衝突した結果，最先端である北西部に向けて DI が大きくなったと解釈できる．またこの陥入方向は日応寺断層が衝上型であることと調和的である．

文

献Google (2021) Google earth.
Ihaka, R. and Gentleman, R. (1996) R: a language for data
analysis and graphics. Journal of Computational andGraphical Statistics, vol.5, pp.299–314.
Izumi, F. and Momma, K. (2007) Three-Dimensional
Visualization in Powder Diffraction. Solid State
Phenomena, vol.130, pp.15–20.
Matheron, G. (1973) The intrinsic random functions andtheir applications, Advances in Applied of Probability,
vol.5, pp.439–468.
Pebesma, E. (2004) Multivariable geostatistics in S: thegstat package. Computers & Geosciences, vol.30,pp.683–691.
今 村 外 治 (1966) 岡 山 県 津 高 町 日 応 寺 よ り Miogypsina
kotoi HANZAWA の発見と岡山市周辺のいわゆる第三系の地質時代．岡山大学地学研究報告, vol.1, pp.1-10．
国土地理院(2021) 基盤地図情報数値標高モデル.
産総研地質調査総合センター(2021) 20 万分の 1 日本シーム
レス地質図 V2
柴田賢(1979) 東中国における花崗岩類の K-Ar 年代． 地質学論集, vol.17, pp.69-72.
鈴木 茂之, 檀原 徹, 田中 元(2003) 吉備高原に分布する第三系のフィッション・トラック年代, 地質学雑誌, vol.112,pp.35-49.
多賀優, 村田守, 草地功(2004) 岡山市北西部, 日近花崗閃
緑岩体の岩石記載および全岩化学組成, 岩石鉱物科学,vol.33, pp.33-45.
濡木輝一, 浅見正雄, 光野千春(1979) 岡山県中・南部の花崗岩類. 地質学論集・中生代酸性火成作用, vol.17, pp.35-46.