japan.quakes {ETAS}

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)