PoR_azi {tectonicr} | R Documentation |
Azimuth Conversion from Geographical to PoR Coordinate Reference System
Description
Transforms azimuths and models the direction of maximum
horizontal stress
\sigma_{Hmax}
in the Euler pole (Pole of Rotation)
coordinate reference system. When type of plate boundary is given, it also
gives the deviation from the theoretically predicted azimuth of
\sigma_{Hmax}
, the circular distance, and the normalized
\chi^2
statistics.
Usage
PoR_azimuth(x, PoR, axial = TRUE)
PoR_shmax(x, PoR, type = c("none", "in", "out", "right", "left"), axial = TRUE)
Arguments
x |
|
PoR |
|
axial |
logical. Whether the azimuth is axial (0-180) or directional (0-360). |
type |
Character. Type of plate boundary (optional). Can be
|
Details
The theoretical azimuth of \sigma_{Hmax}
in the pole of
rotation reference system is
0 (or 180), 45, 90, 135 degrees if the stress is sourced by an
outward, sinistral, inward, or dextral moving plate boundary, respectively.
directions of \sigma_{Hmax}
with respect to the four
plate boundary types.
Value
PoR_azimuth
returns numeric vector of the transformed azimuth in
degrees.
PoR_shmax
returns either a numeric vector of the azimuths in the
transformed coordinate system (in degrees), or a "data.frame"
with
azi.PoR
the transformed azimuths (in degrees),
prd
the predicted azimuths (in degrees),
dev
the deviation between the transformed and the predicted azimuth in degrees (positive for clockwise deviation of observed azimuth wrt. predicted azimuth),
nchisq
the Norm
\chi^2
test statistic, andcdist
the angular distance between the transformed and the predicted azimuth.
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 compute the theoretical direction of
\sigma_{Hmax}
in the geographical reference system.
deviation_shmax()
to compute the deviation of the modeled direction
from the observed direction of \sigma_{Hmax}
.
norm_chisq()
to calculate the normalized \chi^2
statistics. circular_distance()
to calculate the angular distance.
Examples
data("nuvel1")
# North America relative to Pacific plate:
PoR <- subset(nuvel1, nuvel1$plate.rot == "na")
data("san_andreas")
res <- PoR_shmax(san_andreas, PoR, type = "right")
head(res)