Anomalous passive subsidence of deep‐water sedimentary basins: a prearc basin example,southern New Caledonia Trough and Taranaki Basin,New Zealand

Stratigraphic data from petroleum wells and seismic reflection analysis reveal two distinct episodes of subsidence in the southern New Caledonia Trough and deep‐water Taranaki Basin. Tectonic subsidence of ~2.5 km was related to Cretaceous rift faulting and post‐rift thermal subsidence, and ~1.5 km of anomalous passive tectonic subsidence occurred during Cenozoic time. Pure‐shear stretching by factors of up to 2 is estimated for the first phase of subsidence from the exponential decay of post‐rift subsidence. The second subsidence event occured ~40 Ma after rifting ceased, and was not associated with faulting in the upper crust. Eocene subsidence patterns indicate northward tilting of the basin, followed by rapid regional subsidence during the Oligocene and Early Miocene. The resulting basin is 300–500 km wide and over 2000 km long, includes part of Taranaki Basin, and is not easily explained by any classic model of lithosphere deformation or cooling. The spatial scale of the basin, paucity of Cenozoic crustal faulting, and magnitudes of subsidence suggest a regional process that acted from below, probably originating within the upper mantle. This process was likely associated with inception of nearby Australia‐Pacific plate convergence, which ultimately formed the Tonga‐Kermadec subduction zone. Our study demonstrates that shallow‐water environments persisted for longer and their associated sedimentary sequences are hence thicker than would be predicted by any rift basin model that produces such large values of subsidence and an equivalent water depth. We suggest that convective processes within the upper mantle can influence the sedimentary facies distribution and thermal architecture of deep‐water basins, and that not all deep‐water basins are simply the evolved products of the same processes that produce shallow‐water sedimentary basins. This may be particularly true during the inception of subduction zones, and we suggest the term ‘prearc’ basin to describe this tectonic setting.     

The diffuse transition between the Zagros continental collision and the Makran oceanic subduction (Iran): microearthquake seismicity and crustal structure

The nature of the transition between the Zagros intra-continental collision and the Makran oceanic subduction is a matter of debate: either a major fault cutting the whole lithosphere or a more progressive transition associated with a shallow gently dipping fault restricted to the crust. Microearthquake seismicity located around the transition between the transition zone is restricted to the west of the Jaz-Murian depression and the Jiroft fault. No shallow micro-earthquakes seem to be related to the NNW–SSE trending Zendan–Minab–Palami active fault system. Most of the shallow seismicity is related either to the Zagros mountain belt, located in the west, or to the NS trending Sabzevaran–Jiroft fault system, located in the north. The depth of microearthquakes increases northeastwards to an unusually deep value (for the Zagros) of 40 km. Two dominant types of focal mechanisms are observed in this region: low-angle thrust faulting, mostly restricted to the lower crust, and strike-slip at shallow depths, both consistent with NS shortening. The 3-D inversion of P traveltimes suggests a high-velocity body dipping northeastwards to a depth of 25 km. This high-velocity body, probably related to the lower crust, is associated with the deepest earthquakes showing reverse faulting. We propose that the transition between the Zagros collision and the Makran subduction is not a sharp lithospheric-scale transform fault associated with the Zendan–Minab–Palami fault system. Instead it is a progressive transition located in the lower crust. The oblique collision results in partial partitioning between strike-slip and shortening components within the shallow brittle crust because of the weakness of the pre-existing Zendan–Minab–Palami faults.     

Flexural uplift of the Stara Planina range, central Bulgaria

The Stara Planina is an E–W-trending range within the Balkan belt in central Bulgaria. This topographically high mountain range was the site of Mesozoic through early Cenozoic thrusting and convergence, and its high topography is generally thought to have resulted from crustal shortening associated with those events. However, uplift of this belt appears to be much younger than the age of thrusting and correlates instead with the age of Pliocene–Quaternary normal faulting along the southern side of the range. Flexural modelling indicates the morphology of the range is consistent with flexural uplift of footwall rocks during Pliocene–Quaternary displacement on S-dipping normal faults bounding the south side of the mountains, provided that the effective elastic plate thickness of 12  km under the Moesian platform is reduced to about 3  km under the Stara Planina. This small value of elastic plate thickness under the Stara Planina is similar to values observed in the Basin and Range Province of the western United States, and suggests that weakening of the lithosphere is due to heating of the lithosphere during extension, perhaps to the point that large-scale flow of material is possible within the lower crust. Because weakening is observed to affect the Moesian lithosphere for ≈10  km beyond (north of) the surface expression of extension, this study suggests that processes within the uppermost mantle, such as convection, play an active role in the extension process. The results of this study also suggest that much of the topographic relief in thrust belts where convergence is accompanied by coeval extension in the upper plate (or 'back arc'), such as in the Apennines, may be a flexural response to unloading during normal faulting, rather than a direct response to crustal shortening in the thrust belt.     

