the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Influence of crustal mechanical layering on the seismic potential of active faults: insights from the southwestern Valencia Trough (W Mediterranean)
Abstract. We present a structural and seismotectonic analysis of active faults in the southwestern Valencia Trough (western Mediterranean) on the basis of subsurface datasets. In our study, we identify and characterise three major active faults: the Cullera Fault, with a long-term slip rates that vary over time between 0.15 ± 0.1 mm/yr and 0.4 ± 0.1 mm/yr; the oblique Albufera Fault, with a long-term slip rate of 0.2 ± 0.1 mm/yr; and the normal Valencia Fault.
The seismogenic character of the southwestern Valencia Trough is controlled by a mechanically weak layer consisting of Triassic evaporites. This weak layer induces partial to complete decoupling between the suprasalt and subsalt successions, leading to two distinct mechanisms driving fault displacement: tectonic activity and salt withdrawal. A quantitative evolutionary analysis of the Cullera Fault reveals that these two mechanisms alternate over time.
The presence of a mechanically weak layer has implications for seismicity. Earthquakes can nucleate within both sub- and suprasalt successions, with total or partial decoupling influencing rupture propagation. We discuss how these two scenarios lead to different earthquakes and thus impact the seismic hazard of a region. Empirical source-scaling relationships, which are commonly used to estimate the seismogenic potential of active faults, generally assume a homogeneous seismogenic crust. To address this limitation, we propose a methodological approach based on the use of the aspect ratio. Our findings highlight the need to incorporate stratigraphic mechanical layering into seismic hazard assessments, particularly in salt-influenced tectonic settings.
- Preprint
(3020 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2025-1782', Raul Perez-Lopez, 06 Aug 2025
-
AC1: 'Reply on RC1', Ivan Martin-Rojas, 24 Oct 2025
Below is our detailed response to Reviewer 1's comments.
S1-I suggest that the authors consider including the phrase “salt-influenced tectonic setting” in the title of the manuscript. This represents the main novel contribution of the work, along with the parameterization of the three active faults in the Valencia Trough.
We greatly appreciate this suggestion, as it will help clarify the scope of our work for readers. Therefore, if our manuscript is ultimately accepted for publication, we will revise the title as follows: Influence of salt-related mechanical layering on the seismic potential of active faults: insights from the southwestern Valencia Trough (Western Mediterranean)
S2-Section 4. Stratigraphy. I recommend reducing this section, as it does not represent the main objective of the work. Providing such a detailed stratigraphic description may distract the reader from the key focus of the study and could make the narrative less clear. A more concise summary would keep the manuscript streamlined and focused.
That is true. We thank Reviewer 1 for bringing this to our attention. If the manuscript is ultimately accepted for publication, we will follow the reviewer’s suggestion and reduce this section accordingly.
S3-Some figures and graphics are difficult to read, and I have included some recommendations into the attached pdf file. For instance, letters with a grey color, small numbers in the isopach lines, and similar colors I plots for the SR long-term estimation.
Once again, we appreciate Reviewer 1’s recommendations regarding the overall arrangement of the figures. We will follow her/his suggestions to improve the presentation of our figures.
S4- I have some major concerns regarding the comparison between the Valencia Trough and the Zagros Thrust and Belt based solely on the presence of a “mechanically weak layer.” This comparison overlooks several tectonically significant parameters, starting with the fundamentally different tectonic settings: The Zagros is characterized by crustal shortening and active mountain building, whereas the Valencia Trough corresponds to a uniaxial extensional marine basin. In addition, isostatic processes are not considered, and the rupture mechanisms are inherently different. Reverse faulting in the Zagros involves the mobilization of very large rock volumes in the hanging wall, which is not directly comparable to normal faulting in a sedimentary marine basin. A more cautious discussion is required to avoid oversimplifying the analogies between these two geodynamically contrasting regions. Thermal coupling in different depth are also ignored.
S5- You have not demonstrated this claim (similarities between both tectonic settings, Valencia Trough and Zagros region). What is shown is only a convergence in the geometry of the spatial distribution of seismicity, but this does not constitute a “pattern” based on the underlying mechanical behavior of the faults. A true mechanical pattern would require evidence of consistent fault kinematics, stress distribution, and deformation mechanisms, which are not demonstrated in the current version of the manuscript. Please, consider to remove this part, I have no the relevance in your paper. This only distracts from the methodology that authors propose.
We understand Reviewer 1’s concern regarding the analogy presented in our manuscript between the Valencia Trough and the Zagros Mountains. We are fully aware of the significant geodynamic differences between these two regions, as we have already acknowledged in the manuscript. Additionally, Reviewer 2 also expressed concerns about this comparison. The convergence of both reviewers’ opinions suggests that our interpretation may have gone too far. Consequently, we agree with Reviewer 1 and believe it is best to remove this part from our manuscript.
A more relevant regional comparison might be found in central Italy, which has been undergoing extension since the Pliocene. This tectonic regime has produced several active faults that were the seismogenic sources of multiple recent destructive earthquakes, including those in the 2016–2017 seismic sequence. This sequence comprises the 2016 Amatrice Earthquake (Mw 6.0), the 2016 Mt. Visso Earthquake (Mw 5.9), and the 2016 Norcia Earthquake (Mw 6.5) as the more significant events.
The geology of this region includes a significant presence of Triassic salt, which is organized into several layers due to pre-extensional thrusting. Recent studies have shown that some of the main events during the 2016–2017 sequence nucleated within the salt layers at very shallow depths (3.5–4 km for events up to Mw 5.9; Barchi et al., 2021; Barchi & Collettini, 2019). Moreover, some of the seismogenic faults in this region offset the basement-cover interface.
We believe that a comparison between the southwestern Valencia Trough and central Italy could provide valuable insights into the seismotectonic configuration of our study area, as both regions are actively extending and feature a seismogenic crust with mechanically weak layers. Therefore, if our manuscript is accepted for publication, we propose to include a paragraph discussing the analogies between these two regions.
S6- I also recommend removing the estimation of the maximum earthquake magnitude based solely on the total fault length, as the results appear unrealistic (Table 1). This approach tends to overestimate the potential magnitude and is not fully supported by seismotectonic principles. Instead, relying on fault surface area provides a more consistent and physically based estimate, as it is directly related to the seismic energy release (Kanamori and Anderson).
We agree with Reviewer 1 that calculations based solely on fault length tend to overestimate the maximum expected magnitude, especially in a region with the geodynamic characteristics of the Valencia Trough. We also concur that estimations based on fault surface area provide a more reliable constraint on the seismic potential of active faults.
Therefore, as suggested by Reviewer 1, we believe it would be appropriate to include in Table 1 the results of maximum magnitude calculations based on fault surface area. We originally included the length-based calculations to highlight the tendency for overestimation. However, the scaling relationships based on fault length, such as those proposed by Wells and Coppersmith (1994) and Stirling et al. (2002), are widely used in the literature. As such, we consider them a useful reference to frame the discussion on seismic potential.
Rather than removing the length-based estimations, we propose to complement Table 1 with the area-based results.
S7-Conclusions have to be revisited. Some relevant information derived from your results are missed.
We are grateful to Reviewer 1 for pointing out this issue. As noted, some important information is missing from the conclusions section—particularly the seismic characterization of the active faults. If our manuscript is accepted for publication, we will include the missing information to strengthen the conclusions.
Technical suggestions and questions
TSQ1-Page 4; Line 85 You refer “mechanical weak layer” to the Triassic evaporitic formation. My question: Are you referring to ductile behavior under a tectonic stress field? If so, how do you estimate the creep movement associated with ductile salt deformation?
“Mechanical weak layer” is a commonly used term in many studies addressing the mechanical behavior of the crust. It refers to a layer with lower rupture or ultimate strength than the surrounding crust, regardless of whether its behavior is brittle or ductile.
In our context, as mentioned in the manuscript, the Triassic evaporitic succession is considered a “weak layer”, as it acts as the regional detachment. However, we believe that a detailed characterization in terms of ductile or brittle behavior, as well as an estimation of the creep associated with the potential ductile deformation, would require extensive quantitative analyses that go beyond the scope of our study.
Therefore, we consider that a qualitative description of the Triassic layer as a weak layer, supported by robust literature data, is sufficient for the purposes of our work.
TSQ2- Reading the whole manuscript, I still do not fully understand what you mean by the term “mechanically weak layer” regarding the evaporite Traissic layer. Could you please clarify its implication for the seismogenic behavior of the fault system?
Are you suggesting that such a weak layer:
- makes the system more prone to triggering earthquakes?
- would favor medium-sized events without significant tectonic strain accumulation?
- or behaves differently in terms of long-term stress buildup?
A clearer explanation would help in understanding the role of this layer in the seismic cycle.
When we refer to the Triassic succession as a “mechanical weak layer,” our intention is to highlight its influence on the vertical propagation of seismic rupture, as discussed in several parts of the manuscript. However, Reviewer 1’s comments suggest that this point is not clearly presented.
To address this, we will clarify and emphasize our use of the term “mechanical weak layer” from the beginning of the manuscript. Specifically, we will add explanatory sentences in the Introduction, Tectonic Setting, and Discussion sections to explicitly define what we mean by this term and how it relates to the seismic behavior of the study area.
TSQ3- The Cullera Fault appears to show a decrease in slip rate from the Pliocene to the Quaternary (Fig. 8). Could the authors clarify the underlying reason for this apparent behavior? Do you have a conceptual model that could explain this trend? Is it a Nubia plate deceleration?
As Reviewer 1 indicates, the offset of the markers observed in the seismic lines point to a decrease of the slip rate of the Cullera Fault from 0.40 ± 0.1 mm/y for the Pliocene to 0.15± 0.1 mm/y for the Quaternary. The reasons for this decrease could be different. As Reviewer suggests, the complex geodynamic interaction between Nubia and Eurasia plates could be one of these reasons. Also, as we discuss in section 6, the evolution of this salt-influenced fault is controlled by both tectonics and salt withdrawal. So, fault slip rate is the result of a complex interaction between these two mechanisms. Therefore, a decrease of the tectonic component of deformation, as Reviewer 1 suggests, could be one of the reasons for the decrease in the slip rate. But a decrease in the salt-related offset could also explain this variation. The available data permit us to discuss which of these two mechanisms was predominant during the Plio-Quaternary time span. But we believe that to totally explain the decrease in slip rate observed during the Quaternary, it is necessary to have a model that would encompass both kinematic and dynamic data in the region. Our results help to clarify the kinematic situation, but further data is necessary to solve the dynamic part of the problem. However, we believe that, as Reviewer 1 indicates, this is a question that is worth mentioning. Therefore, following that if our manuscript is finally accepted for publication, we will add a brief paragraph discussing the hypothetical reasons that could explain the slip rate decrease.
We thank Reviewer 1 for highlighting this issue. As noted, the offset of markers observed in the seismic lines suggests a decrease in the slip rate of the Cullera Fault—from 0.40 ± 0.1 mm/year during the Pliocene to 0.15 ± 0.1 mm/year during the Quaternary.
Several factors could explain this decrease. As Reviewer 1 suggests, the complex geodynamic interaction between the Nubia and Eurasia plates may be one contributing factor. Additionally, as discussed in Section 6, the evolution of this salt-influenced fault is governed by both tectonics and salt withdrawal. Thus, the fault slip rate reflects a complex interplay between these two mechanisms.
A reduction in the tectonic component of deformation, as proposed by Reviewer 1, could be one explanation. However, a decrease in salt-related displacement could also account for the observed variation. While the available data allow us to discuss which of these mechanisms may have been dominant during the Plio-Quaternary period, we believe that a comprehensive explanation of the slip rate decrease during the Quaternary would require a model integrating both kinematic and dynamic data for the region.
Our results contribute to clarifying the kinematic framework, but additional data are needed to address the dynamic aspects. Nevertheless, we agree with Reviewer 1 that this is an important issue worth mentioning. Therefore, if our manuscript is accepted for publication, we will include a brief paragraph hypothesizing the possible reasons behind the observed decrease in slip rate.
TSQ4- The curves of the Figure 10, Throw depth (T-z) and expansion index plots for the Cullera Fault computed from profiles SGV01-115 and SGV01-117, show a fault throw variation. Could the geometry of these curves be explained by the listric geometry of the Albufera Fault? Variations in fault throw along planes with different dips could potentially generate this type of curve. Could the authors clarify if this interpretation is consistent with their data?
As Reviewer 1 points out, throw along normal listric faults varies with depth due to the curved geometry of the fault plane. Consequently, a passive horizon located at a deeper structural level experiences greater displacement than one situated in the shallower part of the structure.
However, in the case of the Cullera Fault, the horizons analyzed are not passive; they were deposited contemporaneously with fault displacement. This has two main implications:
i) The position of the horizons relative to the fault changes over time, and
ii) Younger markers have been affected by fault activity for a shorter duration than older markers.Therefore, we propose that the upward decrease in fault throw observed in Figure 10 is a result of the differing ages of the horizons—older horizons show greater offset because they have been displaced over a longer time interval.
In any case, Reviewer 1’s comment clearly indicates that the reasons behind this throw decrease are not sufficiently addressed in our current manuscript. To improve clarity, we intend to add a brief discussion on this point.
TSQ5- The accommodation space created by the Albufera Fault, and cited in the line 559, page 30, can be explained by salt- tectonics? Please explain your model.
The mechanisms responsible for the creation of accommodation space along the Cullera Fault are extensively discussed in our manuscript (lines 559–600). However, Reviewer 1’s comment suggests that the connection between the description of the accommodation space (lines 548–558) and the subsequent discussion of our hypothesis for its formation (lines 559–600) may not be sufficiently clear.
To address this, we will add a sentence at the beginning of the discussion paragraph (line 559) to explicitly highlight the relationship between the observed accommodation space and the proposed mechanisms for its development.
TSQ6- Line 641 Tectonic driving mechanism and differences between tectonic driving and sal-withdrawal driving. I have some doubts regarding this point. The “mechanically weak layer” is described at a depth of 2–5 km, but tectonic earthquakes with significant magnitudes (M6–7) typically require rupture depths greater than this. Moreover, the regional tectonic strain is expected to be distributed below the weak layer as well. This raises the question: does the fault release only the strain accumulated in the supra-salt section, or can it also rupture deeper segments (>10 km) where larger tectonic strain is stored? A discussion of this point would clarify the seismogenic potential of the fault system.
In this comment, Reviewer 1 refers to two distinct aspects.
First, Reviewer 1 notes that significant earthquakes (Mw > 6) typically nucleate at depths greater than 5 km. However, recent events in central Italy and France have demonstrated that very shallow earthquakes of moderate magnitude (Mw 5–6) can also occur (Chiaraluce et al., 2017; Godano et al., 2025; Improta et al., 2019). In the southwestern Valencia Trough, the thickness of the supra-salt succession is slightly less than 5 km. Therefore, while the probability of a significant earthquake nucleating within the supra-salt succession is low, it is not zero.
Assuming that such an event could occur, our aim is to evaluate the seismic potential of these (unlikely) supra-salt ruptures. In other words, we focus on describing a characteristic event rather than assessing its likelihood. Reviewer 1’s comment suggests that this distinction was not clearly stated in our manuscript. To address this, we will add a sentence clarifying that events nucleating within the supra-salt succession are considered unlikely.
Secondly, Reviewer 1 asks whether “the fault releases only the strain accumulated in the supra-salt section, or can it also rupture deeper segments.” In our manuscript, we explicitly state that a tectonically driven earthquake could rupture the entire seismogenic crust. This implies that strain accumulates in both the supra-salt and sub-salt successions, as tectonic stress affects the entire crust.
Once again, Reviewer 1’s comment indicates that this point may not have been clearly conveyed. Therefore, if our manuscript is accepted for publication, we will add the necessary clarifying sentences to ensure this aspect is properly explained.
Chiaraluce, L., Di Stefano, R., Tinti, E., Scognamiglio, L., Michele, M., Casarotti, E., Cattaneo, M., De Gori, P., Chiarabba, C., Monachesi, G., Lombardi, A., Valoroso, L., Latorre, D., & Marzorati, S. (2017). The 2016 Central Italy Seismic Sequence: A First Look at the Mainshocks, Aftershocks, and Source Models. Seismological Research Letters, 88(3), 757–771. https://doi.org/10.1785/0220160221
Godano, M., Larroque, C., Delouis, B., Ampuero, J., Arzu, F., Courboulex, F., Deschamps, A., van den Ende, M., Baize, S., & Ritz, J. (2025). Back to the Source: Connecting the Seismological Observations of Le Teil Earthquake (M w 4.9, 2019/11/11, France) to the Local Geology. Journal of Geophysical Research: Solid Earth, 130(9). https://doi.org/10.1029/2025JB031133
Improta, L., Latorre, D., Margheriti, L., Nardi, A., Marchetti, A., Lombardi, A. M., Castello, B., Villani, F., Ciaccio, M. G., Mele, F. M., Moretti, M., Battelli, P., Berardi, M., Castellano, C., Melorio, C., Modica, G., Pirro, M., Rossi, A., Thermes, C., … Di Maro, R. (2019). Multi-segment rupture of the 2016 Amatrice-Visso-Norcia seismic sequence (central Italy) constrained by the first high-quality catalog of Early Aftershocks. Scientific Reports, 9(1), 6921. https://doi.org/10.1038/s41598-019-43393-2
TSQ7- Methodology and the “correction factor” (line 712, page 36). I agree that using the total fault length to estimate earthquake magnitude through common empirical relationships is not the most reliable approach. However, the total fault surface area provides a more reasonable approximation because it is directly related to the seismic energy released during an earthquake (see Kanamori and Anderson). The main concern with introducing a “correction factor” is the potential to generate an artifact in the magnitude estimation, which could lead to biased results if not properly justified or calibrated. On the other hand, why the “aspect radio” factor of the active fault could improve the magnitude estimation? It is not clear in the text.
We agree with Reviewer 1 that using an area-based scaling relationship is a more reliable approach for estimating the seismic potential of a fault. While all scaling relationships are based on extensive empirical observational data, they do not explicitly account for the specific geological and mechanical characteristics of the seismogenic crust, as we mention in our manuscript (see also response to Reviewer 2).
Our goal is to propose a method to refine these scaling relationships by incorporating the mechanical properties of the crust. This is why we suggest using the aspect ratio as a correction factor for the fault area used in magnitude calculations.
Once again, Reviewer 1’s comment indicates that our objective may not have been clearly stated in the manuscript. Therefore, we intend to add clarifying sentences to this paragraph to better explain our approach.
TSQ8- PAY ATTENTION! revise the fault areas of the Table 2. There are some incongruences with the values obtained from the 2D Move analyses and the values indicated in the Table 2 for the Cullera, Albufera and Valencia faults (pages 20, 24 and 27).
As Reviewer 1 points out, there appears to be an inconsistency between the input area used for calculating the maximum expected magnitude and the fault areas reported in Section 5. This apparent discrepancy arises from the use of two different calculation methods.
The values presented in Section 5 are direct measurements obtained from the fault surfaces in our 3D models. In contrast, in Table 2 (Section 8), we chose to compute fault area by multiplying fault length by fault width. This simplified approach was preferred because, in a subsequent step, we apply a weighting using the aspect ratio as a correction factor.
For this weighted calculation, we use the thickness of the supra-salt succession as an input parameter. Therefore, the only feasible way to compute the fault area in this context is by multiplying this thickness by the weighted fault length. To ensure consistency when comparing magnitude estimates derived from both weighted and unweighted areas, we used the same calculation method for both. Also, this calculation would be easily implemented in regions where detailed subsurface data are available.
This methodological detail was not explained in the manuscript. Thanks to Reviewer 1’s careful review, we will add the necessary clarifying sentences to make this point explicit.
TSQ9- Table 2. I don’t understand how a great fault area triggers an earthquake of lower magnitude as appeared in Albufera and Valencia faults). Please, revise carefully TABLE 2.
We thank Reviewer 1 for the careful revision. As noted, there was an error in Table 2, where the Mw values for the Albufera and Valencia faults were inadvertently swapped. This mistake will be corrected in the revised version of the manuscript.
Additionally, I have attached the PDF manuscript with some minor detailed comments posted into the PDF file.
We will address properly the minor comments included by Reviewer 1 in the pdf file.
Citation: https://doi.org/10.5194/egusphere-2025-1782-AC1
-
AC1: 'Reply on RC1', Ivan Martin-Rojas, 24 Oct 2025
-
RC2: 'Comment on egusphere-2025-1782', Anonymous Referee #2, 23 Sep 2025
The authors present a detailed model for stratigraphy, lithology and structure of the southwestern Valencia Trough, based on vintage and modern seismic reflection lines, well data, electric data, gravity data, and taking into account previous studies. The presented reflection lines appear well readable, and the interpretation appears clear or at least reasonable. A crucial feature of the model is a layer of evaporitic Upper Triassic (Keuper) rocks, that behaves mechanically weak and may induce decoupling of the Palaeozoic–Middle Triassic basement from the Upper Triassic –Quaternary units.
The authors identify three major faults in the study area that cut through the salt layer (Albufera, Valencia and Cullera Fault). They appear to be disturbed in their geometry by the Keuper unit in some instances (profiles in Fig. 6a and b), while in other profiles they show a more classic, continuously listric geometry of upper crustal normal faults. The authors estimate 3D fault geometries, slip rates and attempt to distinguish tectonic deformation from salt withdrawal.
So far, this study applies careful analysis and reasoning and is a solid and relevant contribution for understanding a very interesting setting. In my opinion, it is necessary to comment on the seismotectonic implications of the fault geometry. How can we explain the presence and possible interaction of perpendicular normal faults? Could the authors name and compare to other places on Earth where perpendicular normal faults are present? Are there hidden strike slip faults (not easy to identify in sections) necessary to reconcile the geometry (apart from a right-lateral component on the Albufera fault)? Are all faults active at the same time (in particular the Cullera fault is not consistent with available focal mechanisms, and according to throw-depth plots, activity of the Cullera fault has decreased significantly in the Quaternary)?
In contrast, I find the second part of the manuscript (chapters 7 and 8) problematic. Here, the authors shift from the general analysis to the specific focus of mechanical layering and seismic potentials. They distinguish rupture of the entire seismogenic crust or only the suprasalt layers. They hypothesise that the weak layer could limit the width of seismic ruptures, in which case the standard scaling relationships to relate fault length to magnitude no longer apply. This is a valid and interesting question, however in my opinion, some decisions taken by the authors are not justified and major revision is required in this part.
My first concern is the downdip width for suprasalt ruptures. The value of 7.1 km chosen by the authors allows for significant magnitudes to be obtained for suprasalt ruptures with reasonable aspect ratio, however the value should be much smaller according to the profiles shown in Fig. 6. The depth of the Keuper unit at the fault locations is between 1.5 and 3 km, and fault dip is 60 degrees or larger (not 45 degrees, as the authors assume), which rather suggests average downdip width of 3 km or less. Suprasalt ruptures could reach ~15 km^2 and magnitudes slightly above 5, which is interesting, but not relevant to define the seismic potential.
Furthermore, it is not clear if ruptures in the Valencian Trough do actually nucleate in the thin suprasalt layer. The comparison with the mechanical layering in Iran is not helpful here, as the tectonic regime is compressive and the fault geometry is different, the salt layer is significantly deeper, and yet the relevant earthquakes in the Zagros belt do nucleate in the subsalt layer, NOT in the suprasalt layer (i.e., they may be contained below the salt layer, or occasionally involving both, sub and suprasalt units).
Unfortunately, available depth estimates are insufficient to resolve this issue. Those are offshore earthquakes, so we have a lack of short distance observations, and the fact that the catalogue reaches back to 1950 (Fig. 2) does not make it any better… In consequence, the disclaimer in the manuscript “depths assigned to these earthquakes present high uncertainties” should lead us to acknowledge that there is no evidence for significant suprasalt earthquakes in the study area. If we decide to believe the depths in Fig. 2, it even appears as if the suprasalt layer was depleted of earthquakes.
More fundamentally, the presented sections show only limited evidence for mechanical and seismogenic decoupling between subsalt and suprasalt units. Among the shown profiles the Cullera fault in Fig. 6a, b displays a listric profile that may be flattening out in the Keuper, supporting the presence of such a mechanical detachment horizon. However, other faults cross the Keuper with more or less planar geometry, not suggesting any decoupling. In general, the salt horizon is discontinuous under some of the profiles and varies in thickness, and some reasoning or comparison to other settings should be provided to strengthen the mechanical model.
Minor comments:
Introduction: line 75 “these relationships do not consider the influence of potential heterogeneities within the seismogenic crust”: these relationships are empirical, which means they DO consider the real influence of all heterogeneities that were present around the seismic ruptures that were evaluated to build them. The point is rather that the mean regression we use averages out those heterogeneities…
Fig. 1: Could you put labels (dates) on the focal mechanisms?
Fig 2, depth histogram: Looks like depth 0 is excluded (which may mean explosions or fixed depth), please state explicitly. Y-axis annotations are not equidistant and difficult to associate, please check.
I’d suggest rename section 3, “Methods”, to “Data”, because that’s clearly the main focus here.
Fig. 3: Could it make sense to distinguish modern lines and vintage lines in the figure? (from your description in the text, I guess I can, but if you use different line style it would be easier).
On p. 20 you say the Cullera Fault is a planar fault (inside the basement) without significant tilting of the basement-top horizon, however this is in contradiction with the interpretation in figure 6e, where you draw a listric fault that would flatten out near 5 km and shows tilted units in the hanging wall.
In Section 5.3, you claim to “define for the first time the Valencia Fault”, but I wonder if the Valencia Fault could be equivalent to the fault labelled as Central-Western Cabo Cullera Fault in Fig. 1 / Qafi?
In chapters 6 and 7, some explanations are quite long and complicated and some of the information appears repeatedly. Please check writing of this part.
Citation: https://doi.org/10.5194/egusphere-2025-1782-RC2 -
AC2: 'Reply on RC2', Ivan Martin-Rojas, 24 Oct 2025
The authors present a detailed model for stratigraphy, lithology and structure of the southwestern Valencia Trough, based on vintage and modern seismic reflection lines, well data, electric data, gravity data, and taking into account previous studies. The presented reflection lines appear well readable, and the interpretation appears clear or at least reasonable. A crucial feature of the model is a layer of evaporitic Upper Triassic (Keuper) rocks, that behaves mechanically weak and may induce decoupling of the Palaeozoic–Middle Triassic basement from the Upper Triassic –Quaternary units.
The authors identify three major faults in the study area that cut through the salt layer (Albufera, Valencia and Cullera Fault). They appear to be disturbed in their geometry by the Keuper unit in some instances (profiles in Fig. 6a and b), while in other profiles they show a more classic, continuously listric geometry of upper crustal normal faults. The authors estimate 3D fault geometries, slip rates and attempt to distinguish tectonic deformation from salt withdrawal.
So far, this study applies careful analysis and reasoning and is a solid and relevant contribution for understanding a very interesting setting. In my opinion, it is necessary to comment on the seismotectonic implications of the fault geometry. How can we explain the presence and possible interaction of perpendicular normal faults? Could the authors name and compare to other places on Earth where perpendicular normal faults are present? Are there hidden strike slip faults (not easy to identify in sections) necessary to reconcile the geometry (apart from a right-lateral component on the Albufera fault)? Are all faults active at the same time (in particular the Cullera fault is not consistent with available focal mechanisms, and according to throw-depth plots, activity of the Cullera fault has decreased significantly in the Quaternary)?
There appears to be some misunderstanding regarding the kinematics of the faults. In lines 449–465, we present data supporting the interpretation of the Albufera Fault as an oblique structure, with both right-lateral strike-slip and normal components.
However, since the lateral component cannot be observed in the seismic profiles, this aspect may not have been clearly conveyed in the manuscript. To address this, in the revised version we will explicitly state that the Albufera Fault is an oblique structure, characterized by a dominant strike-slip component and a secondary normal offset.
Additionally, we will include kinematic indicators—arrows for strike-slip motion and ticks for normal faulting—on the fault traces in all relevant figures to improve clarity.
In contrast, I find the second part of the manuscript (chapters 7 and 8) problematic. Here, the authors shift from the general analysis to the specific focus of mechanical layering and seismic potentials. They distinguish rupture of the entire seismogenic crust or only the suprasalt layers. They hypothesise that the weak layer could limit the width of seismic ruptures, in which case the standard scaling relationships to relate fault length to magnitude no longer apply. This is a valid and interesting question, however in my opinion, some decisions taken by the authors are not justified and major revision is required in this part.
My first concern is the downdip width for suprasalt ruptures. The value of 7.1 km chosen by the authors allows for significant magnitudes to be obtained for suprasalt ruptures with reasonable aspect ratio, however the value should be much smaller according to the profiles shown in Fig. 6. The depth of the Keuper unit at the fault locations is between 1.5 and 3 km, and fault dip is 60 degrees or larger (not 45 degrees, as the authors assume), which rather suggests average downdip width of 3 km or less. Suprasalt ruptures could reach ~15 km^2 and magnitudes slightly above 5, which is interesting, but not relevant to define the seismic potential.
We agree with Reviewer 2 that some of our calculations for fault downdip width should be revisited. In the initial version of our manuscript, we used a simplified approach by assuming a uniform thickness of the supra-salt succession across the entire study area. However, this represents an oversimplification of our own data.
As previously mentioned, the Albufera Fault exhibits oblique kinematics, with a dominant strike-slip component and a minor normal component. As a result, the thickness of the supra-salt succession in the downthrown block of the Albufera Fault is thinner than in the other two active structures. Similarly, for the Valencia Fault, assuming a 5 km thickness for the supra-salt succession is also an overestimation. Consequently, the downdip width and the resulting estimates of seismogenic potential should be recalculated for these two faults.
Regarding the Cullera Fault, the total thickness of the supra-salt succession at its location is approximately 4.8 km (see Fig. 6a, b, and c). We used 5 km for rounding purposes, but if the manuscript is accepted for publication, we can easily revise our calculations using 4.8 km instead. That said, we would like to emphasize that the 200 m difference between 4.8 and 5 km falls within the margin of error for depth conversion and implies a level of precision that exceeds the accuracy supported by the data.
Regarding the apparent discrepancy between the fault dip used in our calculations (45º) and the dip angle suggested by Reviewer 2 (60º) we would like to emphasize that seismic sections shown in figure 6 are vertically exaggerated. This exaggeration is indicated in the figure by the difference between the vertical and horizontal graphic scales. It is true that we did not clearly stated this fact in the figure caption. Therefore, we will add an annotation indicating this exaggeration and a sentence to caption of figure 6 to avoid confusion.
Furthermore, it is not clear if ruptures in the Valencian Trough do actually nucleate in the thin suprasalt layer. The comparison with the mechanical layering in Iran is not helpful here, as the tectonic regime is compressive and the fault geometry is different, the salt layer is significantly deeper, and yet the relevant earthquakes in the Zagros belt do nucleate in the subsalt layer, NOT in the suprasalt layer (i.e., they may be contained below the salt layer, or occasionally involving both, sub and suprasalt units).
Both reviewers expressed concerns about our comparison between the Zagros Belt and the study area, due to differences in tectonic regime and the depth of the mechanically weak layer. We acknowledge that this indicates the comparison is not sufficiently supported in our manuscript.
Therefore, as mentioned in our response to Reviewer 1, we believe it is best to remove this section from the manuscript. As previously noted, we consider that a comparison with central Italy would be more relevant and meaningful for our discussion.
Unfortunately, available depth estimates are insufficient to resolve this issue. Those are offshore earthquakes, so we have a lack of short distance observations, and the fact that the catalogue reaches back to 1950 (Fig. 2) does not make it any better… In consequence, the disclaimer in the manuscript “depths assigned to these earthquakes present high uncertainties” should lead us to acknowledge that there is no evidence for significant suprasalt earthquakes in the study area. If we decide to believe the depths in Fig. 2, it even appears as if the suprasalt layer was depleted of earthquakes.
We agree with Reviewer 2 that there is no direct evidence of significant supra-salt earthquakes in our study area. However, this region is characterized by low to moderate seismicity, with expected recurrence intervals for significant events likely spanning hundreds to thousands of years.
Therefore, we believe that the absence of supra-salt earthquakes in the instrumental record should not be interpreted as evidence that such shallow events cannot occur—absence of evidence is not evidence of absence. That said, we acknowledge that the probability of occurrence for these supra-salt earthquakes is low.
It seems this point was not sufficiently emphasized in the initial version of our manuscript. As mentioned in our response to Reviewer 1, we intend to add a sentence explicitly stating that events nucleated within the supra-salt succession are considered unlikely.
More fundamentally, the presented sections show only limited evidence for mechanical and seismogenic decoupling between subsalt and suprasalt units. Among the shown profiles the Cullera fault in Fig. 6a, b displays a listric profile that may be flattening out in the Keuper, supporting the presence of such a mechanical detachment horizon. However, other faults cross the Keuper with more or less planar geometry, not suggesting any decoupling. In general, the salt horizon is discontinuous under some of the profiles and varies in thickness, and some reasoning or comparison to other settings should be provided to strengthen the mechanical model.
As Reviewer 2 notes, one of the key pieces of evidence for mechanical decoupling between the supra- and sub-salt units is the listric geometry of faults that flatten within the mechanically weak layer. This is clearly observed in the case of the Cullera Fault, as shown in Figure 6 and highlighted by Reviewer 2.
The same applies to the Valencia Fault. In Figure 6a, the fault trace is represented by the eastern termination of the Mesozoic–Quaternary succession, which displays a curved (listric) geometry that also flattens within the Keuper evaporitic unit. Similarly, the Albufera Fault, as shown in Figures 6d and 6e, exhibits a listric geometry. Although the quality of the seismic profiles makes it difficult to identify the detachment horizon with certainty, we hypothesize that the Albufera Fault also flattens within the Keuper unit, which is a well-established regional detachment level (Guimera, 2018; Guimerà & Álvaro, 1990; Muñoz Martín & Vicente Muñoz, 1988).
Reviewer 2’s comment highlights the need to better address and discuss the mechanical decoupling between the supra- and sub-salt successions in our manuscript. We will revise the relevant sections to clarify this point using the observations and references mentioned above.
Guimera, J. (2018). Structure of an intraplate fold-and-thrust belt: The Iberian chain. A synthesis. Geologica Acta, 16(4), 427 – 438. https://doi.org/10.1344/GeologicaActa2018.16.4.6
Guimera, J., & Alvaro, M. (1990). Structure et evolution de la compression alpine dans la Chaine iberique et la Chaine cotiere catalane (Espagne). Bulletin de La Société Géologique de France, VI(2), 339–348. https://doi.org/10.2113/gssgfbull.VI.2.339
Muñoz Martín, A., & Vicente Muñoz, G. De. (1988). Cuantificación del acortamiento alpino y estructura en profundidad del extremo sur-occidental de la Cordillera Ibérica (sierras de Altomira y Bascuñana). Revista de la Sociedad Geológica de España. http://hdl.handle.net/20.500.14352/57932
Minor comments:
Introduction: line 75 “these relationships do not consider the influence of potential heterogeneities within the seismogenic crust”: these relationships are empirical, which means they DO consider the real influence of all heterogeneities that were present around the seismic ruptures that were evaluated to build them. The point is rather that the mean regression we use averages out those heterogeneities…
Reviewer 2 is right. We will rephrase this sentence following his/her suggestion.
Fig. 1: Could you put labels (dates) on the focal mechanisms?
We will add the labels.
Fig 2, depth histogram: Looks like depth 0 is excluded (which may mean explosions or fixed depth), please state explicitly. Y-axis annotations are not equidistant and difficult to associate, please check.
We will correct the figure and the figure caption following Reviewer 2 indications.
I’d suggest rename section 3, “Methods”, to “Data”, because that’s clearly the main focus here.
We will rename this section following the reviewer comment.
Fig. 3: Could it make sense to distinguish modern lines and vintage lines in the figure? (from your description in the text, I guess I can, but if you use different line style it would be easier).
We will follow this suggestion, using different line styles for vintage and modern seismic lines.
On p. 20 you say the Cullera Fault is a planar fault (inside the basement) without significant tilting of the basement-top horizon, however this is in contradiction with the interpretation in figure 6e, where you draw a listric fault that would flatten out near 5 km and shows tilted units in the hanging wall.
It should be noted that seismic profile in figure 6e (and 6d) is not as deep as those in figures 6a-6c. Therefore, the basement cannot be observed in profiles 6d and 6e. Actually, as we mentioned above, fault traces depicted in figures 6d and 6e flattens out in the suprasalt succession. We will add a mention to this point in the figure caption to avoid confusion.
In Section 5.3, you claim to “define for the first time the Valencia Fault”, but I wonder if the Valencia Fault could be equivalent to the fault labelled as Central-Western Cabo Cullera Fault in Fig. 1 / Qafi?
We agree with Reviewer 2 that our Valencia Fault may correspond to the Western Cabo Cullera Fault shown in Figure 1. However, since we present a new structural interpretation for the entire study area, we believe that retaining the name "Western Cabo Cullera" could lead to confusion. This is because that name has previously been associated with other faults in the sector (e.g., Central Cabo Cullera, Eastern Cabo Cullera, etc.).
For this reason, we chose to assign the name "Valencia Fault" to this structure. Before making this decision, we consulted the author who originally defined the Western Cabo Cullera Fault, and he agreed with our naming choice (Perea, pers. comm.).
Nevertheless, we acknowledge that the potential correlation between our Valencia Fault and the Western Cabo Cullera Fault should be mentioned. Therefore, we will include a sentence in the manuscript to clarify this point.
In chapters 6 and 7, some explanations are quite long and complicated and some of the information appears repeatedly. Please check writing of this part.
Citation: https://doi.org/10.5194/egusphere-2025-1782-AC2
-
AC2: 'Reply on RC2', Ivan Martin-Rojas, 24 Oct 2025
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 674 | 83 | 23 | 780 | 20 | 34 |
- HTML: 674
- PDF: 83
- XML: 23
- Total: 780
- BibTeX: 20
- EndNote: 34
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
Dear Editor of NHESS and authors,
Here you can find my comments regarding the work entitled:I
Influence of crustal mechanical layering on the seismic potential of active faults: insights from the southwestern Valencia Trough (W Mediterranean)
By: Martin-Rojas, Iván; Ramos, Adriá; De Ruig, Menno; Medina-Cascales, Iván; Santamaría-Pérez, Eva; Alfaro, Pedro
Overall comments:
This manuscript presents an excellent study describing active faulting with seismogenic potential in the Valencia Trough, eastern Spain. The analysis of 2D seismic profiles is particularly valuable, given both the total profile length and the extensive dataset processed. Despite the limitations posed by the vintage seismic profiles, the authors achieve an impeccable interpretation using 2D Move software. In this regard, the work represents a significant step forward in the characterization of potential seismic sources that could affect the densely populated eastern coast of Spain. The application of 2D software for estimating fault areas in active tectonic settings offers a realistic approach to seismic hazard assessment, particularly in offshore regions where constraining the 3D geometry of faults remains challenging.
Some doubts arise regarding the methodology used to estimate the maximum earthquake magnitude based on fault lengths, as well as the consistency of the terminology used throughout the manuscript. I think that authors have to explain better the methodology and the “correction factor” to avoid misinterpretations. Please refer to my specific comments above for details.
Greater doubts are in the comparison between the Valencia Trough and the Zagros region. Please read above in main suggestions.
Main suggestions
S1-I suggest that the authors consider including the phrase “salt-influenced tectonic setting” in the title of the manuscript. This represents the main novel contribution of the work, along with the parameterization of the three active faults in the Valencia Trough.
S2-Section 4. Stratigraphy. I recommend reducing this section, as it does not represent the main objective of the work. Providing such a detailed stratigraphic description may distract the reader from the key focus of the study and could make the narrative less clear. A more concise summary would keep the manuscript streamlined and focused.
S3-Some figures and graphics are difficult to read, and I have included some recommendations into the attached pdf file. For instance, letters with a grey color, small numbers in the isopach lines, and similar colors I plots for the SR long-term estimation.
S4- I have some major concerns regarding the comparison between the Valencia Trough and the Zagros Thrust and Belt based solely on the presence of a “mechanically weak layer.” This comparison overlooks several tectonically significant parameters, starting with the fundamentally different tectonic settings: The Zagros is characterized by crustal shortening and active mountain building, whereas the Valencia Trough corresponds to a uniaxial extensional marine basin. In addition, isostatic processes are not considered, and the rupture mechanisms are inherently different. Reverse faulting in the Zagros involves the mobilization of very large rock volumes in the hanging wall, which is not directly comparable to normal faulting in a sedimentary marine basin. A more cautious discussion is required to avoid oversimplifying the analogies between these two geodynamically contrasting regions. Thermal coupling in different depth are also ignored.
S5- You have not demonstrated this claim (similarities between both tectonic settings, Valencia Trough and Zagros region). What is shown is only a convergence in the geometry of the spatial distribution of seismicity, but this does not constitute a “pattern” based on the underlying mechanical behavior of the faults. A true mechanical pattern would require evidence of consistent fault kinematics, stress distribution, and deformation mechanisms, which are not demonstrated in the current version of the manuscript. Please, consider to remove this part, I have no the relevance in your paper. This only distracts from the methodology that authors propose.
S6- I also recommend removing the estimation of the maximum earthquake magnitude based solely on the total fault length, as the results appear unrealistic (Table 1). This approach tends to overestimate the potential magnitude and is not fully supported by seismotectonic principles. Instead, relying on fault surface area provides a more consistent and physically based estimate, as it is directly related to the seismic energy release (Kanamori and Anderson).
S7-Conclusions have to be revisited. Some relevant information derived from your results are missed.
Technical suggestions and questions
TSQ1-Page 4; Line 85 You refer “mechanical weak layer” to the Triassic evaporitic formation. My question: Are you referring to ductile behavior under a tectonic stress field? If so, how do you estimate the creep movement associated with ductile salt deformation?
TSQ2- Reading the whole manuscript, I still do not fully understand what you mean by the term “mechanically weak layer” regarding the evaporite Traissic layer. Could you please clarify its implication for the seismogenic behavior of the fault system?
Are you suggesting that such a weak layer:
A clearer explanation would help in understanding the role of this layer in the seismic cycle.
TSQ3- The Cullera Fault appears to show a decrease in slip rate from the Pliocene to the Quaternary (Fig. 8). Could the authors clarify the underlying reason for this apparent behavior? Do you have a conceptual model that could explain this trend? Is it a Nubia plate deceleration?
TSQ4- The curves of the Figure 10, Throw depth (T-z) and expansion index plots for the Cullera Fault computed from profiles SGV01-115 and SGV01-117, show a fault throw variation. Could the geometry of these curves be explained by the listric geometry of the Albufera Fault? Variations in fault throw along planes with different dips could potentially generate this type of curve. Could the authors clarify if this interpretation is consistent with their data?
TSQ5- The accommodation space created by the Albufera Fault, and cited in the line 559, page 30, can be explained by salt- tectonics? Please explain your model.
TSQ6- Line 641 Tectonic driving mechanism and differences between tectonic driving and sal-withdrawal driving. I have some doubts regarding this point. The “mechanically weak layer” is described at a depth of 2–5 km, but tectonic earthquakes with significant magnitudes (M6–7) typically require rupture depths greater than this. Moreover, the regional tectonic strain is expected to be distributed below the weak layer as well. This raises the question: does the fault release only the strain accumulated in the supra-salt section, or can it also rupture deeper segments (>10 km) where larger tectonic strain is stored? A discussion of this point would clarify the seismogenic potential of the fault system.
TSQ7- Methodology and the “correction factor” (line 712, page 36). I agree that using the total fault length to estimate earthquake magnitude through common empirical relationships is not the most reliable approach. However, the total fault surface area provides a more reasonable approximation because it is directly related to the seismic energy released during an earthquake (see Kanamori and Anderson). The main concern with introducing a “correction factor” is the potential to generate an artifact in the magnitude estimation, which could lead to biased results if not properly justified or calibrated. On the other hand, why the “aspect radio” factor of the active fault could improve the magnitude estimation? It is not clear in the text.
TSQ8- PAY ATTENTION! revise the fault areas of the Table 2. There are some incongruences with the values obtained from the 2D Move analyses and the values indicated in the Table 2 for the Cullera, Albufera and Valencia faults (pages 20, 24 and 27).
TSQ9- Table 2. I don’t understand how a great fault area triggers an earthquake of lower magnitude as appeared in Albufera and Valencia faults). Please, revise carefully TABLE 2.
Additionally, I have attached the PDF manuscript with some minor detailed comments posted into the PDF file.
OVERALL RECOMMENDATION: MAJOR REVISIONS.
Thank you very much for your work, excellent analysis and results. I hope that authors improve the manuscript.
Thanks also for the editor for the consideration as a reviewer.
Dr Raul Pérez-López
Instituto Geológico y Minero de España – IGME CSIC
Email: r.perez@igme.es