LoginSign Up
← Back to all sparks
T

tectonicr

INFRA · APIS
Velocity0.0

Analyzing the Orientation of Maximum Horizontal Stress

github.com/tobiste/tectonicr ↗

Stress-field analysis tracks the 2025 World Stress Map and keeps adding plotting shortcuts.

geophysicsstress-fieldcircular-statisticsplate-tectonicsggplot2
◆Current state
tectonicr analyses the orientation of maximum horizontal stress, comparing observed azimuths against those predicted by plate motion in a plate-of-rotation reference frame. Version 0.4.7 added support for the 2025 World Stress Map release and made the version selectable at download; 0.4.8 followed with new summary statistics for axial data, ggplot2 geoms for plotting directions, and a faster kappa estimator.
◆Where it's heading
The package is maturing along two axes at once. On the statistics side the circular methods keep getting corrected and accelerated, with a major fix to est.kappa() for doubling angles before treating them as directional data and a faster approximation replacing the maximum likelihood path by default. On the usability side the additions are convenience wrappers, data2PoR() and the geom_azimuth pair, which is what a package does once its core methods stop moving. A sign convention for deviation was also flipped, which quietly changes existing results.
◆Prediction
Expect the convenience layer to keep growing around a stable core, and the next substantive release to track whatever the World Stress Map publishes after 2025.

◆Recent moves

  1. 8mo ago

    ggplot2 azimuth geoms and axial summary statistics added

    The largest feature release in the window, and mostly convenience: geom_azimuth() and geom_azimuthpoint() for plotting directions as bars, data2PoR() for adding reference-frame coordinates in one call, ortensor2d() for axial summary statistics, and a weighting() helper. Underneath, est.kappa() switches to a faster approximation and the deviation sign convention changes so counterclockwise deviations are positive.

    View source ↗
  2. 1y ago

    World Stress Map 2025 supported with selectable version

    The package's reference dataset moves forward: download_WSM() defaults to the 2025 World Stress Map with a version argument to pin older ones, and parse_quality() encodes the 2025 treatment of X-ranked data as 90-degree uncertainty. A major fix to est.kappa(), which was not doubling angles before treating the data as directional, matters more than the data update for anyone who reported a concentration parameter.

    View source ↗
  3. 1y ago

    Spatial interpolation performance improved

    A one-line release claiming a performance boost for stress2grid() and related functions. Gridding is the expensive step in this package's workflow, so the target is the right one even without numbers attached.

    View source ↗
  4. 2y ago

    Weighting powers added to spatial interpolation

    Interpolation gains weighting powers in stress2grid(), giving users control over how strongly nearby observations dominate the grid. Small, but it is the parameter that determines how smooth a published stress map looks.

    View source ↗
  5. 2y ago

    Bootstrap dispersion statistics ahead of CRAN submission

    The pre-CRAN release, adding bootstrap statistics for circular dispersion, multiple-angle input across the dispersion functions, and spatial distribution of dispersion. It sets up the circular statistics surface that the later releases spend their time correcting and speeding up.

    View source ↗