スポンサーサイト

上記の広告は1ヶ月以上更新のないブログに表示されています。
新しい記事を書く事で広告が消せます。

Gutenberg-Richter則

確率過程の理論的関心も高い、地震学分野におけるグーテンベルク-リヒター則(Gutenberg-Righter law)を、日本周辺における地震のオープンデータで確認する。
次のサイトからデータを得た。

地震加速度情報の検索ページ

1997年3月6日からデータが始まっている。その時点から現在までの16年分のデータを借りる。
地震発生時刻、震源、マグニチュードに注目する。

Gutenberg-Righter則を再現するには、震源が「日本周辺」であるデータに限って集計に入れる必要があり、例えばインド洋で発生した地震でも日本に揺れを及ぼす場合があるが、そういうデータを入れてはならない。

# Gutenberg-Righter law
data <- read.table("EarthquakeJapan.dat",header=F)
data <- data[data[,4]>=3.0,] #M3.0以上に限定
deletelist <- read.table("deletelist.dat",header=F)
data <- data[is.element(data[,3],deletelist[,1])==FALSE,]
temp <- as.matrix(data[,4])
temp <- as.numeric(temp)
plotdatay <- hist(temp,breaks=seq(3,9.1,0.1))[[2]]
plotdatax <- hist(temp,breaks=seq(3,9.1,0.1))[[1]]
plotdatax <- plotdatax[c(1:(length(plotdatax)-1))]
plotdatay <- plotdatay[c(1:47)]
plotdatax <- plotdatax[c(1:47)]
pdf("Gutenberg_Righter_law_1.pdf", height=6, width=8, pointsize=14)
plot(plotdatax, log(plotdatay/length(data[,1])),
xlim=c(3,9), ylim=c(-10,-2),
xlab="magnitude", ylab="log(N/Ntot)")
result <- lsfit(plotdatax,log(plotdatay/length(data[,1])))
abline(result,col="#DB0005", lwd=0.5)
dev.off()

result$coef
Intercept X
3.115510 -1.601156

Gutenberg_Righter_law_1_1.png



国外のデータを含めてしまうと、M7程度の大きな地震を余計に集計する結果となり、結果に異常が出る。

# Gutenberg-Righter law (default data)
data <- read.table("EarthquakeJapan.dat",header=F)
data <- data[data[,4]>=3.0,] #M3.0以上に限定
# deletelist <- read.table("deletelist.dat",header=F)
# data <- data[is.element(data[,3],deletelist[,1])==FALSE,]
temp <- as.matrix(data[,4])
temp <- as.numeric(temp)
plotdatay <- hist(temp,breaks=seq(3,9.1,0.1))[[2]]
plotdatax <- hist(temp,breaks=seq(3,9.1,0.1))[[1]]
plotdatax <- plotdatax[c(1:(length(plotdatax)-1))]
plotdatay <- plotdatay[c(1:47)]
plotdatax <- plotdatax[c(1:47)]
pdf("Gutenberg_Righter_law_2.pdf", height=6, width=8, pointsize=14)
plot(plotdatax, log(plotdatay/length(data[,1])),
xlim=c(3,9), ylim=c(-10,-2),
xlab="magnitude", ylab="log(N/Ntot)")
result <- lsfit(plotdatax,log(plotdatay/length(data[,1])))
abline(result,col="#DB0005", lwd=0.5)
dev.off()

Gutenberg_Righter_law_2_1.png



ちなみに消去したデータは以下:

台湾付近
東シナ海
ニューギニア付近
中東
朝鮮半島南部
?
中国西部
カムチャツカ半島付近
ニュージーランド付近
南極付近
南米西部
インドネシア付近
ウラジオストク付近
中米
北米西部
地中海
インド洋
フィリピン付近
中央アジア
インド付近
マリアナ諸島
グアム付近
南太平洋
中国東北部
アリューシャン列島
南大西洋
シベリア南部
北東太平洋
インドシナ半島付近
アフリカ東部
シベリア付近
オーストラリア付近
中国中部

+chronos
スポンサーサイト
Profile

zoa

Author:zoa
Metropolis (1927)

Calendar
11 | 2012/12 | 01
- - - - - - 1
2 3 4 5 6 7 8
9 10 11 12 13 14 15
16 17 18 19 20 21 22
23 24 25 26 27 28 29
30 31 - - - - -
Labels
Admin
Previous Posts
Recent Comments
Recent Trackback
Archive
このページのトップへ
上記広告は1ヶ月以上更新のないブログに表示されています。新しい記事を書くことで広告を消せます。