tectonicr is a free and open-source R package for modeling and analyzing the direction of the maximum horizontal stress (SHmax) based on the empirical link between the direction of intraplate stress and the direction of the relative motion of neighboring plates (Wdowinski, 1998; Stephan et al., 2023). The following methods are available:
PoR_shmax()
gives the predicted stress field adjacent to a
plate boundary, calculated using the relative plate motion of the
neighboring plates using the function model_shmax()
. The
goodness-of-fit can be statistically tested by
e.g. norm_chisq()
, circular_dispersion()
,rayleigh_test()
, and
confidence_interval()
.distance_from_pb()
gives the distance between the stress
data point and the plate boundary measured along the stress
trajectories.eulerpole_paths()
generates an sf
object
containing spatial information that is suitable to plot with, for
instance, ggplot()
.equivalent_rotation()
transfers a set
of plate motion parameters into the relative plate motions among the
given plates.PoR_stress2grid()
uses distance, method, and
quality-weighted mean direction of stress data without being affected by
angular distortions.rose()
shows the frequencies
of the orientations in polar coordinatesYou must have R installed on your system (see http://r-project.org). To install tectonicr from CRAN, type the following code at the R command line prompt:
install.packages("tectonicr")
The most recent development version of tectonicr is available from Github and can be installed on your system as follows:
# install.packages("remotes") # install if needed
remotes::install_github('tobiste/tectonicr')
library('tectonicr')
Using the stress measurements from the San Andreas Fault - Gulf of California example, a quick analysis and test against the right-lateral transform plate boundary can be achieved by:
# load example data
data(san_andreas)
data(nuvel1)
data(plates)
stress_analysis(
san_andreas,
PoR = equivalent_rotation(nuvel1, "na", "pa"),
type = "right",
pb = subset(plates, pair == "na-pa")
)
The code produces the underlying results for Fig. 7 from Stephan et al. (2023)
The detailed documentation can be found at https://tobiste.github.io/tectonicr/articles/tectonicr.html
Tobias Stephan
When referencing this package, please cite
Stephan, T., Enkelmann, E., and Kroner, U. (2023). Analyzing the horizontal orientation of the crustal stress adjacent to plate boundaries. Scientific Reports, 13(1). DOI: 10.1038/s41598-023-42433-2.
and the package DOI: 10.5281/zenodo.8372508.
GPL-3.0 License