Western Corinth Rift Research Articles

An Atypical Shallow Mw 5.3, 2021 Earthquake in the Western Corinth Rift (Greece)

AbstractModerate‐to‐large earthquakes in rifts may occur on leading boundary faults or inner antithetic faults. Here we show a rare case of the 2020–2021 seismic sequence in the Corinth rift, that culminated in the shallow rupture of the antithetic fault, neither preceded nor followed by the leading fault rupture. The hypocenter of the largest shock (Mw 5.3 of 17 February 2021) was located at ∼8 km depth. However, seismic waveform data, supported by satellite‐geodetic and tide gauge measurements, pointed to rupture at shallow depth (∼3 km), where no earthquakes were previously observed. We show that the earthquake most probably ruptured two orthogonal, conjugate fault segments: a weak nucleation phase occurred in the microseismically highly active sub‐horizontal detachment layer, followed – a few seconds later – by a larger, shallow moment release on a high‐angle, south‐dipping normal fault. The latter is the Mornos offshore fault, antithetic to the leading, north‐dipping Psathopyrgos fault. Our study presents the first instrumental/observational evidence of a very shallow Mw 5+ event in this rift – and one of the few reported worldwide. The depth limit of the main shallow slip patch coincides with the expected crossing of the Mornos fault with the Psathopyrgos fault, stressing the importance of fault segmentation and rooting inherited from the rift history. This unusual shallow slip in a depth range with little background seismicity and few aftershocks needs to be further investigated by dynamic modeling as a possible prototype of hazardous events in rift environments.

Dynamics of microseismicity and its relationship with the active structures in the western Corinth Rift (Greece)

We analyse the complete earthquake archive of the western Corinth Rift using both cross-correlations between pairs of event waveforms and accurate differential traveltimes observed at common stations, in order to identify small-scale fault structures at depth. The waveform database was generated by the dense Corinth Rift Laboratory network and includes about 205 000 events between 2000 and 2015. Half of them are accurately relocated using double-difference techniques. The novelty of this relocated catalogue is the integration of the recent westernmost earthquakes due to the extension of the network in 2010 to the western extremity of the Corinth Rift and the consideration of the whole database over more than 15 yr. The total relocated seismicity exhibits well-defined clusters at the root of the main normal faults mainly between 5 and 10 km depth in the middle of the gulf and illuminates thin active structure planes dipping north about 20° under the northern coast. Some seismicity is observed in the footwall of the main active faults, along the West and East Helike faults. We also built a multiplet database based on waveform similarity taking into account cross-correlation coefficients weighted by signal-to-noise ratios. Short-term multiplets are concentrated in the middle of the gulf along the Kamarai fault system, in a 1–2 km thick layer at 6–8 km depth, interpreted as a highly fractured geological layer. They are often associated to slow seismic migration velocities occurring in this zone during strong swarm episodes and are thus likely to be triggered by pore pressure variations. On the other hand, most long-term and regular multiplets are located deeper (7–10 km), under the northern coast, within a layer less than 0.3 km thick. They occur at the border of nearly planar structures with low seismicity rate, which we identify as fault planes, and they may be explained by aseismic slip on the fault surface around them. This supports the existence of an immature structure growing downdip towards the north at the base of the active geological layer, which possibly connects to the ductile middle crust around 15 km depth, as suggested by the occurrence of deeper events in the continuity of the 1995–fault plane. The different migration velocities (from 0.05 km d−1 to several km d−1) highlighted during the western 2014–swarms indicate that both pore pressure and creep diffusion are operating in the fault zone. The fast migrations observed in the Psathopyrgos fault zone, where a slow slip event was detected by dilatometers in 2002, compare with that for creeping faults. To the west, from spatial distribution of events, we show that the Rion–Patras fault connecting the western extremity of the Corinth Rift fault system to the Patras Rift, is dipping around 60° north–west with a rake angle of −115°. Finally, we identified two new areas within the central active zone which may correspond to large scale, locked asperities on active fault surfaces, similar in size to the main asperity broken during the 1995, MW 6.3, Aigion earthquake.