Shear-wave anisotropy and the stress field from borehole recordings at 2.5 km depth at Cajon Pass

53 local earthquakes recorded at 2.5 km depth in the Cajon Pass scientific borehole are analysed for shear-wave splitting. The time delays between the split shear waves can be positively identified for 32 of the events. Modelling these observations of polarizations and time delays using genetic algorithms suggests that the anisotropic structure near Cajon Pass has orthorhombic symmetry. The polarization of the shear waves and the inferred strike of the stress-aligned fluid-filled intergranular microcracks and pores suggests that the maximum horizontal compressional stress direction is approximately N13°W. This is consistent with previous results from earthquake source mechanisms and the right-lateral strike-slip motion on the nearby San Andreas Fault, but not with stresses measured within the uppermost 3 km of the borehole. This study suggests that the San Andreas Fault is driven by deeper tectonic stresses and the present understanding of a weak and frictionless San Andreas Fault may need to be modified. The active secondary faulting and folding close to the fault are probably driven by the relatively shallow stress as measured in the 3.5 km deep borehole.     

Evidence of recent surface faulting and surface rupture in the Fore-Alps of Veneto and Trentino (NE Italy)

The Fore-Alps of western Veneto and Trentino regions belong to the central Southern Alps (NE Italy), in which there is little evidence of very fresh surface ruptures or surface faulting. This does not seem to match historical data about earthquakes, some of which have been very intense. The strong influence of the inherited structures makes it difficult to detect a direct link between morphotectonic features and present-day stress fields.In the present study, four areas (Orsara, Scandole, Naole and Soran) with surface faulting and surface rupture features were examined, and models of morphotectonic evolution are discussed. In the Lessini Mountains, the Orsara graben and Scandole ridge show examples of surface faulting and surface rupture, respectively, reactivating Paleogene normal faults and fractures. Within the Orsara graben, rocky bluffs displace the previous morphological features. The bluffs are some decimetres to some metres high and are practically devoid of evidence of either physical or chemical weathering; on the slopes above them are steep areas which may be interpreted as the remnants of previous strongly weathered bluffs. The Scandole ridge has many trenches, some with rocky walls, which may be the result of several episodes of morphotectonic rupture.In the Giudicarie Belt, the Naole and Mt. Soran surface faulting landforms are details within large frontal culmination walls of Neogene thrusts. The Naole ridge corresponds to the southeastern sector of Monte Baldo. Here, inside a fault angle valley, a sinuous scarp originating from surface faulting marks the base of the fault scarp slope. Ridge splitting is the expression of the backward migration of separation niches due to slope tectonics, also evidenced on the slope by several terrace-like features and by a lower belt of very thick slope breccias. On Mt. Soran, in the Gruppo di Brenta massif, the surface faulting scarp faces uphill, giving rise to a trench-like feature. Downvalley of the scarp, there is the niche of a large landslide dated to 3 kyears B.P.All these landforms are consistent with slope tectonic movements caused by intense earthquakes. Whereas the morphostructures in the Lessini Mountains are the result of responses by sensitive structures, the Naole and Mt. Soran features express the evolution of frontal culmination walls of thrusts, with clear evidence of present-day tectonic activity.On the basis of the weathering of the scarps and associated features, the relative seismotectonic episodes probably occurred between the Bronze and Middle Ages.     

Shear-wave polarizations in the Peter the First Range indicating crack-induced anisotropy in a thrust-fault regime

