地球内部物理学のソースを表示

[[File:Slice earth.svg|thumb|right|240px|地球の内部構造]]
'''地球内部物理学'''（ちきゅうないぶぶつりがく）は[[地球]]内部を研究対象とする[[自然科学]]である。[[地球物理学]]の一分野に属する。

地球内部を直接掘削して調査することは困難を伴い、これまで地下10km内外を掘削したに過ぎず、内部を探求する方法は主に地球内部を透過する[[地震波]]の研究に依るところが大きい。

== 地球の慣性モーメント ==
[[密度]]が完全に均一な[[半径]] ''a''、[[質量]] ''M'' の[[球]]の[[慣性モーメント]] ''C'' は以下の式で表される。
: ''C'' = 0.4 ''Ma''<sup>2</sup>
一方で、地球の形は回転楕円体で近似され、さらに[[月]]や[[人工衛星]]の軌道の解析、および地球の[[歳差運動]]周期から三軸不等としてそれぞれ以下の慣性モーメントが求められている。これらは約 0.33 ''Ma''<sup>2</sup> と、密度が均一と仮定した場合より小さいことから中心に向かうにつれ密度が増大していることが判る。また均質な物質であっても地球内部では深度、圧力の増大に伴い圧縮され、密度は徐々に増大し均一とはならない。ここで ''C''軸がほぼ地球の自転軸となる。
: ''C'' = 0.330701 ± 0.000002 ''Ma''<sup>2</sup>
: ''A'' = 0.329615 ± 0.000002 ''Ma''<sup>2</sup>
: ''B'' = 0.329622 ± 0.000002 ''Ma''<sup>2</sup>

== 地震波の伝播速度 ==
[[地震波]]の地球内部の伝播速度は、物質の弾性定数、[[密度]] <math>\rho</math>、および非圧縮率 <math>K</math> に依存する。これらの剛性率および非圧縮率は体積変化およびねじれ変形に対する力学物性定数である<ref name="shimadzu">島津康男 『物理科学選書 地球内部物理学』 裳華房、1966年</ref>。

体積変化を伴う[[疎密波]]である縦波の ''P''波の伝播速度は以下の式で表される。ここで <math>\mu</math> および <math>\lambda</math> は[[弾性]]に関するラメの定数、特に <math>\mu</math> は[[剛性率]]と呼ばれる。
: <math>V_{\mbox{p}} = \sqrt{\frac{\lambda + 2\mu}{\rho}} = \sqrt{\frac{K + \frac{4}{3} \mu}{\rho}}</math>
また、体積変化を伴わない横波の ''S''波の伝播速度は以下のようになる。
: <math>V_{\mbox{s}} = \sqrt{\frac{\mu}{\rho}}</math>

[[File:Refraction of P-wave.PNG|thumb|right|240px|地震波の伝播経路]]
地震波の走時曲線からマントルの地震波速度分布を計算するために、震央距離 <math>\Delta</math> と、波線パラメーター <math>p = \frac{\mbox{d}T}{\mbox{d}\Delta}</math> および速度分布 <math>V(r)</math> の関係式であるHerglotz-Wiechert法が用いられる<ref name="Kanamori2">[[金森博雄]] 『岩波地球科学選書 地震の物理』岩波書店、1991年</ref>。
: <math>\Delta = 2p \int_{r_m}^{r_0} \frac{\mbox{d}r}{\sqrt{(r/\bar V)^2 - p^2}}</math>
この[[アーベル群|アーベル]]の[[積分方程式]]の解は走時解析に用いられ以下の式で表される。ここで <math>\bar V</math> は走時曲線の傾きから得られる見かけの[[角速度]]である。
: <math>\int_{0}^{\Delta_1} \cosh^{-1} \frac{\bar V \Delta_1}{\bar V \Delta} \mbox{d}\Delta = \pi \ln \left( \frac{r_0}{r_1} \right)</math>