Methodology for earthquake rupture rate estimates of fault networks: example for the western Corinth rift, Greece

Abstract. Modeling the seismic potential of active faults is a fundamental step of probabilistic seismic hazard assessment (PSHA). An accurate estimation of the rate of earthquakes on the faults is necessary in order to obtain the probability of exceedance of a given ground motion. Most PSHA studies consider faults as independent structures and neglect the possibility of multiple faults or fault segments rupturing simultaneously (fault-to-fault, FtF, ruptures). The Uniform California Earthquake Rupture Forecast version 3 (UCERF-3) model takes into account this possibility by considering a system-level approach rather than an individual-fault-level approach using the geological, seismological and geodetical information to invert the earthquake rates. In many places of the world seismological and geodetical information along fault networks is often not well constrained. There is therefore a need to propose a methodology relying on geological information alone to compute earthquake rates of the faults in the network. In the proposed methodology, a simple distance criteria is used to define FtF ruptures and consider single faults or FtF ruptures as an aleatory uncertainty, similarly to UCERF-3. Rates of earthquakes on faults are then computed following two constraints: the magnitude frequency distribution (MFD) of earthquakes in the fault system as a whole must follow an a priori chosen shape and the rate of earthquakes on each fault is determined by the specific slip rate of each segment depending on the possible FtF ruptures. The modeled earthquake rates are then compared to the available independent data (geodetical, seismological and paleoseismological data) in order to weight different hypothesis explored in a logic tree.The methodology is tested on the western Corinth rift (WCR), Greece, where recent advancements have been made in the understanding of the geological slip rates of the complex network of normal faults which are accommodating the ∼ 15 mm yr−1 north–south extension. Modeling results show that geological, seismological and paleoseismological rates of earthquakes cannot be reconciled with only single-fault-rupture scenarios and require hypothesizing a large spectrum of possible FtF rupture sets. In order to fit the imposed regional Gutenberg–Richter (GR) MFD target, some of the slip along certain faults needs to be accommodated either with interseismic creep or as post-seismic processes. Furthermore, computed individual faults' MFDs differ depending on the position of each fault in the system and the possible FtF ruptures associated with the fault. Finally, a comparison of modeled earthquake rupture rates with those deduced from the regional and local earthquake catalog statistics and local paleoseismological data indicates a better fit with the FtF rupture set constructed with a distance criteria based on 5 km rather than 3 km, suggesting a high connectivity of faults in the WCR fault system.

Ambient noise tomography of the western Corinth Rift, Greece

Three years of continuous waveform data recorded at 22 stations from the Corinth Rift Laboratory and the Hellenic Unified Seismological Network are used to perform an ambient noise surface-wave tomography of the western Corinth Rift. All available vertical component time-series were cross-correlated to extract empirical Rayleigh-wave Green's functions. Group velocity dispersion curves were measured for each station-pair by applying frequency-time analysis and then inverted to build 2D group velocity maps between 1 and 6 s period. Finally, we locally inverted these velocity maps using a neighborhood algorithm to assess the 3D shear-velocity model of the shallow crustal structure of the western Corinth Rift. Across all studied periods the southern coast of the Corinth Gulf is generally imaged as a region of lower velocities compared to the northern coast. At periods up to 3 s, the spatial variation of the group velocities is correlated with the surface geology of the area. Lower velocities are observed in areas where mostly Plio-Quaternary syn-rift sediments are present, such as off-shore regions of the rift, the Mornos delta and the largest part of the southern coast. Higher velocities are observed in pre-rift basement structures which are dominated mostly by carbonates. At periods above 3 s, where Rayleigh-waves begin to sense deeper structures below the sediments within the underlying basement, our study highlights the presence of a distinct zone of lower velocities across the southern part of the rift with an elongation in the WNW-ESE direction. The interpretation of this low velocity includes two arguments, the present-day active tectonic regime and the possible involvement of fluids circulation processes at depth within a highly fractured upper crust in the vicinity of the major faults zones. In general, the results demonstrate good agreement with the major geological and tectonic features of the area, as well as with previous local earthquake tomography studies and support the assumption of fluid circulations at depth. This work intends to be the base for further investigations towards the study of the Corinth Rift structure using long-time series of ambient noise data.

