Skip to contents

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

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). doi:10.1038/s41598-023-42433-2 .

See also

model_shmax() to calculate the theoretical direction of \(\sigma_{Hmax}\).

Author

Tobias Stephan

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)
#>     dev.gc    dev.sc dev.ld.cw dev.ld.ccw
#> 1 42.45436 -47.54564  87.45436  -2.545636