deviation_shmax function

Deviation of Observed and Predicted Directions of Maximum Horizontal Stress

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

deviation_shmax(prd, obs)

Arguments

• prd: data.frame containing the modeled azimuths of $SHmax$, i.e. the return object from model_shmax()
• obs: Numeric vector containing the observed azimuth of $SHmax$, same length as prd

Returns

An object of class data.frame

• dev.gc: Deviation of observed stress from modeled $SHmax$ following great circles
• dev.sc: Small circles
• dev.ld.cw: Clockwise loxodromes
• dev.ld.ccw: Counter-clockwise loxodromes

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)

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). tools:::Rd_expr_doi("10.1038/s41598-023-42433-2") .

See Also

model_shmax() to calculate the theoretical direction of $SHmax$.

Author(s)

Tobias Stephan