Interaction of an antecedent fluvial system with early normal fault growth: Implications for syn‐rift stratigraphy, western Corinth rift (Greece)

AbstractContinental ‘overfilled’ conditions during rift initiation are conventionally explained as due to low creation of accommodation compared with sediment supply. Alternatively, sediment supply can be relatively high from the onset of rifting due to an antecedent drainage system. The alluvial Lower Group of the western Plio–Pleistocene Corinth rift is used to investigate the interaction of fluvial sedimentation with early rifting. This rift was obliquely superimposed on the Hellenide mountain belt from which it inherited a significant palaeorelief. Detailed sedimentary logging and mapping of the well‐exposed syn‐rift succession document the facies distributions, palaeocurrents and stratigraphic architecture. Magnetostratigraphy and biostratigraphy are used to date and correlate the alluvial succession across and between fault blocks. From 3·2 to 1·8 Ma, a transverse low sinuosity braided river system flowed north/north‐east to east across east–west‐striking active fault blocks (4 to 7 km in width). Deposits evolved downstream from coarse alluvial conglomerates to fine‐grained lacustrine deposits over 15 to 30 km. The length scale of facies belts is much greater than, and thus not directly controlled by, the width of the fault blocks. At its termination, the distributive river system built small, stacked deltas into a shallow lake margin. The presence of a major antecedent drainage system is supported by: (i) a single major sediment entry point; (ii) persistence of a main channel belt axis; (iii) downstream fining at the scale of the rift basin. The zones of maximum subsidence on individual faults are aligned with the persistent fluvial axis, suggesting that sediment supply influenced normal fault growth. Instead of low accommodation rate during the early rift phase, this study proposes that facies progradation can be controlled by continuous and high sediment supply from antecedent rivers.

A CROSS-CORRELATION TECHNIQUE FOR RELOCATION OF SEISMICITY IN THE WESTERN CORINTH RIFT

Local seismological networks provide data that allow the location of microearthquakes which otherwise would be dismissed due to low magnitudes and low signal-to-noise ratios of their seismic signals. The Corinth Rift Laboratory (CRL) network, installed in the western Corinth rift, has been providing digital waveform data since 2000. In this work, a semi-automatic picking technique has been applied which exploits the similarity between waveforms of events that have occurred in approximately the same area of an active fault. Similarity is measured by the crosscorrelation maxi-mum of full signals. Events with similar waveforms are grouped in multiplet clusters using the nearest-neighbour linkage algorithm. Manually located events act as masters, while automatically located events of each multiplet cluster act as slaves. By cross-correlating the P-wave or S-wave segments of a master event with the corresponding segments of each of its slave events, after appropriately aligning their offsets, the measured time-lag at the cross-correlation maximum can be subtracted from the arrival-time of the slave event. After the correction of the arrival-times, a double-difference technique is applied to the modified catalogue to further improve the locations of clusters and distinguish the active seismogenic structures in the tectonically complex Western Corinth rift.

Sedimentary impacts of recent moderate earthquakes from the shelves to the basin floor in the western Gulf of Corinth

