deviation_shmax: Deviation of Observed and Predicted Directions of Maximum...
In tectonicr: Analyzing the Orientation of Maximum Horizontal Stress

deviation_shmax

R Documentation

Deviation of Observed and Predicted Directions of Maximum Horizontal Stress

Description

Calculate the angular difference between the observed and modeled direction of maximum horizontal stresses (\sigma_{Hmax}) along great circles, small circles, and loxodromes of the relative plate motion's Euler pole

Usage

deviation_shmax(prd, obs)

Arguments

`prd`	`data.frame` containing the modeled azimuths of `\sigma_{Hmax}`, i.e. the return object from `model_shmax()`
`obs`	Numeric vector containing the observed azimuth of `\sigma_{Hmax}`, same length as `prd`

Value

An object of class data.frame

dev.gc: Deviation of observed stress from modeled \sigma_{Hmax} following great circles
dev.sc: Small circles
dev.ld.cw: Clockwise loxodromes
dev.ld.ccw: Counter-clockwise loxodromes

Author(s)

Tobias Stephan

References

Stephan, T., Enkelmann, E., and Kroner, U. "Analyzing the horizontal orientation of the crustal stress adjacent to plate boundaries". Sci Rep 13. 15590 (2023). \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1038/s41598-023-42433-2")}.

Examples

data("nuvel1")
# North America relative to Pacific plate:
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")

# the point where we want to model the SHmax direction:
point <- data.frame(lat = 45, lon = 20)

prd <- model_shmax(point, PoR)
deviation_shmax(prd, obs = 90)

tectonicr documentation built on Sept. 11, 2024, 6:05 p.m.