japan.quakes: A Japanese Earthquake Catalog
In ETAS: Modeling Earthquake Data Using 'ETAS' Model

japan.quakes

R Documentation

A Japanese Earthquake Catalog

Description

A data frame with 13724 rows and 6 columns giving occurrence date, time, longitude, latitude, magnitude and depth of shallow earthquakes (depth < 100 km) occurred since 1926-01-08 till 2007-12-29 in Japan and its vicinity (128-145E and 27-45N). Only earthquakes with magnitude greater than or equal to 4.5 are included.

Usage

data(japan.quakes)

Format

An object of class "data.frame" containing the following columns:

date Occurrence date in the format "yyyy/mm/dd"
time Occurrence time in the format "hh:mm:ss.ss"
long Latitude of epicenter in decimal degrees
lat Latitude of epicenter in decimal degrees
mag Magnitude of each earthquake by JMA (Japan Meteorological Agency)
depth Depth of each earthquake

Source

Data are retrieved from the Japan Meteorological Agency (JMA) by Jiancang Zhuang and acoppany the Fortran code at http://bemlar.ism.ac.jp/zhuang/software.html.

References

Zhuang J, Ogata Y, Vere-Jones D (2002). Stochastic Declustering of Space-Time Earthquake Occurrences. Journal of the American Statistical Association, 97(458), 369–380. doi: 10.1198/016214502760046925.

Zhuang J, Ogata Y, Vere-Jones D (2006). Diagnostic Analysis of Space-Time Branching Processes for Earthquakes. In Case Studies in Spatial Point Process Modeling, pp. 275–292. Springer Nature. doi: 10.1007/0-387-31144-0_15.

Examples

  # set the target polygonal study region
  jpoly <- list(long=c(134.0, 137.9, 143.1, 144.9, 147.8,
      137.8, 137.4, 135.1, 130.6), lat=c(31.9, 33.0, 33.2,
      35.2, 41.3, 44.2, 40.2, 38.0, 35.4))

  # creat an earthquake catalog
  japan.cat <- catalog(japan.quakes, study.start="1953-05-26",
      study.end="1990-01-08", region.poly=jpoly, mag.threshold=4.5)

  ## Not run: 
  plot(japan.cat)
## End(Not run)

  param0 <- c(0.592844590, 0.204288231, 0.022692883, 1.495169224,
  1.109752319, 0.001175925, 1.860044210, 1.041549634)

  ## Not run: 
  nthreads <- parallel::detectCores()
  japan.fit <- etas(japan.cat, param0, nthreads=nthreads)
## End(Not run)

ETAS documentation built on Nov. 28, 2022, 5:23 p.m.