In seismically active areas, long term records of large earthquakes are indispensable to constrain reccurence patterns of large earthquakes. In the western Corinth Rift, one of the most active areas in Europe in terms of seismicity, data about ancient earthquakes are still insufficient, despite historical records covering the last two millenia and several studies in onshore paleoseimology. In this paper, we test the use of offshore sediments from the Gulf of Corinth to identify sediment failures and tsunamis that have been triggered by historical earthquakes. Two shelves (40–100m deep), one sub-basin (180m) and the basin floor (330m) have been sampled by short gravity cores. The cores were analyzed in order to identify and characterize event deposits. The age control has been provided by 137Cs and 210Pb activity measurements showing that the cores represent 2 to 4 centuries of sedimentation. In each site, sandy event deposits are interbeded in the muddy, hemipelagic sedimentation. The age of event deposits has been compared to the record of historical earthquakes using new and published macroseismic data. This comparison shows temporal coincidence of some event deposits and documented earthquakes with a macroseismic intensity≥VII in the area, e.g. in 1861CE, 1888CE and 1909CE. In nearshore, shallow-water settings, the record of event deposits does not exactly fit with the historical record of large earthquakes because too few event deposits are present. This may be due to the absence of sediment failures or to a lower preservation of the deposits in such settings. In the deepest site, in the basin floor, the correspondence is better: a sandy turbidite probably corresponds to each large earthquake since 1850CE, except one aseismic sediment density flow that occurred at the end of the 20th century. Surprisingly, the Ms=6.2, June 15, 1995 Aigion earthquake is only possibly recorded in one nearshore site on the Aigion Shelf, in the form of a tsunami back-wash flow deposit. This study showed that moderate earthquakes (M 5.8–6.5) can significantly impact marine sediments. Regarding the evaluation of seismic hazard in the area, the basin floor is proposed as a promising site for long term paleoseismology in the Gulf of Corinth, while shallower settings need to be considered more carefully.

Rift migration and lateral propagation: evolution of normal faults and sediment-routing systems of the western Corinth rift (Greece)

Abstract The active Corinth rift records hanging-wall migration of faulting and slip-rate acceleration. The rift initiated at approximately 5–4 Ma, and older parts are well exposed in the northern Peloponnese. A new correlation of chrono- and lithostratigraphy and structure across the onland central to westernmost rift with offshore data reveals westward rift propagation, as well as northward fault migration. Northward fault migration ended first in the east, with the stabilization of major north-dipping faults that now bound the Gulf. The basin then propagated to the WNW in two stages, each involving the initiation of a new fault that propagated east to SE to link to the stable fault system. Extension rates accelerated in distinct steps as the rift opened to the west. The youngest faults in the westernmost rift are associated with high seismicity and highest geodetic extension due to rapid fault growth and linkage at depth. The early synrift succession infilled substantial inherited palaeo-relief. Antecedent rivers established vigorous sediment-routing systems that controlled facies distribution throughout rifting, albeit with drainage reorganization during fault-migration events. Multiple deepening events recorded in the stratigraphy can be due to lateral rift propagation. The transition from rift initiation to rift climax is, therefore, diachronous along the rift axis.

Bayesian inversion of seismic spectral ratio for source scaling: Application to a persistent multiplet in the western Corinth rift

AbstractWe propose a method to precisely estimate earthquake source parameters as magnitude, size of rupture, stress drop, and coseismic slip, and their uncertainties. This method, which relies on a Bayesian approach, allows the determination of the scalar seismic moment, corner frequency (fc), and their associated uncertainties, by inverting ratios between seismic displacement spectra of nearby located earthquakes. We apply this method to a large earthquake multiplet (56 events) located under the northern coast of the Corinth gulf at 8 km depth. This multiplet is regularly active between 2001 and 2007. Results show fcP/fcS ratios globally between 1.0 and 1.5 which is compatible with the values predicted by Madariaga's circular rupture model. In detail, six earthquakes, however, exhibit corner frequency variations as a function of the station azimuth compatible with linear rupture propagation. Magnitude ranges 1.08 and 2.80 with a b value of 1.04. Source rupture length globally ranges between 40 and 170 m for stress drop between 1 and 100 MPa. We show that the number of ruptures and the cumulated coseismic slip are maximal at the center of the multiplet: this suggests that Multiplet‐866 could be seen as a weak seismogenic patch surrounded by a locked fault. However, the large value of the maximum coseismic slip cumulated over the period 2000–2008 (10 cm) rather suggests creep allowing rapid stress reloading and repeated earthquakes with short delays. We therefore propose that Multiplet‐866 is surrounded by a heterogeneous fault surface with both locked and creeping areas.

The 2003–2004 seismic swarm in the western Corinth rift: Evidence for a multiscale pore pressure diffusion process along a permeable fault system