Summary. Three-component seismograms of small local earthquakes recorded in the Peter the First Range of mountains near Garm, Tadzhikistan SSR, display shear-wave splitting similar to that previously observed near the North Anatolian Fault in Turkey. The Peter the First Range is in a region of compressional tectonics, whereas the North Anatolian Fault is a comparatively simple strike-slip fault. Detailed analysis of the Turkish records suggests that the splitting is diagnostic of crack-induced anisotropy caused by vertical microcracks aligned parallel to the direction of maximum compression. Preliminary examination of paper records from Garm shows that most shear waves arriving within the shear-wave window display shear-wave splitting, and that the polarizations of leading shear-waves are consistently aligned in a NE/SW direction. The area is complicated and the tectonics are not well-understood, but the NE/SW direction is approximately perpendicular to the compressional axis in many of the fault-plane mechanisms of the earthquakes. These earthquakes are usually at depths between 5 and 12 km, although there are some deeper events nearby.
Parallel shear-wave polarizations, such as those observed, are expected to indicate the strike of nearly vertical parallel microcracks, which would be aligned parallel to the direction of maximum compression. Thus the shear-wave polarizations in the Peter the First Range indicate that the directions of principal stress are reversed in the rock above the earthquake foci where thrust faulting is taking place.     

Geodetically determined strain across the southern end of the Tonga-Kermadec-Hikurangi subduction zone

Summary. A method of simultaneous reduction is presented for determining strain rates from multiple triangulation surveys where common triangulation stations have been used, but the angles of the old survey have not necessarily been reobserved. This method is applied to triangulation in the northern South Island, at the southern end of the Tonga-Kermadec-Hikurangi subduction zone. From a profile of shear strain across the Indian-Pacific plate boundary, the displacement of the Indian plate relative to the Pacific is calculated to be 54 ± 9 mm yr⁻¹ at an azimuth of 84°± 10°, in remarkable agreement with the motion predicted by global plate tectonic models. Most of this motion occurs within a 150 km wide zone bounded on the east by the Hikurangi Trough. Within this zone the motion is partitioned: near the Hikurangi Trough no slip is occurring at the upper surface of the subducting Pacific plate (the subduction thrust) and motion is predominantly thrusting normal to the trough axis: to the west is a region of predominantly dextral strike slip faulting. This pattern is consistent with Fitch's model of oblique subduction. To the south of the profile, a change is observed in the azimuth of the faulting along a line which marks the southern extent of the subduction slab, indicating the end of the partitioned motion.     

Basement faulting and the focal depths of the larger earthquakes in the Zagros mountains (Iran)

Summary. The present day seismicity of the Zagros seems to occur on high angle reverse faults distributed across the whole width of the belt. It does not indicate activity on a single inclined thrust surface and there do not seem to have been any well located events at intermediate depths. Modelling of the long period teleseismic body waves of seven large earthquakes presented here shows their focal depths to be in the range 8–15 km. This is thought to indicate faulting in the uppermost basement beneath the sedimentary cover, though the absence of published seismic refraction results renders the sediment thickness uncertain. Faulting of this type and depth may occur on inherited normal faults which have subsequently been reactivated as thrusts. Such reactivation allows considerable shortening to take place in the early stages of continental collision without the subduction or excessive thickening of continental crust.     

A microearthquake survey at the junction of the East Pacific Rise and the Wilkes (9°S) fracture zone

Summary. An ocean bottom seismograph survey of the junction of the East Pacific Rise and the Wilkes fracture zone detected only three microearthquakes beneath the rise crest during seven days of recording. In contrast, during the same period 41 events were detected on the fracture zone, all at distances greater than 10 km from the junction. These results suggest that near the rise crest the thin crust can support sufficient stress only to generate infrequent small earthquakes and that most faulting may take place by aseismic slip. At 10 km from the rise axis part of the crest has become competent enough to support stress, resulting in earthquakes probably at depths of up to 5 km below the sea-bed. Gear 5-waves on the seismometer records indicate that a magma chamber, if it exists near the junction, is less than 10 km across.     

Source parameters of large historical (1918–1962) earthquakes, South Island, New Zealand

We present the results of body waveform modelling studies for 17 earthquakes of M _w ≥5.7 occurring in the South Island, New Zealand region between 1918 and 1962, including the 1929 M _s = 7.8 Buller earthquake, the largest earthquake to have occurred in the South Island this century. These studies confirm the concept of slip partitioning in the northern South Island between strike-slip faulting in southwestern Marlborough and reverse and strike-slip faulting in the Buller region, but indicate that the zone of reverse faulting is quite localized. In the central South Island, all historical earthquakes appear to be associated with strike-slip faulting, although recent (post-1991) reverse faulting events suggest that slip partitioning also occurs within this region. The difference between historical and recent seismicity in the central South Island may also reflect stress readjustment occurring in response to the 1717 ad rupture along the Alpine fault. Within the Fiordland region (southwestern South Island) none of the historical earthquakes appears to have occurred along the Australian/Pacific plate interface, but rather they are associated with complex deformation of the subducting plate as well as with deformation of the upper (Pacific) plate. Two earthquakes in the Puysegur Bank region south of the South Island suggest that strike-slip deformation east of the Puysegur Trench is playing a major role in the tectonics of the region.     

