Τhe February – March 2019 Seismic Swarm Offshore North Lefkada Island, Greece: Microseismicity Analysis and Geodynamic Implications

A quite energetic seismic excitation consisting of one main and additional three distinctive earthquake clusters that occurred in the transition area between the Kefalonia Transform Fault Zone (KTFZ) and the continental collision between Adriatic and Aegean microplates, is thoroughly studied after high–precision aftershocks’ relocation. The activated fault segments are in an area where historical and instrumental data have never claimed the occurrence of a catastrophic (M>6.0) earthquake. The relocated seismicity initially defines an activated structure extending from the northern segment of the Lefkada branch of KTFZ with the same NNE–SSW orientation and dextral strike slip faulting and then keeping the same sense of motion its strike becomes NE–SW and its dip direction NW. This provides unprecedented information on the link between the KTFZ and the Collision front and sheds more light on the regional geodynamics. The earthquake catalog, which is especially compiled for this study, starts one year before the occurrence of the Mw5.4 mainshock, and adequately provides the proper data source for investigating the temporal variation of the b– value, which might be used for discriminating foreshock and aftershock behavior.


Introduction
A moderate magnitude Mw5.4, earthquake occurred on the 5 th of February 2019 in the offshore area north of Lefkada Island, strongly felt in the Lefkada city and the onshore continental area to the east, with no major damage or injuries reported. In the area an adequate number of moderate earthquakes have occurred in the last decades with no one having magnitude larger than M6.0. The seismic excitation started with an Mw4.2 shock that occurred on January 15 th , 2019, to the northeast of the main shock which was followed by a very productive aftershock sequence with more than 250 located shocks by 15 March 2019.
The activated area is located at the boundary between the Kefalonia Transform Fault Zone (KTFZ) to the south and the Adriatic-Eurasian Collision to the north. The continental collision is expressed by a belt of thrust faulting with NE-SW direction of the axis of maximum compression. It runs along the eastern coastline of Adriatic Sea and terminates just north of Lefkada Island. The KTFZ is a major dextral strike slip fault zone that frequently accommodates strong earthquakes, clustered in space and time probably due to the stress transfer between the fault segments comprised in the fault system [1]. The northernmost fault segment of KTFZ was activated in 2003 with an Mw6.2 main shock and its rich aftershock sequence. The accurate location of these aftershocks provided for the first time the indication that the major fault segments bound the western coastlines of the Island, and the contemporaneously activated secondary fault segments lie onshore [2]. No earthquakes with M>6.0 have occurred during the instrumental era along the collision front in the vicinity of the 2 of 18 activated area. Smaller magnitude earthquakes are also sparse in comparison with the remarkably active KTFZ.
Geophysical and geodetic studies have provided a sufficient knowledge of the seismotectonic properties of the area. Identifying in more detail the mechanism of transition along active boundaries from the one plate to its neighbor, remains a major scientific challenge. Tracing the kinematic and dynamic variability along the active boundary can provide important clues for the transition area. Microseismicity is mainly concentrated along the KTFZ and at the same time manifests the swarm type seismogenesis, as a preferential style of strain release (Figure 1). This latter property becomes the tool to decipher the localized strain pattern, since strong earthquakes occur rarely and not all in optimally oriented faults, a fact that results in appreciable gaps in seismicity catalogs both in time and space. The investigation of this activity is challenging to shed light on the seismogenic setting and the complex geodynamics of the area, prerequisite for any seismic hazard assessment study. Its location is critical to decipher how these secondary structures play a role in accommodating strain at the borderline of the KTFZ, an area where also reverse focal mechanisms were determined for moderate magnitude earthquakes that occurred in the last few decades. The adequacy of proper data provided by an earthquake catalog compiled for the targets of this study allowed the monitoring of temporal variations of the b-value, which from both laboratory experiments [3,4] and real data [5] seems to be linked with the stress state in a region in a way that higher differential stress corresponds to lower values. We then extend our investigations to the temporal variations in values aiming to detect whether they can act as an indicator for the main shock occurrence or not.

