地形補正

火山など、地形の起伏がある場所での地殻変動モデリングにおいては、地形の影響が無視できない。有限要素法などの数値モデリングツールでは地形を入れるのは簡単だが、解析的なモデリングの場合、だいたい平面的な地球を仮定しているので、どうにかして地形の影響を取りこむ必要がある。

簡便な方法としては、Reference Elevation法（一番よく合うように、計算面の標高を決める）やVariable Depth法（観測点毎に、ソースの深さを変える）があるが、いろいろ問題を抱えており(Williams and Wadge, 1998）、摂動近似(Williams and Wadge, 2000)のほうがよりよいようだ（ただし急峻な地形には適用できない）。

Williams and Wadge (2000)のAppendix Bはftpでダウンロードできると書いてあるが、現在ftpサイトはなくなっているようだ。AGUの担当者にメールをすると送ってくれる。著者のWilliamsさん（現在GNS Science, New Zealand)にコンタクトすると、地形補正のプログラム(Fortran 77）を送ってくれた。多謝。これ、あらかじめいろいろな力源のルーチンが入っていてとっても便利（点楕円体も入っている！）。自分で力源を追加するのも比較的簡単。

もう15年前の論文だし、やっぱりこのレベルの近似はできるようになっておかないとなぁ。