AbstractMicroseismic multiplets occurring in the western Corinth rift, Greece, during a large swarm are analyzed to retrieve their spatiotemporal characteristics. These multiplets activated small subfaults at depth (∼7 km), up to 1 km long, at the root of two parallel active normal faults. The swarm migrates westward nearly horizontally over 10 km at an average velocity of 50 m/d with a diffusivity of 0.5 m2 s−1. It successively activates the Aigion fault, a relay zone in its hanging wall, and the Fassouleika fault. Within each multiplet, hypocenters also migrate with diffusivities ranging from 0.001 to 0.4 m2 s−1. The largest internal diffusivities appear at the core of the layer defined by the clusters. These results are interpreted as a hydroshear process caused by pore pressure migration within permeable corridors resulting from the intersection of the major faults with a brittle geological layer inherited from the Hellenic nappe stack.

Seismotectonic analysis of the 2013 seismic sequence at the western Corinth Rift

This study focuses on a series of small intraplate earthquakes that took place during May–August 2013 on the southwestern coast of the Corinth Rift (Central Greece), a few km southeast of Aigion city. The Corinth Rift is one of the most seismically active parts of the Mediterranean. We analyzed more than 1500 events with 0.4≤ML≤3.7, the major part of which was recorded by a dense local network. The seismicity is densely clustered in a volume of dimensions ∼4×2×6km3, aligned in a N110° direction and at depths ranging between 6 and 12km. Precisely relocated hypocenters and reliably constrained focal mechanisms indicate north dipping planar faults with an average dip of ∼60°. Stress inversion of focal mechanisms implies that the dominant local stress field is extensional in a N5° direction, in good agreement with geodetic observations. The swarm evolved in two phases, with a spatiotemporal migration of epicenters from the eastern toward the western part of the rupture zone, while slip distribution appears homogeneous over the eastern part and strongly inhomogeneous in the western part. These two phases also produced different results in scaling relations such as the Gutenberg–Richter law, the Modified Omori Formula and the Epidemic Type of Aftershock Sequence model. Similar results from other studies have been reported and correlated with a fluid driven mechanism, however further research is required to strengthen this hypothesis for the purposes of this study.

Reassessment of the rifting process in the Western Corinth Rift from relocated seismicity

The seismic activity in the western part of the Corinth Rift (Greece) over the period 2000– 2007, monitored by a dense network of three-component stations, is analysed in terms of multiplets and high precision relocation using double difference techniques. This detailed analysis provides new insights into the geometry of faults at depth, the nature and the structure of the active zone at 6–8 km depth previously interpreted as a possible detachment, and more generally into the rifting process. The seismicity exhibits a complex structure, strongly varying along the rift axis. The detailed picture of the seismic zone below the rift indicates that its shallower part (at depths of 6–8 km) is 1–1.5 km thick with a complex microstructure, and that its deeper part (at depths of 9–12 km) gently dipping to the north (10–20°) is 0.1–0.3 km thick with a microstructure consistent with the general slope of the structure. Although the nature of this seismic zone remains an open question, the presence of seismicity beneath the main active area, the strong variability of the structure along the rift over short distances and the complex microstructure of the shallower part revealed by the multiplet analysis are arguments against the hypothesis of a mature detachment under the rift: this active zone more likely represents a layer of diffuse deformation. The geometry of the mapped active faults is not well defined at depth, as no seismicity is observed between 0 and 4 km, except for the Aigion Fault rooting in the seismic layer at 6 km depth with a dip of 60°. A distinct cloud of seismicity may be associated with the antithetic Kalithea Fault, on which the 1909 Fokis earthquake (M s = 6.3) may have occurred. The link between the 1995 rupture (M s = 6.2) and the faults known at the surface has been better constrained, as the relocated seismicity favours a rupture on an offshore, blind fault dipping at 30° , rather than on the deeper part of the East Helike Fault. Consequently, the 1995 event is expected to have decreased the Coulomb stress on the East Helike Fault. To explain these seismic observations along with the geodetic observations, a new mechanical model for the rifting process in this region is proposed, involving non-elastic, mostly aseismic uniform NS opening below the rift axis, coupled with the downward and northward growth of a yet immature detachment: the reported GPS rates would mainly result from this deep, silent source, and the seismicity would reveal the detachment position, not yet connected to the ductile lower crust. In such a model, the strong fluctuations of microseismicity would result from small strain instabilities, undetected by continuous GPS and possibly related to pore pressure transients.

