Pith. sign in

REVIEW 4 major objections 5 minor 5 references

Detailing the stress pattern in the area of central Ionian Islands

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Inverting 485 focal mechanisms shows the Kefalonia Transform Fault Zone is under dextral strike-slip stress, with a shift from extension in the north to transpression toward the south.

desk verdict A solid, incremental stress-inversion study with a plausible but not yet proven north-to-south R trend; worth peer review after adding R uncertainties and a binning-robustness check. read the letter →

arxiv 2506.04777 v1 pith:RFX7NXFD submitted 2025-06-05 physics.geo-ph

V. Anagnostou , E. Papadimitriou , V. Karakostas , T. Bäck This is my paper
classification physics.geo-ph
keywords stressinversionfocalmechanismsKefaloniaTransformFaultZoneIonianIslandsrelativemagnitudestrike-slipfaultingSHmaxseismotectonics
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

Using 485 earthquake focal mechanisms divided into six subareas that follow the five mapped segments of the Kefalonia Transform Fault Zone plus its extensional stepover, this paper inverts for the crustal stress tensor in each subarea and finds a consistent dextral strike-slip regime. In every subarea the maximum compression axis (σ1) is sub-horizontal and coincides with the maximum horizontal stress (SHmax) striking SW–NE, while the intermediate axis σ2 plunges steeply. The relative stress magnitude R increases from 0.41–0.47 in the northern (Lefkada and stepover) subareas to 0.71–0.76 in the southern (Kefalonia) subareas, which the authors read as a transition from a more extensional to a more transpressional character without leaving the strike-slip field. The value of the result is that it details the third-order stress field of the most seismically active part of the Aegean, providing the background needed for stress-triggering and seismic-hazard studies of a region that repeatedly hosts Mw>6 earthquakes.

What carries the argument

The analysis runs on the SATSI damped least-squares stress inversion, implemented through the MSATSI software package, which takes focal mechanism data and returns the orientation of the three principal stress axes and a relative stress magnitude R for each grid point. The R value expresses whether the intermediate principal stress σ2 lies closer in magnitude to the most compressive or the least compressive principal stress, so that small R signals extension and large R signals transpression. A bootstrap resampling procedure, run at least twenty times the number of data per grid point, defines 95% confidence regions for the axis orientations, and the grid itself is imposed by the five mapped fault segments of the KTFZ plus the extensional stepover zone, which together also define the six representative focal mechanisms reported in the paper.

What would settle it

Re-invert the same 485 focal mechanisms with a binning scheme that is not tied to the five fault segments — for example a regular geographic grid or a cluster analysis of mechanism similarity — and check whether the north-to-south rise in R (0.41–0.47 to 0.71–0.76) and the SW–NE SHmax orientation survive; if the gradient or the orientation changes, the reported stress pattern is an artifact of the imposed segmentation. A complementary test would compare the inverted SHmax azimuths with independent borehole breakout or GPS strain-rate directions in the Ionian Islands, since a systematic mismatch would refute the claim that σ1 coincides with SHmax.

Watch

Extended reading notes

Core claim

The paper establishes that the contemporary crustal stress of the central Ionian Islands is dominated by dextral strike-slip faulting of the Kefalonia Transform Fault Zone. At all six grid points the optimal stress tensors show a shallow σ1 trending 242°–249° that coincides with SHmax, a steeply plunging σ2 (62°–84°), and a near-horizontal σ3, with bootstrap uncertainties at 95% confidence. The relative stress magnitude R is 0.41 and 0.45 at the stepover and southern Lefkada grid points, and 0.47 at the northern Lefkada grid point, whereas the three Kefalonia-grid-point values rise to 0.71, 0.75, and 0.76. Because R near 0 indicates that σ2 is close to σ3 and R near 1 indicates that σ2 is close to σ1, the paper interprets this southward increase as a transition from a more extensional regime in the north to a more transpressional regime in the south, while the overall geometry of the principal axes confirms strike-slip as the dominant faulting style throughout.

Load-bearing premise

The inversion assumes that all focal mechanisms inside each of the six grid cells sample one uniform stress tensor, and it fixes the cell boundaries to the five mapped fault segments plus the stepover, so if a cell actually contains more than one stress regime the damped solution returns a smoothed average that may not represent any real subsurface state.

Editorial extensions

