スラブ内部の相転移に伴う応力場と深発地震の発生の場との関連性

板体内相变应力场与深部地震发生位置的关系

• 批准号: 09740356
• 负责人: 吉岡 祥一
• 金额: $ 1.41万
• 依托单位: Kyushu University
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for Encouragement of Young Scientists (A)
• 财政年份: 1997
• 资助国家: 日本
• 起止时间: 1997 至 1998
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/en/grant/KAKENHI-PROJECT-09740356/
• 关键词: オリビン; 相転移; 正・負の浮力; スラブ; 温度分布; 有限要素法; 応力場; 深発地震; 物性パラメター

本研究では、物性値の空間分布を考慮して、スラブ内の応力場の数値シミュレーションを行った。スラブ構成物質をオリビン(Mg0.89Fe0.11)2SiO4と仮定し、差分(FD)法を用いてスラブ内の温度分布を計算し、各タイムステップごとに物性値の空間分布を反映させた。スラブ内部の応力場に影響を与える要素として、スラブと周囲のマントルの密度差により生じる正・負の浮力、沈み込むスラブの上面・下面に働く剪断力、スラブ先端がマントルから受ける抵抗力、熱応力及び不均質体積変化により生じる応力を考慮し、有限要素法(FEM)を用いて弾性応力解析を行った。正・負の浮力に起因する応力は、α相内のdown.dip tension、γ相内に卓越するdown-dip compression、γ→Pv+Mw相転移境界付近のdown-dip tensionで特徴づけられ、最大剪断応力はオリビン準安定相境界及びγ相内深部に集中し30MPa程度を示した。剪断力に起因する応力場は、α相内では圧縮場が、γ相内ではdown-dip tensionが現れ、最大剪断応力の大きさは、数MPaと非常に小さかった。スラブ先端が受ける抵抗力に起因する応力場では、down-dip compressionが卓越し、最大剪断応力はスラブ先端部で最大値を示した。熱応力及び不均質体積変化により生じる応力は、各相転移境界で応力集中が見られた。深発地震発生数と解析結果のdown-dip方向の応力を比較した結果、正・負の浮力による応力が深発地震のピークと一致し、700km以深では最大剪断応力がゼロに近づき、サイスミシティが見られないことがうまく説明されることにより、考慮した5つの要素の中で最も良い一致を示すことがわかった。また、正・負の応力に起因する応力場とハーバード大学のグループによるCMT解との比較を試み、これについてもスラブと周囲のマントルの密度差による正・負の浮力が地震学的観測データを最もうまく説明する支配的な要素であると結論づけた。

This study で は, physical numerical の spatial distribution を consider し て, ス ラ ブ within の 応 force field の the numerical シ ミ ュ レ ー シ ョ ン を line っ た. ス ラ ブ composition を オ リ ビ ン (Mg0.89 Fe0.11) 2 sio4 と 仮 し, finite difference method (FD) を with い て ス ラ ブ を の temperature distribution in computing し, each タ イ ム ス テ ッ プ ご と に property numerical の spatial distribution を reflect さ せ た. ス ラ ブ internal の 応 field に effects を and え る elements と し て, ス ラ ブ と weeks 囲 の マ ン ト ル の density difference に よ り raw じ る is negative の buoyancy, shen み 込 む ス ラ ブ の above, the following に 働 く shearing force, ス ラ ブ apex が マ ン ト ル か ら by け る resistance, thermal 応 force and volume - び heterogeneous に よ り raw じ る 応 を consider し, finite element method ( FEM)を resolve を rows った using the 応 elastic 応 force. Positive and negative <s:1> buoyancy に are caused by する応 force <e:1>, <s:1> down.dip tension in the α phase, に excellent するdown-dip compression in the γ phase, and 転 displacement of the γ→Pv+Mw phase to approach <s:1> down-dip Recogniton で, 徴 づ け ら れ, maximum shearing force 応 は オ リ ビ ン quasi stable phase boundary and び gamma phase within the deep に concentration degree し 30 mpa を shown し た. Cause shearing force に す る 応 force field は, alpha で は が 圧 shrinkage field, gamma phase で は down - dip recogniton が れ now, maximum shearing 応 の potent き さ は, small number of MPa と very に さ か っ た. ス ラ ブ apex が by け る に resistance causes す る 応 force field で は, down - dip compression will が し excellence, maximum shearing force 応 は ス ラ ブ apex of で numerical を し indicated the largest た. The thermal 応 force and び heterogeneous volume change によ biomass じる応 force <e:1>, as well as the phase 転 transition boundary で応 force concentration が can be found in られた. Deep 発 earthquake 発 born several analytical results の と down - dip direction の 応 force を compare し た results, positive, negative の buoyancy に よ る 応 force が deep 発 earthquake の ピ ー ク と consistent し, 700 km to deep で は maximum shearing force 応 が ゼ ロ に nearly づ き, サ イ ス ミ シ テ ィ が see ら れ な い こ と が う ま く illustrate さ れ る こ と に よ り, consider し た 5 つ の Among the equations, the で is the most consistent with the を, which is shown in the す とがわ った った った. ま た, positive, negative の に 応 force causes す る 応 force field と ハ ー バ ー ド university の グ ル ー プ に よ る CMT solution と の is を み, こ れ に つ い て も ス ラ ブ と weeks 囲 の マ ン ト ル の density difference に よ る 観 measurement of positive, negative の buoyancy が seismology デ ー タ を most も う ま く illustrate す る な elements of dominant で あ る と conclusion づ け た.