地震危険度
地震危険度（じしんきけんど、英: seismic risk）とは、ある地点における地震の危険度（リスク）を表すものである。一定期間における最大の地震動、一定期間に一定基準以上の地震動がある確率など、様々な表現方法がある[1]。地図形式のものは地震ハザードマップ(seismic hazard map)とも呼ばれる。
なお、特に日本のように震源が陸地から遠い巨大地震が多いところでは、地震活動の多い少ないが必ずしも地震動の多い少ないには直結しないことから、地震危険度とは別に地震活動度（じしんかつどうど、英: seismic activity）という用語を用いて分けて考えることがある[2]。
更に、地震動だけではなく構造物の被害や損失についても扱う場合があり、1990年代以降多くなってきている[3]。
目次
1歴史と各国の状況
2時間依存と時間非依存
3地点ごとの評価と地震ごとの評価
4主な地震危険度評価
5地震危険度の値の意味と活用方法
6地震危険度評価の問題点
7出典
7.1注釈
7.2脚注
8参考文献
9関連項目
10外部リンク
歴史と各国の状況
近代地震学における初期の地震危険度として、河角広やアリン・コーネル(C. Allin Cornell)によるものが挙げられる。
河角は1951年に、日本の599年から1949年までの歴史地震342個の震源や規模などのデータを基に、一定の震度（震度5,6,7相当の加速度）の地震動が襲う平均間隔、また一定期間（75,100,200年）中の最大加速度を、日本列島の地図上に示して発表している[6]。この地図は「河角マップ」と呼ばれ、1950年に制定された建築基準法下の地震地域系数（1952年決定）に反映されている。しかし、河角マップは16世紀以前の資料が極端に少なく地域的な偏りも大きいほか、算出式にも問題があったことが指摘されている。その後、後藤・亀田(1968)も最大加速度を示した地図を作成し発表している。一方、村松(1966)、金井・鈴木(1968)、表ら(1975)、服部(1976)、尾崎ら(1978)は最大加速度ではなく最大速度の分布を求め発表している。特に金井・鈴木のものは河角と同じく75,100,200年の各期間における最大速度分布を示した比較性の良いもので、「金井マップ」と呼ばれている。
一方、アメリカでは、チャールズ・リヒターが1959年に初めてアメリカ全土を対象として簡易な地震危険度のゾーニングを行っているが、それ以前の1940年代にも簡易なものがあったという(Richter,1958)。これらはアメリカの耐震基準の設定の参考とされていた[7]。この後研究が進められ、コーネルが1968年に確率モデルに基づいた危険度評価を世界で初めて発表した。Milne・Davenport(1969)はカナダにおける歴史地震の資料を用いて最大加速度の分布をポアソン過程とみなして解析した地図を発表している。日本でもその後、ウェスノウスキー(Wesnousky)ら(1984)、島崎ら(1985)、亀田・奥村(1985)が震度や加速度を確率で表現した地図を発表している。また、USGSの協力でベイ・エリア自治体協議会（英語版）(ABAG)が1960年代からサンフランシスコ湾岸の危険度評価の検討を開始し、1980年代には危険度地図の発表に至っている[4]。
更にアメリカでは、USGS・CGS（英語版）・SCEC（英語版）の3者が共同で設立した"Working Group On California Earthquake Probabilities"(WGCEP)がカリフォルニア州の危険度地図の検討を続けている。1990年と1995年には、固有地震の繰り返し発生をモデル化した固有地震モデルを取り入れて長期的な地震発生確率を求め、発表した。Ward(1994)は観測により推定した地殻のひずみの進行率を地震モーメントの解放率に関連付けてモデルに取り入れる方法を提案したが、WGCEPは1995年にこの方式の地図も発表しており、StirlingおよびWesnousky(1998)もこれに続いている[7]。
日本でも1980年代に活断層のデータを考慮する動きが始まった。先述のウェスノウスキー(Wesnousky)ら(1984)、亀田・奥村(1985)の地図は活断層と歴史地震の両方を考慮している。構造物の設計における設計地震動に地震危険度評価を取り入れる動きも、主要構造物を皮切りにしてこの頃から見られるようになった。そのような中、1995年兵庫県南部地震（阪神・淡路大震災）以降、地震予知研究計画が見直された影響で国の方針としても地震危険度の評価に重点が置かれるようになった[3][8]。石川・奥村・亀田(1996)は、同地震において神戸市付近で観測された表面最大加速度(PGA)600-800ガルが再現モデルにおけるPGAの1,000年最大値約460ガルよりも大きいことから、構造物の設計において活断層データを考慮して1,000年以上の期間における想定を行うことで、内陸の活断層地震のような低頻度の大地震の評価効率が向上する可能性を指摘した[9]。他方、吉田・今塚(1998)、長橋・柴野(1999)は地盤による地震動増幅特性を加味した危険度評価を試みたほか、隈元(1999)、損害保険料率算定会(2000)、AnnakaおよびYashiro(2000)、宇賀田(2001)は時間依存モデルの設計を試みている[8]。こうした流れの中で、日本政府の地震調査研究推進本部は2002年に地震動の確率を示した「確率論的地震動予測地図の試作版（地域限定）」を発表、以後何度か改訂を行い、2005年には日本全域を対象とした「全国を概観した地震動予測地図」、2009年には地震動の確率と各断層（固有地震）毎の予想地震動を併せた 「全国地震動予測地図」を発表している[10]。
他の国でも同様の危険度地図が作成されている。旧ソビエト連邦では、Nersesov(1984)やSidorenkoら(1984)などが危険度地図を発表していることが知られているほか、力武の『地震予知 発展と展望』(2001)にはトルコや中国などの例も記載されている[7]。
時間依存と時間非依存
地震の発生に関する確率分布はポアソン分布と仮定して、ポアソン過程により算出する場合が多い。定常的かつランダムに発生している地震（例えば、無数の断層を有する領域内における地震の発生確率）を扱う場合、確率は定常ポアソン過程とグーテンベルグ・リヒターの関係式により表され、時間が経過しても変化しない。一方、発生確率が時とともに変化する地震（例えば、1つの断層や海溝における固有地震の発生確率）を扱う場合は、時間経過を織り込んだ非定常ポアソン過程により表される[1]。前者は時間非依存モデル、後者は時間依存モデルという。
時間依存モデルには、いくつかの手法がある。WGCEPが1995年に発表した評価では、一般的な時間予測モデルに対数正規分布のばらつきを加える手法が用いられた。しかし、このモデルでは、前回の地震からの経過時間があまりに長くなると逆に確率が低下してしまうという問題があった。これを防ぐ手法として、Matthews(1999)はBPT(Brownian Passage Time)分布を用いた評価法を考案した。BPTとは、震源における応力場の擾乱が地震や地殻変動などのブラウン運動により表現できる事に着目して、その擾乱の蓄積により大地震の発生に至るというプロセスをモデル化したものである[11]。地震調査研究推進本部は2001年にこれを用いた手法を開発し、以降の評価で継続的に用いている。
地点ごとの評価と地震ごとの評価
ロジャーズ・クリーク断層帯(Rodgers Creek Fault Zone)およびヘイワード断層帯(Hayward Fault Zone)北部の連動地震(M7.1)で予測される地震動の強さの分布図、ベイ・エリア自治体協議会(ABAG)、2003年発表
カリフォルニア州において30年間にM5を超える地震が発生する確率の分布図、全カリフォルニア地震破壊予測(UCERF 2)、2007年発表
トルコ及びイスラエルとその周辺の地震危険度の分布図、GSHAPによる
地震危険度の評価は、地震調査研究推進本部が用いている用語を引用すると、そのアプローチ方法から確率論的地震ハザード評価（「確率論的地震動予測地図」）とシナリオ型地震動評価（「震源断層を特定した地震動予測地図」）に分けられる。簡単に言えば、前者はある地点における地震動を評価したもので複数の地震の影響を受けるもの、後者はある特定の地震における地震動を評価したものである[12]。
両者は似たような確率値で表現されるが、混同しないよう注意が必要である。確率論的地震ハザード評価においては、ある一定の地震動（最大加速度や震度など）を閾値とした地震動の超過確率が用いられる。シナリオ型地震動評価においては、ある単一の地震そのものが発生する確率が用いられる[13]。これらの確率は、時間依存モデルの場合はハザードカーブという曲線をとる。