Tectono‐sedimentary evolution of the western Corinth rift (Central Greece)

AbstractThe Corinth rift (Greece) is one of the world's most active rifts. The early Plio‐Pleistocene rift is preserved in the northern Peloponnese peninsula, south of the active Corinth rift. Although chronostratigraphic resolution is limited, new structural, stratigraphic and sedimentological data for an area >400 km2 record early rift evolution in three phases separated by distinct episodes of extension rate acceleration and northward fault migration associated with major erosion. Minimum total N–S extension is estimated at 6.4–7.7 km. The earliest asymmetrical, broad rift accommodated slow extension (0.6–1 mm a−1) over >3 Myrs and closed to the west. North‐dipping faults with throws of 1000–2200 m defined narrow blocks (4–7 km) with little footwall relief. A N‐NE flowing antecedent river system infilled significant inherited relief (Lower group). In the earliest Pleistocene, significant fluvial incision coincided with a 15 km northward rift margin migration. Extension rates increased to 2–2.5 mm a−1. The antecedent rivers then built giant Gilbert‐type fan deltas (Middle group) north into a deepening lacustrine/marine basin. N‐dipping, basin margin faults accommodated throws <1500 m. Delta architecture records initiation, growth and death of this fault system over ca. 800 ka. In the Middle Pleistocene, the rift margin again migrated 5 km north. Extension rate increased to 3.4–4.8 mm a−1. This transition may correspond to an unconformity in offshore lithostratigraphy. Middle group deltas were uplifted and incised as new hangingwall deltas built into the Gulf (Upper group). A final increase to present‐day extension rates (11–16 mm a−1) probably occurred in the Holocene. Fault and fault block dimensions did not change significantly with time suggesting control by crustal rheological layering. Extension rate acceleration may be due to strain softening or to regional tectonic factors.

Structure and evolution of the western Corinth Rift, through new field data from the Northern Peloponnesus

Abstract Extensional structures with geometrical and kinematic features analogous to the known Gulf of Corinth faults, are found further to the south of what is considered to be the southern margin of the of Proto-Corinth Gulf, reaching south to the northern flanks of Mt Mainalon. This mountain front is marked by the North Mainalon Fault Zone, which comprises a series of normal fault segments with NNE dips. Assuming a listric or ramp-flat geometry for the North Mainalon Fault Zone, it could flatten at a depth of 6–8 km, underneath Mt Khelmos. Its southern, shallow part has been truncated by NNE- and NNW-trending faults, which may be linked to northward propagation of the east–west extension in the Southern Peloponnesus, causing further uplift in the central and northern Peloponnesus, while its deeper part is still active and may reach further north and sole onto the hypothesized detachment zone beneath the Gulf of Corinth.

First results of the CRLN seismic network in the western Corinth Rift: evidence for old-fault reactivation

Abstract The 12 stations Corinth Rift Laboratory Seismological Network (CRLNET) aims at monitoring the seismicity (Ml>1) in the CRL area and at constraining the geometry of active structures at depth. Two years of microseismicity (2000–2001) recorded by the CRLNET in the Aigion area shows: (1) background seismicity inside the Corinth rift at depth of 4.5–11 km, deepening towards the north and no activity in the upper 4 km of the crust – this seismicity is not clearly related to major faults observed at the surface –; (2) a swarm, 6 km south of the city of Aigion, associated with the Mw=4.2, 8 April 2001 earthquake. This earthquake occurred at 6 km depth, on a SW–NE oriented fault dipping 40° to the northwest and corresponds to normal faulting with a right lateral component of slip. It likely occurred on an old structure reactivated in the present stress field. To cite this article: H. Lyon-Caen et al., C. R. Geoscience 336 (2004).