=== 不連続面 ===
[[File:Earthquake wave paths.svg|thumb|right|240px|地球内部の地震波の伝播経路]]
[[File:Earthquake wave shadow zone.svg|thumb|right|240px|シャドウゾーン]]
地球内部を伝わる地震波である実体波の到達時刻と震央距離との関係である、走時曲線を描くことにより、地球内部の層状構造がわかる。地球内部を構成する物質の剛性率および非圧縮率は深度、圧力の上昇に伴い連続的に増大しその進路はカーブを描き、物質の構成が不連続的に変化する深度で地震波は屈折する。

[[1909年]]に[[アンドリア・モホロビチッチ]]は地下数kmから数十kmのところで ''P''波の伝播速度が不連続的に増大する面を発見し、これは[[地殻]]と[[マントル]]の境界である[[モホロビチッチ不連続面]]と呼ばれるようになった。

[[1906年]]にOldhamが震央距離100°を超える距離で ''P''波が急激に減衰するのを見出し、地球内部の地震波伝播速度の遅い領域が地震波を屈折させ影を形成していると唱えた。この震央距離103°- 143°の領域はシャドウゾーンと呼ばれている。[[1926年]]に[[ベノー・グーテンベルグ]]は深さ約2900kmのところで、''P''波の伝播速度が不連続的に減少し、''S''波が伝播しなくなる面を発見し、これはマントルと[[核 (天体)|核]]の境界である[[グーテンベルク不連続面]]と呼ばれるようになった。''S''波が伝播しないことは核（外核）の剛性率が0、すなわち[[流体]]であることを意味する。