The deep seismicity of the Tyrrhenian Sea

The study reappraises the deep seismicity of the Tyrrhenian Sea. Careful examination of the quality of reported hypocentres shows that the earthquakes define a zone dipping NW, about 200 km along strike, 50 km thick, and reaching a depth of about 500 km. The zone is slightly concave to the NW at a depth of 300 km, but, contrary to many previous reports, is not tightly concave, nor are there significant spatial gaps in the seismicity, which is effectively continuous with depth. Seismicity is, however, concentrated in the depth interval 250–300 km, where the dip of the seismic zone changes from 70° (above 250 km) to a more gentle dip of 45° at greater depths. Seven fault-plane solutions are available for the largest earthquakes in this depth interval, all of them consistent with a P -axis down the dip of the seismic zone, and all of them requiring movement on faults out of the plane of the subducting slab.
Two deep earthquakes near Naples lie well outside the main zone of activity; for one of which a fault-plane solution is available that has a P -axis not aligned with the dip of the seismic zone. The tightly concave slab-geometry favoured by other reports is supported mainly by the location of these events near Naples, which we think may represent deformation in a separate, probably shallower dipping, piece of subducted lithosphere.
The lack of shallow seismicity, and particularly of thrust faulting earthquakes, at the surface projection of the Benioff zone suggests that active subduction has ceased. Estimates of the convergence rate responsible for subduction in the last 10 Myr far exceed the present convergence rate of Africa and Eurasia, suggesting that the subduction was related instead to the stretching and thinning of the crust in the Tyrrhenian Sea.     

The evolution of the Gulf of Corinth (Greece): an aftershock study of the 1981 earthquakes

Summary. A preliminary study of the aftershocks of three earthquakes that occurred near to Corinth (Greece) in 1981 is combined with observations of the morphology and faulting to understand the evolution of the Eastern Gulf of Corinth. The well located aftershocks form a zone 60km long and 20km wide. They do not lie on the main fault planes and are mostly located between the north-dipping faulting on which the first two earthquakes occurred and the south-dipping faulting associated with the third event. A cluster of aftershocks also lies in the footwall of the eastern end of the south-dipping fault of the third event.
Morphologically, it is observed that in the evolution of the Eastern Gulf of Corinth, antithetic faulting apparently predates the appearance of the main faulting at the surface. This evolution can be explained by motion on a deep seated, shallow angle, aseismic normal fault. A model based on such a fault also accounts for the aftershock distribution of the 1981 earthquakes.     

Modelling the effect of atmospheric pressure variations on gravity

Summary. The 1973 Hawaii earthquake occurred north of Hilo, at a depth of 40 to 50km. The location was beneath the east flank of Mauna Kea, a volcano dormant historically, but active within the last 4000 yr. Aftershocks were restricted to a depth of 55–35km. The event and its aftershock sequence are located in an area not normally associated with the seismicity of the Mauna Loa and Kilauea calderas. The earthquake was a double event, the epicentres trending NE-SW. The events were of similar size and faulting mechanism. The fault plane solutions obtained by seismic waveform analysis are a strike-slip fault striking EW and dipping 55° S, the auxiliary plane a NS vertical plane with a faulting plunge of 35°. The axis of maximum compressive stress is aligned with the direction of the gravity gradient associated with the island of Hawaii. The fault plane striking EW parallels a surface feature, the Mauna Kea east rift zone. The earthquakes were clearly not associated with volcanic activity normally associated with Mauna Loa and Kilauea and may indicate a deep seated prelude to a resumption of activity at Mauna Kea.     

Normal faulting, extension and uplift in the outer thrust belt of the central Apennines (Italy): role of the Caramanico fault

