The Hellenic Subduction System (HSS) in the eastern Mediterranean is the oldest active subduction margin on earth. It is a segmented boundary that hosts the continuum of faulting styles over a ~200km range in depth and can generate large earthquakes with high tsunamigenic potential.  The complexity of deformation styles and rates leave key aspects of the system poorly understood. For example, historical records of Mw<8 earthquakes fail to explain the current observed convergence rate (~35mm/year), and recent geodetic measurements suggest that the degree of locking within the system is heterogeneous. The density of geodetic measurements is increasing rapidly, nevertheless, the inherent time lag required to accumulate data that will enable identifying regions that undergo slower (than seismic) deformation transients will necessitate inferences from seismic signals. In this work, we aim to further close the observational gap between heterogeneous deformation styles and rates using the features of seismicity distributions to infer where deformation rates, and by inference, locking, vary most.    To that scope, we will present new results of an enhanced earthquake catalog that we will use to explore the spatio-temporal distribution of seismicity features (e.g., b-value, effective stress drop, seismic-moment-release skewness) to infer variability in deformation rates and loading. Catalog enhancement exploits data from the temporary (EGELADOS) broadband seismometer network that operated between 2005 until 2007 combined with permanent stations leading to a station spacing of ~40 km and covering the entire southern Aegean Sea. We first use the combined network to detect earthquakes using machine learning approaches (EQTransformer, PhaseLink) for detection, phase picking and association. After performing initial locations using NonLinLoc combined with a 1D velocity model and quality control procedure, we enhance the number of small-magnitude detections using a multi-station template-matching approach. Next, we scan the enhanced high-resolution catalog for distinct spatial and temporal patterns of seismicity using unsupervised clustering. We then quantify the clustered seismicity using b-value, effective stress drop, and seismic-moment-release skewness (among other parameters). We will present our clustering results in the context of the variability in slip phenomena related to earthquake-earthquake interactions (e.g., static and dynamic triggering) as well as in the context of external forcing (e.g., aseismic triggering or fluid migration).   The preliminary results that we will present will provide a basis for our more broad-scale study of interplay between seismic and aseismic deformation. In particular, where the latter is gradually becoming increasingly resolvable using GNSS data within the HSS, this work will provide a basis for links with geodetically observed deformation in the future.  
Abstract Small stress changes such as those from sea level fluctuations can be large enough to trigger earthquakes. If small and large earthquakes initiate similarly, high‐resolution catalogs with low detection thresholds are best suited to illuminate such processes. Below the Sea of Marmara section of the North Anatolian Fault, a segment of 150 km is late in its seismic cycle. We generated high‐resolution seismicity catalogs for a hydrothermal region in the eastern Sea of Marmara employing AI‐based and template matching techniques to investigate the link between sea level fluctuations and seismicity over 6 months. All high resolution catalogs show that local seismicity rates are larger during time periods shortly after local minima of sea level, when it is already rising. Local strainmeters indicate that seismicity is promoted when the ratio of differential to areal strain is the largest. The strain changes from sea level variations, on the order of 30–300 nstrain, are sufficient to promote seismicity.
Abstract We investigated the Coulomb stress changes in the active faults surrounding a moderate‐magnitude normal‐faulting earthquake (2009 L'Aquila, Mw 6.3) and the associated variations in the expected ground motion on regional probabilistic seismic hazard maps. We show that the static stress variations can locally increase the seismic hazard by modifying the expected mean recurrence time on neighbouring faults by up to ~290 years, with associated variations in the probability of occurrence of the maximum expected earthquake of up to ~2%. Our findings suggest that the increase in seismic hazard on neighbouring faults following moderate‐magnitude earthquakes is probably not sufficient to necessitate systematic upgrades of regional probabilistic seismic hazard maps, but must be considered to better address and schedule strategies for local‐scale mitigation of seismic risk.
Earthquake stress drop values estimated from ground-motion spectra commonly vary by several orders of magnitude, particularly for small earthquakes (~M < 3). Stress-drop values have been found to vary with faulting style, faulting type (intraplate, interplate), depth, and to exhibit differences between natural and induced earthquakes. Nevertheless, distinguishing uncertainties from real trends across data sets is challenging, in part due to the variation in methodological approaches and observational constraints. However, the proliferation of high-quality, dense seismic data in recent years has shown that at least some of the variability in stress drop values almost certainly reflects diversity in fault strength and geological conditions. Coupling well-constrained observations to a variety of modeling approaches will help uncover what controls earthquake rupture processes, but deconvolving observational constraints from real variation in rupture behavior is key.We present our stress drop estimates from data sets representing a wide range of fault loading conditions and geological environments, from interplate, intraslab and forearc subduction faults, to volcanic, intraplate, and human induced events. Stress-drop values range primarily between 1 – 100 MPa for events that meet the criteria for spectral-ratio analysis.  We present correlations of low relative stress drop values in areas of high seismic attenuation indicative of lower rock strength, and a slight correlation with depth that corresponds to modeled deviatoric stress values. We also show one notable subset of induced events near active injection wells that exhibit stress drop values of ~0.1 MPa and have distinctive low-frequency content. Their spatial distribution, waveform, and source spectral characteristics suggest either slower rupture, lower stress drop values, or a combination of both, and may represent part of the transition between aseismic and seismic slip. We show using a Large-n array that while stress drop values are roughly constant (within 2 orders of magnitude), estimates can vary by roughly 25% when station coverage is limited to 15 stations or less with a maximum azimuthal gap of 90°.  Our findings highlight the importance of using modeling approaches to explore relative influence of fault strength and methodological approaches in stress drop variation. In particular, models that incorporate both frictional and thermoelastic approaches may provide clues to the variability of conditions that can activate faults, both within stable sliding and seismic rupture conditions.        
<p>The month-to-year-long deformation of the Earth&#8217;s crust where active subduction zones terminate is poorly explored. Here we report on a multidisciplinary dataset that captures the synergy of slow-slip events, earthquake swarms and fault-interactions during the ~5 years leading up to the 2018 M<sub>w</sub> 6.9 Zakynthos Earthquake at the western termination of the Hellenic Subduction System (HSS). It appears that this long-lasting preparatory phase initiated due to a slow-slip event that lasted ~4 months and released strain equivalent to a ~M<sub>w</sub> 6.3 earthquake. We propose that the slow-slip event, which is the first to be reported in the HSS, tectonically destabilised the upper 20-40 km of the crust, producing alternating phases of seismic and aseismic deformation, including intense microseismicity (M<4) on neighbouring faults, earthquake swarms in the epicentral area of the M<sub>w</sub> 6.9 earthquake ~1.5 years before the main event, another episode of slow-slip immediately preceding the mainshock and, eventually, the large (M<sub>w </sub>6.9) Zakynthos Earthquake. Tectonic instability in the area is evidenced by a prolonged (~4 years) period of overall suppressed b-values (<1) and strong earthquake interactions on discrete strike-slip, thrust and normal faults. We propose that composite faulting patterns accompanied by alternating (seismic/aseismic) deformation styles may characterise multi-fault subduction-termination zones and may operate over a range of timescales (from individual earthquakes to millions of years).</p>