# Revealing the hidden signature of fault slip history in the morphology of degrading scarps

Regina Holtmann, Rodolphe Cattin, Martine Simoes, Philippe Steer

► To cite this version:

Regina Holtmann, Rodolphe Cattin, Martine Simoes, Philippe Steer. Revealing the hidden signature of fault slip history in the morphology of degrading scarps. Scientific Reports, 2023, 13, pp.3856. 10.1038/s41598-023-30772-z. insu-04021048

HAL Id: insu-04021048

https://insu.hal.science/insu-04021048

Submitted on 9 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Revealing the hidden signature of fault slip history in the morphology of degrading scarps

Regina Holtmann¹,³,⁴, Rodolphe Cattin²,⁴✉, Martine Simoes¹,⁴ & Philippe Steer³,⁴

Active faults accommodate tectonic plate motion through different slip modes, some stable and aseismic, others characterized by the occurrence of large earthquakes after long periods of inactivity. Although the slip mode estimation is of primary importance to improve seismic hazard assessment, this parameter inferred today from geodetic observations needs to be better constrained over many seismic cycles. From an analytical formulation developed for analyzing fault scarp formation and degradation in loosely consolidated material, we show that the final topographic shape generated by one earthquake rupture or by creep (i.e., continuous slip) deviates by as much as 10–20%, despite a similar cumulated slip and a constant diffusion coefficient. This result opens up the theoretical possibility of inverting, not only the cumulated slip or averaged slip rate, but also the number of earthquakes and their sizes from scarp morphologies. This approach is all the more relevant as the number of rupture events is limited. Estimating the fault slip history beyond a dozen earthquakes becomes very difficult as the effect of erosion on scarp morphology prevails. Our modeling also highlights the importance of trade-offs between fault slip history and diffusive processes. An identical topographic profile can be obtained either with a stable fault creep associated with rapid erosion, or a single earthquake rupture followed by slow erosion. These inferences, derived from the simplest possible diffusion model, are likely to be even more pronounced in nature.

Scarps are landscape features forming a steep break in the local or regional topographic slope. They relate to fluvial incision (i.e., terrace risers) or faulting (i.e., fault scarps) and can be of different sizes, from ca. 1–10 m high scarps for terrace risers or recent seismic scarps, and up to > 100 m high scarps for mature faults. Scarps result from a relative vertical movement (river incision in the case of terrace risers, tectonic uplift in the case of fault scarps). They are degraded by surface processes, leading to a decrease in their steepness over time. This latter observation provided the long-lived basis for the relative dating of these landscape features¹ and later for various attempts at quantifying their age. An illustration is provided from field and satellite views of the Kongur Shan extensional fault system²,³ (Fig. 1). Bedrock has been exhumed along the large (> 1000 m high) cumulated fault scarp that forms the mountain front. The overall slope of this cumulated scarp decreases with altitude, in line with the progressively increasing age of bedrock exposure and exhumation³ (Fig. 1, left picture). Along strike, the fault trace is marked by more subtle (1–20 m high) and younger scarps that disrupt recent deposits (Fig. 1, right and left pictures). Based on their offset and age, these scarps have been used to quantify Quaternary fault slip rates². Among these, the older higher scarp disrupting glacio-fluvial deposits (Fig. 1, left) appears less steep than the younger and smaller one within a more recent alluvial fan (Fig. 1, right). Altogether, these various features are clear indications that the fault is active. However, their simple inspection does not permit at the moment to assess how fault slip is accommodated over time, during earthquakes or by creep.

In loosely consolidated material (e.g., fluvial deposits, alluvial fans—as in the case of the most recent scarps in Fig. 1), slope-dependent processes favor the erosion of the higher and steeper parts of a scarp, while deposition of the eroded material occurs downhill, with a flux of sediment overall proportional to the local slope. These processes can be modeled by the progressive diffusion of a steep topographic slope⁴,⁵, providing from there a conceptual framework to constrain the diffusion age $\kappa t$ of such scarps and, in turn, infer their age $t$ or their