The outer Adriatic zones of the central Apennines (Italy) provide good conditions for analysing geometry and kinematics of the earliest normal faults, superposed onto the thrust belt. During the latest stages of thrusting onto the Adriatic foreland (late Pliocene–early Pleistocene), the outermost imbricates of the thrust belt were subjected to normal faulting, coeval with differential uplift. Crosscutting normal faults get younger towards the foreland, thus the easternmost normal faults record the latest stages of fault propagation and growth. The Caramanico fault, on the western flank of Mt. Maiella, is the largest outcropping normal fault of the outer zones. This high‐angle fault (dip > 70°) has cumulative offsets ≤ €4.2 km, and propagated with slip rates of 2.6 mm/year in a short time interval (≤ 1.6 Ma), concomitant with intense uplift of Mt. Maiella. In contrast with normal faults in a more internal position, the Caramanico fault maintains a high‐angle planar geometry, and does not reach the major basal detachment of the thrust belt. Thus the fault did not cause large extensional displacements; its major role was rather to accommodate ongoing components of vertical uplift of the overthickened thrust wedge. Downfaulting of the thrust belt on the western flank of Mt. Maiella represents the youngest end member of the same processes that have operated since 11 Ma in the Tyrrhenian hinterland, where large extensional strains and crustal thinning of the orogenic belt were achieved by long‐lasting activity of listric normal faults detached at lower crustal depths.     

Earthquake processes of the Himalayan collision zone in eastern Nepal and the southern Tibetan Plateau

Focal mechanisms determined from moment tensor inversion and first motion polarities of the Himalayan Nepal Tibet Seismic Experiment (HIMNT) coupled with previously published solutions show the Himalayan continental collision zone near eastern Nepal is deforming by a variety of styles of deformation. These styles include strike-slip, thrust and normal faulting in the upper and lower crust, but mostly strike-slip faulting near or below the crust–mantle boundary (Moho). One normal faulting earthquake from this experiment accommodates east–west extension beneath the Main Himalayan Thrust of the Lesser Himalaya while three upper crustal normal events on the southern Tibetan Plateau are consistent with east–west extension of the Tibetan crust. Strike-slip earthquakes near the Himalayan Moho at depths >60 km also absorb this continental collision. Shallow plunging P -axes and shallow plunging EW trending T -axes, proxies for the predominant strain orientations, show active shearing at focal depths ∼60–90 km beneath the High Himalaya and southern Tibetan Plateau. Beneath the southern Tibetan Plateau the plunge of the P -axes shift from vertical in the upper crust to mostly horizontal near the crust–mantle boundary, indicating that body forces may play larger role at shallower depths than at deeper depths where plate boundary forces may dominate.     

Thrust faulting and crust—upper mantle structure in East Australia

Summary. The mid-crustal earthquake of 1973 March 9 (m_b= 5.5, h ≤ 20 km) located 60 km south-west of Sydney, Australia, provides unambiguous evidence of contemporary thrust faulting in South-eastern Australia — a region of high heat flow and Cenozoic basaltic volcanism. Aftershock locations suggest a steeply dipping fault in the depth range from 8 to 24 km with a lateral extent of about 8 km. The mechanism solution is consistent with a tectonic stress field that is dominated by east—west horizontal compression. A seismic moment of 5.7 ± 10²³± 20 per cent dyne-cm was computed from surface-wave amplitudes. Minimum values of slip and stress drop, 2 cm and 1 bar respectively, were estimated from the moment and a fault size taken' from aftershock locations.
Refinement modelling by a controlled Monte Carlo technique was used to provide unbiased models directly from multimode group velocities. The dispersion of fundamental and higher mode surface waves recorded at the digital high-gain station at Charters Towers, Queensland, and the WWSSN station at Adelaide, South Australia, is satisfied by crust- and upper-mantle models which have neither pronounced S-wave low-velocity zones nor thick high-velocity lids within 140 km of the Earth's surface. These models have subcrustal shear velocities of 4.20–4.32 km/s which are 0.4–0.5 km/s slower than Canadian shield shear velocities (CANSD).     

A lower mantle S-wave triplication and the shear velocity structure of D"

Summary. A lower mantle S-wave triplication detected with short- and long-period WWSSN and CSN recordings indicates a substantial shear velocity discontinuity near 280 km above the core–mantle boundary. The triplication can be observed in rotated SH seismograms from intermediate and deep focus events throughout the distance range from 70° to 95°. Three distinct source region–receiver array combinations that have been investigated in detail demonstrate consistent travel time and relative amplitude behaviour of the triplication, with slight systematic shifts in the triplication indicating up to 40 km variations in the depth of the discontinuity. Modelling of the observations with synthetic seismograms produced with the Cagniard de Hoop and reflectivity methods constrains the shear velocity increase to be 235 ± 0.25 per cent, comparable to upper mantle discontinuities. Short-period observations indicate that the velocity increase may be a sharp first-order discontinuity, or may extend over a transition zone no more than 50 km thick. The shear velocity gradient below the discontinuity, within the D" layer, is not well-constrained by the SH data, but slightly positive or near zero velocity gradients are consistent with the long-period amplitude ratios of ScSH/SH .     

