deviation_shmax {tectonicr}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 (σHmax\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 σHmax\sigma_{Hmax}, i.e. the return object from model_shmax()

obs

Numeric vector containing the observed azimuth of σHmax\sigma_{Hmax}, same length as prd

Value

An object of class data.frame

dev.gc

Deviation of observed stress from modeled σHmax\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). doi:10.1038/s41598-023-42433-2.

See Also

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

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)

[Package tectonicr version 0.3.2 Index]