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 (\sigma_{Hmax}
) along
great circles, small circles, and
loxodromes of the relative plate motion's Euler pole
Usage
deviation_shmax(prd, obs)
Arguments
prd |
|
obs |
Numeric vector containing the observed azimuth of
|
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). doi:10.1038/s41598-023-42433-2.
See Also
model_shmax()
to calculate the theoretical direction of
\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)