¹CNRS, Institut de physique du globe de Paris, Université Paris Cité, 75005 Paris, France. ²CNRS, Géosciences Montpellier, Université de Montpellier, 34000 Montpellier, France. ³CNRS, Géosciences Rennes, Université de Rennes, 35000 Rennes, France. ⁴These authors contributed equally: Regina Holtmann, Rodolphe Cattin, Martine Simoes and Philippe Steer. ✉email: rodolphe.cattin@umontpellier.fr
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Normal fault scarp from the Kongur Shan extensional fault system<sup>2,3</sup>, using satellite and field views. A large (> 1000 m high) cumulated scarp, forming the mountain front, has exhumed bedrock over several million years. The overall slope of this long-term scarp decreases with altitude, as the exposure age of bedrock increases (left picture). Along strike, the normal fault disrupts more recent Quaternary deposits, such as glacio-fluvial sediments (lower left picture) or alluvial fans (lower right picture) and forms more subtle (1–20 m high) scarps. Erosion of scarps by diffusion, as considered hereafter in this study, only applies to these most recent features in loosely non-consolidated materials. Satellite general view from © Google Earth (2023 Maxar Technologies and Landsat/Copernicus), dating from July 9th 2007. Field views taken on October 18th 2007, along the road G314 from Kashgar to Tashgorgan (38° 43' 42.97" N, 75° 0' 59.51" E). The contrast of the original field pictures has been increased using Adobe Photoshop for an easier visibility of the landscape features. The original non-interpreted satellite view and field pictures are provided in supplementary material (Figs. S3 to S5). Figure generated using Adobe Illustrator.</figcaption></figure>

coefficient of diffusion $κ$. This approach has been applied for the dating of terrace risers and associated terrace abandonment and incision<sup>6–8</sup> or to provide chronological constraints on seismic fault scarps assuming that they formed after one single earthquake<sup>1,5,9,10</sup>, even though some complexities to these various approaches on the analysis of diffused scarps have been pointed out: some studies<sup>9</sup> advocate for a decrease in the diffusion coefficient $κ$ with time to reconcile the morphology of young (high $κ$) and old (low $κ$) fault scarps and climate warming and the shift from glacial to inter-glacial are generally suggested to explain possible variations<sup>10</sup>.

Interestingly, in the specific case of fault scarps, most do not result from a single earthquake rupture but from several consecutive ruptures or even creep. Multiple fault ruptures or creep will rejuvenate the overall fault scarp and maintain it steep (or relatively steeper) over time<sup>5,6,9</sup>. This latter point raises the possibility that there could be a hidden morphological signature of whether a fault scarp was built by one or multiple earthquake ruptures or by creep. Active faulting leaves a morphological signature (Fig. 1) that has been amply used in numerous studies to quantify active tectonics<sup>11</sup>. However, whether landscapes keep a record of how active faults slip, either seismically or aseismically, has remained largely unexplored<sup>12,13</sup>. With the increasing availability of high-resolution topographic data (< 1 m resolution), there is a timely opportunity to investigate this issue<sup>14,15</sup>.

Here, we aim to illuminate fundamental aspects of the influence of tectonic slip history, seismic or aseismic, on the morphology of degrading fault scarps. We investigate and explore this idea by considering the simple case of scarps in non-cohesive materials, such as those disrupting alluvial fans or fluvial deposits (see recent scarps pointed out in Fig. 1), where diffusion-like erosion is the primary process degrading the scarps. We first develop a new analytical formulation for scarp evolution, including fault slip by earthquake ruptures or by creep, coupled to linear topographic diffusion, which is considered as the simplest model of scarp degradation<sup>5,8,16</sup>. We then systematically investigate the morphological signature of fault slip history, including the role of the number of ruptures, the coefficient of diffusion, and the time since the onset of fault activity. Last, we revisit with this model the time dependency of diffusion coefficient associated with a previously documented field case in Central Greece to reveal the potential role of fault slip history in scarp morphology.

## Method

**Fault uplift and scarp degradation model.** In our approach, for the sake of simplicity, we consider a very steep fault (dip angle → 90°) for which we neglect the effect of horizontal tectonic advection. Hence, the elevation profile $z$ across the scarp is initially:

$$z(x < 0, t = 0) = 0 \quad \text{and} \quad z(x \ge 0, t = 0) = u, \qquad (1)$$

with $x$ the horizontal distance from the fault trace, positive towards the uplifting block (hanging-wall for a reverse fault, footwall in the case of a normal fault), $t$ the time, and $u$ the uplift leading to the height of the initial scarp.
We consider different fault slip and uplift scenarios, which all reach the same final total uplift: (1) a fault scarp formed during one single earthquake rupture; (2) a fault scarp resulting from multiple identical and periodical earthquake ruptures; and (3) a fault scarp formed by continuous slip on a creeping fault. In all scenarios, the fault scarp is subjected to surface processes, modeled here using the well-established linear diffusion model$^{5,8,16}$.

$$ \frac{\partial z}{\partial t} = \kappa \frac{\partial^2 z}{\partial x^2}, \qquad (2) $$

where $κ$ is a diffusion coefficient, assumed here constant through space and time. This model is generally considered relevant to describe erosion of loosely unconsolidated rocks subjected to creep, bioturbation and rainsplash on low-gradient slopes.

**Analytical solution for a single fault uplift.** The solution to Eq. (2) for the degradation of a topographic step related to a single vertical rupture, as defined by the initial conditions in Eq. (1), is$^5$:

$$ z(x, t) = u \times \left[ 1 + \operatorname{erf}\left(\frac{x}{\sqrt{4\kappa t}}\right) \right] / 2 \qquad (3) $$

where $u$ is the height of the rupture, $t$ is the time elapsed since the rupture, and $erf$ is the error function defined as

$$ \operatorname{erf}(\eta) = \frac{2}{\sqrt{\pi}} \int_{0}^{\eta} e^{-\gamma^2} d\gamma. \qquad (4) $$

Figure 2a illustrates the diffusion of a fault scarp formed after one earthquake rupture (i.e., one uplift event). Erosion and sedimentation modify rapidly and significantly the topography near the fault. Over time, as the scarp profile degrades and flattens, the area affected by surface processes widens, and the scarp degradation rate decreases.

Equation (3) can be expressed in a dimensionless form with the dimensionless parameters $z/u$ and $\eta = x/\sqrt{4\kappa t}$.

$$ \frac{z(x, t)}{u} = [1 + \operatorname{erf}(\eta)]/2. \qquad (5) $$

The resulting dimensionless diffusion profile across the scarp is independent of both the uplift magnitude and the time elapsed since the uplift event (Fig. 2b). In the case of a single earthquake rupture, the final topographic profile that diffused over a total time-lapse $t_{end}$ can be calculated using the initial (and total) vertical offset as $= u_{end}$:

$$ \frac{z(x, t_{end})}{u_{end}} = [1 + \operatorname{erf}(x/\sqrt{4\kappa t_{end}})]/2. \qquad (6) $$

**Analytical solution for multiple fault uplifts or creep.** As the considered diffusion degradation is uni-dimensional and linear, the scarp profiles resulting from multiple events can be constructed through a com-

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Evolution of a degrading fault scarp. (a) Elevation profiles $z(x)$ across a fault scarp assuming a single rupture uplift $u = 2.5$ m and a diffusion coefficient $κ = 10^{-4}$ m$^2$/yr. Blue, orange, and green lines show the degraded scarp profiles after 100, 5000, and 10,000 years, respectively. (b) Same profiles using the dimensionless parameters $\eta$ and $z/u$. Figure generated using the Python's matplotlib library.</figcaption></figure>
bination of diffused incremental topographies. The case of a double-uplift scenario, with diffusion during a total duration of $t_{end}$, is illustrative (See Supplementary information S1 and Fig. S1). Two uplift events, each with a vertical offset of $u_{end}/2$ occur at $t = 0$ and $t = t_{end}/2$. The final topographic profile results from the sum of two diffused incremental topographies: one, associated with the first event at $t = 0$, diffuses throughout the whole model duration $t_{end}$, whereas the other, related to the second event at $t = t_{end}/2$, diffuses during a time-lapse of $t_{end}/2$ only. The initial topographies being similar (i.e., a Heaviside step function with a vertical offset of $u_{end}/2$), the formulation for the resulting profile is:

$$ \frac{z(x, t_{end})}{u_{end}} = \frac{1}{2} \left( \left[ 1 + \operatorname{erf}\left(x/\sqrt{4\kappa t_{end}/2}\right) \right] / 2 + \left[ 1 + \operatorname{erf}\left(x/\sqrt{4\kappa t_{end}}\right) \right] / 2 \right). \quad (7) $$

By extension, the degradation of a scarp formed during a multiple-uplift scenario, with periodical ruptures, is given by:

$$ \frac{z(x, t_{end})}{u_{end}} = \frac{1}{2 \times n} \sum_{i=1}^{n} \left[ 1 + \operatorname{erf}\left(x/\sqrt{4\kappa t_{end} \times \frac{i}{n}}\right) \right], \quad (8) $$

with $n$ the number of uplift events. In the case of a continuous uplift scenario (i.e., creep), $n$ tends to infinity, and each event generates a vertical uplift of $u_{end}/n$. Here, we simulate the creeping behavior using Eq. (8) with $n = 10^4$. This formulation gives similar results to those obtained from Eq. (13) of Hanks⁵ (Fig. 3a), in which the scarp morphology of a creeping fault is the solution of the equation

$$ \frac{\partial z}{\partial t} = \kappa \frac{\partial^2 z}{\partial x^2} + A, \quad (9) $$

where $A = u_{end}/t_{end}$ is the uplift rate (see supplementary material S2 for mathematical formulations).

## Results

**Effect of uplift scenarios.** Using Eq. (8), we first calculate the diffusion profiles associated with the deg radation of a scarp generated by $n$ total identical events, with $n$ ranging from 1 (single rupture) to $10^4$ (creep) (Fig. 3). Although the cumulative uplift equals $u_{end}$ at the final time $t_{end}$ in all scenarios, each uplift scenario pro duces a topographic profile that significantly differs from the others (Fig. 3a). The greater the number of events, the steeper the slope near the fault.

Compared to the single uplift scenario, the deviation between the resulting diffusion profiles is maximal at the dimensionless horizontal distance of $\eta = \pm \eta_{max}$ from the fault (Fig. 3b). The peak-to-peak deviation reaches

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Influence of the number of rupture events on scarp profiles associated with a total uplift of $u_{end}$ and a final time-lapse of $t_{end}$. The thick and dashed black lines give the result for the single-event scenario and a creeping fault, respectively. The color scale is related to the number of events. (a) Calculated dimensionless profiles. Hanks' formulation⁵ is given in Supplementary information S2. (b) Topographic deviation with respect to the single-event scenario. (c) Variation of the maximum deviation with the number of uplift events. Figure generated using the Python's matplotlib library.</figcaption></figure>
up to ca. 20% of $u_{end}$ for the creep scenario. The greater the number $n$ of uplift events, the smaller the difference in deviation associated with an additional uplift (Fig. 3b,c). For example, this deviation is ca. 8%, 16%, 18%, and 19% of $u_{end}$ when $n$ equals 2, 6, 12, and 22, respectively. The parameter $\eta_{max}$ is relatively constant, slightly decreasing from ca. 0.6 to ca. 0.4 when increasing $n$ (see Supplementary material S3).

A maximum deviation of a few percent may at first seem insignificant. However, rates of erosion or sedimentation by diffusion vary significantly with time. They are fast at the beginning due to greater slopes and slow down with the decrease of the topographic curvature⁶. Therefore, the observed deviation, even though modest, could significantly impact the estimate of parameter $\kappa t$ and, from there, the inferred diffusion coefficient or the age of a scarp.

**Synthetic tests.** We examine the potential to identify an uplift history (i.e., single or multiple earthquake ruptures, or creep) by analyzing a given scarp profile. This approach compares a well-defined synthetic profile to profiles produced from various scenarios by varying $n$, $\kappa t$, and topographic noise.

First, we consider a reference model obtained with $n^{ref} = 5$ ruptures and $\kappa t^{ref} = 1 \text{ m}^2$. A normal distributed random noise is added to this profile with a mean of zero and a standard deviation of $\sigma^{ref} = 10^{-3} \times u_{end}$. This noise accounts for the accuracy of topographic data and the natural heterogeneity and variability of scarp profiles. Next, we perform sensitivity tests assuming “true” models similar to the reference model, except for the tested parameter (either $n$, $\kappa t$, or $\sigma$). Finally, we perform the inversion by systematically exploring the model parameters with $1 \le n \le 20$ and $0.1 \le \kappa t \le 10 \text{ m}^2$. We calculate for each tested profile its likelihood $L$ relative to the “true” profile, defined as

$$L(n, \kappa t) = \exp \left( -\frac{1}{N} \sum_{i=1}^{N} \left[ \frac{z_i - z_i^{\text{true}}}{\sigma_i^{\text{true}}} \right]^2 \right) \quad (10)$$

with $N$ the amount of data along the resulting profiles ($N = 5000$ in our approach), $z_i$ the topographic profile associated with $n$ and $\kappa t$, $z_i^{\text{true}}$ the “true” profile, and $\sigma_i^{\text{true}}$ the noise. The likelihood distribution is normalized to the likelihood of the best-fitting profile to compare the different tests. The normalized likelihood ranges from 0 (poor agreement with the “true” profile) to 1 (good agreement).

For the reference model, we obtain good normalized likelihoods ($L > 0.9$) for $0.4 \le \kappa t \le 1.6 \text{ m}^2$ and $4 \le n \le 7$ (Fig. 4a), similar to the values of reference parameters. The obtained likelihood distribution also suggests a strong trade-off between $\kappa t$ and $n$ for models with $n < 10$. Beyond this threshold, the likelihood only depends on $\kappa t$. This result is consistent with our previous observation that the diffusion profile is only weakly affected by the number of events for $n > 10$ (Fig. 3b,c). Additionally, the likelihood of a profile resulting from one or many ruptures can be similar if $\kappa t$ is adjusted (Fig. 4a). This result confirms that even though the effect of $n$ on the resulting profiles remains modest (with a difference of only a few percent; Fig. 3), assuming a scarp formed after a single rupture—instead of 5—can lead to underestimating $\kappa t$ by a factor of two.

We test the influence of changing the number of earthquake ruptures on the likelihood distribution, with $n^{true}$ equal to 3, 5, 10, and 15, while keeping $\kappa t^{true} = 1 \text{ m}^2$ and $\sigma^{true} = 10^{-3} \times u_{end}$ (Fig. 4b–d). In each case, the best likelihood is obtained for—or close to—the “true” values of the parameters. Our results suggest that the larger the parameter $n^{true}$, the more difficult it is to retrieve this parameter from our inversions, and the easier it is to constrain $\kappa t^{true}$. Indeed, a larger $n^{true}$ leads to a wider zone with a good likelihood ($L > 0.9$) in the parameter space.

We also test the influence of varying $\kappa t^{true}$, between 0.2 and $5 \text{ m}^2$, while keeping $n^{true} = 5$ and $\sigma^{true} = 10^{-3} \times u_{end}$ (Fig. 4e–g). In each case, the best likelihood is obtained for—or close to—the “true” parameter values. Yet, we observe a wider zone associated with good likelihood, in the parameter space, for low values of $\kappa t^{true}$. Therefore, the inversion is more likely to confidently retrieve model parameters for large $\kappa t^{true}$, when diffusion has significantly modified the scarp morphology (i.e., old scarp or high diffusion rates).

To assess the impact of topographic noise on our inversion results, we perform modeling with $\sigma^{true}$ ranging between $10^{-4} \times u_{end}$ and $10^{-1} \times u_{end}$ (Fig. 4h–j). Regardless of the noise level, the best likelihood is obtained for—or close to—the “true” parameter values. Our results un-surprisingly show that noise affects the estimated range of good likelihood around the best-fit parameters, with a lower noise level leading to a narrower and better-defined range. Parameter $\sigma^{true} = 10^{-1} \times u_{end}$ gives a mostly uniform likelihood, suggesting that the approach presented here cannot be applied to the analysis of noisy topographic profiles (i.e., low data accuracy or large along-scarp variability).

**Application to the Eliki fault scarp (Greece).** The Gulf of Corinth (Central Greece) is one of the most seismically active regions in Europe¹⁷. Its rapid extension rate of 11–15 mm/year is accommodated along active normal faults bounding the rift¹⁸⁻²⁰. Kokkalas and Koukouvelas⁹ investigated the degradation of normal fault scarps formed after the 1861 (Eliki) and 1981 (Kaparelli) earthquakes, which affected the southern and eastern regions of the Gulf of Corinth (Fig. S2). By considering that the scarps degraded by linear diffusion after these earthquakes, they determined the diffusion coefficient $\kappa$ that best fitted the observed morphologies and applied a similar approach to older scarps re-interpreted from paleoseismological trenches. By combining these results with the ages of the scarps, they propose that $\kappa$ decreases exponentially with time (Fig. 5a). Interestingly, these authors noted that the $\kappa$ of composite scarps (i.e., recording multiple ruptures over a longer time) was lower when considering the overall scarp rather than only the most recent segment formed during the last earthquake, in line with the documented decrease of $\kappa$ with time. They also note that these composite scarps can be erroneously modeled as if they formed during a single earthquake.

Here, we apply our multiple-rupture approach to test the impact of this simplification. To avoid lithological effects that may affect $\kappa$, we focus on the results from the Eliki fault. This fault system is located in the
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Likelihood distribution associated with synthetic tests. The black cross corresponds to the “true” model parameters. (a) Reference model associated with a five-events scenario, a constant $\kappa t = 1 \text{ m}^2$, and a noise level of $10^{-3} \times u_{\text{end}}$. (b–d) Effect of the number of events with $n^{\text{true}}$ ranging from 3 to 15, $\kappa t^{\text{true}} = \kappa t^{\text{ref}}$ and $\sigma^{\text{true}} = \sigma^{\text{ref}}$. (e–g) Influence of the $\kappa t$ parameter, with $\kappa t^{\text{true}}$ varying from 0.2 to 5 m$^2$, $n^{\text{true}} = n^{\text{ref}}$ and $\sigma^{\text{true}} = \sigma^{\text{ref}}$. (h–j) Sensitivity to the topographic noise, with noise level between $10^{-4} \times u_{\text{end}}$ and $10^{-1} \times u_{\text{end}}$, $n^{\text{true}} = n^{\text{ref}}$ and $\kappa t^{\text{true}} = \kappa t^{\text{ref}}$. Figure generated using the Python's `` library.</figcaption></figure>

southwestern Gulf of Corinth. It consists of two north-dipping fault segments, both of ca. 15–20 km in length (Fig. S2). The Eliki sites studied by Kokkalas and Koukouvelas⁹ are located along the eastern segment. The investigated scarps are formed within poorly indurated Plio-Quaternary syn-rift fan delta conglomerates and alluvial gravels, and are therefore consistent with our modeling framework. We do not re-analyze the observations and processing performed on trench logs and scarps profiles in the following. Instead, we consider the mean diffusion coefficient $\kappa \sim 1.3 \text{ m}^2/\text{kyr}$ derived for the 1861 scarp and the one $\kappa'_{\text{eff}} \sim 0.5 \text{ m}^2/\text{kyr}$ obtained in the T13 trench and covering a time-lapse of ca. 8500 years, as derived and determined by Kokkalas and Koukouvelas⁹. Based on the estimated recurrence interval of 450 ± 109 years, these diffusion coefficients can be related to a single event and to 16–25 events, respectively (Fig. 5b).

We consider the diffusion profile associated with a single-earthquake scenario, $\kappa^{\text{true}} = 0.5 \text{ m}^2/\text{kyr}$ (old scarp), and $\sigma^{\text{true}} = 10^{-2} \times u_{\text{end}}$. We then calculate the likelihood of profiles associated with $1 \le n \le 30$ and $10^{-1} \le \kappa' \le 10^1 \text{ m}^2/\text{kyr}$. In the case of $n=16–25$, we obtain good likelihoods ($L > 0.9$) for $0.89 \le \kappa' \le 1.17 \text{ m}^2/\text{kyr}$, i.e., for $\kappa'$ similar to the average $\kappa$ estimated from the recent 1861 scarp (Fig. 5c). A diffusion parameter of $1–1.2 \text{ m}^2/\text{kr}$ constant over time can therefore explain the diffusion profiles of both recent single-event scarps and multiple-event scarps in Eliki. This result suggests that the deduced decrease of the diffusion coefficient with time might be apparent and instead reflect the influence of tectonic slip history on the morphology and shape of the composite scarp.
<figure><img src="image_7.png" /></figure>

**Figure 5.** Influence of the number of events on the diffusivity coefficient $κ$ as inferred from measurements along active faults in Greece. (a) Time variation of $κ$ along the Eliki, Kaparelli, and Skinos faults, as proposed by Kokkalas and Koukouvelas⁹. (b) Relationship between the number of events and the degradation time along the Eliki fault⁹. (c) Likelihood of scarp profiles along the Eliki fault, assuming a single-event scarp that degraded during 8500 years with $κ_{eff}$. ~ 0.5 m²/kyr as the “true model”. The vertical black line and the white square give the effective diffusivity determined in a trench, assuming a single event scenario⁹. The thick vertical dashed line is associated with $κ' = 1.05$ m²/kyr as re-assessed from the calculated likelihood distribution and taking into account the 16–25 events that potentially ruptured this fault over the last 8500 years. Small circles give the diffusivity from field measurements on the fault scarp formed after the 1861 Eliki earthquake. These individual values yield a mean diffusivity of 1.26 m²/kyr for a duration of t = 144 years (large circle and vertical gray line). Figure generated using the Python's matplotlib library.

Discussion

Our results rely on the simplest model of fault scarp evolution, with a block uplift by slip on a vertical fault and degradation by linear diffusion. Even though linear diffusion has been successfully applied to specific field settings²¹, there is ample evidence that hillslope and scarp diffusion may not be linear²²,²³. Mass wasting and non-local sediment transport processes are also ignored while likely playing a significant role after each earthquake rupture. Additionally, we consider only an initial flat topography, even though the contribution of initial topography can be added to our mathematical formulation⁶,²⁴. Given these simplifications, our modeling should not be considered as a faithful description of reality. However, it should be taken as a helpful way to isolate the specific impact of fault slip history on scarp morphology. We recall here that erosion by diffusive processes is only applicable where scarps disrupt loose and unconsolidated material such as alluvial fans or fluvial deposits.

Under these limitations, our results demonstrate that the slip history of a fault, whether it creeps or slips during one or several earthquake ruptures, leaves a potentially detectable morphological signature. For the same diffusion coefficient $κ$ and cumulated slip $u_{end}$, a creeping fault will lead to a steeper scarp than the one formed after a single earthquake rupture, with up to 20% of topographic deviation. In turn, this result implies that it is theoretically possible to invert the diffusion age and the number of uplift events that formed a scarp; however, only when this number is reduced (< 10 uplift events), otherwise only the diffusion age can be robustly inferred. These inferences can only be achieved for a limited topographic noise, significantly lower than ca. 10% of the cumulated scarp height. More specifically, topographic noise needs to be significantly lower than the maximum
deviation expected between topographic profiles resulting from different slip scenarios – i.e., ≪ 20% of the total scarp height (Fig. 3b,c), otherwise the morphological signature of fault slip history remains hidden and undetectable.

The model developed in this paper succeeds in reconciling the apparent paradigm⁹ of young and old scarps associated with high and low apparent $κ$, respectively, by accounting for their contrasting slip history. We note that, due to the principle of superposition, a linear model of diffusion should lead to the least obvious effect of fault slip history on scarp morphology. This effect is, therefore, likely to be enhanced when considering more complex and potentially non-linear models. Since high-resolution (<1 m) topographic data are now becoming extensively available²⁵, our study will help revisit the age-morphology of scarps and, more interestingly, help better constrain the slip history of active faults. Our results imply that there could be more to be extracted from the morphology of landscapes than just the rates of active faulting, as extensively done in numerous studies. Indeed, the kinematics of faulting (seismic ruptures vs. creep) may leave a subtle signature yet to be fully explored.

## Data availability

Data used in this study come from previously published sources. Data are available through the article of Kokkalas and Koukouvelas⁹—https://doi.org/10.1016/j.jog.2005.07.006.

Received: 19 January 2023; Accepted: 28 February 2023

Published online: 08 March 2023

## References

1. Wallace, R. E. Profiles and ages of young fault scarps, north-central Nevada. *Geol. Soc. Am. Bull.* **88**, 1267–1281. https://doi.org/10.1130/0016-7606(1977)88%3c1267:paaoyf%3e2.0.co;2 (1977).

2. Chevalier, M. L., Pan, J., Li, H., Liu, D. & Wang, M. Quantification of both normal and right-lateral late quaternary activity along the Kongur Shan extensional system, Chinese Pamir. *Terra Nova* **27**(5), 379–391. https://doi.org/10.1111/ter.12170 (2015).

3. Arnaud, N. O., Brunel, M., Cantagrel, J. M. & Tapponnier, P. High cooling and denudation rates at Kongur Shan, eastern Pamir (Xinjiang, China) revealed by $^{40}$Ar/$^{39}$Ar alkali feldspar thermochronology. *Tectonics* **12**(6), 1335–1346. https://doi.org/10.1029/93TC00767 (1993).

4. Nash, D. B. Morphologic dating of degraded normal fault scarps. *J. Geol.* **88**, 353–360. https://doi.org/10.1086/628513 (1980).

5. Hanks, T. C., Bucknam, R. C., Lajoie, K. R. & Wallace, R. E. Modification of wave-cut and faulting-controlled landforms. *J. Geophys. Res. Solid Earth* **89**(B7), 5771–5790. https://doi.org/10.1029/JB089iB07p05771 (1984).

6. Andrews, D. J. & Hanks, T. C. Scarp degraded by linear diffusion: Inverse solution for age. *J. Geophys. Res. Solid Earth* **90**(B12), 10193–10208. https://doi.org/10.1029/JB090iB12p10193 (1985).

7. Avouac, J. P. & Peltzer, G. Active tectonics in southern Xinjiang, China: Analysis of terrace riser and normal fault scarp degradation along the Hotan-Qira fault system. *J. Geophys. Res. Solid Earth* **98**(B12), 21773–21807. https://doi.org/10.1029/93JB02172 (1993).

8. Carretier, S., Ritz, J. F., Jackson, J. & Bayasgalan, A. Morphological dating of cumulative reverse fault scarps: Examples from the Gurvan Bogd fault system, Mongolia. *Geophys. J. Int.* **148**(2), 256–277. https://doi.org/10.1046/j.1365-246X.2002.01007.x (2002).

9. Kokkalas, S. & Koukouvelas, I. K. Fault-scarp degradation modeling in central Greece: The Kaparelli and Eliki faults (Gulf of Corinth) as a case study. *J. Geodyn.* **40**(2–3), 200–215. https://doi.org/10.1016/j.jog.2005.07.006 (2005).

10. Peltzer, G. et al. Stable rate of slip along the Karakax section of the Altyn Tagh Fault from observation of interglacial and postglacial offset morphology and surface dating. *J. Geophys. Res. Solid Earth* **125**(5), e2019JB018893. https://doi.org/10.1029/2019JB018893 (2020).

11. Lavé, J. & Avouac, J. P. Active folding of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal. *J. Geophys. Res. Solid Earth* **105**(B3), 5735–5770. https://doi.org/10.1029/1999JB900292 (2000).

12. Garcia-Estève, C., Caniven, Y., Cattin, R., Dominguez, S. & Sylvain, R. Morphotectonic evolution of an alluvial fan: Results of a Joint analog and numerical modeling approach. *Geosciences* **11**(10), 412. https://doi.org/10.3390/geosciences1100412 (2021).

13. Lu, L., Zhou, Y., Zhang, P. & Cheng, X. Modelling fault scarp degradation to determine earthquake history on the Muztagh Ata and Tahman faults in the Chinese Pamir. *Front. Earth Sci.* **10**, 342. https://doi.org/10.3389/feart.2022.838866 (2022).

14. Hodge, M. et al. Evidence from high-resolution topography for multiple earthquakes on high slip-to-length fault scarps: The Bilila-Mtakataka Fault, Malawi. *Tectonics* **39**(2), e2019TC005933. https://doi.org/10.1029/2019TC005933 (2020).

15. Charreau, J. et al. Denudation outpaced by crustal thickening in the eastern Tianshan. *Earth Planet. Sci. Lett.* **479**(179–191), 2017. https://doi.org/10.1016/j.epsl.2017.09.025 (2017).

16. Culling, W. E. H. Analytical theory of erosion. *J. Geol.* **68**(3), 336–344. https://doi.org/10.1086/626663 (1960).

17. Mesimeri, M., Karakostas, V., Papadimitriou, E., Tsaklidis, G. & Jacobs, K. Relocation of recent seismicity and seismotectonic properties in the Gulf of Corinth (Greece). *Geophys. J. Int.* **212**(2), 1123–1142. https://doi.org/10.1093/gji/ggx450 (2018).

18. Armijo, R., Meyer, B. G. C. P., King, G. C. P., Rigo, A. & Papanastassiou, D. Quaternary evolution of the Corinth Rift and its implications for the Late Cenozoic evolution of the Aegean. *Geophys. J. Int.* **126**(1), 11–53. https://doi.org/10.1111/j.1365-246X.1996.tb05264.x (1996).

19. Briole, P. et al. Active deformation of the Corinth rift, Greece: Results from repeated global positioning system surveys between 1990 and 1995. *J. Geophys. Res.* **105**(B11), 25605–25625. https://doi.org/10.1029/2000JB900148 (2000).

20. Chousianitis, K., Ganas, A. & Evangelidis, C. P. Strain and rotation rate patterns of mainland Greece from continuous GPS data and comparison between seismic and geodetic moment release. *J. Geophys. Res. Solid Earth* **120**(5), 3909–3931. https://doi.org/10.1002/2014JB011762 (2015).

21. Pelletier, J. D. et al. Evolution of the Bonneville shoreline scarp in west-central Utah: Comparison of scarp-analysis methods and implications for the diffusion model of hillslope evolution. *Geomorphology* **74**(1–4), 257–270. https://doi.org/10.1016/j.geomorph.2005.08.008 (2006).

22. Roering, J. J., Kirchner, J. W. & Dietrich, W. E. Evidence for nonlinear, diffusive sediment transport on hillslopes and implications for landscape morphology. *Water Resour. Res.* **35**(3), 853–870. https://doi.org/10.1029/1998WR900090 (1999).

23. Tucker, G. E., Hobley, D. E., McCoy, S. W. & Struble, W. T. Modeling the shape and evolution of normal-fault facets. *J. Geophys. Res. Earth Surf.* **125**(3), e2019JF005305. https://doi.org/10.1029/2019JF005305 (2020).

24. Hanks, T. C. & Andrews, D. J. Effect of far-field slope on morphologic dating of scarp like landforms. *J. Geophys. Res. Solid Earth* **94**(B1), 565–573. https://doi.org/10.1029/JB094iB01p00565 (1989).

25. Kersten, T. P. & Lindstaedt, M. Geometric accuracy investigations of terrestrial laser scanner systems in the laboratory and in the field. *Appl. Geom.* **14**(2), 421–434. https://doi.org/10.1007/s12518-022-00442-2 (2022).
## Acknowledgements

This study was supported by the Agence Nationale de la Recherche (France): project TopoExtreme (PI: R.C.—Grant ANR-18-CE01-0017), project SLOWDEF (PI: M.S.—Grant ANR-18-CE31-0008), and IdEx Université de Paris (Grant ANR-18-IDEX-0001). This study has also received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (PI: P.S.—Grant 803721). We thank Christopher Scholz and the two anonymous reviewers for their constructive comments.

## Author contributions

R.C., M.S. and P.S. initiated this study. R.H. and R.C. developed the analytical formulations and conducted the numerical experiments. R.C., M.S. and P.S. analyzed the data associated with the Eliki fault scarp. All authors contributed to the writing of the paper and provided critical feedback.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-30772-z.

Correspondence and requests for materials should be addressed to R.C.

Reprints and permissions information is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023