Rift-phase extensional fabrics of the back-arc Drummond Basin, eastern Australia

The base of the Late Devonian–Early Carboniferous Drummond Basin, a major backarc extensional feature in eastern Australia which formed in response to detachment faulting, is extensively exposed in central Queensland. Here a crystalline basin floor is overlain by the Silver Hills Volcanics, a synrift sequence of predominantly silicic ash flow tuffs and lavas ranging to over 2 km in thickness. Detailed mapping of faults and stratigraphic logging of thickness changes within the Silver Hills Volcanics have allowed the rift-phase structural architecture that accompanied initial subsidence near the basin margin to be resolved. A complex mosaic of block faults with throws of up to 1 km is indicated. Locally developed mosaics may conform to, or depart from, the configuration predicted by the detachment faulting model. Structural fabric of the basement was a critical determinant of the extensional geometry. Distributed shear along pre-existing penetrative planar fabrics is considered to have accommodated hangingwall extension at lower strain rates whereas the propagation of tension fractures and the development of block faults by failure on pre-existing, brittle, basement dislocations facilitated extension at higher strain rates. The detachment fault inferred to lie beneath the extended hangingwall carapace has not been mapped at the surface and is thought to dissipate into a broad zone of distributed shear within basement to the east of the basin. Volcanism coincided with the initiation of extensional movements at which time deep crustal repositories for evolved magma were tapped by extensional fractures. The main extensional faults cutting the basinal succession were not used as conduits for magmatic products which were sourced from the basin margin and from extended hinterland to the east.     

Global mapping of upper mantle reflectors from long-period SS precursors

Long-period precursors to SS resulting from underside reflections off upper mantle discontinuities ( SdS where d is the discontinuity depth) can be used to map the global distribution and depth of these reflectors. We analyse 5,884 long-period seismograms from the Global Digital Seismograph Network (1976-1987, shallow sources, transverse component) in order to identify SdS arrivals. Corrections for velocity dispersion, topography and crustal thickness at the SS bounce point, and lateral variation in mantle velocity are critical for obtaining accurate estimates of discontinuity depths. The 410 and 660 km discontinuities are observed at average depths of 413 and 653 km, and exhibit large-scale coherent patterns of topography with depth variations up to 40 km. These patterns are roughly correlated with recent tomographic models, with fast anomalies in the transition zone associated with highs in the 410 km discontinuity and lows in the 660 km discontinuity, a result consistent with laboratory measurements of Clapeyron slopes for the appropriate phase changes. The best resolved feature in these maps is a trough in the 660 km discontinuity in the northwest Pacific, which appears to be associated with the subduction zones in this region. Amplitude variations in SdS arrivals are not correlated with discontinuity depths and probably result from focusing and defocusing effects along the ray paths. The SdS arrivals suggest the presence of regional reflectors in the upper mantle above 400 km. However, only the strongest of these features are above probable noise levels due to sampling inadequacies.     

Viscoelastic relaxation and long-lasting after-slip following the 1997 Umbria-Marche (Central Italy) earthquakes

We combine Global Positioning System (GPS) measurements with forward modelling of viscoelastic relaxation and after-slip to study the post-seismic deformation of the 1997 Umbria-Marche (Central Apennines) moderate shallow earthquake sequence. Campaign GPS measurements spanning the time period 1999–2003 are depicting a clear post-seismic deformation signal. Our results favour a normal faulting rupture model where most of the slip is located in the lower part of the seismogenic upper crust, consistent with the rupture models obtained from the inversion of strong motion data. The preferred rheological model, obtained from viscoelastic relaxation modelling, consists of an elastic upper crust, underlain by a transition zone with a viscosity of 10¹⁸ Pa s, while the rheology of deeper layers is not relevant for the observed time-span. Shallow fault creep and after-slip at the base of the seismogenic upper crust are the first order processes behind the observed post-seismic deformation. The deep after-slip, below the fault zone at about 8 km depth, acting as a basal shear through localized time-dependent deformation, identifies a rheological discontinuity decoupling the seismogenic upper crust from the low-viscosity transition zone.