schuster: Plot of a phasor walkout for the Schuster's test

In GRTo: Tools for the Analysis of Gutenberg-Richter Distributions of Earthquake Magnitudes

Description

This function plots a phasor walkout for the Schuster's test.

Usage

schuster(finame, title = "Schuster's diagram", color = c("black","red","orange",
"green","cyan","navy"), hd = FALSE, colidye = NULL, colidmo = NULL, colidda = NULL, 
colidho = 1, colidmi = 2, colidma = 3, colidz = NULL, utccor = 0, dayt1 = NULL, 
dayt2 = NULL, ysel1 = NULL, ysel2 = NULL, mosel1 = NULL, mosel2 = NULL, 
dasel1 = NULL, dasel2 = NULL, magsel1 = NULL, magsel2 = NULL, zsel1 = NULL, 
zsel2 = NULL, weekday1 = c("mo","tu","we","th","fr","sa","su"), 
weekday2 = c("mo","tu","we","th","fr","sa","su"))

Arguments

finame	name of the file to be loaded
title	main title that will appear in the phasor walkout plot
color	color of the phasor walkout (black/red/orange/green/cyan/navyblue). The first letter is enough to discriminate the color name
hd	logical. Whether finame contains headers or not
colidye	field number for the year values in the loaded file
colidmo	field number for the month values in the loaded file
colidda	field number for the day values in the loaded file
colidho	field number for the hour values in the loaded file
colidmi	field number for the minute values in the loaded file
colidma	field number for the magnitude values in the loaded file
colidz	field number for the depth values in the loaded file
utccor	correction for local time. This value and values of the start and end hours of nighttime are used to filter the daytime seismic noise from the walkout plot
dayt1	the start hour of nighttime in local time
dayt2	the end hour of nighttime in local time
ysel1	minimum year value for selection (Default is NULL, e. g. no selection)
ysel2	maximum year value for selection (Default is NULL, e. g. no selection)
mosel1	minimum month value for selection (Default is NULL, e. g. no selection)
mosel2	maximum month value for selection (Default is NULL, e. g. no selection)
dasel1	minimum day value for selection (Default is NULL, e. g. no selection)
dasel2	maximum day value for selection (Default is NULL, e. g. no selection)
magsel1	minimum magnitude value for selection (Default is NULL, e.g. no selection)
magsel2	maximum magnitude value for selection (Default is NULL, e.g. no selection)
zsel1	minimum depth value for selection (Default is NULL, e. g. no selection)
zsel2	maximum depth value for selection (Default is NULL, e. g. no selection)
weekday1	first day in the week for selection ("mo" for "Monday", "tu" for "Tuesday", "we" for "Wednesday", "th" for "Thursday", "fr" for "Friday", "sa" for "Saturday", "su" for "Sunday")
weekday2	last day in the week for selection ("mo" for "Monday", "tu" for "Tuesday", "we" for "Wednesday", "th" for "Thursday", "fr" for "Friday", "sa" for "Saturday", "su" for "Sunday")

Details

This function reads earthquake times (hour and minutes values are required) in an input file and plots a phasor walkout into a file with name file\_schu.png (png format file). Selections can be performed based on ranges in hour of the day, year, month, day, magnitude, depth and day of the week (in the same week). The correction for local time is such that GMT time + correction (e. g. utccor) = local time.

Author(s)

Daniel Amorese \& Paul A. Rydelek

References

Rydelek, P. A. and Hass, L. (1994) On Estimating the Amount of Blasts in Seismic Catalogs with Schuster's Method Bulletin of the Seismological Society of America, Vol. 84, No. 4, pp. 1256-1259.

Zurn, W. and Rydelek, P. A. (1996) Revisiting the phasor-walkout method for detailed investigation of harmonic signals in time series Surveys in Geophysics, Vol. 15, No. 4, pp. 409-431.

Examples

schuster(finame=system.file("extdata","IDYLLdeep.data.txt",package="GRTo"), 
title = "Schuster's diagram", hd = FALSE, colidye = 1, color="n", colidmo = 2, 
colidda = 3, colidho = 4, colidmi = 5, colidma = 15, colidz = 9, utccor = -9, 
dayt1 = NULL, dayt2 = NULL, ysel1 = 1983, ysel2 = 1990, mosel1 = NULL, 
mosel2 = NULL, dasel1 = NULL, dasel2 = NULL, magsel1 = NULL, magsel2 = NULL, 
zsel1 = NULL, zsel2 = NULL, weekday1 = NULL, weekday2 = NULL)

Example output

Year range 1981 1998 
1983  1990
Magnitude range 0 3.9 
Depth range 8.01 15 
[1] 573

Schuster's plot in Schuster's diagram_schu.png 

GRTo documentation built on May 2, 2019, 6:46 a.m.