If this is right

  • The uniform SW–NE SHmax orientation with σ2 steeply plunging confirms that all five KTFZ segments plus the stepover are loaded for dextral strike-slip failure, so the derived tensors can serve as the background stress field for Coulomb stress-transfer modeling of the region's Mw>6 earthquakes.
  • The north-to-south rise in R implies that the southern Kefalonia branch operates under more transpressional conditions than the northern Lefkada branch, which may translate into differing recurrence behavior and rupture styles along the zone.
  • The six representative focal mechanisms derived from the best tensors offer an internally consistent source model set for strong-motion and finite-fault modeling of future scenarios in the central Ionian Islands.
  • Because the grid was chosen to match known structures rather than a regular geography, the paper demonstrates that fault-segment-guided binning can resolve third-order stress variations that regional stress maps miss.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A testable next step is to compare the R-value gradient against GPS strain rates: rapid along-strike changes in dilatation or shear strain rate would corroborate the north-to-south extension-to-transpression transition, and their absence would suggest the gradient is an artifact of the six-cell binning.
  • The fixed segmentation assumption could be probed by running the same SATSI inversion with the grid shifted by half a segment length; if the main pattern persists, the stress field is robust, and if not, the mapped segment boundaries are exerting more control than the data.
  • The paper's extension-to-transpression interpretation relies solely on R, but a steeply plunging σ2 with σ1 subhorizontal could be equally well described as a strike-slip regime with a minor normal component in the north; a future study could test which description better predicts aftershock geometries.
  • Applying the same fault-segment-binned inversion workflow to the broader Hellenic Arc, where the plate boundary transitions from subduction to collision, could reveal whether the R gradient seen here is a regional signature of that transition.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper inverts a compiled set of 485 focal mechanisms from the central Ionian Islands using the MSATSI implementation of the damped least-squares stress inversion method of Hardebeck and Michael (2006). The dataset is binned into six subareas that follow the five known segments of the Kefalonia Transform Fault Zone (KTFZ) plus the extensional stepover zone between the Lefkada and Kefalonia branches. For each subarea the authors report the optimal principal stress axes, bootstrap uncertainties for axis orientations, relative stress magnitudes R, and a representative focal mechanism derived from the best tensor. The main conclusions are that the area is uniformly dominated by dextral strike-slip faulting, with σ1 roughly SW-NE and subhorizontal, σ2 steeply plunging, and σ3 subhorizontal; and that R values increase from the northern subareas (0.41–0.47) to the southern Kefalonia subareas (0.71–0.76), interpreted as a transition from more extensional to more transpressional conditions.

Significance. If the central spatial claim holds, the paper would provide a usefully detailed, third-order stress picture of a seismically active transform fault zone, with a clear north-to-south transition in stress regime that could inform Coulomb stress-transfer and hazard studies. Strengths include the use of a standard, well-established inversion method; bootstrap uncertainty estimates for the principal-axis orientations; a careful compilation from published and routine catalogues with duplicate handling; and an explicit statement of the 50% fault-plane probability choice. The paper is, however, a short conference contribution, and the main novelty—the R-value trend—requires additional rigor before it can be considered established. In particular, no uncertainties are given for R, and the binning into fault-segment-based cells is assumed without testing whether a single uniform tensor adequately represents each cell.

major comments (4)
  1. [Exploitable Results; Table 1] The north-to-south R transition (0.41–0.47 vs 0.71–0.76) is the central quantitative result of the paper, yet no uncertainties are reported for any of the R values. The bootstrap resampling used for the principal-axis orientations should be extended to R, and the 95% confidence intervals for R should be reported. Without these, the reader cannot judge whether the apparent contrast between the northern and southern groups is significant or within bootstrap scatter.
  2. [Data and Method; Figure 1] The gridding scheme follows the five KTFZ segments plus the stepover zone, and the SATSI inversion assumes a single uniform stress tensor within each grid cell. The paper provides no test of this homogeneity assumption, for example by comparing results with alternative binning, by performing within-cell jackknife tests, or by reporting the data misfit and model length per cell. Because the R contrast is the main claim, the possibility that the trend is an artifact of where the segment-based boundaries were drawn must be addressed.
  3. [Exploitable Results; Summary and Conclusions] The representative focal mechanisms are constructed directly from the best stress tensor derived from the inversion (Data and Method), and the text then states that all of them point to a strike-slip regime. This is circular: the representative mechanisms cannot serve as independent evidence for the strike-slip interpretation. The strike-slip conclusion should rest on the principal-axis orientations and on the original focal-mechanism data, and the representative mechanisms should be explicitly labeled as derived products.
  4. [Data and Method; Exploitable Results] The text states that the optimal damping parameter is chosen from the trade-off curve between data misfit and model length, but it does not report the chosen damping value, the resulting data misfit, or the model length for any of the six inversions. Reporting these quantities would allow readers to assess the degree of spatial smoothing and the fit quality, which is directly relevant to the reliability of the fine-scale R differences.