Data acquisition and earthquake relocation
Since the occurrence of the 2003 Lefkada Mw6.3 main shock, continuous seismic monitoring was intensified in the area of central Ionian Islands, namely Lefkada and Kefalonia, by firstly installing temporary broadband and short period digital seismological stations, then by upgrading this network into a permanent one [6] (http://geophysics.geo.auth.gr/ss/) and finally incorporating the stations in the Hellenic Unified Seismological Network (HUSN). For the purpose of this study we used some additional stations located in the mainland, at distances not exceeding 90 km, in order to reduce azimuthal gap in the events location.
In order to identify the activated structure of the northern end of the Kefalonia Transform Fault Zone (KTFZ) we compiled a catalog of 662 earthquakes, for the period February 2018 -October 2019. The starting time of the catalog was decided to be a year before the mainshock in order to study any pre-seismic patterns because it was observed that the seismic activity was considerable. The locations were made by manual picking of 3560 P-and 3197 S-phases, using 14 broadband and short period stations of the HUSN. Phase pickings for the earthquakes with M≥2.5 were performed manually by the analysts of the Geophysics Department of the Aristotle University of Thessaloniki (GD -AUTh). For the smaller magnitude earthquakes phase picking was held manually by the authors of the present work, as an ongoing procedure of deeper understanding the seismotectonic and geodynamic properties of the central Ionian islands.
Relocation was performed with the HYPOINVERSE [7] program using the 1D model (Table 1) along with the vP/vS ratio for this region suggested by [8], for the relocation of the 2015 Lefkada aftershock sequence. Station delays were calculated for further refining the locations, following a procedure described in [9], to account for lateral crustal variations that are not considered in a 1D model. The medians of the hypocentral errors derived by HYPOINVERSE are of the order of 2900m for the horizontal and 1500m for the vertical dimension. method (LSQR) [11] with appropriate damping. For the cross correlation differential times the time domain cross correlation technique was realized [12,13]. Waveforms with a duration of 60 seconds, starting at the origin time of each earthquake, were prepared and band-pass filtered at 2-10 Hz with a second order Butterworth filter. Cross correlation was performed for a 1 -second window length starting at the arrival of P and S phases, searching over a lag of 1 second. Differential times with a correlation coefficient of 0.8 or more (CC≥0.8) were kept being used in hypoDD.
In the application of the DD algorithm five sets of iterations with five iterations each were performed. The cross correlation differential times were down weighted by a factor of 100 for the first 10 iterations, to obtain locations from the catalog data [8,14]. The catalog data was downweighted by the same factor for the remaining 15 iterations to allow cross correlation differential times to better define the active structures. For the first five iterations, and considering the catalog data alone, weighting was not applied. The relocation procedure took place for interevent distances less than or equal to 5 km for the first 10 iterations of catalog data and cross correlation times, reaching at 3km for the cross-correlation dataset at the final set of iterations.

Location errors estimation
The location errors calculated by hypoDD are underestimated, usually being a few tens of meters. Thus, a bootstrap method was applied to the final residuals, creating a dataset by drawing samples from the observed residual distribution and adding them to the differential times and then repeating the relocation [15,16]. This procedure was repeated 200 times and then the 95% error ellipsoid per event was computed by the resulting catalogs.
The jackknife method was also applied to evaluate the station distribution effect on the final catalog. The relocation was repeated omitting one station at a time, followed by calculation of the standard deviations per event. The medians calculated for the axis of the ellipsoids produced by the bootstrap method and the medians of the standard deviations calculated from the jackknife method are summarized in Table 2.
The errors produced by both methods show that in the West -East horizontal direction they are much larger than in the North -South direction. This reveals the effect of the station distribution in the location accuracy, as the closer stations are located south of the activated area in Lefkada and Kefalonia islands, in a rather narrow N-S oriented zone covering a small azimuthal range in the E-W direction. In addition, the number of stations is small to the east and north of the study area, in western Greece, and obviously there are no stations to the west, an area covered entirely by sea.

Fault plane solutions determination
A centroid moment tensor for the mainshock of the sequence was published soon by GCMT (https://www.globalcmt.org/CMTsearch.html), suggesting an almost pure strike slip motion with the two nodal planes steeply dipping and striking either N-S or E-W. None of the nodal planes agrees with the epicentral distribution in the area surrounding the main shock. In the framework of this study, the first polarities were exploited with the use of the FPFIT program [17]. The recordings of the HUSN stations for earthquakes with magnitudes ML≥3.0 were used for polarities identification. Since most of the stations are located south of the activated area, providing a poor azimuthal coverage, only impulsive onsets were used for the focal mechanism determination. Eight focal mechanisms were calculated, information on which is provided in Table 3.
The one of the nodal planes determined for the main shock of Mw5.4, listed fourth in Table 3, striking at 227 o , dipping to the northwest at an angle of 45 o , shows an excellent agreement with the epicentral distribution of the relocated seismicity, as this will be discussed in the next section. The significant dextral strike slip component (rake=160 o ) agrees well with the regional stress field pattern. Most of the focal mechanisms determined for the smaller magnitude earthquakes that belong to the same cluster along with the main shock (numbered 1, 2, 3 and 5 in Table 3), are similar to that of the mainshock. The fault plane solutions of the stronger earthquakes that occurred south and west of the main cluster, and impressively all of them following the stronger events occurrence of the main cluster, exhibit differentiation of the first nodal plane strike (10 o and 30 o for the events numbered 6 and 7 in Table 3 that belong to the second cluster, and 18 o for the 8 th earthquake belonging to the southernmost cluster). They also differ from the gross faulting properties of the Lefkada segment of the KTFZ [8,18]. Table 3. Information on the fault plane solutions determined in the present study. The first column gives the serial number of each solution ordered according to the occurrence date and time, shown in the second and third columns, respectively. The relocated epicentral coordinates and focal depths follow in the next three, 5th to 7th and the magnitude in the 8th column. The strike, dip and rake of the 2 nodal planes are shown in columns 8-10 and 11-13, respectively.

Temporal distribution of earthquakes
The seismic excitation attracted our attention from the first days, given that it started on the 13 th of January with an M3.6 earthquake, followed by numerous smaller ones, then an M4.5 shock on the 15 th of January, followed by an M3.7 after an hour and a half (see Table 3 for occurrence dates and times), all of them felt by the inhabitants of the surrounding urban areas. The activity continued for months, exhibiting alteration of active and relatively quiet periods. In order to investigate and analytically describe the temporal behavior of this rather complex seismic activity, the cumulative number and magnitude distribution of the earthquakes as a function of time are investigated and discussed. Figure 2a shows the cumulative number of earthquakes that occurred since February 2018. The total period is divided into four subintervals separated by three turning points where the color in the symbols is changed. The first period spans from February 2018 until 15 January 2019, where a magnitude 4.5 earthquake occurred close to the incoming mainshock. The activity of the first period (green symbols in Figure 2a impressive is the relative quiescence for about 20 days, which is more obvious in Figure 2b, when on February 05 the mainshock with Mw5.4 occurred. One hundred aftershocks followed in two days and the activity presented an exponential like decrease, up to shortly before the occurrence of the M4.0 earthquake on 26 th of February. It seems that its occurrence followed a rejuvenation that could be considered as a foreshock phase, and then its imminent aftershocks created a new short seismicity intensification. Soon after, the activity declined as a function of time, as expected for an aftershock sequence.  Figure 3, where four colors are used to illustrate the epicenters, corresponding to the four different periods as defined in Figure 2. The seismicity is sparse and evenly distributed in the study area during 2018 (green color epicenters in Figure 3). The M4.5 earthquake triggers significant microseismicity close to its epicenter and in an area of several kilometers around (blue colored epicenters in Figure 3), much longer than its source dimensions. The Mw5.4 strongest earthquake in the seismic excitation, triggered aftershocks in a rather large seismic zone having NE-SW orientation (red epicenters in Figure 3) for a short time interval during which the activity was considerably increased (red dots in Figure 2b). In one day, the activity expanded in the entire study area, manifesting a remarkable cluster comparatively far from the main shock, to the southwest of its aftershock area. The earthquake magnitudes in this cluster are all smaller than 4.0.  In less than twenty (20) days of regular aftershocks rate decay, the activity was appraised again and on 26 th of February an M4.0 earthquake followed (yellow star in Figure 2 and Figure 3). This indicates a migration and expansion of the activity to the southern part of the study area. It is of paramount importance then to explore further the sequential occurrence of the stronger members of this activity and the possible triggering through transfer of stress among the activated fault segments, which will be given in a next section of this study.

Spatial distribution
The relocated seismicity allowed refined and more reliable picture than the disperse epicentral distribution provided from the regional catalog as shown in Figure 1, revealing four distinctive clusters along with the disperse activity (Figure 4a). The northernmost cluster is associated with the mainshock epicenter, the location of which is shown with a red star, occupying an area of almost 7km in length, slightly larger than expected for an MW5.4 earthquake, according to well accepted scaling laws [19,20]. The epicentral elongation at an NE-S W direction agrees with one of the nodal planes of the focal mechanisms of the main shock and the four aftershocks that belong to this cluster ( Figure  4a and Table 3).
The activated fault segment with a strike of 220 o dips to NW at an angle of 45 o and is associated with dextral strike slip faulting (rake=160 o ) with a considerable thrust component. This agrees with the dextral strike slip motion of the KTFZ to the south of the study area, and the continental collision to the north. The almost E-W azimuth of the axis of maximum compression also agrees with the regional stress pattern. For a more detailed investigation of the geometry of this fault segment, the relocated epicenters of the specific cluster alone are plotted in Figure 5a strike parallel and strike normal that are shown in Figures 5b and 5c, respectively. It is worth to note here that the focal depths are restricted in depths between 14 and 20 km, with the vast majority between 16 and 18 km, an important result for the identification of the seismogenic layer in the study area. The dip of the fault plane as derived from the fault plane solution is again in perfect agreement with geometry of the activated fault segment as this is revealed by the strike normal cross section in Figure 5c. The next seismicity concentration is located to the south of the previous cluster and encompasses several earthquakes the spatial distribution of which show a barely preferable E-W elongation, in agreement with the strike of the one nodal plane from the 05/02/2019 earthquake of M3.3 (Figure 4a). The depth range is between 10 and 15 km, shallower than and exactly above the activated layer in the main cluster and with a steeper dip angle (Figures 5b and 5c). The faulting type, in addition to the rotation of its strike more to the north as mentioned before, exhibits dextral strike slip motion with a slight normal component. The epicentral distribution delineates a distinctly N-S striking structure, as one proceeds more to the south, where two sub -clusters might be recognized, forming an almost N-S striking activated structure. The northern one, hosting the M4.0 strong event (yellow circle in Figure 4a) and forming a narrow strip of relocated epicenters, and one to its south where the activity appears sparser. The focal depths are now even more shallow, between 5 and 11 km, as the depth range of the relocated seismicity further south to Lefkada [8,21]. Almost N-S striking, dipping to the east dextral strike slip faulting holds for the two focal mechanisms determined for the M4.0 earthquake and one of its foreshocks, again in full agreement with the faulting properties of the Lefkada fault branch of the KTFZ [2].

Coulomb stress changes and possible triggering
The spatial aftershock distribution is known to be affected by the static stress changes caused by the mainshock [22]. The change in Coulomb failure function (ΔCFF) is used to quantify closeness to failure [23]. A positive value of ΔCFF on a fault indicates increased likelihood that this fault will where μ is the coefficient of friction. Δp is expressed by the relation [25]: if undrained conditions are considered during the coseismic phase. B is the Skempton's coefficient (0≤B<1), with experimental values between 0.5 and 0.9 shown to be appropriate for rocks [25], and Δσkk is the trace of the induced stress tensor. Assuming a ductile fault zone enclosed by a homogeneous and isotropic medium, equation (1) becomes: where μ' is the apparent coefficient of friction [23], which contains the effects of pore fluid as well as the material properties of the fault zone. It is estimated by the equation: For this study, the values used to compute ΔCFF are μ=0.75 and Β=0.5 as in [26] and as adopted in other studies in Greece [27,28]. These values result in an apparent coefficient of friction of μ'≈0.4, suggested by [1] for the area of central Ionian Islands after testing different values and found the most appropriate one for this estimation. We calculated the stress changes due to the coseismic slip of the main shock with the causative fault being approximated as a planar area with a length of 7 km and a width of 3km, equal to the width of the seismogenic layer as explained in the previous section. Α seismic moment of Μο=1.68•10 17 Ν•m was adopted from the GCMT solution and with the above source dimensions and a shear modulus of 33 GPa, a mean slip of 0.1745 m was estimated. The Poisson ratio was fixed equal to 0.25.
The Coulomb stress changes were calculated for the faulting type of the main shock (strike=227 o , dip=45 o , rake=160 o ) at a depth of 16.5 km, which is approximately the middle of the seismogenic layer ( Figure 5c) and are shown in Figure 6a along with the relocated aftershocks. It is evident that the vast majority of the off-fault activity is located into stress enhanced areas. Given that the faulting type of the earthquakes encompassing in the other clusters differs than the one of the main shock, and that the activity was more shallow, the stress field is inverted according to the faulting type of the M4.0 shock (strike=18 o , dip=62 o , rake=139 o ) at a depth of 11.5 km (Figure 6b). The earthquakes that occurred south of the main cluster are located inside areas of positive stress changes, an evidence for possible triggering, although with smaller values than in Figure 6a, possibly due to the difference in the calculation depth.

Aftershock sequence analysis in relation with temporal variations of b value
The evolution of an aftershock sequence is among the key factors for operational forecasting, given that the areas already hit by the main shock are continuously under the threat of a strong aftershock or even a higher magnitude earthquake. Monitoring the seismicity is the means to evaluate the imminent seismic hazard and several powerful tools have been developed to this respect. The occurrence of a strong aftershock, or of a larger magnitude earthquake that is afterwards considering as the main shock, along with the occurrence frequency of aftershocks, depends upon the stress residuals and stress redistribution before and after the target earthquake. The direct manifestation of the stress state is the seismicity and one approach to detect fluctuations in its temporal behavior might be through the detection of the temporal variations of the b value. The disadvantages of this procedure are frequently connected with the lack of adequate data before the main shock that could prevent us from a robust estimation of the value, and the overlapping of the coda-waves in the first hours after the main shock that can limit the early aftershocks detection capability [29]. However, the duration of our relocated catalog starting one year before the main shock and careful inspection and manual phase picking for aftershocks, augmented the data sample and diminished the possible uncertainties in the b value estimations. The approach that we apply is described in the next steps and follows the procedure of [30].
A catalog encompassing 662 earthquakes that occurred from the 4 th of February 2018 up to the 10 th of October 2019 consists our data sample for investigating the temporal behavior of b value before and after the main shock that occurred on the 5 th of February 2019, one year after the incipience of the catalog. Firstly, the completeness magnitude is detected through the maximum-curvature (MAXC) method [31], which is insensitive to sample sizes and shows great stability in estimating [32]. The magnitude bin with the highest frequency of events corresponds to the threshold value, which in our case is found to be = 1.3 (Figure 7). The value of the GR law is computed through the Maximum Likelihood Estimation (MLE) method [33]: where ̅ is the mean value of the magnitudes and is the magnitude bin, usually chosen to be equal to 0.1. This relation evidences the dependence of value on the choice of . For the assessment of the temporal variations of the value before and after the main shock, firstly we divided the catalog into two sub-catalogs, before and after the main shock with and sample sizes, respectively. Then, for each dataset we recompute the magnitude of completeness with the same method, since it is more appropriate than the Goodness of Fit method for small sample sizes ( * ≤ 200), which is the case here. A correction factor of +0.2 is considered in order to avoid biased results in the value estimation, since the tends to be underestimated [34,32]. Table 4 summarizes the details of the two datasets. In both cases, the magnitude bin with the highest frequency corresponds to = 1.3, so the completeness magnitude after the addition of the correction factor +0.2 is = 1.5 leading to = 94 events before the main shock and = 226 aftershocks (see Figure S1 for the GR law distribution of both datasets). Then, a sample window of size and is determined for the datasets before and after the main shock, respectively, for which we iteratively estimated the values moving at each iteration one event forward. A minimum number of events for the sample size, ℎ = 50, is defined, which is the threshold for an approximately 15% error in the estimation of value [35]. In this way, we achieved statistically reliable results. We tested different sample sizes in order to examine the dependency of the results on the choice of the free parameters, and (Figure 8). The time-varying evolution of value does not significantly depend on the choice of these parameters in the size ranges (50, 85). We set = 50 and = 70, with 45 and 157 number of samples, respectively. However, we stress that with any other combination the same conclusions would have been reached. In our case, the magnitude completeness is assumed constant, however, it is generally suggested to reevaluate it at each step since it can change in time, due to seismic network modifications or detection limitations. We keep the same initially estimated completeness threshold for two reasons. Firstly, the dataset covers a short period with relatively few events since it concerns the aftershock sequence of a moderate size main shock, which must be taken into account with caution for a robust estimation of the value while evaluating iteratively the completeness threshold (many samples fall under the minimum value ℎ ). Secondly, we verified that the geometry of the seismic network during the study period was stable, meanwhile we are highly confident about the monitoring analysis quality since it was under our surveillance for the entire time interval. Nevertheless, we tested the temporal fluctuations of the values both with constant and moving magnitude completeness threshold and we did not observe significant variations through time ( Figure S2 and S3). In all cases, the main observation of decreasing b values after the 4.5 foreshock (dashed black line in Figure  8 ) compared to the increasing ones just after the main shock is holding true, therefore, in our analysis we continued with a constant threshold that let us include more samples.
Finally, we computed the median of the values for the period of almost one year (from February 2018 until the last event before the foreshock of Mw4.5 as a reference value for the comparison with the period between the foreshock and the main shock that occurred 21 days later, as well as with the entire aftershock period. In Figure 9 the temporal variation of the values with = 0.64 is illustrated. We can see that after the foreshock of Mw4.5 on 15 January 2019 the values decrease rapidly compared to the reference value by approximately 10%. In contrast, a general increase in the values after the main shock occurrence is observed, that exceeds the 10% a day after, continues to increase and remains at a level reaching the 20% or more up to the end of the study period.

Discussion and conclusions
The seismic activity that manifested in 2019 in the offshore area north of the Lefkada Island motivated our investigation for the identification of the geometry and kinematics of the activated fault segments, given that they are located at the junction between the two major active boundaries, the dextral strike slip Kefalonia Transform Fault Zone and the collision between the Adriatic microplate and the Greek mainland. It was expected to unveil the mode of kinematics and fill in the gap for the transition from transform to collision, adding value in the exhaustive relocation and manually picking of arrival phases for as many earthquakes as possible of the recorded ones, and the first polarities of the stronger shocks that are associated with the activated fault segments.
The exact location and geometry of the offshore part of the KTFZ and its termination against the thrust front, was not identified with seismological data up to now. This lack of knowledge about its northern extent is mostly due to the lack of moderate earthquakes at this location during the last years where the seismological network detectability became efficient to provide the data for this analysis. The geometry of the offshore faults identified in this study exhibit mostly strike slip motion do not match the orientation of the Lefkada fault branch. The difference in strike is ~20 o , with the incorporation of a small thrust component in agreement with the adjacency of the compressional stress field regime. The three observed clusters that are traced from the relocated seismicity in this study exhibit spatial properties differentiation and reveal the 3-D local structure, given that the focal depths were relocated with high accuracy. The northern fault segment although traced in already suggested kinematic models [36,37], it is the first time to be identified and investigated with seismological recordings, and in particular, with such fine details. [36] referred to a right lateral strike-slip system that resulted from the collision of western Greece with the Apulian platform, and runs through the western Greece to the KTFZ. The position of the identified activated fault segments, their strike and sense of slip agrees with the active fault mapping and suggested regional tectonics of [38]. A clockwise rotation 2 B discussed of the Ionian Islands and Akarnania block accommodated by major marginal strike-slip zones that appeared segmented along their strike is suggested by [39].
Available fault plane solutions of moderate and strong earthquakes (M≥5.0) in the Lefkada segment of KTFZ and the collision front extends as far in the past as 1973, from various sources ( [40] for the earthquake of 1973/11/04, [41] for the earthquake of 25/02/1994, and the GCMT catalogue). Nevertheless, there is not any available focal mechanism solutions for earthquakes of such magnitude in the study area. Fault plane solutions for smaller magnitude earthquakes in the offshore area north of and close to Lefkada Island can be found in several studies ( [8,42,43] among others). Those solutions are either covering the offshore area that is the southern edge of our study area, thus consistently following the NNE -SSW strike of the main fault zone, or they are not accompanied by a highly accurate focal distribution as in our study, because even in the most recent cases of the 2003 and 2015 earthquakes the seismic activity in the area of the main fault under study was not intense. Figure 10 shows the identified active structures and a kinematic sketch of the study area. Τhe focal mechanisms determined in the present study (compression quadrants shown in red, and numbering according to Table 3) along with the focal mechanisms of earthquakes with M≥5.0 already published (green compression quadrants). Their origin time (YYΥΥ/MM/DD) is shown on the top of the beach balls. The thick gray lines denote the major active boundaries, namely the KTFZ where the dextral strike slip motion is shown by the antiparallel yellow arrows and the collision front by the sawtooth line. The prolongation of the KTFZ in the study area, as identified and documented in the present study, is colored in red. The inset map shows the axes of maximum stress, P-axes, as they have been taken from the fault plane solutions shown in the main part of the figure. The thick brown curved arrows follow the P-axis rotation in agreement with the strike slip motion along the KTFZ and the counterclockwise rotation of Adria microplate. The local fault system sheds more light to the transition from the dextral strike slip motion along the KTFZ to the thrust faulting along the collision front. The thrust component present in the fault plane solutions determined for the stronger earthquakes of the 2019 swarm, put in evidence the influence of the compressive stress field to the strike slip faults that are either striking along the KTFZ (the southern ones) or "kinking" (the northern one) for compensating the kinematic direction prevalent at this location. The slight but gradually changing rotation in the azimuth of the maximum stress axes, from south (green arrows in Figure 10) to the northernmost edge of the 2019 activated structures (red arrows) support the hypothesis of a transpression movement in a transition zone between the transform and collision zones.
Preprints (www.preprints.org) | NOT PEER-REVIEWED | Posted: 25 May 2020 doi:10.20944/preprints202005.0412.v1 Figure 10. Fault plane solutions of 16 earthquakes projected as equal area lower hemisphere projections, where compression quadrants are in green for earthquakes of M>5.0, that occurred before the earthquake sequence in study and in red for earthquakes that were determined and belong to the swarm investigated in this study. Thick grey lines denote the active boundaries of the KTFZ and the Apulian collision front. Red lines depict the identified fault segment offshore Lefkada Thin black line represents the Ionian thrust. Yellow arrows indicate relative plate motion. The Inset shows the Paxis inferred from the focal mechanisms. Black arrows illustrate the characteristic direction of P -axis, NE -SW for the earthquakes along the KTFZ and N -E for those Examination of the b value changes show lower values before and higher after the mainshock occurrence. A considerable decrease of the b-value (from 0.6 to 0.5) was observed immediately after the occurrence of the M4.5 foreshock on January 15, 2019 and remained low until the mainshock occurrence on February 05, 2019. Then, sharply increased to the values before the M4.5 foreshock and gradually became higher approaching the value of 0.9.
Supplementary Materials: The following are available online at www.mdpi.com/xxx/s1, Figure S1: Frequency -Magnitude Distribution (FMD of pre-mainshock and post-mainshock datasets, Figure S2: Time-series of the bvalue of different sample sizes for the pre -mainshock dataset, Figure S3: Same as in Figure S2 for the postmainshock dataset.