[[1936年]]に[[インゲ・レーマン]]はシャドウゾーンに到達する回折波の存在から核の中に速度が7 - 9%程速くなる領域を見出し、剛性を持った固体からなると考えられ、外核と内核の境界である[[レーマン不連続面]]と呼ばれるようになった<ref>Lehmann, I., P' , Publ. Bur. Centr. Seismol. ''Int. Trav. Sci. Ser. A'', '''14''', 87, (1936).</ref>。

地球内部を伝わる実体波は地表面で反射し、''P''波が1回反射すると ''PP''あるいは ''PS''となり、後者はモード変換が含まれる。2回反射した波は ''PPP''と表記される。また外核の表面で反射した地震波は ''PcP''あるいは ''ScS''と表記し、核および地表で反射を繰り返した ''PcPPcP''あるいは ''ScSScS'' (''ScS''<sub>2</sub>)などの地震波も存在する。

また、外核を伝播する ''P''波は ''K''と表記され、外核を通過し再びマントルに出た地震波は ''PKP''と表記される。モード変換により、''PKS''、''SKS''および ''SKP''というフェーズも存在する。内核を通過する地震波は ''PKIKP''などと表記され、内核表面で反射し外核を通過する地震波は ''PKiKP''などと表記される<ref name="Terry">T.レイ・T.C.ウォレス 『地震学 上巻』 柳谷俊訳、古今書院、2002年</ref>。

=== モンテカルロ法による速度分布 ===
でたらめに選んだ多数のモデルで ''P''波および ''S''波の走時、地球の質量、慣性モーメントが説明できるか否かのテストを繰り返し、深さの関数としての ''P''波および ''S''波の速度分布、および密度分布の範囲が求められた。

このような手法を[[モンテカルロ法]]と呼び、その範囲は古典的方法により求められた速度分布および密度分布に一致する<ref name="Rikitake">[[力武常次]] 『固体地球科学入門』 共立出版、1994年</ref>。

=== インバージョン ===
''P''波および ''S''波の走時曲線、表面波の分散、地球震動スペクトルなどの観測結果から、剛性率および密度分布などの地球モデルのパラメーターを直接求める方法をインバージョンと呼ぶ。

''m'' 個の物理量の観測データ ''O''<sub>''j''</sub> (''j'' = 1, 2, … ''m'') と ''n'' 個の地球モデルパラメーター ''P''<sub>''i'' </sub>(''i'' = 1, 2, … ''n'') の場合、''P''<sub>''i''</sub> を与えると、
: ''C''<sub>''j''</sub> = ''Fi''(''P''<sub>''i''</sub>) , &nbsp; &nbsp;''j'' = 1, 2, … ''m''
によりこのモデルに期待される物理量 ''C''<sub>''j''</sub> が計算される。この計算値と観測値との差 ''O''<sub>''j''</sub> - ''C''<sub>''j''</sub> が極小と成るような[[データセット]]を選ぶことにより最適化されたパラメーターを求める。モンテカルロ法もインバージョンの一種である<ref name="Rikitake" />。

観測値の数 ''n''、モデルパラメーターの数 ''m'' である[[正方行列]]では連立方程式の数と未知数の数は等しく ''m'' = ''n'' のとき、行列 '''G''' に対する逆行列 '''G<sup>-1</sup>''' が存在すると仮定すれば、
: '''G<sup>-1</sup>''' d = ''m''
の式により ''m'' について解くことができる。

== 地球内部のパラメーター ==
[[File:Radial seismic velocity PREM.GIF|thumb|right|320px|地震波の伝播速度分布 ''V''<sub>p</sub>(''r''), ''V''<sub>s</sub>(''r'')]]
過去には和達ら<ref>Wadati ''et al''. (1933,34)</ref>、[[ハロルド・ジェフリーズ]]<ref>Jeffreys, H., Geophysical discussion : Observatory, '''59''', 267-269 (1936).</ref><ref>Jeffreys, H., The times of P, S and SKS and the velocities of P and S : ''M.N.-R.A.S., Geophys. Suppl.'', '''4''', 498-533 (1939).</ref><ref>Jeffreys, H., The times of core waves (2nd papre) : ''M.N.R.A.S., Geophys. Suppl.'','''4''', 594-615 (1939)</ref>および[[ベノー・グーテンベルグ]]<ref>B. Gutenberg (1958), Caustics Produced by Waves Through the Earth's Core, ''Geophysical Journal of the Royal Astronomical Society'', '''1''', 238–248.</ref>により ''P''波および ''S''波の深度に対する速度分布が求められた。[[:en:Keith Edward Bullen|Bullen]](1936)は地殻をA層、マントルをB層, C層, D層、外核をE層, F層、内核をG層に分け、速度分布、密度分布を求めている<ref>Bullen K. E. (1936): The variation of density and the ellipticities of strata of equal density within the Earth. ''MNRAS Geophys. Suppl''. '''3''', 395-401.</ref>。Jordan and [[:en:Don L. Anderson|Anderson]](1974)も速度分布および密度分布を求めている<ref>T. H. Jordan and D. L. Anderson : ''Geophys. J.'', '''36''', 411, (1974).</ref>。これらは何れも古典的方法によるものであった。

地球内部の密度、剛性率、圧力、地震波の伝播速度などの物性の分布を深さの関数として表現した地球モデルの一つにPREM (Preliminary Reference Earth Model)がある<ref name="PREM">Dziewonski, A. M., and D. L. Anderson, Preliminary Reference Earth Model, ''Phys. Earth Planet. Inter.'', '''25''', 297-456, 1981.</ref><ref name="PREM2">[http://geophysics.ou.edu/solid_earth/prem.html Preliminary Reference Earth Model PREM]</ref>。

PREM作成には膨大な数の実体波の走時、自由震動の固有周波数、表面波の減衰などの結果が用いられた。このPREMは海洋と大陸の[[リソスフェア]]について個別のモデルが求められているが、球対称から外れる不均質な分布を三次元的に表現するには至っていない。

=== 密度分布 ===
[[File:RadialDensityPREM.jpg|thumb|right|320px|地球内部の密度分布 ρ(''r'')]]
地球の平均密度は、[[キャベンディッシュの実験]]により得られる[[万有引力定数]]から[[地球質量]]を算出し、体積を用いてその平均密度ρ = 5.5 × 10<sup>3</sup> kg m<sup>-3</sup> が求められている。しかし、地球内部の密度は均一ではなく、慣性モーメントおよび、地球表面を造る[[岩石]]の密度ρ = 2.6-3.0 × 10<sup>3</sup> kg m<sup>-3</sup> から深部はより高密度な物質であることが窺われる。

地球内部のマントルなどは固体からなるが、地球サイズでみれば全体を液体と見做すことが可能で[[静水圧平衡]]が成立していると仮定される。
: <math>\frac{\partial P}{\partial \rho} = \frac{K}{\rho} = \phi(r)</math>

: <math>\frac{\mbox{d} \rho}{\mbox{d} r} = \frac{1}{\phi} \frac{\mbox{d} P}{\mbox{d} r}</math>

一方、弾性論的には以下の関係式が成立する。
: <math>\frac{K}{\rho} = V_{\mbox{p}}^2 - \frac{4}{3}V_{\mbox{s}}^2 = \phi(r)</math>
これらの式から一定組成の部分におけるある深度の密度変化が求まる。この式はAdams-Williamsonの式と呼ばれる<ref>L. H. Admas and E. D. Williamson : ''J. Wash. Acad. Sci.'', '''13''', 413, (1923).</ref>。''G'' は万有引力定数、''M'' は半径 ''r'' 内の質量である。
:: <math>\frac{\mbox{d} \rho}{\mbox{d} r} = \frac{-GM \rho}{r^2 \left( V_{\mbox{p}}^2 - \frac{4}{3}V_{\mbox{s}}^2 \right)}</math>
: <math>M(r) = M(R_0) + 4 \pi \int_{R_0}^{r} \rho r^2 \mbox{d}r</math>

マントル、外核の境界など不連続面における密度の差は直接これらの式から求めることは出来ないが、適当な境界条件を与えて連立微分方程式を解き、地球の質量や慣性モーメントが束縛条件となり最も適切と考えられる値 ''V''<sub>p</sub>(''r''), ''V''<sub>s</sub>(''r''), ρ(''r''), ''P''(''r''), ''g''(''r'') が定められる<ref name="Rikitake" /><ref name="jiten">[[宇津徳治]]、嶋悦三、吉井敏尅、山科健一郎 『地震の事典』 朝倉書店、2001年</ref>。

マントル上部、モホロビチッチ不連続面における密度はρ = 3.4 × 10<sup>3</sup> kg m<sup>-3</sup>、グーテンベルク不連続面ではマントル底部の密度はρ = 5.6 × 10<sup>3</sup> kg m<sup>-3</sup>、外核上部の密度はρ = 9.9 × 10<sup>3</sup> kg m<sup>-3</sup> と推定されている。また地球中心部の密度はρ = 13.1 × 10<sup>3</sup> kg m<sup>-3</sup>であり、外核も 9.9-12.2 × 10<sup>3</sup> kg m<sup>-3</sup> と高密度となっているが、これを常圧における密度に換算すると 7.0 × 10<sup>3</sup> kg m<sup>-3</sup> となり、[[鉄]]-[[ニッケル]]合金から予想される密度よりやや小さくなる。これは軽元素の混入を示唆している。

またマントルは均質というわけではなく、深度220km、400km、670kmにそれぞれ不連続面が存在すると仮定され、これらの不連続面は橄欖岩質物質の高圧における[[相転移]]が原因とされる。しかし実際にはマントル物質は純粋な[[ケイ酸マグネシウム]]から成るわけでなくPREMのモデルのようにシャープな不連続面ではなく、遷移帯と呼ばれる。

地球内部の高圧下では物質が高度に圧縮され、以下のような[[状態方程式 (化学)|状態方程式]]が成立する。ここで ''P'' は圧力、 ''K<sub>0</sub>'' は常圧における等温体積弾性率、ρは圧力''P''における密度、ρ<sub>0</sub> は常圧における非圧縮状態の密度である。この式は有限歪弾性論から導出される<ref>F. D. Murnagham : ''Amer. J. Math.'', '''59''', 235, (1937).</ref>。
: ''P'' = 3/2'' K<sub>0</sub>'' [(ρ/ρ<sub>0</sub>)<sup>7/3</sup> - (ρ/ρ<sub>0</sub>)<sup>5/3</sup>]

=== 圧力分布 ===
ある地点に於ける圧力 <math>P(r)</math> は、その上部にある物質の密度 <math>\rho</math>、[[重力加速度]] <math>g</math>、高さ <math>h</math> の積で表される。
: <math>P = \rho g h</math>
ある深度における圧力は地表からその深度までの間の密度と重力加速度との積を高さで積分して求められる。ここで <math>M</math> は半径 <math>r</math> 内の全質量である。
: <math>\frac{\mbox{d} P}{\mbox{d} r} = -\rho g = -\rho \frac{r M}{r^2}</math>
: <math>\frac{\mbox{d} M}{\mbox{d} r} = 4\pi \rho r^2</math>
このようにして求められた地球内部の圧力は深度2890mのマントル-外核境界（グーテンベルク不連続面）において1.36 × 10<sup>11</sup> [[パスカル (単位)|Pa]] (134万[[気圧]])、地球中心で3.64 × 10<sup>11</sup> Pa (359万気圧)である<ref name="PREM" />。

地球内部で物質が如何なる状態となるかを研究するのは高圧実験の分野であり、地球中心部の圧力を実験室で再現することは困難を伴う。[[ダイヤモンドアンビル]]などの実験装置により、漸く地球中心部の圧力を再現することが可能となった。また、[[衝撃波]]実験により物質の高圧状態を研究することも可能である<ref name="Rikitake" />。

=== 重力加速度分布 ===
[[File:EarthGravityPREM.svg|thumb|right|320px|地球内部の重力加速度分布 ''g''(''r'')]]
地球を完全な球体と見做した場合、ある深度における重力加速度 <math>g(r)</math> はその中心からの距離までの球状の天体と見做した場合の重力加速度となり、その深度より上部の物質は存在していないと仮定した場合と等しくなる。

従って、ある深度における重力加速度は、中心からその深度における密度分布に基づく半径内の全質量と[[万有引力|万有引力の法則]]から導かれる式により求められる<ref name="jiten" />。重力加速度はマントル内部では大きな変化は無く、グーテンベルク不連続面で最大値 10.68 m s<sup>-2</sup> となり、核内部では中心に向かって減少し、中心で0となる。

=== 温度分布 ===
地球内部は重力による断熱圧縮により[[内部エネルギー]]が増大し、深くなるにつれ断熱的に温度が上昇する<ref name="Rikitake" />。

均質系の断熱変化では以下の式が成立する。ここで<math>S</math>は[[エントロピー]]、<math>T</math>は[[絶対温度]]、<math>P</math>は[[圧力]]である。
: <math>\mbox{d}S = \left( \frac{\partial S}{\partial T} \right)_P \mbox{d}T + \left( \frac{\partial S}{\partial P} \right)_T \mbox{d}P = 0</math>
静水圧平衡の式からある深度における断熱勾配が求められる。ここで<math>\alpha</math>は[[熱膨張率|体積膨張係数]]、<math>C_\mbox{p}</math>は[[定圧比熱]]である。
: <math>\frac{\mbox{d} T}{\mbox{d} r} = - \frac{\alpha g T}{C_\mbox{p}}</math>
弾性論的に導かれる地震波伝播速度と密度との関係式から、断熱圧縮率<math>\chi_s</math>は ''P''波および ''S''波の伝播速度と以下の関係が成立する。
: <math>\frac{1}{\rho \chi_s} = V_{\mbox{p}}^2 - \frac{4}{3}V_{\mbox{s}}^2</math>

[[物質]]の[[融点]]は圧力の増大に伴い上昇し、マントルの2点の深度 <math>a</math> および <math>d</math> における融点の間には以下の関係がある<ref>R. J. Uffen : ''Trans. Am. Geophys. Union'', '''33''', 893, (1952).</ref>。ここで <math>\chi</math> は圧縮率である。
: <math>\frac{\tau_d}{\tau_a} = \frac{(\chi \rho)_a}{(\chi \rho)_d}</math>
核を構成する物質については半実験式として以下の式が得られている<ref>J. J. Gilvarry : ''J. Atmosph. Terrest. Phys.'', '''10''', 84, (1957).</ref>。ここで <math>\tau_0</math> は常圧における融点、<math>\alpha</math> は内部圧力に関する定数、<math>c</math> も定数である。
: <math>\frac{P}{\alpha} = \left( \frac{\tau}{\tau_0} \right)^c - 1</math>
マントルおよび内核は固体と考えられ、外核は液体である事実から凡その温度分布が求められ、地球中心部の温度は凡そ 6000 - 7000 [[ケルビン|K]] と推定されている。

また、地球を構成する物質中に含まれる[[放射性同位体]]の壊変により熱エネルギーが生産される。地球内部で発熱源として重要な放射性同位体は[[ウラン238]] (<sup>238</sup>U)、[[ウラン235]] (<sup>235</sup>U)、[[トリウム232]] (<sup>232</sup>Th)、および[[カリウム40]] (<sup>40</sup>K)である。これらは何れも[[花崗岩]]に濃縮されやすく、これらの熱源は何れも地球表面に集中している。地熱の元となる地殻熱流量の数十%は地殻内の熱源により賄われているという計算になる<ref name="Rikitake" />。

マントルは[[橄欖岩]]質の物質からなるが、時間的に緩慢な運動に対しては[[ニュートン流体]]として扱うことが可能である。この長時間のスケールで見ると流体として振舞うマントルは[[対流]]による熱輸送が可能となる。

=== 地磁気 ===
{{main|地磁気}}

== 地球の内部構造 ==
地球内部を構成する物質は地震波の速度分布から求められた、密度、剛性に当てはまり、且つ[[太陽系の元素組成|宇宙を構成する元素]]として多量に存在するものが主な材料として推定される。地球を構成する物質の平均[[原子量]] ''A''と密度ρおよび縦波の伝播速度 ''V''<sub>p</sub> との間にはBirchの法則と呼ばれる経験式が知られる<ref name="shimadzu" /><ref>Birch, F. (1961): CompositionComposition of the Earth's Mantle. ''Geophysical Journal of the Royal Astronomical Society'', '''4''', 295–311.</ref>。
: ''V''<sub>p</sub> = a + bρ + c(''A'' - 21) , &nbsp; &nbsp;(a = -1.87, b = 3.05, c = -0.80)

=== 地殻 ===
{{main|地殻}}
地殻の構造は近地地震および[[人工地震]]の走時からモホロビチッチ不連続面の深さの分布が求められている。人工地震では爆破時刻と位置を正確にコントロールできる点が長所である<ref name="Rikitake" />。モホロビチッチ不連続面より上部は ''P''波の伝播速度は6.0&nbsp;km s<sup>-1</sup>、下部は7.9&nbsp;km s<sup>-1</sup>前後となっている。

=== マントル ===
{{main|マントル}}
上部マントルは地震波の伝播速度や密度の条件を満たすものとして橄欖岩質の物質から成ると推定され、[[玄武岩]]の高圧相と考えられる[[エクロジャイト]]も構成する物質の候補とされる。

深さ70 - 200 kmの低速度層は[[アセノスフェア]](athenosphere)と呼ばれ剛性率がやや低く、その上部はリソスフェア(lithosphere)と呼ばれる。[[プレートテクトニクス]]において、プレートとは剛性率の高いリソスフェア（地殻 + マントル上部）を指す。

深さ650 kmから720km付近で[[橄欖石]] (Mg,Fe)<sub>2</sub>SiO<sub>4</sub> は[[ペロブスカイト構造]]に相転移し、下部マントルはペロブスカイト型(Mg,Fe)SiO<sub>3</sub> および[[塩化ナトリウム|岩塩]]型 (Mg,Fe)O から構成され、さらに[[スピネル]]型 (Mg,Fe)<sub>2</sub>SiO<sub>4</sub> および[[スティショバイト]]型の[[二酸化ケイ素]] SiO<sub>2</sub> が加わると推定される。

=== 外核 ===
{{main|核 (天体)}}
外核は液体であり、鉄-ニッケル合金に、[[硫黄]]、[[酸素]]、[[水素]]などの軽元素が合計で10%程度加わっているとされる。核-マントル境界面を決めるために、境界上側の反射波 ''PcP'' あるいは下側の反射波 ''PKKP'' および透過波 ''PKP''、''SKS'' あるいは変換波 ''PKIKP'' が用いられ、その境界面の起伏は数km以下であるとされている。外核は均質と考えられ、密度および ''P''波伝播速度は滑らかに増大する。

=== 内核 ===
内核は固体と考えられるが、[[ポアソン比]]は0.44程度と液体の0.5に近く、''S''波伝播速度は小さく 3&nbsp;km s<sup>-1</sup> 程度で剛性率が低いことを示唆する。

内核は地球の自転軸方向は赤道面内よりも地震波を速く伝えるため、異方性が存在することが判っている。この軸対称の走時異常は内核の対流に起因し、鉄の結晶軸を特定の方向に揃える役割を果たしていると考えられている<ref>Morelli, A., A. M. Dziewonski and J. H. Woodhouse, Anisotropy of the inner core inferred from ''PKIKP'' travel times, ''Geophys. Res Lett.'', '''13''', 1545-1548, 1986.</ref>。

== 地球内部の三次元構造 ==
PREMは地球内部の物性が深さのみの関数であり、水平方向は均質とした球対称のモデルであるが、現実には横方向にも不均質であり地震波の伝播速度は ''P''波は1%以内、''S''波は10%以内の不均一性が見られる。その地球内部の不均一性を三次元的にマッピングする方法を[[地震波トモグラフィー]]と呼ぶ。

地球を三次元的に細かい領域に分割し、ある領域の地震波伝播速度の観測値と、均一と仮定した地球も出るとの差をとることにより、三次元構造が空間的な関数として解明される。

地震波トモグラフィーから海洋プレートの沈み込みや、中央海嶺におけるマントルの上昇流が見出されている<ref>{{PDFlink|[http://docsrv.godac.jp/MSV2_DATA/12/be92_01.pdf 地球内部変動研究センター]}} 海底下から見えてきた巨大地震</ref>。

また、重力異常からも凡その事が推定され、プレート沈み込み地域（低温地域）で重力が大きく、湧き出しの地域（高温地域）で小さい傾向がある<ref name="Rikitake" />。

== 地球振動 ==
[[巨大地震]]の発生などにより、数日間にわたって地球の自由震動が観測される。この震動モードは[[球面調和関数]]を用いて表現され、3つの基本モードがある。
* バルーンモード - 半径が伸縮する
* フットボールモード - 半径方向に回転楕円体状に歪む
* 捩れモード - 球の表面に平行なせん断運動を伴う

地球表面を伝わる表面波は原理的にこれらの震動モードの重ね合わせとして理解され、周期の長い表面波は地球震動として扱うほうが便利である<ref name="Rikitake" />。

地球の震動モードは[[球座標]]系を用いて解析され、均一で自転していない球状の液体の[[波動方程式]]は以下のようになる<ref name="Terry" />。
: <math>\frac{1}{r^2} \frac{\partial}{\partial r} \left( r^2 \frac{\partial S}{\partial r} \right) + \frac{1}{r^2 \sin\theta} \frac{\partial}{\partial \theta} \left( \sin\theta \frac{\partial S}{\partial \theta} \right) + \frac{1}{r^2 \sin^2\theta} \frac{\partial^2 S}{\partial \phi^2} = \frac{1}{c^2} \frac{\partial^2 S}{\partial t^2}</math>
この波動方程式の解は以下のような[[ベッセル関数|球ベッセル関数]]で表される。
: <math>j_l(x) = x^l \left( \frac{-1}{x} \frac{\mbox{d}}{\mbox{d}x} \right) \frac{\sin x}{x}</math>

== 参考文献 ==
{{脚注ヘルプ}}
{{reflist}}

== 関連項目 ==
* [[プレートテクトニクス]]

{{地球}}
{{Earthquake}}

{{DEFAULTSORT:ちきゆうないふふつりかく}}
[[Category:地球物理学]]
[[Category:地震学]]