model_shmax function

Theoretical Direction of Maximum Horizontal Stress in the geographical reference system.

Models the direction of maximum horizontal stress $SHmax$ along great circles, small circles, and loxodromes at a given point or points according to the relative plate motion in the geographical coordinate reference system.

model_shmax(df, euler)

Arguments

• df: data.frame containing the coordinates of the point(s) (lat, lon).

• euler: "data.frame" or object of class "euler.pole"

  containing the geographical coordinates of the Euler pole

Returns

data.frame

• gc: Azimuth of modeled $SHmax$ following great circles
• sc: Small circles
• ld.cw: Clockwise loxodromes
• ld.ccw: Counter-clockwise loxodromes

Details

$SHmax$ following great circles is the (initial) bearing between the given point and the pole of relative plate motion. $SHmax$ along small circles, clockwise, and counter-clockwise loxodromes is 90$degree$, +45$degree$, and 135$degree$

(-45$degree$) to this great circle bearing, respectively.

Examples

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

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

model_shmax(point, euler)

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

deviation_shmax() to compute the deviation of the modeled direction from the observed direction of $SHmax$. PoR_shmax() to calculate the azimuth of $SHmax$

in the pole of rotation reference system.

Author(s)

Tobias Stephan