63
課題 3
重力加速度の測定を通してみた地球
Section 3.1
はじめに
リンゴが木から落ちるのを観察して，ニュートンが万有引力の法則を思いついた話は有名である．質 量を持つものはお互い，まるで恋人どうしが寄り添うように，引き付けあう力が働くのである．サクラ の花びらがひらひら散るのも，熟した柿が落ちるのも，体重の増減に一喜一憂するのも，人工衛星ひま わりが地球の周りをはなれずにいるのも，ブラックホールに星が飲み込まれていくのも，重力（万有引 力）が支配する世界の出来事である．地球上にある物体の場合，強く引っ張りあう相手は地球そのもの で，その引力の方向はおおよそ地球の中心に向かう．ある地点で物体に働く重力はその質量に比例し， その物体が自由落下するときの加速度は質量によらないので，重力を表現するために重力加速度（g） をもちいる．ただし地球は自転をしていて，地球の自転軸に直角外向きの遠心力が働くため，測定され る重力加速度はほぼ地球からの引力と遠心力との合力になる（図 3.1）．さらに，重力加速度は地球中心 図 3.1: 地球の層構造: 重力加速度は地球内部の主として，金属核（内核・外核）からの引力と地球の自 転による遠心力とのバランスで決まる．また，地球は図のように層構造を呈していて，中心部分から， 固体金属の内核，液体金属の外核，柔らかな岩石からなるマントルと，硬い岩石からなる地殻によって 構成されている． からの距離の２乗に反比例するので，高い山に行けば重力加速度は小さくなる．したがって，地球上の重力加速度は測定場所（緯度や高度）によって値が違いうるのである．例えば，体重を軽くみせるため には赤道上の高い山に行けばよい．ただし，物質の重さを測る秤には各地の重力加速度の値を補正する 方法が記述されているので，簡単に体重をごまかすことはできない（補正値：北海道で g=9.805 m/s2， 沖縄県で g=9.791 m/s2）．また，高度 36,000 km 上空で地球の周りを一日一回まわり続けている衛星は， まさに遠心力と地球からの引力とがつりあった位置に存在しており，静止してみえる． 地球は層状構造をしていて，その最深部に密度の高い金属 (内核・外核)1が存在し，その周りを高温 で柔らかくなった岩石 (マントル)2が取り巻いている．さらに地球表層部は卵の殻のように岩石 (地殻)3 がその周りを取り囲み，その上で生命が育まれている．地球上の重力加速度の値が“ ほぼ ”9.8 m/s2であ ることは，地球を構成する物質のうち，主として密度の高い金属核によって決まっている．しかし重力 加速度を実測すると，上にも述べたとおり，地球の遠心力，測定場所の高さや局所的な地殻の構造（密 度変化）から影響を受け，場所により異なる値を示す．逆に，この特性を用いてわれわれが直接覗くこ とができない地球内部の密度構造を推定することができる． 重力加速度を測定するには，ニュートンの発想を利用して落体の加速度を直接観察する方法や，振り 子の等時性を利用して間接的に重力加速度を求める方法がある（図 3.2）．これらの方法では，重力加 図 3.2: さまざまな重力加速度の測定方法: a) りんごを落下させて，その時間から推定する方法（ニュー トンの発想）．b) 振り子の周期の等時性を利用して推定する方法（ガリレオの発想）． 速度の小数点以下 5 桁目の分解能で重力加速度の違いを検出できる．最近では，超伝導現象を利用した 超伝導重力計が開発され，15 桁目の分解能でその違いを検出できるようになってきている． 地球表面はヒマラヤ山脈のような高い山もあれば，日本海溝4のような深い海もあって，大変複雑な 形をしており，高さを定義するための基準面が必要である．一体どのように高さの基準面は決められて いるのだろう？この高さの基準面として，海水面を長期間平均して得られる平均海水面を仮想的に陸地 1金属核の平均密度は 10.0 × 103∼14.0 × 103kg/m3 2マントルの平均密度は 3.3 × 103∼5.8 × 103kg/m3 3地殻の平均密度は約 2.65 × 103kg/m3 4日本列島太平洋岸沖に南北に延びる水深 6000 m 以上のくぼみ地形で，海洋プレートが沈み込んでいる場所である．
3.1. はじめに 65 0 0 0 0 0 0 -50 50 50 100 50 100 150 -50 50 0 50 100 100 100 50 100 0 0 0 0 0 0 0 0 0 -50 -50 -50 -50 -100 -150 0 0 50 100 150 0 0 図 3.3: 日本付近の重力加速度異常図： 数字（単位は 10−5m/s2）の正負が重力加速度異常の正負に対応 し，その後の数字がその大きさである．一般にマントルの物質は地殻の岩石よりも密度が大きいことか ら，地殻の薄いところではマントル物質の影響を受け，補正後の重力加速度の異常は正になる．逆に地 殻の厚いところでは負になる． にも延長して地球表面全体を滑らかに覆った曲面，ジオイド5を採用し，このジオイド面からの高さに よって高度が表現されている．また，このジオイド面をもっともよく近似する回転楕円体を考えたとき， その楕円体の各緯度での計算された重力加速度の値を正規重力6といい，各緯度での標準的な重力加速 度を定義する．重力加速度を測定し，正規重力と比較する場合には，注意が必要である．測定された重 力加速度はジオイド面とは異なる有限の高度での測定値であるため，地表面とジオイド面との間にある 物体・空間の影響を取り除く補正作業が必要である（重力の補正7）．その補正後にもかかわらず，正 規重力との間に差が生じる場合，その差を重力異常と呼ぶ．この差は測定地点より内層の地球物質密度 の指標であり，地球内部の密度構造の推定に用いられる．図 3.3 に日本付近の重力加速度異常図を示す． この図は正規重力との差を等高線で表現したものであり，灰色の部分は重力異常が負の部分である． 5ジオイド：章末の参考 1 を参照 6正規重力：章末の参考 2 を参照 7重力の補正：章末の参考 3 を参照
この課題では，仙台における重力加速度を実 際に測定してみて，東北大学川内北キャンパス 学生実験棟の地下の密度について考察すること を目的としている． 重力加速度を測定するために，図 3.4 のよう なケーターによって改良された可逆振り子（ケー ター振り子）と呼ばれる振り子を使う． おもり B おもり C 9 ∆ 6 ? おもり C′ X X X y 副尺 H H H H Y h おもり A Q Q Q kXX調節ねじ X X y ▽9支点 S △XyXXX支点 S′ l0 6 ? lS 6 ? ld 6 ? lu 6 ? r 重心 G lS′ 6 ? 図 3.4 ケーター振り子
Section 3.2
実験装置の概要と使用方法
ケーター振り子は，図 3.4 のように支点 S と S′があり，おもり A が下の状態（順方向と呼ぶ）とお もり A が上の状態（逆方向と呼ぶ）に下げた時の周期が測定できるようになった剛体振り子である．剛 体振り子を使って重力加速度を求めるには一般に振り子の慣性モーメントを知る必要があるが，ケー ター振り子の実験ではその必要がない．これは，剛体振り子を順方向と逆方向に振らせた時の周期がわ かれば，慣性モーメントによらない次のような関係式が成り立つことを利用するからである8． 8章末の 3.5.4 参考４：ケーター振り子の理論的考察を参照．3.2. 実験装置の概要と使用方法 67 計算式 ケーター振り子の支点 S (S′) と重心の距離を lS(lS′)，支点間の距離を l0(= lS+ lS′)，振 り子の順方向（逆方向）の周期を Td(Tu)，重力加速度を g とすると，振り子の振れ角 α（ラジアン）が十分小さい時に (2π)2 g = [T2 d + T 2 u 2l0 + T 2 d − T 2 u 2(lS − lS′) ]( 1−1 8α 2+ O(α4)) (3.1) 重心 (lS の値) を正確に求めるのは難しいが，順方向と逆方向に下げた時の周期が等しくなるような 状況では，関係 (3.1) から lS − lS′を含む項が落ちて g= (2π) 2l T2 ( 1+1 8α 2) (T = T d = Tu) (3.2) と書き換えられる．周期 T 及び支点間の距離 l（l0の温度補正後の値）は，比較的簡単によい精度で 測定できるので，この関係を用いると非常に高い精度で重力加速度を測定できる．
●● 3.2.1 実験装置 ●●
本実験で使う装置は，ケーター振り子および光電スイッチ付きストップウォッチである．各部名称等 は，図を参照せよ．ケーター振り子は，図 3.4 のように可動なおもり A, B, C-C′があり，それらの位置 を適当に調整することにより順方向と逆方向に下げた時の周期が等しくできるように設計されている． 注意）本実験では，支点間の距離 l0（正確には，エッジ間の距離）の実測は行わない．l0はマイクロ メーターにより測定し，その測定値と誤差および測定時の温度が各装置毎に与えてあるのでその値と振 り子の熱膨張率θ(ただし，θ = 0.000019 deg−1とする) 及び当日の実験室の温度 TR(◦C) から l を求める． l= l0(1+ (TR− 15)θ) (3.3)●● 3.2.2 振り子とその振らせ方 ●●
振り子は，支点のエッジが壁に据え付けてある台の滑らかな面にまたがるように下げる．振り子を振 らせても，台に他の部分が接触しないように滑らかな面のほぼ真ん中に置く．この時，エッジが台など とぶつかって欠けることの無いように慎重におこなう．振り子は，光電スイッチの光源とセンサーの間 を通過する．振り子の前後の揺れ をできるだけ小さくするように注意して振らせる．∩ スイッチ @ @@R 9 9 9 9 9 9. 図 3.5 ストップウォッチ
●● 3.2.3 ストップウォッチの使い方 ●●
ストップウォッチは， 100 µsec ( 1 10000秒) 単位で測定 できるデジタル式時計で最大 99.9999 sec まで表示で きる (図 3.5)．時計自身の誤差は無視して考えてよい． ストップウォッチのスタート・ストップは，装置に つながっている光電スイッチにより行われる．時計本 体にあるボタン式スイッチを押すとリセットされ，ス イッチを離すと次に振り子が光電スイッチを横切った 時スタートする．時計は 10 周期後に自動的にストッ プするように設定されている．測定された 10 周期分 の時間を 1/10 にすることで，10 µsec（1/100000 秒） 単位で周期が求められる．●● 3.2.4 光電スイッチの設置 ●●
光電スイッチは，発光ダイオードの（赤外） 光源とセンサーからなり，その間を振り子が通 過して光が遮られたり当たったりすることで， スイッチのオン・オフが繰り返される (図 3.6)． ストップウォッチはその回数と時間を計測して いる．光電スイッチを設置する時は，振り子の 剣先部分の高さが順方向と逆方向で異なるので 注意を要する．測定を速やかに行うため，セン サーの高さは，振り子を 順方向に下げた状態 で スタンドを 台座 I に置いて 調整し，逆方向の 時は台座 I と II の両方を重ねて調整する． : 光 源 > 7 振り子の先端 -定規 @@R b ) センサー 図 3.6 光電スイッチSection 3.3
実験
本実験の目的は，Td= Tuとなる条件を実験で見いだし，その交点座標と理論式から学生実験棟での 重力加速度の絶対値を求めることである．そのために，下記実験 1，2 を通して，重力加速度の値を求 めてみよう．3.3. 実験 69 （全般的注意） • ふりこの振れ幅は 片側 3.0 cm(±1 mm が望ましい) とすること． • 周期の測定は，１回の測定について振り子を止めずに連続して５回おこない， その平均値を測定値とする． • データはすべて各自がその都度ノートに記録する（まとめて後で書き写すこ とはしないこと）． • ふりこを逆方向（または順方向）に設置する前に，TA または教員に許可を 取ること． • C −C’の副尺の読みは，77 cm，79 cm，81 cm，83 cm，85 cm を使う（後述）．
実験 1 予備実験：T
d= T
uとなる交点の予測 A
目的 予備実験では，Td= Tuとなる，おもり C− C’のおよその位置を予測するための実験を行う． 手順概略 1. Tdの値について，おもり C− C’の副尺の読みが，77 cm の場合と，85 cm の場合のそれぞれで求 める． 2. Tuの値について，同様に，おもり C− C’の副尺の読みが，77 cm の場合と，85 cm の場合のそれ ぞれで求める． 3. グラフを描く．縦軸を「振り子の周期」，横軸を「C− C’の位置」とする．測定した Tdと Tuの値 （各２点ずつ）をグラフにプロットし，Tdと Tuのそれぞれの値を結ぶ直線を定規で描き，Td= Tu となる T を目測で読み取る（以下，「交点予測 A」と呼ぶ）．図 3.7 を参照せよ． 4. 「交点予測 A」の読み取り値をノートに記録する．実験 2 T
d= T
uとなる交点の予測 B
目的 実験１によって，およその交点の予測値 A が得られた．次に，交点の予測精度を高めるため，お もり C− C’の副尺の読みが 79 cm, 81 cm, 83 cm のそれぞれの場合における Tdと Tuをそれぞれ追加し て計測して，Td= Tuとなる交点をより精度よく求めてみる． 手順概略 1. 逆方向（または順方向）について，C− C’の副尺の読みが，79 cm, 81 cm, 83 cm の場合の周期を それぞれ測定する．測定値を予測するために，測定ごとに５回分の測定値の平均値をグラフにプ ロットしてゆくことが望ましい． 2. 順方向（または逆方向）について，C− C’の副尺の読みが，79 cm, 81 cm, 83 cm の場合の周期を それぞれ測定する．測定値を予測するために，測定ごとに５回分の平均値をグラフにプロットし てゆくことが望ましい． 3. 実験中に配布される「データ報告用紙」に，必要事項を記入して，TA に提出する． 4. 「実験１ 予備実験」で作成した同じグラフ用紙上に，既に描いた直線を消さずに，新たに測定 した結果と合わせて新しい直線を Td, Tuとも定規で描き，新たに目測で交点を予測する（「交点 予測 B」とよぶ）． 5. 交点予測 B をノートに記録する．グラフは，実験１と実験２の区別がつくように工夫すること．3.3. 実験 71 実験２測定記録例 （表中の d= 順方向 (down) u= 逆方向 (up)） 測定 C-C′の位置 (cm) 回数 周期 (sec) 1 d77.0 1 2.00468 2 2.00471 3 2.00469 4 2.00470 5 2.00467 平均 2.00469 2 d79.0 1 2.00560 2 2.00545 3 2.00552 : : : 平均 2.00552 3 d81.0 1 1.79956× 2 2.00650 3 2.00654 4 2.00652 : : : 平均 2.00652 : : : 測定 C-C′の位置 (cm) 回数 周期 (sec) 1 u85.0 1 2.00976 2 2.00971 3 2.00974 : : : 平均 2.00974 2 u83.0 1 1.52235× 2 2.00721 3 1.90720 4 1.90722 : : : 平均 2.00721 3 u81.0 1 2.00444 : : : : 実験中の注意 • l0を補正するのに必要なので，当日の室 温をノートに記録する． • データ報告用紙は退出前に回収するので， TA が求めた最小二乗法の式（順方向，逆 方向）をノートに書き写すこと（すべて の桁をそのまま書き写すこと）． • グラフの縦軸は，できればグラフ用紙に 2 cm= 0.001 秒で描くこと．縦軸の最小 値と最大値は，実験１の結果から各自判 断せよ． 2.0010 2.0020 2.0030 2.0040 2.0050 sec 77.0 79.0 81.0 83.0 85.0 C-C′の位置 - cm r r r r !!!! !!!! !!!! !!! 図 3.7 実験２の交点測定のグラフ例
実験 3 予測 B の検証（時間が余った班のみ）
目的 予測を検証する．時間が余った場合，交点予測 B のおもり C−C’の副尺の位置で，実際に Td= Tu となるのか，検証してみよう（「検証値 A」とよぶ）． 手順概略 1. 交点予測Ｂのおもり C− C’の値に，副尺の位置を合わせる． 2. その C− C’の副尺の位置のまま，順方向と逆方向のそれぞれの Td，Tuの値を求める． 3. ここで得た Tdと Tuを引き算（∆T = |Td− Tu|）して，予測の精度を確かめよ． 実験 3 の注意 • 交点予測 B は目測で得た値なので，Tdと Tuがグラフ上で重なる程度の精度があれば十分である． • ここでいう「検証値 A」は，「C− C′の副尺の値，実測した Td，実測した Tu，∆T」の 4 種類の値 をセットとする．◆◆ レポートに最低限必要な項目・事項 ◆◆
レポート作成に際し, テキスト p. 10 の「レポートに関して」を参考にすること．この課題では以下の 指針も参考にする． 目的 • 重力加速度を測ることで何を知ろうとするのかを述べる． 原理 • ケーター振り子で重力加速度が測定できる原理に力点を置く（例えば振り子の等時性）．授業で は振り子や装置の使い方を詳しく述べているが，それらはここで期待する「原理」ではない． • 式を使うので，そのパラメータが何を示すのかわかるような図も含める． 方法 • 実験 1 と 2 を行なう目的に触れながら書く． • 振り子の各寸法（単位に注意），使用した振り子の番号，実験日の室温は，「方法」のなかの「実 験条件」にあたる．必ず記述すること．3.3. 実験 73 結果 • 実験データとして，以下の情報を提示すること． 1. 実験 1，実験 2，（実験 3）で測定したすべてのデータ（平均値も含む）．データは１周期に直 したものを示すこと． 2. 実験 1 と 2 で作成したグラフ（グラフ用紙に手書き．縦・横軸の項目と単位も忘れない．グ ラフの直線は定規で引く）．一つのグラフに，C−C′が 77，79，81，83，85 cm の場合の Td， Tuの値をプロットする． 3. 「実験 2 交点の予測 B」のデータを元に TA が算出した順・逆方向の最小二乗法の式．な お，自分自身で再計算してもよい（データ・式ともに小数点第５位までで行う事） 4. 交点予測 A，交点予測 B の値． • 実験中に気づいた点も含めて記載すると良い． 考察 • この課題では，問題を解答することで考察の流れが理解できるようになっているので，Section 3.4 の「問題 1」，「問題 2」，「問題 3」のすべてに解答すること．この順番で，記述すること． • ここでは解答に至る過程の妥当性を評価する．計算する問題については，第三者が検証できるよ うに計算過程を明記すること． • 単に計算するのではなく，計算結果から「何が言えるのか」を解説すること． • 単位変換などは細心の注意を払うこと．随所で単位をふること． • 実験条件，結果，考察の過程で気がついた事を書きたい場合，「考察」のなかに別途「実験全般に 関する考察」と小見出しをつけ，そこで述べる．かならず「根拠となるデータ」などを引用しつ つ「論拠」を明示しながら書くこと． 結論 • 何を知るためにこの実験を行ったのか，原理や方法にも簡単に触れ，得られた結果と考察（問題 1∼3）によってわかったことを簡潔にまとめる．
Section 3.4
問題
●● 3.4.1 問題 1 ●●
この設問では，予測精度の高い重力加速度を求めるために必要なことを学び，誤差を理解する． 1. 重力加速度を測定するのに必要な Td = Tuを知るために，予備実験（実験１）で交点があること を確かめ（交点予測 A），本実験（実験２）で交点の予測精度を高めた（交点予測 B）．しかし， これらの交点はいずれも目測で読み取った Td = Tuの値であり，「結果の再現性」で不安が残る． そこで，最小二乗法で計算した直線近似式（実験ノートに記載済み）を連立し，交点の座標を求 めよ（「交点予測 C」と呼ぶ）．座標のうち，Td= Tuとなる T の値は，小数点第 5 位まで書くこ と．また，計算過程を書くこと． 2. 「交点予測 C」が正しく計算されているとした場合，「交点予測 A」，「交点予測 B」はどれくらい 誤差があるのだろうか？「交点予測 C」と「交点予測 A」，「交点予測 B」の差をそれぞれ計算し， その値に基づいてズレを評価せよ．なお，交点予測 B は，C− C’が 77, 79, 81, 83, 85 cm の場合 の Td，Tuの値を用いて求めること． 3. 交点予測 C が正しく計算されているとして，式（3.2）に交点予測 C の T の値を代入して，重力 加速度 g( m/s2) の最確値を計算せよ．計算過程を示すとともに，g は小数点第 4 位まで記せ．な お，π = 3.1415 とする．（本来は，有効数字を考える必要がある（付録 A.2 有効数字，p. 271））． 4. 測定にはいくつかの誤差が考えられる．過去の調査から東北大学でもちいている実験装置では， 支点間の距離 l0の測定誤差がもっとも影響を与えることが分かっている．そこで，ふりこの支点 間の距離 l0の測定誤差分による重力加速度の最小値 g(−) と最大値 g(+) を求め，次のような形式 で重力加速度とその誤差を表せ（単位も明記）．計算過程も簡潔に示せ．誤差に関して，付録 A.3 誤差，A.4 誤算の伝播（p. 272∼274）と Section 3.5 参考５，６，７を参照せよ． g= 9.8003 ± 0.0002 m/s2 (3.4) この式は，真の重力加速度は，9.8001 m/s2∼ 9.8005 m/s2の間にある，という意味と同等である． 5. 「結果の再現性」を重視して，交点予測 C で重力加速度を求めた．交点予測 A や交点予測 B に もとづく T の値は利用出来ないほど誤差が伝播するのだろうか？交点予測 A と交点予測 B のそ れぞれについて，重力加速度の最確値をもとめよ（問題１-３と同じように計算する）．その値が， 問題１-４（一つ前の小問のこと）でもとめた誤差の範囲に収まるかどうか，差を計算して定量的 に記述せよ．「誤差の伝播」を詳しく知りたい人は付録 A.4（p. 274）を参照．3.4. 問題 75
●● 3.4.2 問題 2 ●●
この設問では，測定値と正規重力値とを比較するための事前準備として，測定値に補正を施し，補正 する意義を学ぶ． 重力は地下の密度だけでなく，標高（万有引力の法則で考えると，物体間の距離に相当する）によっ ても変化する．そのため標高を統一して比べないと，地下の密度を議論できない．そこで海抜 0 m（ジ オイド面）で期待される平均的な重力値（これが正規重力で，緯度だけから計算で求められる）と比較 を行う．東北大学自然科学総合実験棟の緯度は北緯 38 度 15 分 37 秒なので，Section 3.5 の参考２から 正規重力値 (gN) は 9.80016 m/s2である．正規重力値と比較するために，井戸を掘って海抜 0 m で重力 測定できればよいのだが，それは無理なので，本実験で得られた重力加速度に２種類の補正を施してジ オイド面上での値を得る必要がある（Section 3.5 の参考 3 を参照）．正規重力値と，標高による重力加 速度の変化を補正した値に差が生じるとき，その差を「重力異常値（∆g）」と呼ぶ．この差は，測定地 点近傍の地下の密度分布を反映している．重力異常値を求めるための考え方の流れを追ってみよう． 1. 測定値は標高 h= 67.8 m の丘の上での値である．まず，重力加速度の測定値（g）に，高度によ り減少した分を足すため，1 m 当たり 3.086 × 10−6 m/s2加味する補正値（g f:フリーエア補正） を次式（gf = 3.086 × 10−6h）より求めよ．つぎに海抜 0 m と測定地点の間に存在する岩石から の引力を差し引くため，その値をブーゲー補正 （gB = 2πGρh）により求めよ．地殻の平均密度 ρ = 2.650 × 103 kg/m3と万有引力定数 G= 6.672 × 10−11m3kg−1s−2を用いよ．それぞれの単位も 書くこと． 2. 正規重力値（gN）と上記２つの補正（gf，gB）の大きさから重力異常値を次式より計算せよ． ∆g = (g + gf− gB)− gN (3.5)●● 3.4.3 問題 3 ●●
この設問では，「誤差」と「有意差」の違いを学ぶ． 1. 問題 2-2 で得た重力異常値が有意と判断するためには，その絶対値が問題 1-4 で検討した「誤差」 より大きいことを確かめる必要がある．そこで，問題 2-2 で求めた重力異常値と，問題 1-4 で計 算した式（3.4）の誤差とを比較し，定量的に記述せよ．そのうえで，有意な重力異常値が得られ ているか，理由を明示しつつ論ぜよ． 2. 有意な重力異常値が得られている場合，重力異常を議論することが可能である．議論可能な場合， 自然科学総合実験棟の地下深部にある岩石密度が，地殻の平均密度（ρ = 2.650 × 103kg/m3）よ り大きいか，あるいは小さいかを推定してみよ．推定にあたり，その判断となった根拠も明記せ よ．なお，重力異常を含むのは測定値 g である．従って「重力異常を小さくするにはブーゲー補 正のρ が大きければ（あるいは，小さければ）良い」といった議論は誤りである．Section 3.5
参考
●● 3.5.1 参考１：ジオイド ●●
われわれは高度をあらわすとき，「この場所での高度は海抜 200 m である」のように海水面を基準にし ている．地球表層を覆っている流体（海水）は自由に移動できるため，この海水面を求めることによっ て，「地球の形」そのものを決定できる．海水の表面は，波や干満などにより絶えず変動しているが，長 い年月の平均をとれば世界中を覆う滑らかな曲面になる．これを平均海水面と呼び，この平均海水面を 陸地にも仮想的に延長し，地球表面を流体で覆ったときの曲面を“ ジオイド ”とよぶ．山の高さなど各 地の高さは，このジオイドからの高さであらわされる．当然ジオイドはその性質上，重力の方向に垂直 で，人工衛星もこのジオイド面と平行に高度数万 m を飛行している．またジオイドにもっともよく合 うように決められた回転楕円体を地球楕円体と呼び，この楕円体上での重力加速度を基準として正規重 力を決定している．●● 3.5.2 参考２：正規重力 ●●
地球の重力加速度はだいたい地球の引力と自転の遠心力の和であり，緯度と高度で重力加速度の値は ほぼ決まる．したがって，ジオイド面上のある地点での重力加速度の値は，回転楕円体の扁平率・赤道 半径・極半径が決まれば自動的に計算によって決定されることになる．この正規重力 (gN) を与える国 際正規重力式は，1984 年に改訂され， gN= 9.7803267714( 1+ 0.00193185138639 sin2λ √ 1− 0.00669437999013 sin2λ ) (3.6) である．(λ は緯度をラジアン単位であらわしたもの）●● 3.5.3 参考３：重力の補正 ●●
• フリーエア補正 重力加速度は地球の中心から遠ざかると小さくなる．測定値をジオイド面上の 値に計算しなおす作業がフリーエア補正である．高度 h の点で測定した重力加速度を g(= g0− gf)， ジオイド上での値を g0，ジオイド面を半径 R の球とすれば， g g0 = R2 (R+ h)2 ( ．．．mg 0= G mM R2 , mg = G mM (R+ h)2) (3.7) となる (m: 地球上での物質の質量，M: 地球の質量，G: 万有引力定数）．したがって，測定値 g を ジオイド面上の重力値 g0に引き戻すと，g0=g+3.086×10−6h となる．すなわち，フリーエア補正 の値は gf = 3.086×10−6h で与えられる． • ブーゲー補正 フリーエア補正ではジオイドと測定点との間に物体が何もない空間を考えていた が，実際は地殻岩石が存在するので，その岩石の質量による影響を補正する必要がある．ジオイ3.5. 参考 77 ドと測定点との間に密度ρ，厚さ h の無限に広い平板を敷き詰めたとして，その分の引力を差し 引くことによってこの補正を行なうことにする．密度ρ の無限平板による引力 gBは， gB= 2πGρh (3.8) である． • 重力異常 測定値 g にフリーエア補正（高度補正）とブーゲー補正（密度補正）を施し，ジオイ ド面上での重力値に補正した値は，その地点の地下の密度が平均的ならば，その緯度の正規重力 gNと等しいはずである．しかし，実際は地下の密度が平均からずれていて，補正後の重力値は正 規重力と一致しないことが多い．この補正された重力値から正規重力の値を差し引いた値∆g ∆g = (g + 3.086 × 10−6h− 2πGρh) − g N (3.9) を重力異常という．この異常値は測定点の地下に存在する物質の密度を反映している．∆g > 0 な ら高密度な物質，∆g < 0 なら低密度な物質の存在を示唆している．このことは石油やダイヤモン ドの探査に利用されたり，火山噴火予測のための常時観測にも利用される場合がある．
●● 3.5.4 参考４：ケーター振り子の理論的考察 ●●
実験で見たように，ケーター振り子は支点 S と S′があり順方向（おもり A が下の状態）と逆方向 （おもり A が上の状態）に下げた時の周期が測定できるようになった剛体振り子である．順方向に下げ た時，振り子の質量を M，慣性モーメントを I0，支点と重心の距離を lS，振り子の振れの角度をθ と すると，運動方程式は， d2 dt2θ(t) = −ω 2 0sinθ(t) (3.10) で与えられる．ここに ω2 0= MglS I0+ Ml2S (3.11) である．振り子の振れθ が十分小さいとして三角関数をテーラー展開すると運動方程式は d2 dt2θ(t) = −ω 2 0 ( θ(t) −1 6θ(t) 3+ O(θ5)) (3.12) この運動方程式を振れ角α が小さいとして α4以上を無視する近似で解くと（このような方法を摂動法 という）その近似解は θ(t) = α(sinωt − 1 192α 2 sin 3ωt + O(α4)) (3.13) ω = ω0 ( 1− 1 16α 2+ O(α4)) (3.14) で与えられることがわかる．よって順方向に下げた時の振り子の周期は Td= 2π √ I0+ Ml2S MglS ( 1+ 1 16α 2− O(α4)) (3.15)この関係式をつかって重力加速度を求めるには慣性モーメントを知る必要があるがケーター振り子のよ うな複雑な物体では，慣性モーメント I0を高い精度で求めることは非常に難しい．ところが，ケーター 振り子では２つの周期 Tdと Tuが測定できる．Tuに関しては， Tu= 2π √ I0+ Ml2S′ MglS′ ( 1+ 1 16α 2− O(α4)) (3.16) が同様にして導かれる．この２つの式 (3.15) と (3.16) を組み合わせて慣性モーメント I0を消去すると 定理の関係式 (3.1) が得られる．
●● 3.5.5 参考５：誤差評価の考え方と各種補正の意味 ●●
本実験において，重力加速度は (3.2) 式 g=(2π) 2l T2 ( 1+1 8α 2) (3.2) に，温度補正後の支点間の距離 l，周期 T ，振れ角α を入れることで与えられる．これらの値はいずれ も誤差を含み，それがわかれば次に述べる方法により，最終的に g がもつ誤差を求められる．しかし T やα の誤差による影響は，l に比べてはるかに小さいのと，計算が煩雑なので，問題 1-4 では l の誤差だ け考えればよいことにする（順方向と逆方向で振り子長が変わることによるα の変化すら，l に比べる と影響が小さい）．この場合，次の参考６で述べる偏微分を行わなくとも，支点間の距離として l− δl と l+ δl を (3.2) 式に代入することで，g の誤差を見積もることができる（もちろん，参考６で変数を一 つとすれば同じ結果が得られる）． ここでいう l の誤差（δl）とは，机の上にある表で，l0= 999.14 ± 0.02 mm などと書かれているうちの 0.02 mm のことであり，(3.3) 式の温度補正のことではない．金属の熱膨張率（1◦C の違いで長さが変化 する割合）は正確に求められており，それを用いることで任意の温度での長さを正確に予測できる．例 えば 15◦C の時に比べると，22◦C では+0.000019 × (22 − 15) = 1.000133 倍の長さになる（999.14 mm だったのが，999.27 mm になる）．この変化を誤差と混同してはいけない．なお，±0.02 mm の部分も 熱膨張の影響を受けるが，桁数を考えるとその変化は無視して構わない（0.02 × 1.000133 を計算して 見よ）． 測定値を正規重力と比較可能にするために行うフリーエア補正とブーゲー補正の大きさは，測定地点 の標高に比例する．つまり，海抜 0 m で測定した場合の補正量は 0 で，標高が高いほど補正量は大きく なる．従って，「補正量が小さいので重力異常を議論できる」だとか，「補正が必要だったので重力異常 は議論できない」などと考えるのは全くの誤りである（この理屈だと，重力異常を議論できるのは海抜 0 m の場所だけになる）． 測定値に補正を施した重力値と正規重力の差が，見積もった誤差の範囲内に入らないとしたら，それ は測定精度が悪いせいではない（求めた誤差こそが，精度の悪さの見積もりなのだから）．地下の物質 に原因があって，そこでの重力が平均と異なると考えねばならない．これが重力異常の意味である．3.5. 参考 79