minor comments (5)
  1. [Data and Method, Eq. (1)] Equation (1) is referred to but the defining expression for R is not actually displayed in the text; please include the formula.
  2. [Figure 1; Table 1] Figure 1 shows the subareas as rectangles, but the grid-point coordinates (0,0), (0,1), (0,2), (1,2), (1,3), (2,3) used in Table 1 are not labeled on the map, making it difficult to link the table rows to the geographic subareas.
  3. [Data and Method] The text says each subarea comprises at least 20 focal mechanisms, but the actual number of mechanisms per grid point is not reported; including this information in Table 1 or the text would be useful for assessing robustness.
  4. [Exploitable Results] The phrase 'maximum compression axis is equal to the SHmax' should be qualified: for a strike-slip regime with σ2 vertical, SHmax coincides with the azimuth of σ1, but the equality is a regime-dependent relationship rather than an identity; a brief clarification would avoid ambiguity.
  5. [Abstract and Background] The degree symbols in 'striking from 12ο to 40ο' are rendered inconsistently; use a consistent degree symbol.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the stress inversion is a data-driven inference; the representative focal mechanisms are explicitly constructed from the inversion output, and the self-citations to prior KTFZ segmentation are not load-bearing for the central claim.

full rationale

The paper's derivation chain is: compile 485 focal mechanisms -> bin by KTFZ segments (a methodological choice, not a fitted result) -> invert each subset with SATSI/MSATSI -> obtain principal axes and R -> interpret as strike-slip with north-south R transition. Each step is an inference from the focal mechanism data. The representative focal mechanisms are constructed from the best-fitting tensor ('we construct a representative focal mechanism for each grid point, according to the best stress tensor derived from our methodology'), so their strike-slip character is a restatement of the inversion, not an independent confirmation; however, the paper does not use them as load-bearing evidence for the stress regime. The stepover zone and five-segment geometry are cited to prior work by overlapping authors (Karakostas et al. 2015; Kourouklas et al. 2023), but those citations only define the binning scheme; the R values and axis orientations are estimated from the focal mechanisms themselves. No uniqueness theorem is invoked, no parameter is fitted to a subset and then 'predicted' for a closely related quantity, and no known result is renamed. The north-to-south R trend is a direct output of the inversion; whether a single uniform tensor per cell is adequate is a modeling assumption, not a circularity. Thus the paper is self-contained against external benchmark data (published focal mechanisms), and no circular step meets the quoting standard.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The central claim rests on the standard Wallace-Bott assumption, uniform stress per grid cell, and the adopted fault segmentation. Two hand-set parameters (50% fault-plane probability and damping parameter) and a minimum-count threshold shape the inversion; R and axes are outputs, not external benchmarks. No new physical entities are introduced.

free parameters (3)
  • optimal damping parameter = not reported
    MSATSI selects it from the trade-off curve between data misfit and model length; it controls smoothing between grid points and affects the recovered stress tensors.
  • fault-plane probability = 50%
    Authors set the probability that each input focal mechanism plane is the true fault plane to the typical 50%, which affects inversion weights.
  • minimum focal mechanisms per grid point = 20
    Threshold chosen to ensure robustness; affects which grid points are included and how the dataset is binned.
assumptions (4)
  • domain assumption Stress is uniform within each gridded subarea.
    SATSI solves for one smoothed tensor per grid point; mixed regimes inside a cell would be averaged. Invoked in the gridding step of Data and Method.
  • domain assumption Focal mechanisms are samples of the regional deviatoric stress via the Wallace-Bott hypothesis.
    Standard assumption of all focal-mechanism stress inversions; local fault heterogeneity and pore pressure changes are not modeled. Invoked throughout Data and Method.
  • domain assumption The five-segment KTFZ geometry and stepover location used for binning are correct.
    Grid points follow Kourouklas et al. (2023) and Karakostas et al. (2015); incorrect geometry would misassign mechanisms to stress cells. Invoked in Data and Method and Figure 1.
  • standard math A deviatoric stress tensor is fully described by three orthogonal principal axes plus the R ratio.
    Used in Eq. (1) and MSATSI output; standard linear algebra.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Detailing the stress pattern in the area of central Ionian Islands." pith.science (2026). https://pith.science/paper/RFX7NXFD

@misc{pith2026250604777,
  author       = {Pith},
  title        = {Pith review of: Detailing the stress pattern in the area of central Ionian Islands},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RFX7NXFD}},
  note         = {Machine review of arXiv:2506.04777}
}
read the original abstract

The Kefalonia Transform Fault Zone (KTFZ) is the most seismically active area in the Mediterranean and consists of two major branches, the Lefkada fault segment to the north and the Kefalonia fault segment to the south. KTFZ acts as an active boundary between the subduction zone of the remnants of the oceanic lithosphere of the Eastern Mediterranean that subducts under the Aegean microplate to the south and the continental collision between the Eurasian plate and the Adriatic Microplate to the north. The tectonic activity in the region is reflected in the rapid crustal deformation rates of the region and subsequently the frequent occurrence of strong earthquakes (Mw>6.0) that occurred during both the historical and instrumental era of seismology. Those strong earthquakes and their temporal distribution can be explained due to stress transfer between closely located fault segments (Papadimitriou, 2002) and as such, studying those stress interactions is an integral part of understanding the long-term tectonic loading in the region.

Figures

Figures reproduced from arXiv: 2506.04777 by the authors.

Figure 1
Figure 1. Map of the study area of the central Ionian Islands along with the major fault segments of the KTFZ and the extensional stepover area, denoted with red lines and enclosed by rectangles that correspond to the subareas or grid points containing the focal mechanisms used in the inversion. The focal mechanisms are shown as equal area lower hemisphere projections [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

  1. [1]

    Three-body structure of $^{19}$B: Finite-range effects in two-neutron halo nuclei

    Hardebeck, J. L., & Michael, A. J . (2006). Damped regional‐scale stress inversions: Methodology and examples for southern California and the Coalinga aftershock sequence. Journal of Geophysical Research: Solid Earth , 111(B11). https://doi.org/10.1029/2005JB004144 Hergert, T., & Heidbach, O. (2011). Geomechanical model of the Marmara Sea region —II. 3 -D...

    work page Pith review arXiv 2006

  2. [295]

    https://doi.org/10.3390/geosciences13100295 Martínez‐Garzón, P., Kwiatek, G., Ickrath, M., & Bohnhoff, M. (2014). MSATSI: A MATLAB package for stress inversion combining solid classic methodology, a new simplified user‐handling, and a visualization tool. Seismological Research Letters, 85(4), 896-904. https://doi.org/10.1785/0220130189 Papadimitriou, E. E...

    work page doi:10.3390/geosciences13100295 2014

  3. [1983]

    https://doi.org/10.1007/BF00880738 Zoback, M

    Pure and Applied Geophysics, 123(3), 388–397. https://doi.org/10.1007/BF00880738 Zoback, M. L. (1992), First- and second-order patterns of stress in the lithosphere: The World Stress Map Project , J. Geophys. Res., 97(B8), 11703–11728, doi:10.1029/92JB00132 Data Sources Benetatos, C., Kiratzi, A., Papazachos, C., & Karakaisis, G. (2004). Focal mechanisms ...

    work page doi:10.1007/bf00880738 1992

  4. [2002]

    https://doi.org/10.1016/J.PEPI.2004.03.008 Pondrelli, S., Salimbeni, S., Morelli, A., Ekström, G., & Boschi, E

    Physics of the Earth and Planetary Interiors , 145(1–4), 127 –147. https://doi.org/10.1016/J.PEPI.2004.03.008 Pondrelli, S., Salimbeni, S., Morelli, A., Ekström, G., & Boschi, E. (2007). European –Mediterranean Regional Centroid Moment Tensor catalog: Solutions for years 2003 and

    work page doi:10.1016/j.pepi.2004.03.008 2007

  5. [2004]

    https://doi.org/10.1016/J.PEPI.2007.05.004 Pondrelli, S., Salimbeni, S., Morelli, A., Ekström, G., Postpischl, L., Vannucci, G., & Bo schi, E

    Physics of the Earth and Planetary Interiors , 164(1–2), 90–112. https://doi.org/10.1016/J.PEPI.2007.05.004 Pondrelli, S., Salimbeni, S., Morelli, A., Ekström, G., Postpischl, L., Vannucci, G., & Bo schi, E. (2011). European – Mediterranean Regional Centroid Moment Tensor catalog: Solutions for 2005 –2008. Physics of the Earth and Planetary Interiors, 185...

    work page doi:10.1016/j.pepi.2007.05.004 2011

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.