1 Geophysical Journal International - Supplementary material for

2

3 The GPS velocity field of the Aegean. New observations, contribution of the earthquakes,

4 crustal blocks model.

5

6 Briole, P., Ganas, A., Elias, P & Dimitrov, D.

7

8 SUMMARY

9

10 The analysis of the secular component of the velocity field of the Aegean presented, discussed and

11 modelled in the main text, requires the accurate determination of the transient part of the velocity

12 field. This transient part is dominated by the coseismic displacements produced by the major

13 earthquakes that occurred in the area during the analysed time window. Another significant

14 component is due to the postseismic relaxations associated with those earthquakes. In this

15 supplementary material we review the coseismic and postseismic displacements induced by the

16 crustal earthquakes of magnitude Mw ≥ 5.3 during the period 2000-2020. In addition, several GPS

17 stations have their time series disrupted by other sources of transients with different origins, that we

18 also present and discuss. Once the transient velocity field is estimated it can be removed from the

19 total velocity field to extract what can be considered as the secular velocity field. We discuss the

20 gradients of this secular velocity field as they are measured along three sections originated at the

21 Euler pole of rotation Anatolia-Eurasia.

22

23

24 S1. INTRODUCTION

25

1 26 During the period 2000-2020 several strong earthquakes occurred in the Aegean region, with a

27 concentration in western Greece around and near the Ionian Islands: Lefkada 2003, Movri 2008,

28 Cephalonia doublet 2014, Lefkada 2015, and Zakynthos 2018. The strongest events of the period

29 are the Methoni earthquake of February 14, 2008 and the Samothraki earthquake of May 24, 2014.

30 Those large events generated coseismic displacements observable in a large number of GPS

31 stations, up to distances larger than 200 km. The data from the Bulgarian and Turkish stations are

32 essential for the correct analysis of the Samothraki 2014 and Kos 2017 earthquakes.

33 Tab. 1 reports twenty-one earthquakes (plotted in Fig. 1 and Fig. S1) with minimum magnitude

34 5.3, having induced a movement observed at least at one of the analysed GPS stations. A previous

35 inventory of displacements induced by earthquakes can be found in Ganas et al. (2018a). In two

36 cases (Zakynthos 2006 and Gulpinar 2017) we are not considering a single earthquake but a swarm.

37 The formulas of Okada (1992) show that for a shallow earthquake of magnitude Mw ≥ 5.3, a station

38 located in a radius of 5 km from the epicentre can in most cases detect a deformation signal large

39 enough to be separated from the noise. The largest population of analysed events in 2019 is not only

40 due to the occurrence of several moderate earthquakes that year but also to the higher density of the

41 array in the recent years. The period we studied includes also the volcanic unrest of Santorini in

42 2011-2012 (Foumelis et al. 2013; Parks et al. 2015). The deformation produced by the analysed

43 earthquakes comprises also in most cases a postseismic relaxation which can represent a significant

44 percentage of the coseismic deformation. Fig. S2 shows the example of the coseismic and

45 postseismic signals of the north component of the LEMN station located near the epicentre of the

46 Samothraki 2014 earthquake. In the following we model the postseismic signals by assuming that

47 they obey the negative exponential law of a Maxwell material (see Thatcher & Rundle 1984).

48 Fig. S3 shows the locations of the sites mentioned in this supplementary material. Tab. S2 gives

49 details on the toponymy of the sites used in the main text and in this supplementary material.

50

2 51

52 S2. GPS DATA PROCESSING

53

54 First, the RINEX (Receiver Independent Exchange, v2.11) files were checked for antenna type,

55 receiver type, and a priori point coordinates. Tab. S1 indicates the observation period of each

56 station. Fig. S4 shows the histogram of duration of the time series, the average being 6.5 years. Fig.

57 S5 shows the percentage of daily data completeness of the time series, the average value being 85%.

58 All a priori coordinates are the result of a precise preliminary calculation that we carried out

59 beforehand. Their accuracy is within a few centimetres, which is better that the needs of the

60 processing software in terms of initial conditions. The reference height in the files and for

61 calculations is the Antenna Reference Point (ARP) height of the antenna. For stations with log-

62 sheets providing information on changes in equipment (receiver, antenna) that could have led to

63 jumps, we took this information into account from the beginning of the calculations. For some of

64 the other stations, we found jumps during the examination of the time series, and we corrected

65 them. We processed the data with Gipsy-Oasis software (https://gipsy-oasis.jpl.nasa.gov) version

66 6.4. We used the flinnR orbits, 15° minimum elevation and ambiguity resolution, but no parameter

67 of lateral variation of the troposphere and no a priori tropospheric model. These last settings are

68 intended to keep all tropospheric residuals together (stratified and laterally variable components) so

69 that they can be used directly to correct the InSAR interferograms produced on the area. By

70 comparing our results with those published by the Nevada Geodetic Laboratory (NGL,

71 http://geodesy.unr.edu) we found that the dispersion of our daily solutions is ~30% higher than that

72 of the NGL. This small increase has no impact on the average velocity of the stations, even with

73 only one year of time series. Fig. S6 shows that the difference between the velocities estimated by

74 the two calculations is below 0.05 mm yr-1.

75

3 76

77 S3. COSEISMIC AND POSTSEISMIC DEFORMATIONS

78

79 This section supports the eponymous section 3 of the main text, providing a detailed analysis of

80 each seismic event. The events of the Tab. 1 are analysed one by one. The parameters of the

81 modelled fault that are used to correct the velocity field are listed in Tab. 2. For a number of

82 earthquakes, we have used models already published by us. For the others we re-analysed the data

83 and produced a model, this is mentioned in Tab. 2.

84 Our modelled faults are simple, with the assumption of homogeneous slip on a single fault plane.

85 Some articles published for certain earthquakes present more refined slip models characterized by

86 heterogeneous slip and geometries that may not be planar. These models are of interest, particularly

87 when they are built in synergy with the constraints given by seismology, notably an accurate

88 aftershocks location. However, they also have their drawbacks. The first of them is that their

89 complexity makes it easy to obtain ad hoc solutions fitting well the GPS displacements thanks to

90 the insufficient density of the GPS stations. The other reason for us to use simple models is that we

91 do not need highly refined coseismic models to extract a precise velocity field and to extract and

92 model the postseismic signals. Indeed, with the exception of the rare cases of GPS stations located

93 very close to the sources, both simple and complex models produce theoretical displacements

94 almost equal in medium (~20 km) or far field (> 50 km). This is because, once the depth of the fault

95 and its angular parameters are defined, all well constrained by seismology, as well as the seismic

96 moment, the detail of the slip on the fault has minor, if not undetectable, impact on the theoretical

97 displacements expected at a medium and large distances.

98 The formalism we use to model coseismic and postseismic deformations is that of Okada (1992).

99 It assumes a homogeneous elastic half-space. The earthquake is represented by a dislocation

100 occurring homogeneously in a rectangle originally centred at the hypocentre of the earthquake, and

4 101 then tuned by inverting the geodetic data. We model the postseismic deformations with this same

102 formalism. We therefore assume that the postseismic slip has a nature comparable to coseismic one,

103 thus affecting a fault or a narrow fault zone, but in a continuous manner instead of stick slip. Under

104 these simple assumptions, we can produce models that fit the data at first order. Even if, in some

105 cases, these models may not reflect completely the exact physics of the processes, their capacity to

106 map with fidelity the postseismic movements is good and they therefore allow a good separation

107 between the transient effects of the earthquakes and the secular part of the deformation.

108 We used the inverse6 code (Briole et al. 1986; Briole 2017). The GPS observations used as input

109 to model the fault parameters are those already published for the earthquakes of Movri 2008,

110 Efpalion 2010, Kefalonia 2014, Lefkada 2015, Gulpinar 2017, Kos 2017, Zakynthos 2018, and

111 those calculated here for the earthquakes of Methoni 2018, Samothraki 2014, Lesvos 2017, and all

112 other small earthquakes.

113 We found several cases in which it was needed to combine two exponentials laws to fit the

114 postseismic relaxation, e.g. in the case of the Samothraki 2014 earthquake (Fig. S2), one short

115 relaxation time (65 days) and one long relaxation time (1295 days). The need to involve two time

116 constants arises also for the earthquakes of Methoni 2008, Lefkada 2015 and Zakynthos 2018. In

117 the other earthquakes with magnitude below Mw = 6.5, the amplitude of the postseismic signal is

118 less and the dual process, if it exists, cannot be constrained by GPS, so we modelled the postseismic

119 signal with a single relaxation time.

120

121

122 S3.1 The Mw = 6.4 Skyros earthquake of July 26, 2001

123

124 This earthquake was studied with geodesy by Hollenstein et al. (2008) and Müller (2011). Using

125 the GPS displacements reported by the latter at the campaign points NSKR and CG38 and the

5 126 permanent station DION, and taking into account the seismological and structural constraints given

127 by Ganas et al. (2005), we calculated the location of the best fitting fault plane (Tab. 2; Tab. S4).

128 This plane is located 10 km south of the cloud of aftershocks with an upper edge of the fault not

129 very shallow (~5 km). The accuracy of its north-south location is around 5 km thanks to the

130 constraints of the nearby NSKR station. This suggests a blind left lateral fault located beneath the

131 south of the island rather than aligned with its northern coast as suggested by Ganas et al. (2005).

132 The antithetic plane, aligned with the overall trend of the branches of the North Anatolian fault, fits

133 slightly better the geodetic data, but it has been rejected in the past studies because it does not fit

134 with the shape of the aftershocks which is elongated in the north-west/south-east direction.

135

136

137 S3.2 The Mw = 6.2 Lefkada earthquake of August 14, 2003

138

139 We used the fault parameters determined by Ilieva et al. (2016). There were some permanent

140 stations in the near field at that time (Müller 2011) but their data are not available for reprocessing.

141 Among the station belonging to our data set the closest are those of the western Gulf of Corinth

142 where a very small signal is observed at EYPA only.

143

144

145 S3.3 The Zakynthos swarm of April 2006

146

147 The six stronger events of this swarm occurred between April 3 and April 19. Their focal

148 mechanisms are very similar. The focal mechanism in Tab. 1 was obtained by averaging the G-CMT

149 determinations, after removing for each angle the two extreme values, i.e. averaging four values.

150 Müller (2011) reports a small ground motion at the permanent station KERI. Our best fitting model

6 151 (Tab. 2) is obtained by using this unique data (Tab. S4) and inverting only for the three coordinates

152 of the upper edge of the fault, the three angles being those of the composite focal mechanism and

153 the size of the fault and amount of slip being constrained by the seismic moment tensor using the

154 formulas of Briole (2020). This earthquake does not affect any of the stations of our catalogue. The

155 modelled reverse fault is almost coplanar with that of the 2018 Zakynthos earthquake.

156

157

158 S3.4 The Mw = 6.9 Methoni earthquake of February 14, 2008

159

160 This major event on the Hellenic megathrust was studied be several teams and an exhaustive

161 review and detailed geodetic analysis can be found in Howell et al. (2017) along with plots of

162 coordinates time series. Using their GPS observations and our observations at the PYL1 station, we

163 computed a model making the assumption of single fault plane. Indeed, we found that the GPS data

164 available at the seven GPS stations (Tab. S4), located relatively far from the fault, were not

165 sufficiently discriminant to constrain simultaneously the parameters of the main shock and the main

166 Mw = 6.2 aftershock of 12:08 UT on February 14, even by making use of the seismological

167 constraints. Therefore, and especially for the purpose of this work which does not requires the

168 maximum details of the fault slip but a robust determination, we preferred to estimate a single

169 average fault merging in the best way the two events. The best fitting model (Tab. 2) gives a mean

170 scatter between observations and model of 2.2 mm in east and 1.2 mm in north, instead of 2.3 and

171 0.9 mm respectively with a dual faults model. When using the fault parameters of Howell et al.

172 (2017) we found larger scatters of 3.7 mm in east and 6.0 mm in north. Although only a limited

173 number of GPS stations are located sufficiently close to be used for modelling, this large earthquake

174 propagated small offsets of a few millimetres at the stations located in the north Peloponnese in the

175 Corinth rift area, and even beyond.

7 176 For the ten first days of postseismic period we found a good fit to the data by enlarging slightly

177 the initial fault, by 4 km in length and 3 km in width at the bottom of the fault, with no additional

178 slip in the central part already ruptured. The mean scatter of the residuals is 2.4 mm in east and 2.2

179 mm in north.

180 For the long-term postseismic we first estimated the time constant by using the various available

181 time series. We found a mean value of 449 ± 70 days. This value is well constrained by the signal at

182 the GPS station PYL1, see also the Fig. 3 of Durand et al. (2014). The GPS coverage in azimuth

183 and distance did not allow us resolving a robust variable slip model. We thus calculated a simple

184 model of homogeneous slip on a planar rectangular fault. We found the best fitting plate located at

185 shallow depth south and west of the coseismic fault, instead of a deeper fault north of the coseismic

186 one for Howell et al. (2017). The long lateral extension towards east of the fault is mostly the

187 consequence, in the inversion, of the azimuth of the postseismic vectors at PYL1 and METH. The

188 mean scatter of the residuals is 7.9 mm in east and 7.5 mm in north. Despite the difference between

189 our model and the model of Howell et al. (2017), the predicted displacements of both are close,

190 especially at the stations located in the northern Peloponnese where they are small, not usable for

191 modelling, but important to be considered to derive a correct estimation of the secular component of

192 the velocity field.

193

194

195 S3.5 The Mw = 6.4 Movri earthquake of June 8, 2008

196

197 We used the parameters of Serpetsidaki et al. (2014). The time series of RLSO (Fig. S8) shows the

198 existence of a postseismic signal not mentioned previously representing 32% of the east component

199 and 14% of the north component. We could estimate its amplitude and time constant (103 ± 35

200 days). There is a postseismic signal also at the station PAT1, located near the EUREF station PAT0,

8 201 with amplitude of -7.3, 3.4 and -5.2 mm yr-1 in east, north and vertical respectively, and a time

202 constant consistent with the one of RLSO. But PAT0 does not record this signal which, therefore,

203 might be a local disturbance of PAT1 only. Having reliable postseismic observations at RLSO only

204 we did not try to model the postseismic deformation of the Movri earthquake.

205

206

207 S3.6 The two Mw = 5.3 Efpalion earthquakes of January 18 and 22, 2010

208

209 Those two events were modelled by Elias (2013) using a stack of interferograms. The offsets

210 recorded at the GPS station EYPA clearly located in the hanging wall give an accurate tie to the

211 interferograms. While InSAR encompassed the two events, the GPS data of EYPA allows to

212 estimating approximately the relative location of the two fault planes, the first one being to the west

213 and the second to the east of the activated area. Our preferred fault plane is the north-dipping for the

214 two events. This scenario, also suggested by Ganas et al. (2013), fits slightly better the data, yet

215 other publications suggest the south-dipping solution. In additional to the clear signal at the GPS

216 station EYPA, used by Elias (2013) for his model, we found a small signal at the station LAMB

217 located on the opposite shore of the Gulf of Corinth.

218

219

220 S3.7 The Mw = 6.4 western Crete earthquake of October 12, 2013

221

222 This large earthquake that occurred close the western coast of Crete (Papadimitriou et al. 2016)

223 produced a signal slightly visible at TUC2: -5.2 ± 2 mm in east, -4.2 ± 2 mm in north, -0.2 ± 4 mm

224 in vertical. The various available moment tensor solutions (see EMSC web portal) indicate a

225 relatively wide range of depths, from 17 to 50 km, and a seismic moment ranging between 6.77 1018

9 226 N m (and depth 15 km) according to SCARDEC (http://scardec.projects.sismo.ipgp.fr; Vallée &

227 Douet 2016), and 1.7 1019 N m (and depth 44 km) according to G-CMT. With data recorded at only

228 one station, we cannot provide strong constraints but we observe that the fit is better with both the

229 depth and the seismic moment in their lower range.

230

231

232 S3.8 The Mw = 6.1 and 6.0 Cephalonia earthquakes of January 26 and February 3, 2014

233

234 We use the models of Briole et al. (2015). The time series at the station VLSM, 15 km east of the

235 epicentre, do not show significant postseismic deformation for those earthquakes. However, a

236 postseismic signal is recorded at the station KIPO, located in the close vicinity of the ruptured

237 faults, with 85 ± 15 days of relaxation time, and amplitude of -7.5, 9.6 and -7.9 mm yr-1 in the east,

238 north and vertical directions respectively. Because of the very shallow character of those

239 earthquakes, especially the second one, we have not tried to model this postseismic motion of KIPO

240 with an elastic model. We can note however that the motion is rotated towards the north-west

241 compared to the coseismic that is almost purely towards north, and the station is subsiding while it

242 was uplifted during the two earthquakes.

243

244

245 S3.9 The Mw = 6.9 Samothraki earthquake of May 24, 2014

246

247 This earthquake was studied by Kiratzi et al. (2016) and using geodesy by Saltogianni et al.

248 (2015) and Konca et al. (2018). Despite the length of the fault, the modelling of the coseismic

249 deformations can be done with a single rectangular plane and homogeneous slip on a fault located

250 between the surface and a depth of 8.5 km. Such a simple model fits ~95% of the signal. See the

10 251 observed and modelled values in Tab. S4.

252 For the postseismic we find a good fit with a model that extends the initial fault 4.5 km deeper,

253 until a depth of 13 km, and 8 km towards the west, i.e. offshore to the north-west of the island of

254 Lemnos. The relaxation time is 1295 ± 200 days. It is not very well constrained because its duration

255 is commensurable with the entire duration of the time series of the stations that sample it, e.g.

256 LEMN (Fig. S2). This uncertainty does not affect the value of 4.5 km found for the downward

257 extension of the slip during the relaxation process. Below this second layer we speculate that there

258 is a continuous creep between the two lithospheric blocks located on either side of the North

259 Anatolian fault. This locking depth of 13 km (8.5 km coseismic + 4.5 km postseismic) is close to

260 that found by Konstantinou (2017) for the same area (16-17 km). The analysis of the spatial

261 distribution of the interseismic GPS velocities that is often used to estimate locking depths is

262 difficult to be used here because of the insufficient number of islands located in appropriate

263 locations on both sides of the fault. Further east, in the Marmara Sea the reported locking depth is

264 between 10 and 15 km (e.g. Schmittbuhl et al. 2016; Bohnhoff et al. 2016; Uchida et al. 2019) thus

265 comparable to the one we find from the modelling of the deformation of this earthquake.

266

267

268 S3.10 The Mw = 6.5 Lefkada earthquake of November 11, 2015

269

270 There are several studies involving the geodetic data GPS and/or InSAR, e.g. Melgar et al. (2017),

271 Bie et al. (2017), Ganas et al. (2016). We use the model of the latter. The fault plane is well

272 constrained by the efficient combination of InSAR and GPS. The time series show the existence of

273 two time constants for the postseismic signal. The first one, with 138 days of time constant, is

274 observed only at the nearby stations. We can model it with a small extension of the fault plane all

275 around the coseismic one, with an increase of 11% of the seismic moment. The second part of the

11 276 postseismic signal is different and can be observed and quantified a greater distance from Lefkada,

277 until Agrinion and the centre of Cephalonia (Tab. S4). In the near field, with the exception of PONT

278 and SPAN, the stations located on Lefkada do not have sufficient preseismic data to be used for

279 postseismic extraction since their secular velocity is not yet known. The structure of the postseismic

280 signal imposes a fault dipping at low angle, 15 ± 5° in our best fitting model. The upper tip of this

281 postseismic plane is located beneath Lefkada in the area of the coseismic event, yet difficult to

282 constrain accurately within ± 5 km in horizontal. In our model we forced the postseismic fault to be

283 geometrically in the continuity of the coseismic one. However, we did not put other constrains and

284 in particular the azimuth of the postseismic slip was estimated. Although the residuals are large at

285 some stations, there are many stations involved in a large range of azimuths so the parameters of the

286 modelled fault are robust within the error bars. We find a significant component of thrusting with an

287 azimuth N54 ± 10° for the slip vector.

288

289

290 S3.11 The Gulpinar seismic swarm of February 2017

291

292 We use the geodetic model of Ganas et al. (2018b). This earthquake swarm, well constrained by

293 InSAR, has almost no impact on the GPS stations of Greece except, very slightly, for the point

294 AGPA located in the northern part of Lesvos.

295

296

297 S3.12 The Mw = 6.3 Lesvos earthquake of June 12, 2017

298

299 Because of its offshore location, this earthquake is not well covered by GPS and not covered by

300 InSAR as it occurred too far away from the south shore of Lesvos. GPS does not allow resolving

12 301 the ambiguity between south and north-dipping solution. The south-dipping is preferred by Kiratzi

302 (2018), Papadimitriou et al. (2018), Chousianitis & Konca (2018) based on the seismological data

303 and the geomorphology of sea bottom. It was also proposed by Briole et al. (2018) based on the

304 analysis of the GPS and InSAR data. In Tab. 2 we assume the south-dipping plane. For the purpose

305 of removing the coseismic signal at the surrounding GPS stations and for assessing the secular

306 velocity field, the two planes are in any case equivalent. A postseismic signal is observed at the four

307 stations of Lesvos with two time constants, one of 4 ± 1 days, and one of 140 ± 20 days.

308

309

310 S3.13 The Mw = 6.6 Kos earthquake of July 20, 2017

311

312 The combination of GPS and InSAR allows to resolving the fault plane ambiguity, with the north-

313 dipping plane being largely preferred, mostly thanks to the structure of the fringes in the Karaada

314 Islet and inland further north (Karasözen et al. 2018; Ganas et al. 2019; Konca et al. 2019). There

315 are very few points of the Greek network nearby and this earthquake affects very marginally the

316 overall velocity field of Greece. A slight postseismic relaxation is observed far away at RODO, with

317 time constant 220 ± 50 days, but not at KATC which has more noise than RODO. There is no

318 evidence of postseismic motion at TILO. The postseismic time series of the nearby Turkish stations,

319 as well as those of the HEPOS station 086A in the Kos Island, could provide further constraints.

320

321

322 S3.14 The Mw = 6.8 Zakynthos earthquake of October 25, 2018

323

324 This earthquake was studied by combining geodetic and seismological data by Chousianitis &

325 Konca (2019), Cirella et al. (2020), and Ganas et al. (2020), and by using mainly seismological data

13 326 by Sokos et al. (2020). Here we use the model of Ganas et al. (2020) which is the simplest and we

327 believe the most robust available. It is based on GPS only but we checked that the synthetic

328 interferograms it produces fit with the coseismic interferograms calculated inland Zakynthos. There

329 is a fast postseismic described by Ganas et al. (2020), with 21 ± 5 days of time constant, that they

330 modelled with the connection of the upper edge of the coseismic fault (modelled at a depth of 3 km

331 beneath the seafloor) with the surface. The slip on this shallow fault in the days following the

332 earthquake adds 4.45 1018 N m (i.e. 18%) to the coseismic seismic moment. Over longer time period

333 there is a postseismic signal that can be observed in a much wider area of radius larger than 100 km.

334 To the north, in particular, the dense array of the Corinth rift laboratory (http://crlab.eu) documents

335 well the postseismic signal and its attenuation with distance. We could produce a well constrained

336 model (Tab. 2). In that model the stations located in the near field, in particular AMAL and RLSO,

337 play an important role to constrain the depth and the dip angle. The best fitting model is found with

338 an elongated low angle plane that could correspond to the subduction interface, with a total slip of

339 90 mm and a relaxation time of 183 ± 25 days, much shorter than that of Methoni 2008.

340 This postseismic plane for the 2018 earthquake appears as a narrow band elongated downwards

341 (78 km). This is well constrained by the data. We tested various geometries, and this one is clearly

342 preferred by the modelling, strongly constrained in this respect by the postseismic vectors at

343 AMAL, PYRG, ZAKU, and ZAKY for the length of the fault trace (26 km), and by the stations in

344 the north and north-west of the Gulf of Corinth for the downward elongation. The reverse slip

345 amplitude, 90 mm, is well constrained, in a consistent manner, by all stations. The modelled plane

346 begins in its upper part at the junction with the lower end of the coseismic plan and extends

347 northwards to end at a depth of 30 km towards the southern end of the rupture zone of the Movri

348 earthquake of 2008. In terms of connection between the subduction plane and the base of the crustal

349 faults, this geometry is consistent with the schematic description made by Serpetsidaki et al. (2014)

350 in their conclusive Fig. 13. The azimuth of the slip is N40°E which is rotated 21° clockwise with

14 351 respect to the N19°E relative motion of the Peloponnese with respect to Africa mentioned by Ganas

352 et al. 2020a (see their Tab. 7 and Fig. 9), using the average vector of AMAL, RLSO and PYRG.

353

354

355 S3.15 The Mw = 5.4 Preveza earthquake of February 5, 2019

356

357 We found a small signal in the time series of LEFK, LEUK, SPAN, LEFK (Tab. S4, and see

358 additional KML file) that is very likely to correspond to that earthquake. The observed average

359 displacement of the triplet of nearby stations LEFK, LEUK and LFKD is 2.1 and -2.4 mm in east

360 and north respectively, while the modelled values, using the parameters of the G-CMT, are 1.6 and

361 -1.9, thus close to the observations. This offshore earthquake was studied in detail by Kostoglou et

362 al. (2020) who located the main shock slightly further north. Such location gives a lower quality of

363 fit to the GPS data.

364

365

366 S3.16 The Mw = 5.3 Eratini earthquake of March 30, 2019

367

368 This moderate event occurred at the east termination of the rupture of the Mw = 6.2 earthquake of

369 June 15, 1995 (Briole et al. 2000), and west of the Mw = 5.9 of November 18, 1992 (Hatzfeld et al.

370 1996). It is an unusual event in the Gulf of Corinth area because of the small number of aftershocks.

371 A coseismic signal is visible and can be quantified at PSAR and VALI. For the modelling we make

372 the assumption of north-dipping low angle plane for consistency with the 1992 and 1995 planes, but

373 the south-dipping one would fit the GPS data equally well. InSAR shows no deformation inland.

374

375

15 376 S3.17 The Mw = 5.3 Magoula earthquake of July 19, 2019

377

378 This earthquake awakened in the Athenians the memory of the Mw = 5.9 earthquake of September

379 7, 1999 (Kontoes et al. 2000) as it occurred in the western Parnitha region of Attica. ATHI register a

380 small signal visible in its time series (Fig. S9), of 1.8 ± 1 mm in east, -2.8 ± 1 mm in north, -0.7 ± 3

381 mm in up. The fault strikes N106°E, dips to the south, and is located 10 km west of the one of 1999.

382 It may be coplanar to the 1999 one according to Kapetanidis et al. (2020) and Kouskouna et al.

383 (2019). It should be noted that the above authors indicate a magnitude Mw = 5.1 while the global

384 networks (e.g. G-CMT, GFZ) indicate Mw = 5.3.

385

386

387 S3.18 The Mw = 5.7 Kanallaki earthquake of March 21, 2020

388

389 This significant earthquake produced ground deformation very well measure by ascending and

390 descending Sentinel-1 SAR interferometry, which allowed to producing a robust model for the fault

391 plane (Valkaniotis et al. 2020). The nearby GPS station of GARD contains a small signal but we

392 cannot yet measure it accurately. To this end, like for other moderate events, there is the need of

393 long enough postseismic measurements, at least one year. Nevertheless, as the InSAR coverage is

394 good, we do not expect GPS to bring additional constraints for the coseismic field.

395

396

397 S3.19 Other earthquakes

398

399 The Mw = 7.4 Izmit earthquake of August 17, 1999 and the Mw = 5.9 Athens earthquake of

400 September 7, 1999 are visible in the time series of ISTA and DION respectively. For Izmit the time

16 401 series allows to retrieve the two relaxation times already described in literature (Ergintav et al.

402 2002) and estimate their values: 1150 ± 200 days for the long one and 282 ± 35 days for the short

403 one. The long relaxation time is comparable to that found for the Samothraki 2014 earthquake. This

404 suggests comparable viscosity of the underlying material involved in the relaxation. For Athens

405 1999 the displacement we find at DION is consistent with the one found by Müller (2011). From

406 InSAR data it was determined that the total seismic moment calculated in terms of afterslip is 49%

407 of the seismic moment released by the mainshock (Atzori et al. 2008), however the GPS data of

408 DION do not show this postseismic motion but the coseismic displacement only.

409

410

411 S4. OTHER TRANSIENTS

412

413 We identified specific events in the time series of several stations. The most noteworthy are listed

414 below. Some of these events are temporally correlated with earthquakes. Others are not. For some

415 events we propose an interpretation. For others we could not find yet a satisfactory interpretation.

416

417

418 S4. 1 A slow earthquake in May 2018 below Zakynthos?

419

420 In Fig. S10 and Fig. S11 we see, in the time series of the stations AMAL, PYRG and TROP,

421 located in the western Peloponnese, a signal that affects the horizontal and vertical components of

422 the stations. The amplitude of this signal is low but sufficient to be detected, especially because its

423 effect is permanent along the east and vertical axis, and partly recovered before the event of October

424 2018 along the north axis. By assuming a simple instantaneous step in the time series, and

425 minimizing the residuals, we find that the date of the discontinuity is May 10, 2018 ± 3 days. We

17 426 tried to model the displacements assuming a fault plane coplanar with the postseismic of Zakynthos

427 2018. The best fitting solution is obtained with a plane of 24 km length and 13 km width and 90

428 mm of right lateral slip. The azimuth of slip is imposed by the direction of the transient vectors

429 which are close to perpendicular from the azimuth of the place convergence there estimated at

430 N19°E by Ganas et al. (2020). The complete parameters are at the bottom of Tab. 2. Tab. S5 gives

431 the observed and modelled displacements at AMAL, PYRG and TROP as well as the predicted

432 displacements at the nearby stations KOPA, ZAKY, STRF, and GARG where, by visual inspection,

433 we could not see an offset in the time series. The geodetic moment of this modelled event is 78 10 16

434 N m which corresponds to a magnitude Mw ~ 5.8.

435

436

437 S4.2 The unrest of Santorini 2010-2012

438

439 From 2010 to 2012 a seismic crisis occurred at shallow depth within the caldera of Santorini

440 accompanied with ground deformation (Foumelis et al. 2013; Parks et al. 2015). All GPS time

441 series available at Santorini are biased by this unrest. The station NOMI, which documents

442 continuously the unrest period, shows (Fig. S12) that the unrest period can be quantified precisely

443 with approximately one-week accuracy from February 12, 2011 to May 3, 2012. Fig. S13 (top)

444 shows the east component of the eleven available stations. During the unrest the data of NOMI and

445 the other data show an almost linear trend. We therefore modelled the unrest period with a linear

446 anomaly lasting this 446 day interval. Fig. S13 (bottom) shows the time series after applying this

447 correction. The anomalies of velocities during the unrest are listed in Tab. S5. The secular

448 velocities, corrected for the effect of the unrest, are included with the others in Tab. S1. For the

449 whole island, the average secular ITRF2014 velocity of the eleven stations of Santorini, calculated

450 by weighting each velocity with its uncertainty is 7.4, -15.5 and 0.7 mm yr-1 in east, north and up

18 451 respectively. This is the value we recommend for local geodetic works made for the analysis of the

452 volcanic activity. We found no significant change of the average velocity between before and after

453 the unrest. As a whole, the island of Santorini is uplifting by 0.7 mm yr-1. This global uplift exists

454 before and after the unrest. We modelled the anomalous slopes observed during the unrest with a

455 simple model of Mogi (1958). The best fitting solution is with a centre of inflation located at

456 25.386° E, 36.428° N and 3.7 km depth, and an increase of volume of the source of 0.27 m3 s-1,

457 which corresponds to a yearly uplift on top of the reservoir of 189 mm yr-1 and a total volume of

458 12.6 106 m3 for the whole 446 days of unrest. The mean scatter of the residuals is 8.4 mm yr-1 which

459 is 18.4% of the mean scatter of the initial data. We tried also a model of vertical dyke of 2 x 2 km

460 and 3.15 m opening but we could not reduce the variance of our residuals below 35% of the initial

461 values (i.e. two times higher than the Mogi model). For the modelling we removed the vertical

462 transient velocities of MKMN and DSLN which were found to be outliers. For the post-unrest

463 period we analysed the relative velocities obtained by subtracting the average secular velocity of

464 Santorini indicated above. Those relative velocities are listed in Tab. S5 along with the modelled

465 velocities obtained using again a Mogi source. This best fitting source for the post-unrest period

466 corresponds to a slight deflation on top of the unrest source of -3.1 mm yr -1. However, the mean

467 residual is 1.2 mm yr-1 thus close to the amplitude of the data. This local deflation can be observed

468 by InSAR (Parks et al. 2015) better than with GPS. The local character of this post-unrest deflation

469 does not hide the global continuation of the uplift motion of the island. Two stations located to the

470 east of the island, WNRY and THIR, show a trend very different from the one predicted by the

471 deflation source. This suggests the occurrence of local deformation around those stations.

472

473

474 S4.3 Earthquake-triggered subsidence of basins or deltas

475

19 476 After the Zakynthos earthquake of October 25, 2018, a subsidence is recorded (Fig. S14) at the

477 GPS stations VALI and XILI located in the deltas of the rivers Selinountas and Mornos in the

478 western Gulf of Corinth. We modelled the signals with decreasing exponentials which proves to fit

479 well the data. The time constant is 36 days for VALI and 16 days for XILI and the amplitudes -10.2

480 mm and -9.6 mm respectively. We looked at the borehole water level data, acquired by the Institute

481 of Geology and Mineral Exploration of Greece (IGME), for possible variations correlated with

482 these transients, but we did not find any. However, the IGME data are not from permanent sensors

483 but from episodic observations made once or twice a year. They are therefore not appropriate to

484 detect rapid changes unless those changes have a permanent signature afterwards, which appears

485 not to be the case. We note also that, while the pre-seismic trend of XILI is subsidence at a rate of

486 6.2 mm yr-1, the pre-seismic trend of VALI is uplift at a rate of 2.4 mm yr -1. We explain this

487 peculiarity of VALI by the activity of the nearby off-shore normal north-dipping Aigion fault that

488 tends, on the long-term, to raise the Selinountas delta. A modelling of the activity of this fault was

489 made by Elias & Briole (2018). Fig. S15 shows that there is no significant displacement recorded by

490 InSAR at a broad scale of a few kilometres around the GPS stations VALI and XILI. Therefore, the

491 signals recorded by GPS are likely to correspond to local subsidence (or acceleration of the

492 subsidence in the case of XILI) in the vicinity of the buildings that host the GPS antennas.

493

494

495 S4.4 Movements related to aquifers not associated with earthquakes

496

497 Land subsidence occurs in several places in the Thessaly plain, central Greece, because of the

498 water pumping for agriculture and it can be observed by InSAR (Foumelis et al. 2016). The most

499 flagrant case is that of the STEF station in Thessaly, whose temporal variations were studied by

500 Argyrakis et al. (2020). Fig. S16 shows that, despite the very high subsidence rate of STEF, this

20 501 subsidence is missed by P-SBAS. This suggests that the deformation captured by the GPS station

502 STEF is local with its origin at very shallow depth.

503 There are other cases such as the KAL1 and KAL2 stations located in the Gallikos River delta

504 west of Thessaloniki. This area is well known for the considerable vertical movements that affected

505 it in the 1990s and 2000s (Raucoules et al. 2008). Subsequently the vertical velocities decreased

506 significantly as shown by the time series of the coordinates of the above-mentioned stations.

507 On January 12, 2019 a ground deformation of -1.9, 0.2, 0.0 mm in east, north and up respectively

508 is detected (Fig. S17) at the station TRIP (Tripoli) in the central Peloponnese. The nearby station

509 TRIU located 240 m to the east detects a change of -1.8, -0.1 and 0.0 mm. In the following days the

510 deformation continues. It can be modelled with an exponential of time constant 17 days and the

511 following final amplitudes of deformation: -5.2, -9.5, 29.7 mm at TRIP, -4.1, -7.9, 28.3 mm at

512 TRIU. No anomalous signal is detected at the other GPS stations around TRIP in that period

513 (ARGU, LEON, SPRT, TROP, VASS) so the origin of this uplift cannot be at the depth of the

514 subduction interface (~45 km beneath TRIP according to Halpaap et al. (2018)) otherwise the signal

515 would affect also those stations, even with lower amplitude. Assuming now a local source of

516 deformation, a Mogi model fits well the horizontal and vertical displacements of both stations with

517 a source located at 22.3721°E, 37.5146°N, i.e. 750 m north of TRIP, buried at 2.2 km beneath the

518 town. The increase of volume needed to produce the observed uplift is 0.66 106 m3. However, none

519 of the eleven ascending and descending Sentinel-1 interferograms we calculated on the area for the

520 period under consideration (e.g. the one shown in Fig. S18) show the deformation detected by GPS.

521 This invalidates the Mogi model of local deformation. Alternatively, there might be the possibility

522 of a signal at intermediate wavelength that would affect the karst of the entire plain of Tripoli.

523

524

525 S4.5 Response of the remote aquifer of Kalymnos Island to the Lesvos earthquake of 2017

21 526

527 Immediately after the Mw = 6.3 Lesvos earthquake of June 12, 2017, the two GPS stations located

528 on the island of Kalymnos show a deformation transient with relaxation time 16 days (Fig. S19, Fig.

529 S20). The amplitude of the displacement is 1mm towards east on June 12 followed by 3.4 mm in the

530 azimuth N330°E for KALU and 0.3 mm towards east followed by 3.5 mm in the azimuth N3°E for

531 KALY. There is no vertical anomaly. The distance between Kalymnos and Lesvos (~230 km) is too

532 large for it to be an elastic response to coseismic dislocation. This obvious signal, perfectly

533 correlated at the two stations, may correspond to the response of the shallow aquifer of Kalymnos to

534 the shaking produced by the Lesvos earthquake. The island hosts an important aquifer located

535 beneath the valley of Vathi, north of the valley of Kalymnos. This aquifer has been studied and

536 mapped in detail by Mantoglou (2003).

537

538

539 S4.6 Landslides

540

541 Very few of the GPS stations have been found to be on landslides. However, this is the case of the

542 KRIN station (Fig. S21) located in the mountains overlooking Aigion. In Fig. S22 the anomalous

543 vectors are plotted along with line of sight velocities obtained from a multitemporal analysis of

544 Sentinel-1 data acquired on the ascending track 175 from January 1, 2017 to May 15, 2020. This

545 multi-temporal interferogram made with the TRE-ALTAMIRA FASTVEL processor on the

546 Geohazards Thematic Exploitation Platform (GEP-TEP, https://geohazards-tep.eu), as well as

547 another made with the CNR-IREA P-SBAS processor on the GEP-TEP (Berardino et al. 2002; Casu

548 et al. 2014; see Fig. S23) allow to map the elongation of the landslide and to quantify its velocity in

549 the line of sight of the InSAR observation. The GPS station KRIN is located inside the sliding area,

550 near its southern edge. The velocity from GPS is consistent with the one from InSAR. The slope of

22 551 the terrain and the slope of the vertical/horizontal velocity are consistent. The velocity is stable with

552 small fluctuations (less than 10%) in the east component, not correlated with the seasons. We do not

553 see any anomaly of velocity at the time of the Zakynthos 2018 earthquake.

554

555

556 S4.7 Abnormal movements correlated with earthquakes and possibly due to local instability of

557 the GPS antenna

558

559 The suspicion of such a case exists for the STRA station (Chalkidiki), which is located inside an

560 industrial excavation site. The displacements recorded there, in correspondence with the Samothraki

561 earthquake of May 24, 2014, cannot be explained by the Okada model which nevertheless explains

562 all other observations very well. This is particularly the case for the postseismic part. We believe

563 that there is an anomaly at this station related to the particular situation of the antenna in the

564 excavation site.

565 In some other cases we have identified in the time series possible clues of local movement of the

566 building supporting the GPS antenna. This is the case of GEYB located on the Gefyra building

567 north of the Rio-Antirio bridge, the building itself being built on thick embankments. We also have

568 suspicions of instability of the high pillar (over 3m) of TRIZ (Trizonia) or the shelter holding ROD3

569 (Ano Rodini). In those three cases the instabilities remain low enough not to reduce too much the

570 estimation of secular velocities, especially at TRIZ where the time series is almost twenty years

571 long. There are however a few stations where the signal is so disturbed that it becomes unusable, for

572 example the eastern component of PYL1 or all components of the POLY station.

573

574

575 S5. ANALYSIS OF VELOCITIES PROFILES

23 576

577 Here we present the velocity profiles made along the three red lines plotted in Fig. 7. Those lines

578 are passing through the Euler pole Anatolia-Eurasia of Le Pichon & Kreemer (2010). The names of

579 the lines are given by large cities located near their northward terminations. Profiles give a

580 representation of the velocity gradients that is easier to quantify visually and analyse, with the

581 disadvantage of concealing the rotations. Fig. S25, show the vectors amplitudes (in the radial and

582 tangential directions) of the GPS stations located at less than 55 km of the “Istanbul profile” of Fig.

583 7. Fig. S26 shows in a map the points used in the “Sandanski profile” and Fig. S27 the vectors

584 amplitudes. Fig. S28 shows in a map the points used in the “Tirana profile” and Fig. S29 the vectors

585 amplitudes. We also show other velocities profiles that are discussed in the main text.

586

587

588 S5.1 Istanbul profile

589

590 In Fig. S25 the tangential velocity of the southern GPS stations of MEGI, FINI, ANTL, DNZL

591 and DEIR is consistent with the rotation of Anatolia around the pole of Le Pichon & Kreemer

592 (2010), yet our best fit is obtained with a rotation rate of 1.56° Ma-1, 19% larger than their value.

593 However, on average, those five stations are showing a small component of southward motion

594 which indicates that the plate is no longer rigid at this longitude. North of DEIR the amplitude of

595 the tangential velocity starts decreasing because we cross the southern branch of the North

596 Anatolian fault. At BAND, which is close to the main segment of the North Anatolian fault, the

597 perpendicular velocity is 19.8 mm yr-1. North of the North Anatolian fault, the vectors are almost

598 zero except for the station YSST which is close to the fault and sensitive to its seismic cycle.

599

600

24 601 S5.2 Sandanski profile

602

603 This profile (Fig. S26, Fig. S27) corresponds to points close to a line, aligned with the AT-EU pole

604 of Le Pichon & Kreemer (2010) connecting approximately Sandanski (Bulgaria) to Rodos. The line

605 is passing close to the islands of Rodos, Kos, Samos, Ikaria, Lesvos, Lemnos and continuing

606 towards the western border of Bulgaria.

607 From north to south, we see that there is minor deformation at KUST, SOFI and YUND with

608 respect to LOVE and the first large step of deformation occurs north of Sandanski, possibly on the

609 latitude range of the Krupnik normal fault which is oriented east-west. Then, there is no large

610 change until the latitude of Vrasna-Kavala where another step is clear in the data, possibly in

611 correspondence to crustal extension across the Volvi (Mygdonia) graben. Then, there is a third

612 change shown by the two stations of Neos Marmaras in the central peninsula of Chalkidiki. The

613 major discontinuity in the profile is, logically, the one between Chalkidiki and the island of Lemnos,

614 corresponding to the location of the main branch of the North Anatolian fault. The radial data show

615 that Lemnos and Lesvos present no anomaly with respect to their expected velocity in a rigid block

616 model. This is also the case of IZMI, AYD1, MUGL and FETH further south-east.

617 Because of the vicinity to the North Anatolian fault the vectors of the three stations of Lemnos

618 contain a component of strain accumulation. The grey dots show the values after correction (3.3, 3.2

619 and 2.7 mm yr-1 for LIMN, LEMN and MOUD respectively), thus corresponding to the velocity at

620 the fault axis. The correction was calculated with the formula v(x) = V/π tg-1 x/D with V fault

621 velocity, D locking depth and x distance to the fault. The average value after correction is 21.4 mm

622 yr-1 which is close to the value measured at BAND (mentioned in section 5.1 earlier and Fig. S25).

623 Lemnos and Lesvos, both show no anomaly of their radial component (open circles), but present a

624 deficit of tangential motion where we see a step twice larger between Lesvos and Chios (8 mm yr -1)

625 than between Lemnos and Lesvos (4 mm yr-1).

25 626

627

628 S5.3 Tirana profile

629

630 This profile (Fig. S28, Fig. S29) is passing close to Karpathos, Astypalea, Naxos, and continues

631 along the gulf of Evia and through Thessaly until Kastoria and Ohrid. Along this profile we can see

632 clearly the transition zone between the Anatolian and the European domains. From the south-east to

633 the north-west, the tangential velocities decrease mostly from the south of KYMI (Kymi, eastern

634 Evia) to around GREV (Grevena), thus in a band of ~240 km, with a variation more abrupt near

635 Kymi than near Grevena. Then there is a remaining 3 mm yr-1 to be absorbed between GREV and

636 stable Europe. The profile shows that the rate is much slower there. Along the radial axis we see

637 that, for the stations located in the 240 km band, there is a small component of radial motion

638 towards the pole (corresponding to extension in the domain) decreasing from ~6 mm yr-1 around

639 KYMI to ~3 mm yr-1 around GREV.

640

641

642 S5.4 Other profiles

643

644 Figs. S30, S31 and S32 show three profiles made respectively across the south Peloponnese (east-

645 west profile), Crete, and central Macedonia (north-south profile). Those three profiles are discussed

646 in the main text.

647

648

649 S6. MISCELLANOUS

650

26 651 There are three figures shown in this supplementary material and not discussed here but mainly in

652 the main text: Fig. S24 shows the difference between the time series of coordinates of the EUREF

653 station PAT0 with and without correction of the effect of the earthquakes, Fig. S33 shows the

654 probability density function of the vertical velocities for the whole set of 322 GPS stations and for

655 the subset of 177 stations with velocity uncertainties below 1 mm yr -1. Fig. S34 shows the ten

656 blocks model with the available focal mechanism of earthquakes for the period 2003-2018.

657 658

659 ATTACHED FILES

660

661 KML file with the location of the GPS points (briole-et-al_GPS-points.kml)

662

663 KML file with the location of the modelled fault of the studied earthquakes (briole-et-al_analysed-

664 earthquakes.kml)

665

666 KML file with the boundaries of the blocks model (briole-et-al_blocks-model.kml)

667

668 Tab. S1 file in text format with fields are separated by semicolons (briole-et-al_tabS1.csv)

669

670 Tab. S3 file in text format with fields are separated by semicolons (briole-et-al_tabS3.csv)

671

672 Tab. S7 file in text format with fields are separated by semicolons (briole-et-al_tabS7.csv)

673

674 Folder with time series of the GPS coordinates in CSV format (briole-et-al_GPS-time-series). The

675 seven columns contain the following data: 1: date, 2,3,4: daily east, north and vertical UTM34

676 coordinates (in m) converted from the ITRF2014 coordinates, 5,6,7: residual daily east, north and

27 677 vertical time coordinates, in mm, obtained after correction of the instrumental jumps, coseismic and

678 postseismic signals, and annual fluctuation; there is a shift of 20, 60 and 120 mm in east, north and

679 vertical to make plots easier to visualise. By removing those residuals to the original time series, the

680 reader can separate and visualise the model, and therefore see the jumps and the post-seismic

681 signals as estimated by the model.

682

683 Folder with the programme HelVel2020 and documentation (briole-et-al_helvel2020-model)

684

685

686 REFERENCES

687

688 Altamimi, Z., Rebischung, P., Metivier, L. & Collilieux, X., 2016. ITRF2014: A new release of the

689 International Terrestrial Reference Frame modeling nonlinear station motions, Journal of

690 Geophysical Research-Solid Earth, 121, https://doi.org/10.1002/2016JB013098

691 Argyrakis, P., Ganas, A., Valkaniotis, S., Tsioumas, V., Sagias, N. & Psiloglou, B., 2020.

692 Anthropogenically induced subsidence in Thessaly, central Greece: new evidence from GNSS

693 data, Natural Hazards, 102, 1, 179-200, https://doi.org/10.1007/s11069-020-03917-w

694 Atzori, S., Manunta, M., Fornaro, G., Ganas, A. & Salvi, S., 2008. Postseismic displacement of the

695 1999 Athens earthquake retrieved by the Differential Interferometry by Synthetic Aperture Radar

696 time series, Journal of Geophysical Research-Solid Earth, 113,

697 https://doi.org/10.1029/2007JB005504

698 Berardino, P., Fornaro, G., Lanari, R & Sansosti, E., 2002. A new algorithm for surface deformation

699 monitoring based on small baseline differential SAR interferograms, IEEE Transactions on

700 Geoscience and Remote Sensing, 40(11), 2375-2383, https://doi.org/10.1109/TGRS.2002.803792

701 Bie, L.D., Gonzalez, P.J. & Rietbrock, A., 2017. Slip distribution of the 2015 Lefkada earthquake

28 702 and its implications for fault segmentation, Geophysical Journal International, 210(1), 420-427,

703 https://doi.org/10.1093/gji/ggx171

704 Bohnhoff, M., Martínez-Garzón, P., Bulut, F., Stierle, E. & Ben-Zion, Y., 2016. Maximum

705 earthquake magnitudes along different sections of the North Anatolian fault zone,

706 Tectonophysics, 674, 147-165, https://doi.org/10.1016/j.tecto.2016.02.028

707 Briole, P., De Natale, G., Gaulon, R., Pingue, F. & Scarpa, R., 1986. Inversion of geodetic data and

708 seismicity associated with the Friuli earthquake sequence (1976-1977), Annales Geophysicae,

709 4(B4), 481-492.

710 Briole, P., Rigo, A., Lyon-Caen, H., Ruegg, J.C., Papazissi, K., Mitsakaki, C., Balodimou, A., Veis,

711 G., Hatzfeld, D., Deschamps, A., 2000. Active deformation of the Corinth rift, Greece: Results

712 from repeated Global Positioning System surveys between 1990 and 1995, Journal of

713 Geophysical Research-Solid Earth, 105(B11), 25605-25625,

714 https://doi.org/10.1029/2000JB900148

715 Briole, P., Elias, P., Parcharidis, I., Bignami, C., Benekos, G., Samsonov, S., Kyriakopoulos, C.,

716 Stramondo, S., Chamot-Rooke, N., Drakatou, M.L. & Drakatos, G., 2015. The seismic sequence

717 of January-February 2014 at Cephalonia Island (Greece): Constraints from SAR interferometry

718 and GPS, Geophysical Journal International, 203(3), 1528-1540,

719 https://doi.org/10.1093/gji/ggv353

720 Briole, P., 2017. Modelling of earthquake slip by inversion of GNSS and InSAR data assuming

721 homogeneous elastic medium, Zenodo. http://doi.org/10.5281/zenodo.1098399

722 Briole, P., Ganas, A., Karastathis, V., Elias, P., Mouzakiotis, P., Valkaniotis, S., Foumelis, M.,

723 Moshou, A., Kourkouli, P., Roger, M., Chouliaras, G., Argyrakis, P., Tselentis, G.A. &

724 Parcharidis, I., 2018. The June 12, 2017 M6.3 Lesvos offshore earthquake sequence (Aegean

725 Sea, Greece): fault model and ground deformation from seismic and geodetic observations,

726 Geophysical Research Abstracts, 20, EGU2018-18189, EGU General Assembly 2018,

29 727 https://meetingorganizer.copernicus.org/EGU2018/EGU2018-18189.pdf

728 Briole, P., 2020. Monitor (estimates seismic fault parameters from available centroid moment

729 tensors, version 3.2). Zenodo. http://doi.org/10.5281/zenodo.4295872

730 Casu, F., Elefante, S., Imperatore, P., Zinno, I., Manunta, M., De Luca, C. & Lanari, R., 2014.

731 SBAS-DInSAR Parallel Processing for Deformation Time-Series Computation, IEEE Journal of

732 Selected Topics in Applied Earth Observations and Remote Sensing, 7(8), 3285-3296,

733 https://doi.org/10.1109/JSTARS.2014.2322671

734 Chousianitis, K. & Konca, A.O., 2018. Coseismic slip distribution of the 12 June 2017 Mw = 6.3

735 Lesvos earthquake and imparted static stress changes to the neighbouring crust, Journal of

736 Geophysical Research-Solid Earth, 123(10), 8926-8936, https://doi.org/10.1029/2018JB015950

737 Chousianitis, K. & Konca, A.O., 2019. Intraslab deformation and rupture of the entire subducting

738 crust during the 25 October 2018 Mw 6.8 Zakynthos earthquake, Geophysical Research Letters,

739 46(24), 14358-14367, https://doi.org/10.1029/2019GL085845

740 Cirella, A., Romano, F., Avallone, A., Piatanesi, A., Briole, P., Ganas, A., Theodoulidis, N.,

741 Chousianitis, K., Volpe, M., Bozionellos, G., Selvaggi, G. & Lorito, S., 2020. The 2018 M w 6.8

742 Zakynthos (Ionian Sea, Greece) earthquake: seismic source and local tsunami characterization,

743 Geophysical Journal International, 221(2), 1043–1054, https://doi.org/10.1093/gji/ggaa053

744 Durand, V., Bouchon, M., Floyd, M.A., Theodulidis, N., Marsan, D., Karabulut, H. & Schmittbuhl,

745 J., 2014. Observation of the spread of slow deformation in Greece following the breakup of the

746 slab, Geophysical Research Letters, 41(20), 7129-7134, https://doi.org/10.1002/2014GL061408

747 Elias, P. 2013. Ground deformation observed in the western Corinth rift (Greece) by means of SAR

748 interferometry, PhD Thesis, Ecole Normale Supérieure of Paris and University of Patras,

749 https://tel.archives-ouvertes.fr/tel-00839348

750 Elias, P. & Briole, P., 2018. Ground Deformations in the Corinth Rift, Greece, Investigated Through

751 the Means of SAR Multitemporal Interferometry, Geochemistry Geophysics Geosystems, 19(12),

30 752 4836-4857, https://doi.org/10.1029/2018GC007574

753 Ergintav, S., Burgmann, R., McClusky, S., Çakmak, R., Reilinger, R.E., Lenk, O., Barka, A. &

754 Özener, H., 2002. Postseismic deformation near the Izmit earthquake (17 August 1999, M7.5)

755 rupture zone, Bulletin of the Seismological Society of America, 92(1), 194–207,

756 https://doi.org/10.1785/0120000836

757 Foumelis, M., Trasatti, E., Papageorgiou, E., Stramondo, S. & Parcharidis, I., 2013. Monitoring

758 Santorini volcano (Greece) breathing from space, Geophysical Journal International, 193(1),

759 161-170, https://doi.org/10.1093/gji/ggs135

760 Foumelis, M., Papageorgiou, E. & Stamatopoulos, C., 2016. Episodic ground deformation signals in

761 Thessaly Plain (Greece) revealed by data mining of SAR interferometry time series,

762 International Journal of Remote Sensing, 37(16), 3696-3711,

763 https://doi.org/10.1080/01431161.2016.1201233

764 Ganas, A., Drakatos, G., Pavlides, S.B., Stavrakakis, G.N., Ziazia, M., Sokos, E. & Karastathis,

765 V.K., 2005. The 2001 Mw = 6.4 Skyros earthquake, conjugate strike-slip faulting and spatial

766 variation in stress within the central Aegean Sea, Journal of Geodynamics, 39(1), 61-77,

767 https://doi.org/10.1016/j.jog.2004.09.001

768 Ganas, A., Chousianitis, K., Batsi, E., Kolligri, M., Agalos, A., Chouliaras, G. & Makropoulos, K.,

769 2013. The January 2010 Efpalion earthquakes (Gulf of Corinth, Central Greece): earthquake

770 interactions and blind normal faulting. Journal of Seismology, 17, 465–484(2013),

771 https://doi.org/10.1007/s10950-012-9331-6

772 Ganas, A., Elias, P., Bozionelos, G., Papathanassiou, G., Avallone, A., Papastergios, A., Valkaniotis,

773 S., Parcharidis, I., Briole, P., 2016. Coseismic deformation, field observations and seismic fault

774 of the 17 November 2015 M = 6.5, Lefkada Island, Greece earthquake, Tectonophysics, 687,

775 210-222, https://doi.org/10.1016/j.tecto.2016.08.012

776 Ganas, A., Andritsou, N., Kosma, C., Argyrakis, P., Tsironi, V. & Drakatos, G., 2018a. A 20-yr

31 777 database (1997–2017) of coseismic displacements from GPS recordings in the Aegean area and

778 their scaling with Mw and hypocentral distance. Bull. Geol. Soc. Greece, 52, 98–130,

779 http://doi:10.12681/bgsg.18070

780 Ganas, A., Kourkouli, P., Briole, P., Moshou, A., Elias, P. & Parcharidis, I., 2018b. Coseismic

781 Displacements from Moderate-Size Earthquakes Mapped by Sentinel-1 Differential

782 Interferometry: The Case of February 2017 Gulpinar Earthquake Sequence (Biga Peninsula,

783 Turkey), Remote Sensing, 10(7), 1089, https://doi.org/10.3390/rs10071089

784 Ganas, A., Elias, P., Kapetanidis, V., Valkaniotis, S., Briole, P., Kassaras, I., Argyrakis, P.,

785 Barberopoulou, A. & Moshou, A., 2019. The July 20, 2017 M6.6 Kos Earthquake: Seismic and

786 Geodetic Evidence for an Active North-Dipping Normal Fault at the Western End of the Gulf of

787 Gökova (SE Aegean Sea), Pure and Applied Geophysics, 176(10), 4177-4211,

788 https://doi.org/10.1007/s00024-019-02154-y

789 Ganas A., Briole, P., Bozionelos, G., Barberopoulou, A., Elias, P., Tsironi, V., Valkaniotis, S.,

790 Moshou, A. & Mintourakis, I., 2020. The 25 October 2018 Mw = 6.7 Zakynthos earthquake

791 (Ionian Sea, Greece): a low-angle fault model based on GNSS data, relocated seismicity, small

792 tsunami and implications for the seismic hazard in the west Hellenic Arc, Journal of

793 Geodynamics, 137, https://doi.org/10.1016/j.jog.2020.101731

794 Halpaap, F., Rondenay, S. & Ottemoller, L., 2018. Seismicity, deformation, and metamorphism in

795 the western Hellenic subduction zone: new constraints from tomography. Journal of Geophysical

796 Research-Solid Earth, 123(4), 3000-3026, https://doi.org/10.1002/2017JB015154

797 Hatzfeld, D., Kementzetzidou, D., Karakostas, V., Ziazia, M., Nothard, S., Diagourtas, D.,

798 Deschamps, A., Karakaisis, G., Papadimitriou, P., Scordilis, M., Smith, R., Voulgaris, N., Kiratzi,

799 A., Makropoulos, K., Bouin, M.P. & Bernard, P., 1996. The Galaxidi earthquake of 18 November

800 1992: A possible asperity within the normal fault system of the Gulf of Corinth (Greece),

801 Bulletin of the Seismological Society of America , 86(6), 1987–1991

32 802 Hollenstein, C., Muller, M., Geiger, A. & Kahle, H.G., 2008. GPS-Derived Coseismic

803 Displacements Associated with the 2001 Skyros and 2003 Lefkada Earthquakes in Greece,

804 Bulletin of the Seismological Society of America, 98(1), 149-161,

805 https://doi.org/10.1785/0120060259

806 Howell, A., Palamartchouk, K., Papanikolaou, X., Paradissis, D., Raptakis, C., Copley, A., England,

807 P. & Jackson, J., 2017. The 2008 Methoni earthquake sequence: the relationship between the

808 earthquake cycle on the subduction interface and coastal uplift in SW Greece, Geophysical

809 Journal International, 208(3), 1592-1610, https://doi.org/10.1093/gji/ggw462

810 Ilieva, M., Briole, P., Ganas, A., Dimitrov, D., Elias, P., Mouratidis, A. & Charara, R., 2016. Fault

811 plane modelling of the 2003 August 14 Lefkada Island (Greece) earthquake based on the analysis

812 of ENVISAT SAR interferograms, Tectonophysics, 693, 47-65,

813 https://doi.org/10.1016/j.tecto.2016.10.021

814 Kapetanidis, V., Karakonstantis, A., Papadimitriou, P., Pavlou, K., Spingos, I., Kaviris, G. &

815 Voulgaris, N., 2020. The 19 July 2019 earthquake in Athens, Greece: a delayed major aftershock

816 of the 1999 Mw = 6.0 event, or the activation of a different structure?, Journal of Geodynamics,

817 https://doi.org/10.1016/j.jog.2020.101766

818 Karasözen, E., Nissen, E., Büyükakpınar, P., Cambaz, M.D., Kahraman, M., Kalkan Ertan, E.,

819 Abgarmi, B., Bergman, E., Ghods, A. & Özacar, A.A., 2018. The 2017 July 20 Mw 6.6 Bodrum–

820 Kos earthquake illuminates active faulting in the Gulf of Gökova, SW Turkey, Geophysical

821 Journal International, 214(1), 185-199, https://doi.org/10.1093/gji/ggy114

822 Kiratzi, A., Tsakiroudi, E., Benetatos, C. & Karakaisis, G., 2016. The 24 May 2014 M w 6.8

823 earthquake (North Aegean Trough): spatiotemporal evolution, source and slip model from

824 teleseismic data, Physics and Chemistry of the Earth, 95, 85-100,

825 https://doi.org/10.1016/j.pce.2016.08.003

826 Kiratzi, A., 2018. The 12 June 2017 Mw 6.3 Lesvos Island (Aegean Sea) earthquake: Slip model and

33 827 directivity estimated with finite-fault inversion, Tectonophysics, 724, 1-10,

828 https://doi.org/10.1016/j.tecto.2018.01.003

829 Konca, A.O., Cetin, S., Karabulut, H., Reilinger, R., Dogan, U., Ergintav, S., Cakir, Z. & Tari, E.,

830 2018. The 2014, Mw 6.9 North Aegean earthquake: seismic and geodetic evidence for coseismic

831 slip on persistent asperities, Geophysical Journal International, 213(2), 1113-1120,

832 https://doi.org/10.1093/gji/ggy049

833 Konca, A.O., Guvercin, S.E., Ozarpaci, S., Ozdemir, A., Funning, G.J., Dogan, U., Ergintav, S.,

834 Floyd, M., Karabulut, H. & Reilinger, R., 2019. Slip distribution of the 2017 Mw 6.6 Bodrum-

835 Kos earthquake: resolving the ambiguity of fault geometry, Geophysical Journal International,

836 219(2), 911-923, https://doi.org/10.1093/gji/ggz332

837 Kontoes, C. Elias P., Sykioti O., Briole P., Remy D., Sachpazi, M., Veis, G. & Kotsis, I., 2000.

838 Displacement field and fault model for the September 7, 1999 Athens Earthquake inferred from

839 ERS2 Satellite radar interferometry, Geophysical Research Letters, 27(24), 3989-3992,

840 https://doi.org/10.1029/2000GL008510

841 Konstantinou, K.I., 2017. Accurate relocation of seismicity along the North Aegean Trough and its

842 relation to active tectonics, Tectonophysics, 717, 372-382,

843 https://doi.org/10.1016/j.tecto.2017.08.021

844 Kouskouna, V., Ganas, A., Kleanthi, M., Kassaras, I., Sakellariou, N., Sakkas, G., Tsironi, V.,

845 Karamitros, I., Tavoularis, N., Valkaniotis, S. & Bossu, R., 2019. The Athens 19/07/2019 Mw 5.1

846 earthquake, preliminary macroseismic report, EMSC technical report,

847 https://doi.org/10.13140/RG.2.2.32971.34089

848 Kostoglou, A., Karakostas, V., Bountzis, P. & Papadimitriou, E, 2020. Τhe February – March 2019

849 seismic swarm offshore North Lefkada Island, Greece: microseismicity analysis and geodynamic

850 implications, Applied Sciences, 10(13), https://doi.org/10.3390/app10134491

851 Le Pichon, X. & Kreemer, C., 2010. The Miocene-to-Present Kinematic Evolution of the Eastern

34 852 Mediterranean and Middle East and Its Implications for Dynamics, Annual Review of Earth and

853 Planetary Sciences, 38, 323-351, https://doi.org/10.1146/annurev-earth-040809-152419

854 Mantoglou, A., 2003. Pumping management of coastal aquifers using analytical models of saltwater

855 intrusion, Water Resources Research, 39(12), 1335, https://doi.org/10.1029/2002WR001891

856 Melgar, D., Ganas, A., Geng, J., Liang, C., Fielding, E. & Kassaras, I., 2017. Source characteristics

857 of the 2015 Mw 6.5 Lefkada, Greece, strike-slip earthquake. Journal of Geophysical Research-

858 Solid Earth, 122(3), 2260-2273, http://doi.org/10.1002/2016jb013452

859 Mogi, K., 1958. Relations between the eruptions of various volcanoes andthe deformations of the

860 ground surface around them, Bull. Earthquake Res. Inst. Univ. Tokyo, 36, 99– 134.

861 Müller, M.D, 2011. Analysis of long-term GPS observations in Greece (1993-2009) and

862 geodynamic implications for the Eastern Mediterranean, PhD Thesis, 19796, ETH Zürich,

863 https://doi.org/10.3929/ethz-a-006740688

864 Okada Y., 1992. Internal deformation due to shear and tensile faults in a half space, Bull. Seism.

865 Soc. America, 82, 1018–1040.

866 Papadimitriou, E., Karakostas, V., Mesimeri, M. & Vallianatos, F., 2016. The Mw 6.7 12 October

867 2013 western Hellenic arc main shock and its aftershock sequence: implications for the slab

868 properties, International Journal of Earth Sciences, 105(7), 2149-2160,

869 https://doi.org/10.1007/s00531-016-1294-3

870 Papadimitriou, P., Kassaras, I., Kaviris, G., Tselentis, G.A., Voulgaris, N., Lekkas, E., Chouliaras,

871 G., Evangelidis, C., Pavlou, K., Kapetanidis, V., Karakonstantis, A., Kazantzidou-Firtinidou, D.,

872 Fountoulakis, I., Millas, C., Spingos, I., Aspiotis, T., Moumoulidou, A., Skourtsos, E., Antoniou,

873 V., Andreadakis, E., Mavroulis, S. & Kleanthi, M., 2018. The 12th June 2017 Mw = 6.3 Lesvos

874 earthquake from detailed seismological observations, Journal of Geodynamics, 115, 23-42,

875 https://doi.org/10.1016/j.jog.2018.01.009

876 Parks, M.M., Moore, J.D.P., Papanikolaou, X., Biggs, J., Mather, T.A., Pyle, D.M., Raptakis, C.,

35 877 Paradissis, D., Hooper, A., Parsons, B. & Nomikou, P., 2015. From quiescence to unrest: 20

878 years of satellite geodetic measurements at Santorini volcano, Greece, Journal of Geophysical

879 Research-Solid Earth, 120(2), 1309-1328, https://doi.org/10.1002/2014JB011540

880 Raucoules, D., Parcharidis, I., Feurer, D., Novalli, F., Ferretti, A., Carnec, C., Lagios, E., Sakkas, V.,

881 Le Mouelic, S., Cooksley, G. & Hosford, S., 2008. Ground deformation detection of the greater

882 area of Thessaloniki (Northern Greece) using radar interferometry techniques, Natural Hazards

883 and Earth System Sciences, 8(4), 779-788, https://doi.org/10.5194/nhess-8-779-2008

884 Saltogianni, V., Gianniou, M., Taymaz, T., Yolsal-Cevikbilen, S. & Stiros, S., 2015. Fault slip

885 source models for the 2014 Mw 6.9 Samothraki-Gokceada earthquake (North Aegean Trough)

886 combining geodetic and seismological observations, Journal of Geophysical Research-Solid

887 Earth, 120(12), 8610-8622, https://doi.org/10.1002/2015JB012052

888 Serpetsidaki, A., Elias, P., Ilieva, M., Bernard, P., Briole, P., Deschamps, A., Lambotte, S., Lyon-

889 Caen, H., Sokos, E. & Tselentis, G.A., 2014. New constraints from seismology and geodesy on

890 the Mw = 6.4 2008 Movri (Greece) earthquake: evidence for a growing strike-slip fault system,

891 Geophysical Journal International, 198(3), 1373–1386, https://doi.org/10.1093/gji/ggu212

892 Schmittbuhl, J., Karabulut, H., Lengline, O. & Bouchon, M., 2016. Seismicity distribution and

893 locking depth along the Main Marmara Fault, Turkey, Geochemistry Geophysics Geosystems,

894 17(3), 954-965, https://doi.org/10.1002/2015GC006120

895 Sokos, E., Gallovic, F., Evangelidis, C.P., Serpetsidaki, A., Plicka, V., Kostelecky, J. & Zahradnik,

896 J., 2020. The 2018 Mw 6.8 Zakynthos, Greece, Earthquake: Dominant Strike-Slip Faulting near

897 Subducting Slab, Seismological Research Letters, 91(2), 721-732,

898 https://doi.org/10.1785/0220190169

899 Thatcher, W. & Rundle, J.B., 1984. A viscoelastic coupling model for the cyclic deformation due to

900 periodically repeated earthquakes at subduction zones, Journal of Geophysical Research-Solid

901 Earth, 89(NB9), 7631-7640, https://doi.org/10.1029/JB089iB09p07631

36 902 Uchida, N., Kalafat, D., Pinar, A. &Yamamoto, Y., 2019. Repeating earthquakes and interplate

903 coupling along the western part of the North Anatolian Fault, Tectonophysics, 769,

904 https://doi.org/10.1016/j.tecto.2019.228185

905 Valkaniotis, S., Briole, P., Ganas, A., Elias, P., Tsironi, V., Fokaefs, A., Partheniou, H. & Paschos, P.,

906 2020. The Mw = 5.7 Kanallaki earthquake of March 21, 2020 in Epirus, Greece: reverse fault

907 model and seismotectonic implications for Apulia-Eurasia collision, Geosciences, 10, 454,

908 https:/doi.org/10.3390/geosciences10110454

909 Vallée, M. & Douet, V., 2016. A new database of Source Time Functions (STFs) extracted from the

910 SCARDEC method, Physics of the Earth and Planetary Interior, 257, 149-157,

911 https://doi.org/10.1016/j.pepi.2016.05.012

912

913

914 TABLES

915

916 Table S1. Coordinates and velocities in the ITRF2014 (Altamimi et al. 2016) of the 329 points.

917

Code Long. Lat. Name Owner First Last vE vN vU σE σN σU

epoch epoch (mm yr-1) (mm yr-1) ABEL 21.215 38.923 Abelaki CNRS 1/11/15 24/1/20 15.1 8.4 -1.4 0.4 0.4 1.5 AFY0 23.434 40.101 Afytos Metrica 18/12/15 20/2/20 22.3 5.1 -1.2 0.4 0.4 0.9 AFYT 23.435 40.097 Afytos Metrica 17/7/12 16/12/15 23.1 6.1 -0.9 0.4 0.4 0.9 AGNI 25.718 35.192 Agios Nikolaos Uranus 28/11/13 25/2/20 8.4 -13.6 0.4 0.7 0.5 1.2 AGPA 26.269 39.246 Agia Paraskevi Uranus 28/11/13 7/3/20 5.6 1 0.2 0.4 0.4 1.0 AGRI 21.409 38.624 Agrinion Metrica 28/9/10 12/2/19 13.3 8.3 -1.1 0.3 0.3 0.6 AGRU 21.408 38.623 Agrinion Uranus 28/11/13 13/2/20 13.3 8.3 -0.5 0.3 0.4 0.9 AGTH 23.123 37.939 Agii Theodori Uranus 28/11/13 25/2/20 8.5 -10.3 0.2 0.8 1.0 2.1 AIGI 22.073 38.242 Aigion CNRS 11/12/15 10/2/20 10 -9.4 3.5 0.3 0.4 1.1 AIGU 23.440 37.734 Aigina Uranus 24/11/13 23/2/20 8.6 -10.2 1.4 0.3 0.4 0.9 AKYR 24.913 34.981 Antikythira COMET 10/3/04 8/2/08 8.7 -14.1 -0.1 0.4 0.5 1.2 ALE1 25.873 40.846 Alexandroupoli Uranus 15/12/13 14/6/16 23 13.2 0.4 0.6 0.8 1.7 ALEU 25.906 40.854 Alexandroupoli Uranus 17/6/16 14/3/20 22.6 11.7 0 0.4 0.4 1.1 ALEX 25.853 40.849 Alexandroupoli Metrica 1/9/11 8/3/20 22.8 12.3 0.7 0.2 0.3 0.6 ALXA 22.444 40.628 Alexandreia Uranus 28/11/13 15/3/20 24.6 8.9 -6.9 0.3 0.3 0.8 AMAL 21.355 37.796 Amaliada Uranus 28/11/13 12/11/19 7.8 -6.9 3.6 0.3 0.4 0.8 AMFI 21.150 38.892 Amfilochia CNRS 28/9/15 10/2/20 15.7 9 3.6 0.4 0.5 2.9

37 ANAV 23.904 37.733 Anavyssos Metrica 18/1/07 8/5/12 8.6 -11.4 0.3 0.4 0.4 1.0 ANDR 24.737 37.886 Andros Metrica 26/1/11 8/3/20 8.2 -11.8 -1.1 0.4 0.6 1.2 ANIK 22.282 36.826 Agios Nikolaos NOA 4/7/16 23/12/19 4.7 -12.8 1 0.6 0.6 1.8 ANIS 23.760 35.332 Anisaraki Uranus 25/11/13 9/3/20 8.2 -10.9 -0.1 0.5 0.5 1.5 ANKY 23.301 35.867 Antikythira INGV 20/6/17 21/2/20 7.5 -12.6 -1.2 0.6 0.5 1.0 ANOA 21.590 38.138 Ano Alissos CNRS 15/1/15 13/2/20 10.1 -7.3 -0.9 0.5 0.5 1.1 ANOC 21.925 38.594 Ano Chora CNRS 3/10/14 31/12/19 12.7 4.5 -4.5 0.3 0.4 1.2 ANOP 24.096 35.218 Anopoli COMET 25/7/03 8/2/08 7.8 -11.1 1.9 0.5 0.5 1.2 ANTL 30.666 36.888 Antalya Turkey 16/1/09 31/10/15 11.7 8.2 -1.1 0.2 0.2 0.6 ARGO 20.492 38.169 Argostoli Metrica 22/3/16 8/3/20 16.5 4.8 -0.9 0.4 0.4 1.0 ARGU 22.693 37.594 Argos Uranus 28/11/13 13/2/20 7.6 -12.5 0.6 0.4 0.4 1.0 ARKI 23.034 38.755 Arkitsa NTUA 20/6/01 4/11/06 12.7 -5.1 -0.5 0.6 0.5 1.4 ARKL 25.269 35.134 Arkalochori Metrica 30/9/14 8/3/20 8.5 -13 -0.7 0.5 0.4 0.9 ARSA 21.817 38.301 Arsacheio CNRS 12/6/14 10/2/20 8.4 -3.4 -0.2 0.4 0.5 0.9 ART1 20.933 39.153 Arta Metrica 24/8/11 12/3/13 17.3 9.1 -1.4 0.6 0.6 1.3 ARTA 20.992 39.159 Arta Metrica 22/4/13 12/2/20 18.7 10.8 -1.2 0.3 0.3 0.9 ARTU 21.000 39.157 Arta Uranus 7/8/15 19/2/20 17.1 9.9 -2.1 0.4 0.4 1.0 ASPR 22.631 40.899 Aspro Uranus 28/11/13 9/3/20 24.6 8.5 -0.2 0.3 0.3 0.8 ASSO 20.542 38.378 Assos CNRS 11/12/15 18/2/20 17.8 4.6 1.1 0.5 0.5 1.2 ASTU 26.349 36.549 Astypalea Uranus 19/7/15 7/3/20 10.6 -14.8 -0.7 0.4 0.5 1.2 ASTY 26.353 36.545 Astypalea Metrica 19/10/12 9/3/20 10.6 -14.6 0.5 0.3 0.3 0.7 ATAL 22.999 38.653 Atalanti NOA 28/3/09 3/3/20 12.8 -5.2 -0.6 0.2 0.3 0.7 ATHI 23.699 37.935 Athina Uranus 28/11/13 23/2/20 8.9 -11 -0.1 0.3 0.4 1.0 ATRS 22.835 37.244 Ano Tiros COMET 17/5/05 13/2/16 7.1 -12.9 -0.3 0.3 0.4 1.1 AUT1 23.004 40.567 Thessaloniki AUTH 31/3/05 11/3/20 24.9 7.5 -1.2 0.2 0.2 0.5 AYD1 27.838 37.841 Aydin Turkey 16/1/09 31/10/15 4.1 -3.4 -1.3 0.2 0.2 0.8 AYVL 26.686 39.311 Ayvalik Turkey 16/1/09 31/10/15 4.9 1.5 -0.6 0.2 0.2 0.5 BAJR 20.074 42.358 Bajram AlbPos 28/11/09 7/6/19 23.1 15 -0.8 0.4 0.4 1.1 BALK 27.894 39.639 Balikesir Turkey 16/1/09 31/10/15 3.3 6.8 -0.6 0.2 0.2 0.6 BAND 27.997 40.331 Bandirma Turkey 16/1/09 9/1/14 4.8 7.9 2.9 0.2 0.3 0.9 BELE 23.133 39.295 Koropi Metrica 17/10/17 15/3/20 18.8 0.3 -1.1 0.5 0.6 1.3 BERA 19.956 40.703 Berat AlbPos 28/11/09 20/11/19 22.6 14.1 -0.7 0.2 0.2 0.7 BURS 29.015 40.214 Bursa Turkey 16/1/09 31/10/15 2.5 9.5 -1.7 0.2 0.3 0.7 CANA 26.414 40.111 Canakkale Turkey 25/1/09 31/10/15 6.6 4.9 0.5 0.2 0.2 0.6 CESM 26.373 38.304 Cesme Turkey 16/1/09 31/10/15 5.6 -8.1 -1.4 0.2 0.2 0.6 CHA1 24.018 35.511 Chania Uranus 22/7/05 5/8/18 9.1 -11.1 0.5 0.4 0.5 1.2 CHA2 24.035 35.483 Chania Uranus 7/3/18 18/3/20 10.7 -4.4 12.1 0.5 1.2 1.9 CHIO 26.127 38.368 Chios Metrica 5/9/12 9/3/20 6.1 -8.8 -0.9 0.3 0.3 0.7 CHI1 26.142 38.365 Chios Uranus 22/7/15 23/9/16 5.9 -8.4 5.6 0.8 0.7 1.4 CHIU 26.136 38.367 Chios Uranus 25/9/16 8/3/20 6.2 -9.5 -1.6 0.5 0.4 1.0 DATC 27.692 36.709 Datca Turkey 1/1/09 31/10/15 11.7 -13.3 -4.2 0.2 0.2 0.9 DEIR 28.648 39.035 Demirci Turkey 16/1/09 16/6/15 1.5 5.3 -0.3 0.2 0.2 0.6 DIDI 27.269 37.372 Didim Turkey 16/1/09 31/10/15 6.4 -10.9 -1.8 0.2 0.2 0.5 DION 23.933 38.079 Dionyssos NTUA 30/9/97 5/10/06 9 -10.8 0.2 0.5 0.7 1.8 DNZL 29.092 37.762 Denizli Turkey 16/1/09 31/10/15 4.1 1.6 3.3 0.2 0.4 0.9 DRAM 24.069 41.166 Drama Metrica 30/9/14 9/3/20 23.7 12.3 -0.9 0.3 0.6 0.8 DRAN 20.575 38.684 Dragano CNRS 23/11/15 17/6/18 17.3 8.6 0.8 0.6 0.5 1.4 DRAU 24.137 41.155 Drama Uranus 22/7/15 7/3/20 23.3 12 0.4 0.3 0.4 0.9 DSLN 25.389 36.461 Santorini NTUA 15/1/12 13/1/15 6 -19.8 -0.7 0.4 0.5 1.1 DUTH 24.917 41.140 Xanthi AUTH 18/9/08 11/3/20 23.6 11.9 0.3 0.2 0.3 0.6 DYNG 23.932 38.079 Dionissos CNES 11/5/11 31/12/19 8.9 -10.7 0.3 0.3 0.4 0.8 EDES 22.051 40.804 Edessa Metrica 12/6/11 20/2/20 25.6 7.2 -1.9 0.3 0.3 1.4 EGIO 22.075 38.255 Aigion NOA 6/7/16 21/2/20 10.5 -7.6 -0.8 0.6 0.9 1.7 ELAS 22.206 39.892 Elassona Metrica 13/4/18 9/3/20 22.7 6.1 1.3 0.5 0.6 1.4 EPID 23.157 37.638 Epidauros Uranus 19/4/16 13/2/20 8.2 -12.6 -1.7 0.4 0.4 1.0 EXAN 20.619 38.754 Exanthia CNRS 25/11/15 2/5/20 19.4 5.1 -0.7 0.3 0.3 1.5

38 EYPA 21.928 38.427 Efpalion CNRS 14/10/01 31/10/19 11.4 2.7 -1.4 0.2 0.2 0.6 FARS 22.384 39.296 Farsala Uranus 17/7/15 1/12/18 18.7 2.1 0.9 0.4 0.5 1.1 FETH 29.124 36.626 Fethiye Turkey 16/1/09 17/10/15 11 -2.5 -2.9 0.2 0.3 0.6 FINI 30.146 36.302 Finike Turkey 16/1/09 31/10/15 14.8 8 -1.1 0.2 0.3 0.7 FISK 20.577 38.460 Fiskardo CNRS 25/11/15 18/2/20 17.4 3.3 -0.8 0.4 0.4 1.0 FLOR 21.402 40.780 Florina Metrica 2/5/14 12/2/20 24 9.6 -0.1 0.4 0.4 1.3 FLOU 21.456 40.804 Florina Uranus 22/7/15 8/3/20 24.3 10.3 0.1 0.3 0.4 1.0 FREN 33.911 35.037 Frenaros Uranus 22/7/15 26/2/20 17.5 19.9 -1.1 0.3 0.4 0.9 GALA 22.392 38.375 Galaxidi CNRS 30/9/14 31/12/19 10.8 -4 -2.3 0.4 0.5 1.1 GARD 20.571 39.351 Gardiki Uranus 22/7/15 19/2/20 21.1 13.2 -0.8 0.4 0.5 1.1 GARG 21.640 37.066 Gargalianoi Uranus 22/7/15 31/12/19 2.2 -12.1 0.2 0.4 0.4 1.0 GEYB 21.766 38.334 Gefyra Building CNRS 28/9/11 10/2/20 11.9 2.9 -3.8 0.3 0.4 0.7 GODE 23.734 41.566 Goce Delchev BAS 2/5/14 13/9/17 24.9 10.4 2.9 0.4 0.5 1.2 GOUM 22.447 40.947 Gumenissa Metrica 19/5/15 13/2/20 24.3 8.9 0.7 0.3 0.4 0.9 GREU 21.431 40.085 Grevena Uranus 22/7/15 17/7/19 23.3 9.2 -0.7 0.4 0.5 1.3 GREV 21.438 40.085 Grevena Metrica 18/7/13 3/2/20 23.9 9.8 0.4 0.2 0.3 0.8 GYTH 22.566 36.761 Gythio Uranus 22/7/15 1/6/17 3.9 -11.3 -1.3 0.8 0.8 2.1 HALK 23.603 38.465 Halkida Metrica 7/9/10 9/3/20 10.6 -8.5 0 0.2 0.3 0.6 HARC 29.153 39.678 Harmancik Turkey 16/1/09 31/10/15 1.9 9.3 -1 0.3 0.2 0.5 HERA 25.142 35.324 Heraklion Metrica 7/9/10 10/3/20 8.7 -13 -1.6 0.2 0.3 0.6 HGOU 20.268 39.493 Igoumenitsa Uranus 22/7/15 22/2/20 21.3 14.2 -1.7 0.4 0.4 0.9 HRA1 25.191 35.314 Iraklio Uranus 1/1/16 19/2/18 9.3 -9.6 0.6 0.4 0.6 1.1 IDI0 24.890 35.289 Anogia TEIC 20/2/13 31/1/20 8.7 -12.3 -0.5 0.5 0.5 1.4 IERA 25.797 35.053 Ierapetra Metrica 30/6/10 9/3/20 9 -13.6 -0.5 0.4 0.4 1.1 IERI 23.879 40.400 Ierissos Uranus 1/1/16 9/3/20 23.7 7.7 -0.2 0.3 0.4 0.8 IGOU 20.263 39.510 Igoumenitsa Metrica 2/5/14 13/2/20 21.2 13.6 0.1 0.4 0.3 1.0 IKAR 26.224 37.628 Ikaria Metrica 28/6/18 10/3/20 6.3 -11.2 1.4 0.7 0.8 1.9 IKAU 26.273 37.605 Ikaria Uranus 22/7/15 24/2/20 7.2 -11.7 -1.8 0.6 0.7 1.8 IOAN 20.851 39.664 Ioannina Metrica 21/12/10 13/2/20 21.7 11 -0.5 0.3 0.3 0.7 IOAU 20.854 39.655 Ioannina Uranus 22/7/15 20/2/20 22 10.2 0.4 0.5 0.4 1.0 IPSA 26.380 40.918 Ipsala Turkey 16/1/09 31/5/15 22.1 13 0.2 0.2 0.2 0.6 ISTA 29.019 41.104 Istanbul BKG 11/12/99 21/3/20 22.6 13.7 -0.2 0.1 0.1 0.4 ISTI 23.152 38.956 Istiaia Uranus 22/3/13 21/2/20 15.6 -1.8 -0.8 0.3 0.3 0.7 ISTN 28.832 40.991 Istanbul Turkey 16/1/09 31/10/15 22.5 14.4 -0.9 0.2 0.2 0.6 ITEA 22.430 38.431 Itea Metrica 6/10/11 14/2/20 10.9 -2 -1.8 0.3 0.3 0.9 ITEU 22.427 38.434 Itea Uranus 22/7/15 13/2/20 10.6 -2.1 -1.3 0.4 0.5 1.1 IZMI 27.082 38.395 Izmir Turkey 10/8/08 23/1/20 3.1 -3.3 -0.3 0.2 0.2 0.5 KAL1 22.857 40.642 Kalohori Metrica 22/11/13 13/2/20 24.7 6.9 -1.8 0.3 0.4 1.2 KAL2 22.848 40.643 Kalohori Metrica 26/11/13 14/2/20 23.9 7.2 -3.1 0.3 0.3 0.9 KALA 22.102 38.031 Kalavrita CNRS 10/10/14 10/2/20 8.6 -10 -0.7 0.3 0.4 1.2 KALM 22.111 37.034 Kalamata Metrica 3/11/10 6/2/20 4.4 -12.1 -0.4 0.3 0.3 0.7 KALU 26.962 36.962 Kalymnos Uranus 22/7/15 8/2/20 9.7 -11.9 -0.7 0.5 0.5 1.1 KALY 26.976 36.956 Kalymnos Metrica 24/1/13 10/3/20 9 -12.8 0.1 0.3 0.3 0.9 KARB 28.683 41.347 Karaburun Turkey 16/1/09 31/10/15 23.6 13.9 -0.2 0.2 0.2 0.6 KARD 22.082 39.470 Karditsa Uranus 14/9/17 28/2/20 20.4 6.1 -1.2 0.4 0.5 1.1 KARP 21.803 38.912 Karpenissi Metrica 9/6/11 14/2/20 16.2 5.8 -0.8 0.3 0.3 0.9 KARU 27.211 35.508 Karpathos Uranus 22/7/15 27/2/20 13.6 -16.2 -0.3 0.6 0.4 1.2 KARY 24.419 38.016 Karystos Uranus 22/7/15 9/3/20 8.9 -11.2 -0.1 0.5 0.7 1.2 KASI 19.936 39.746 Kassiopi NOA 1/4/07 6/5/17 20.3 14 -0.7 0.3 0.3 0.6 KAST 21.263 40.523 Kastoria Metrica 24/8/11 14/2/20 24.3 10.8 -1.5 0.2 0.3 0.7 KATC 27.781 35.951 Katavia NOA 31/8/12 16/12/19 17.6 -14.2 -1 0.5 0.4 1.1 KATE 22.521 40.270 Katerini Metrica 18/12/09 18/2/20 24.5 6.5 -2 0.2 0.3 0.7 KATU 22.507 40.272 Katerini Uranus 22/7/15 27/2/20 24 6.5 0.7 0.4 0.4 1.1 KAVA 24.399 40.937 Kavala Uranus 22/7/15 9/3/20 23.8 11.3 0.2 0.3 0.4 1.1 KEFA 20.438 38.196 Kefalonia Uranus 22/7/15 12/3/20 17.4 7.4 -0.7 0.4 0.4 1.0 KERA 25.348 36.418 Therasia Unavco 5/5/06 19/1/16 9.6 -14.1 0 0.4 0.3 0.8

39 KERK 19.873 39.494 Kerkyra Metrica 9/12/15 18/2/20 21.1 14.7 -0.7 0.3 0.3 1.1 KERT 23.979 37.806 Keratea Uranus 22/7/15 9/3/20 8.2 -11.1 -0.9 0.4 0.5 1.1 KERU 19.915 39.619 Kerkyra Uranus 22/7/15 20/2/20 22.1 14.5 -0.4 0.4 0.4 0.9 KERY 22.384 36.493 Kerya COMET 14/3/04 7/2/08 6.1 -12.3 0.3 0.4 0.4 1.2 KIAT 22.750 38.014 Kiato Uranus 22/7/15 13/2/20 8.8 -11.2 -1.7 0.4 0.5 1.1 KIKA 27.672 39.106 Kirkagac Turkey 16/1/09 31/10/15 2.8 2.8 -0.1 0.2 0.2 0.7 KILK 22.878 40.995 Kilkis Uranus 22/7/05 4/9/19 23.8 9.5 -0.5 0.3 0.4 1.0 KIPO 20.348 38.203 Kipouria NOA 6/9/10 23/8/18 17.8 9.8 1.5 0.5 0.5 1.0 KIRL 27.218 41.738 Kirlareli Turkey 16/1/09 31/10/15 23.5 13.7 -0.4 0.2 0.2 0.6 KISM 23.654 35.495 Kissamos Metrica 15/2/19 10/3/20 9.5 -11.8 1.6 0.9 0.9 2.2 KITH 23.015 36.275 Kythira COMET 13/3/04 8/2/08 7 -11.8 3.2 0.4 0.4 1.1 KLAM 22.111 37.030 Kalamata Uranus 26/7/15 31/12/19 4.4 -12.1 -1.1 0.5 0.5 1.0 KLOK 22.014 39.565 Klokotos INGV 17/7/08 4/3/20 21.1 6.4 -0.3 0.2 0.2 0.6 KOM1 25.396 41.122 Komotini Metrica 22/7/16 10/3/20 22.9 12.1 0 0.3 0.4 1.0 KOMO 25.409 41.120 Komotini Uranus 22/7/15 7/3/20 23.6 12.2 0.2 0.4 0.6 1.2 KOPA 21.818 37.290 Kopanaki Metrica 2/5/14 9/2/20 4.1 -11.9 0.2 0.3 0.4 1.0 KORI 22.931 37.942 Korinthos Metrica 30/7/08 21/2/20 8.3 -11.6 0.5 0.2 0.2 0.6 KORO 21.956 36.798 Koroni NOA 20/4/18 5/2/20 3.5 -11.1 1.4 0.6 0.7 1.5 KORU 21.913 36.810 Koroni Uranus 22/7/15 13/2/20 4 -11.2 -0.6 0.4 0.4 0.9 KOUN 22.046 38.209 Kounina CNRS 23/5/02 10/2/20 10.1 -10.1 -1 0.2 0.2 0.7 KOZA 21.785 40.299 Kozani Uranus 22/7/15 22/2/20 23.7 9 0.8 0.3 0.4 0.9 KRDI 21.923 39.366 Karditsa Metrica 8/6/12 19/2/20 19.9 6.7 -1.2 0.2 0.3 0.8 KRIN 21.960 38.189 Krini CNRS 4/10/14 26/9/19 111.6 -31.3 -19.7 1.3 0.6 1.3 KRP1 27.120 35.481 Karpathos Metrica 5/2/16 11/3/20 12.5 -16.3 0.9 0.5 0.4 1.0 KRPS 27.161 35.547 Karpathos Metrica 14/8/12 31/10/15 12.5 -16.5 0.8 0.6 0.5 1.4 KRUM 25.652 41.473 Kroumovgrad BAS 15/9/11 12/3/20 23.5 13.4 0.1 0.3 0.3 0.8 KRYO 22.618 37.972 Krioneri NTUA 12/5/05 30/6/16 8.4 -11.7 0.6 0.3 0.4 1.1 KTCH 21.247 38.412 Katochi CUP 9/1/14 31/12/19 14.6 -2.4 -0.6 0.3 0.3 0.9 KTHA 23.063 36.256 Kythira INGV 18/6/17 22/10/18 7.8 -12.8 -2.7 0.9 1.0 2.3 KTIM 22.120 39.134 Ktimeni Metrica 5/5/10 8/5/12 17.5 3.9 0.1 0.8 0.9 2.4 KUKE 20.417 42.077 Kukes AlbPos 28/11/09 7/6/19 24 13.7 -0.4 0.3 0.2 0.8 KUST 22.713 42.284 Kustendil BAS 3/6/07 26/10/14 23.8 12.3 0.8 0.3 0.3 1.0 KYMI 24.105 38.635 Kymi Metrica 6/11/14 18/2/20 10.9 -10 -3.2 0.4 0.4 1.1 LAMA 22.432 38.898 Lamia Metrica 15/12/11 20/2/20 16 2.1 0.6 0.2 0.3 0.7 LAMB 21.973 38.320 Lambiri CNRS 22/7/09 10/2/20 11.6 -6.1 1.7 0.3 0.4 1.0 LAMI 22.403 38.890 Lamia Uranus 22/7/15 20/2/20 16.5 2 1.9 0.4 0.5 1.1 LARI 22.400 39.637 Larissa Uranus 22/7/15 22/2/20 20.5 4.7 -0.7 0.4 0.4 0.9 LARM 22.388 39.614 Larissa Metrica 19/1/11 1/6/19 21.3 5.3 -0.3 0.2 0.3 0.7 LEFK 20.695 38.830 Lefkada Metrica 20/5/18 11/3/20 20 0.5 2.8 0.6 0.7 1.5 LEMN 25.181 39.897 Lemnos NOA 10/7/07 28/2/20 9.9 1.6 0 0.2 0.2 0.5 LEON 22.861 37.163 Leonidio Uranus 22/7/15 12/2/20 6.7 -12.7 -0.9 0.4 0.4 1.1 LEPE 21.290 38.710 Lepenou CUP 27/12/14 5/2/17 13.5 8.7 -3.8 0.6 0.7 1.5 LERO 26.855 37.136 Leros Uranus 22/7/15 8/3/20 7.5 -12.1 -0.7 0.4 0.4 1.0 LESU 26.449 39.031 Lesvos Uranus 22/7/15 7/3/20 6.3 0.4 0.1 0.4 0.5 0.9 LESV 26.554 39.100 Lesvos Metrica 14/10/14 12/3/20 5.4 0.9 -0.3 0.3 0.3 1.0 LEUK 20.708 38.830 Lefkada Uranus 22/7/15 19/2/20 19.4 1.4 -2.7 0.5 0.4 1.0 LFKD 20.712 38.835 Lefkas Metrica 3/12/14 20/2/20 20.8 2 -2 0.4 0.4 0.9 LIBR 20.314 41.179 Librazhd AlbPos 28/11/09 21/11/19 22 13.3 0 0.4 0.3 1.2 LIDO 22.201 38.529 Lidoriki CNRS 13/10/01 10/2/20 12.1 1 0.8 0.2 0.3 0.7 LIMN 25.061 39.875 Limnos Uranus 22/7/15 7/3/20 9.2 3 -2.4 0.4 0.3 0.8 LIVA 22.865 38.440 Livadia Uranus 9/9/16 18/3/20 10.9 -4.6 1.3 0.5 0.6 1.2 LOVE 24.708 43.137 Lovech BAS 3/11/09 19/2/18 23.5 13.5 -1.2 0.2 0.2 0.6 LYGO 23.038 37.613 Lygourio Metrica 8/10/14 31/12/19 7.8 -12.6 -1.4 0.3 0.4 0.9 MANT 23.480 38.797 Mantoudi Uranus 22/7/15 20/2/20 13.5 -5.1 -0.2 0.3 0.5 0.9 MAOR 21.390 38.186 Mavro Oros CRL 15/10/98 28/9/15 10.3 -5.7 0 1.1 1.0 2.5 MARM 23.780 40.094 Marmaras Metrica 31/1/17 18/2/20 22.3 5.8 -0.9 0.5 0.5 1.6

40 MARU 23.783 40.096 Marmaras Uranus 22/7/15 8/3/20 23.5 5.8 0.5 0.4 0.4 0.9 MEGI 29.593 36.151 Kastellorizo Metrica 2/5/14 31/8/16 12.1 3.3 0.9 0.6 0.7 1.6 MENA 23.386 37.561 Methana COMET 24/6/04 8/2/08 7.2 -13.1 0.2 0.3 0.4 1.2 MESA 21.586 38.485 Mesarista CNRS 1/10/14 10/2/20 16.1 7.5 -0.2 0.3 0.4 1.0 MESO 21.475 38.366 Mesologhi CNRS 2/10/14 10/2/20 13.7 2.8 -2.7 0.3 0.4 1.1 MET0 23.762 38.065 Metamorfosi Metrica 27/10/06 21/2/20 9 -11 0.3 0.2 0.2 0.6 METH 21.705 36.826 Methoni COMET 20/3/03 19/7/09 4.2 -11.5 -0.2 0.3 0.3 0.9 MILO 24.431 36.742 Milos Metrica 20/7/11 21/2/20 7.9 -12.3 -0.8 0.2 0.3 0.6 MKMN 25.401 36.408 Nea Kammeni NTUA 3/1/12 24/1/15 6.2 -14.8 -3.5 0.4 0.4 1.2 MLOS 24.519 36.747 Milos COMET 29/6/04 8/2/08 7.5 -12.3 0.3 0.3 0.4 1.1 MOIR 24.878 35.054 Mires Uranus 23/10/12 5/12/19 8.8 -13.1 -2.3 0.3 0.3 0.7 MOLA 22.851 36.806 Molaoi Metrica 8/11/10 10/2/20 6.6 -13.1 -0.3 0.3 0.3 0.9 MOLU 22.852 36.804 Molaoi Uranus 23/7/15 13/2/20 6.4 -12.9 -0.6 0.4 0.4 1.0 MOUD 25.269 39.874 Moudros Metrica 24/5/17 23/2/20 9.6 0.4 0.2 0.4 0.5 1.0 MOZI 25.422 36.456 Mouzaki Bay Unavco 28/8/11 8/12/14 4.6 -18 0.1 0.5 0.5 0.9 MTSV 21.183 39.771 Metsovo Uranus 25/9/18 5/4/20 21.5 8.8 3.8 0.7 0.8 2.6 MUGL 28.364 37.216 Mugla Turkey 16/1/09 31/10/15 6.5 -5.1 -0.2 0.3 0.3 0.6 MURG 21.554 39.739 Mourgani Metrica 30/3/16 23/2/20 21.7 7.9 -0.4 0.4 0.4 1.1 MURT 24.403 35.199 Murthios Uranus 22/7/15 8/3/20 8.3 -12.5 1.6 0.5 0.9 1.6 MYKN 25.329 37.442 Mykonos Metrica 24/5/16 23/2/20 7.4 -11.9 -1.8 0.4 0.4 1.1 NAFP 21.781 38.370 Nafpaktos UPAT 20/3/14 31/12/19 13.6 3 -0.9 0.3 0.5 1.0 NAOU 22.071 40.630 Naoussa Uranus 22/7/15 23/2/20 24.7 7.9 0.2 0.4 0.4 1.0 NAXO 25.381 37.098 Naxos Metrica 17/6/11 23/2/20 7.6 -11.8 -1 0.3 0.3 0.7 NAXU 25.377 37.102 Naxos Uranus 22/7/15 24/2/20 7.7 -12 0.4 0.4 0.4 0.9 NEAB 23.060 36.509 Neapoli NOA 28/6/12 9/2/20 7.3 -12.7 -0.4 0.4 0.3 1.0 NEAP 25.610 35.261 Neapoli COMET 7/3/02 30/10/10 7.8 -14.2 -0.9 0.4 0.4 1.0 NEAU 23.058 36.511 Neapoli Uranus 22/7/15 10/3/20 6.5 -11.7 -0.8 0.6 0.4 1.5 NEOS 23.647 41.098 Neo Souli Uranus 22/7/15 7/3/20 24 11.2 -1.3 0.4 0.5 1.0 NOA1 23.864 38.047 Penteli NOA 10/4/06 6/3/20 8.6 -10.5 0.2 0.2 0.3 0.7 NOMI 25.429 36.422 Santorini Unavco 9/5/06 19/1/16 7.3 -15.7 1.6 0.3 0.3 0.7 NVRK 23.870 41.336 Nevrokopi NOA 13/7/07 6/7/13 24 11.2 0 1.2 0.8 2.1 ORES 26.546 41.414 Orestiada Metrica 2/9/11 4/3/20 23.4 13.1 0 0.2 0.3 0.6 OREU 26.529 41.520 Orestiada Uranus 22/7/15 7/3/20 22.9 13.9 -0.3 0.3 0.4 0.9 ORID 20.794 41.127 Ohrid IGS 24/7/00 18/2/20 24.3 11.8 0.2 0.2 0.2 0.5 OROP 25.488 35.198 Oropedio Uranus 22/7/15 10/3/20 9.1 -12.9 -6.5 0.5 1.3 2.0 PALC 23.641 35.239 Paleochora Metrica 22/10/12 22/2/19 8 -11.6 1 0.4 0.6 1.0 PARO 25.154 37.086 Paros Uranus 22/7/15 10/3/20 7.8 -11.8 0.9 0.4 0.4 1.0 PAT0 21.787 38.284 Patras CNRS 29/9/06 31/12/19 8.1 -5.8 -0.3 0.3 0.3 0.6 PAT1 21.787 38.283 Patras CNRS 1/1/08 11/7/09 8.1 -5.8 -0.3 0.6 0.7 1.9 PAT2 21.764 38.266 Patra2 Uranus 1/1/17 18/3/20 5.9 -5.6 0 0.5 0.6 1.2 PAT3 21.722 38.228 Patra3 Uranus 24/8/17 18/3/20 7.4 -6.1 -6 0.5 0.6 1.5 PATR 21.733 38.241 Patras Metrica 18/1/11 24/2/20 8.7 -6.4 -0.5 0.3 0.3 0.7 PAXO 20.164 39.211 Paxoi Metrica 2/5/14 24/2/20 24.3 16.4 -0.9 0.4 0.4 0.9 PAZA 24.340 42.231 Pazardhik BAS 7/5/07 28/6/19 23.6 12.9 -0.6 0.2 0.2 0.6 PERI 24.137 38.438 Perivolia Uranus 22/7/15 20/2/20 9.8 -10.1 -0.5 0.3 0.4 0.9 PERM 20.353 40.233 Permet AlbPos 28/11/09 6/6/62018 21.7 12.9 0.5 0.3 0.4 0.9 PESH 20.426 41.685 Peshkopi AlbPos 28/11/09 30/9/18 23.4 12.7 0.2 0.2 0.2 0.9 PIRG 22.589 39.919 Pirgetos Uranus 22/7/15 23/2/20 22.4 6 0.1 0.4 0.4 1.0 PKMN 25.385 36.396 Palea Kammeni Unavco 2/7/08 22/2/16 8.6 -13.6 0.1 0.3 0.3 0.7 POLI 23.438 40.376 Poligiros Uranus 22/7/15 8/3/20 23.8 6.3 0.2 0.3 0.4 0.9 POLY 23.444 40.373 Polygiros NTUA 11/9/01 7/12/04 23.1 4 -0.1 1.5 1.2 1.7 PONT 20.585 38.619 Ponti NOA 16/2/07 21/2/20 20.8 7.9 -0.9 0.3 0.3 0.8 PRKV 26.265 39.246 Agia Paraskevi NOA 30/6/07 6/3/20 5.4 0.3 0.2 0.3 0.3 0.7 PROV 23.387 41.068 Provatas Metrica 20/2/13 25/2/20 24 10 -1.2 0.3 0.3 0.7 PSAR 22.184 38.322 Psaromita CNRS 11/10/01 10/2/20 9.4 -0.9 -0.8 0.2 0.3 0.7 PSAT 21.871 38.329 Psathopirgos CNRS 12/6/14 10/2/20 10.6 -0.8 -3.4 0.5 0.6 1.0

41 PTKG 22.490 36.421 Porto Kagio INGV 30/6/18 9/2/20 4.5 -11.9 0 0.9 0.6 1.4 PTOL 21.678 40.511 Ptolemaida Metrica 29/6/11 24/2/20 24.3 8.6 0.2 0.2 0.2 0.7 PVOG 21.523 38.614 Paravola CUP 20/11/13 31/12/19 13.3 7.2 -1.5 0.4 0.4 0.8 PYL1 21.742 36.896 Pylos CNRS 23/6/07 16/7/12 1.7 -11.4 2.6 0.8 0.5 1.0 PYLO 21.695 36.914 Pylos NOA 24/8/11 9/2/20 2.5 -11.8 -1.4 0.3 0.3 0.7 PYRG 21.462 37.679 Pyrgos Metrica 24/8/11 12/3/20 5.1 -11.3 -0.5 0.2 0.3 0.6 RAFI 24.009 38.019 Rafina Uranus 11/9/17 18/3/20 8.7 -11.6 -1.4 0.5 0.5 1.3 RETH 24.613 35.388 Rethimno Metrica 9/6/11 4/3/20 8.2 -12.5 -0.6 0.2 0.3 0.6 RETS 21.429 38.429 Kato Retsina CUP 18/11/14 19/12/19 13.5 4.4 0.5 0.4 0.5 1.0 RIBA 25.343 36.454 Port Riba Unavco 27/8/11 22/2/16 9 -15.3 0.4 0.4 0.4 0.9 RGNI 21.239 38.582 Rigani CUP 28/1/15 31/12/19 16.5 -2 0 0.4 0.5 1.2 RLSO 21.465 38.056 Riolos NOA 29/7/16 31/12/19 10.3 -7.9 0 0.2 0.2 0.6 ROD3 21.892 38.308 Ano Rodini CNRS 14/3/14 10/2/20 11.3 -8.8 0.2 0.3 0.4 0.9 RODO 28.162 36.293 Rodos Metrica 29/4/11 4/3/20 15.5 -11 0.8 0.3 0.3 0.7 ROZH 24.741 41.696 Rozhen BAS 20/6/05 4/4/17 23.4 12.3 0.5 0.2 0.2 0.6 SALH 28.124 38.483 Salihli Turkey 16/1/09 31/10/15 0.6 1.8 -6.2 0.2 0.2 0.7 SAMO 26.705 37.793 Karlovasi Metrica 28/10/11 25/2/20 6.4 -9.3 0.1 0.3 0.3 0.8 SAMU 26.973 37.758 Samos Uranus 22/7/15 8/3/20 7.1 -8 1.1 0.4 0.5 0.9 SAND 23.268 41.551 Sandanski BAS 20/5/07 17/5/12 23.4 11.3 0.1 0.3 0.3 1.1 SANT 25.423 36.434 Santorini NKUA 3/2/12 3/3/20 5.5 -15.8 0.3 0.3 0.3 0.6 SARY 27.916 41.443 Saray Turkey 16/1/09 31/10/15 23.4 13.8 -0.5 0.2 0.2 0.6 SERR 23.037 41.286 Ano Poroia Metrica 21/1/12 26/2/20 24.1 9.9 -0.5 0.2 0.4 0.9 SIFN 24.736 36.961 Sifnos Uranus 22/7/15 24/2/20 7.8 -12 -0.7 0.4 0.5 0.9 SISS 20.659 38.101 Sission NKUA 20/3/14 12/11/18 16.2 2.3 3.3 0.5 0.5 1.5 SITI 26.106 35.205 Sitia Metrica 2/5/14 27/2/20 10.2 -12.6 -5.4 0.3 0.4 0.9 SIVA 24.812 35.018 Sivas TEIC 18/5/14 6/1/17 9.7 -14 -0.9 0.5 0.6 1.3 SKOP 23.728 39.122 Skopelos Metrica 5/8/10 27/4/16 17 -6.6 -1.6 0.3 0.2 0.7 SKOU 23.731 39.119 Skopelos Uranus 22/7/15 4/8/19 14.8 -5.4 -1.7 0.4 0.5 0.9 SKP1 23.733 39.117 Skopelos Metrica 28/4/16 27/2/20 14.3 -5.6 -0.2 0.5 0.5 1.1 SKYR 24.565 38.904 Skyros NOA 29/12/12 6/3/20 11.4 -9.2 0.3 0.3 0.3 0.7 SMLN 24.712 41.575 Smolyan BAS 18/9/17 3/3/20 23.6 11.9 0.5 0.5 0.6 1.3 SNDA 23.277 41.561 Sandanski BAS 18/9/17 3/3/20 23.4 10.2 0.7 0.4 0.5 1.2 SNTR 25.357 36.358 Santorini COMET 26/6/04 26/11/14 7.8 -13.6 1.4 0.3 0.3 0.8 SOFI 23.395 42.556 Sofia BAS 4/6/97 21/3/20 23.5 12.7 -0.4 0.1 0.1 0.4 SPAN 20.674 38.781 Spanochori NOA 22/5/07 31/12/19 21.4 4.7 -0.8 0.3 0.2 0.6 SPAR 22.424 37.078 Sparti Uranus 22/7/15 13/2/20 5.5 -12.1 -0.3 0.4 0.4 1.1 SPET 23.155 37.267 Spetses Metrica 10/12/08 27/2/20 6.9 -12.6 -0.3 0.2 0.2 0.6 SPRT 22.372 37.046 Sparti COMET 21/7/06 8/2/08 6.5 -10.5 4.7 0.6 0.6 2.2 STEF 22.743 39.464 Stefanovikio NTUA 14/9/12 4/9/17 11.7 -0.6 -116 0.5 0.4 3.0 STRA 23.791 40.520 Stratoni Metrica 26/10/10 26/2/20 24.2 7.1 -0.8 0.3 0.3 0.8 STRF 21.016 37.245 Strofades INGV 5/10/16 6/11/18 3.3 -6 1.7 0.4 0.4 1.1 SVI1 26.213 41.759 Svilengrad BAS 19/9/17 2/3/20 22.6 13 0.5 0.4 0.5 1.6 SYR1 24.939 37.440 Ermoupoli Metrica 28/6/16 2/3/20 7.8 -12.3 -0.3 0.4 0.4 1.0 SYRO 24.945 37.446 Syros Uranus 5/3/16 24/2/20 7.6 -12 -1.7 0.4 0.5 1.0 TEI1 23.551 41.074 Serres TEIS 26/11/07 8/7/09 23.7 11.1 -2 0.6 0.6 1.6 TEIS 23.555 41.074 Serres TEIS 24/6/10 28/8/12 24.3 10.6 -4.8 0.6 0.8 1.6 TEKR 27.496 40.958 Tekirdag Turkey 16/1/09 31/10/15 22.7 11.6 -0.6 0.2 0.2 0.6 THAS 24.535 40.714 Thassos Uranus 3/1/16 7/3/20 23.1 10.6 0.3 0.4 0.4 0.9 THES 22.929 40.504 Thessaloniki Uranus 1/1/16 8/3/20 24 6.6 -3.4 0.4 0.4 0.9 THI1 23.317 38.316 Thiva Metrica 18/1/07 28/3/11 9.2 -7.2 0.5 0.4 0.4 1.0 THIR 25.452 36.385 Thira NOA 28/12/12 9/12/15 8.2 -16.3 3.8 0.5 0.5 1.2 THIU 23.319 38.320 Thiva Uranus 1/1/16 20/2/20 9.8 -8 -0.9 0.4 0.4 1.1 THIV 23.318 38.317 Thiva Metrica 9/6/11 29/2/20 9.5 -8 1.1 0.2 0.3 0.7 THS1 22.959 40.627 Thessaloniki Metrica 25/2/12 2/3/20 24.3 7.9 0 0.2 0.3 0.7 TILO 27.394 36.380 Tilos Unavco 25/8/12 27/7/15 13.5 -15.3 -0.9 0.4 0.4 1.0 TINU 25.169 37.539 Tinos Uranus 21/6/16 2/7/17 7.6 -12.1 0 1.0 1.8 3.5

42 TRIK 21.798 39.555 Trikala Uranus 17/6/17 4/4/20 20 7.4 -3.8 0.6 0.8 1.8 TRIP 22.369 37.508 Tripoli Metrica 4/9/11 6/2/20 7.6 -10.6 -4.4 0.3 0.4 1.3 TRIU 22.367 37.509 Tripoli Uranus 1/1/16 10/1/20 7.3 -12.2 -2.8 0.4 0.4 1.1 TRIZ 22.073 38.365 Trizonia CNRS 17/5/02 10/2/20 8.5 0.6 -2.1 0.2 0.3 0.6 TROP 21.959 37.730 Tropaia Uranus 1/1/16 10/1/20 6.5 -9.9 0.7 0.4 0.4 1.0 TUC2 24.071 35.533 Chania TEIC 15/10/04 31/12/19 8.4 -11.5 -0.7 0.2 0.2 0.5 VALI 22.135 38.234 Valimitika CNRS 6/6/15 6/5/20 10.2 -8.8 2.3 0.3 0.5 0.9 VALY 21.984 37.164 Valyra CUP 7/12/11 5/2/20 3.6 -12.4 -0.7 0.5 0.4 1.0 VAM0 24.199 35.406 Vamos TEIC 27/2/13 5/6/15 8.1 -11.1 -2.6 0.8 1.0 2.3 VASS 21.899 37.430 Vasses COMET 11/3/04 8/2/08 5 -12.5 1.3 0.4 0.4 1.3 VERI 22.240 40.496 Veria Metrica 25/11/10 2/3/20 24.2 7.5 1.4 0.3 0.3 0.9 VLSM 20.589 38.177 Valsamata NOA 26/9/06 31/12/19 17.1 5.1 -0.6 0.2 0.3 0.6 VOLO 22.887 39.353 Volos Metrica 7/7/10 3/3/20 19.7 2.7 -0.2 0.3 0.3 1.1 VOLU 22.761 39.181 Volos Uranus 1/1/16 22/2/20 19 1.7 0.1 0.3 0.4 1.0 VOUR 21.185 39.371 Vourgareli Uranus 22/7/15 19/8/19 20 9.9 0 0.4 0.5 1.2 VRAS 23.704 40.712 Vrasna Uranus 1/1/16 7/3/20 23.3 8.4 0.9 0.4 0.5 1.1 WNRY 25.437 36.387 Winary NTUA 1/2/12 31/12/15 8.2 -17.1 4.3 0.5 0.5 1.1 XALK 23.607 38.470 Xalkida Uranus 19/10/16 20/2/20 11.7 -9 -0.7 0.7 0.5 1.9 XANT 24.858 41.115 Xanthi Uranus 31/3/16 7/3/20 23.3 10.5 -0.4 0.4 0.5 1.0 XILI 21.912 38.385 Xiliadou CNRS 8/5/14 10/2/20 11.1 0.6 -6.2 0.3 0.4 1.0 XRSO 23.533 35.311 Xrissokellaria COMET 16/3/03 8/2/08 7.7 -11.5 0.9 0.6 0.4 1.1 YENC 27.242 39.936 Yenice Turkey 16/1/09 31/10/15 5.5 4.3 -1 0.2 0.2 0.6 YSST 28.991 40.866 Yassi Ada Turkey 2/11/11 23/6/15 20.5 15 1.5 0.2 0.3 0.9 YUND 23.854 42.064 Yundola BAS 3/6/07 26/10/10 24.2 12.7 1.6 0.5 0.5 1.6 ZAGO 23.102 39.446 Zagora Uranus 31/3/18 18/3/20 22.7 5.7 3.1 0.6 0.6 1.6 ZAKU 20.889 37.779 Zakynthos Uranus 1/1/16 13/2/20 12.3 0 1.4 0.4 0.4 1.0 ZAK1 20.878 37.781 Zakynthos Metrica 11/6/11 25/7/13 11.6 -0.2 -3.5 0.5 0.5 1.2 ZAKY 20.885 37.779 Zakynthos Metrica 8/8/13 15/3/20 12.2 0.2 1.7 0.3 0.3 0.8 ZKRO 26.217 35.115 Zakros TEIC 18/5/14 9/3/20 9.9 -13.9 -0.2 0.6 0.4 1.2 028Α 20.845 37.729 Laganas HEPOS 1/1/15 10/4/18 10.3 0.2 -0.1 0.9 0.8 1.9 029Α 21.649 37.485 Zacharo HEPOS 1/1/15 2/4/18 5.7 -11.2 0.1 0.7 0.7 1.7 030Α 21.269 37.909 Andravida HEPOS 1/1/15 9/3/18 9.7 -5.6 1.3 0.7 0.7 1.6 918

919

920 Table S2. Glossary and coordinates of site names (Fig. 2 and Fig. S3) used in the text, figures and

921 tables. In the column “Text”, M means main document and S supplementary material.

922 Toponym Long. Lat. Text Figure Table Comment Aegean Sea 25.26 39.03 M, S 2, S3 Sea in the eastern Mediterranean Aetolia- 21.37 38.7 M, S 2 Region of central Greece Acarnania Agrinion 21.4 38.62 M, S 2 5, S1 City in Aetolia-Acarnania Town, known for the earthquake of June Aigion 22.08 38.25 M, S 2, S3 4, 6 1995 Albania 20.26 41.15 M 1, 2, 10 Country Alexandroupolis 25.87 40.84 M 2 Town in north-eastern Greece Alonissos 23.89 39.19 M 2 6 Island of the Sporades archipelago Ambracian Gulf 21.04 38.95 M 2 4, 6 Known also as Gulf of Arta Amfilochia 21.16 38.86 M, S 2 S1 Town in Aetolia-Acarnania Island of the eastern Cyclades, known for Amorgos 25.88 36.84 M 2 4, 6 its earthquake of 1956 Andros 24.87 37.83 M, S 2 4, S1 Island in the Cyclades archipelago Antirio 21.76 38.33 S S3 Town north of Patras

43 Province of south-eastern Italy and Apulian 17.1 40.79 M tectonic platform part of the African plate Arta 20.99 39.16 M, S 2 S1 Town in Epirus Astypalea 26.35 36.54 M, S 2, S3 S1 Island in the Dodecanese archipelago Athens 23.72 37.98 S S3 Capital city of Greece Athos Mnt. 24.22 40.25 M 2 4 Mountain range of Chalkidiki Attica 23.85 38.04 S Region in central Greece Bulgaria 25.48 42 M, S 1, 2, S3 Country Central 23.19 40.6 M, S 2, S32 5 Region in northern Greece Macedonia Island in the Ionian Sea, host place of the Cephalonia 20.56 38.17 M, S 2, 5, S3 1, 2, S4 earthquakes of 2014 Chalkidiki 23.88 40.08 S S3 Peninsula in northern Greece Chios 26.13 38.36 M, S 2, S3 S1 Island in the northern Aegean Sea Corfu 19.88 39.6 M 2 Island in the Ionian Sea Corinth 22.93 37.93 M, S 1, 2, S3, S15 4, 6 City in the north-eastern Peloponnese Crete 24.76 35.24 M, S 2, 9, S3, S31 1, 2 The largest island in southern Greece 1, 2, S1, Village, known for the twin earthquakes Efpalion 21.92 38.41 S 2, S3 S4 of 2010 Epirus 20.76 39.57 M 2 4, 5 Region in north-western Greece Village of the northern coast of the Gulf Eratini 22.22 38.36 M, S 2, S3 1, 2, S4 of Corinth Europe 15.25 54.52 M, S 7, 9 S6 Continent Evia 23.85 38.52 M, S 2, 13, S3 4, 6 Island in central Greece Florina 21.21 40.78 M, S 2 4, S1 City in northern Greece Gallikos 22.89 40.87 S S3 River of northern Greece and its delta Greece 21.82 39.07 M, S 1, 2, S3 Country Town, known for the earthquake of May Grevena 21.42 40.08 M, S 2, S3 S1 1995 Gulf of Arta 21.02 38.95 M 2 Known also as Ambracian Gulf Gulf of Corinth 22.57 38.15 M, S 1, 2, 10, S3, S15 4, 6 Gulf in central Greece Gulf of Evia 23.62 38.54 M 2, 13 4 Gulf in central Greece Gulf of Gökova 27.85 36.9 M 2 4 Gulf in south-western Turkey Gulf of Patras 21.69 38.28 M, S 2 4, 6 Gulf in western Greece Village of western Turkey, known for its Gulpinar 26.11 39.53 M, S 2, S3 1, 2, S4 earthquake swarm of 2017 Ierapetra 25.74 35.01 M, S 2 S1 City in south-eastern Crete Ikaria 26.11 37.59 M, S 2, S3 4, S1 Island in eastern Aegean Ioannina 20.85 39.66 M, S 2 S1 Town in western Greece Ionian Sea 20.27 38.78 M, S 2, S3 5 Sea between Greece and Italy Istanbul 28.97 41 M, S 2, 7, S25, S27 S1 Major town of Turkey Town known in Turkey for the destructive Izmit 29.94 40.76 S earthquake of 1999 Kalamata 22.11 37.03 M, S 2 S1 City in Peloponnese Kalavrita 22.11 28.03 M, S 2 S1 Town in Peloponnese Kalymnos 26.97 36.95 M, S 2, S3, S19, S20 S1 Island of the south-eastern Aegean Sea Kanallaki 20.6 39.23 M, S 2, S3 1, 2, 4, S4 Town, known for its earthquake of 2020 Islet of the eastern Aegean Sea opposite Karaada Islet 27.46 36.97 S S3 to the city of Bodrum Karpathos 27.14 35.55 M, S 2, S3, S28, S29 S1 Island in the Dodecanese archipelago Kastoria 21.26 40.51 S S3 S1 Town in northern Greece Village west of Agrinion, known for its Katouna 21.11 38.78 M 2 4, 6 fault Kavala 24.41 40.93 S S3 S1 City in northern Greece Kiato 22.75 38.01 M, S 2 6, S1 Town in Peloponnese Killini 21.14 37.92 M 2 Town in Peloponnese 1, 2, 3, 4, Island of the Dodecanese archipelago, Kos 27.08 36.79 M, S 2, 4, 5, S3 S4 known for its earthquake of 2017

44 Krini 21.96 38.19 S S21, S22, S23 S1 Village in Peloponnese Krupnik 23.11 41.85 S S3 Village in Bulgaria Kymi 24.1 38.63 S S3 S1 Town in Evia Larissa 22.41 39.63 S S3, S16 S1 Capital of Thessaly Region Lefka Mnt. 24.03 35.29 M 2 Mountain range in western Crete 1, 2, 3, 4, Island in Ionian Sea, known for its Lefkada 20.67 38.77 M, S 2, 4, 5, S3 6, S1, S4 earthquakes of 2003 and 2015 Lemnos 25.14 39.91 M, S 2, 7, S3 6, S1 Island of the north Aegean Sea Leros 26.84 37.14 M, S 2 S1 Island in Dodecanese archipelago 2, 4, 5, S3, S19, 1, 2, 3, Island of the northern Aegean Sea, known Lesvos 26.27 39.26 M, S S34 S1, S4 for its earthquake of 2017 Suburban town of Athens, known for its Magoula 23.52 38.07 M, S 2, S3, S9 1, 2, S4 earthquake of 2019 Marmara Sea 28.11 40.4 M, S 2, 7, S3 Sea in Turkey Marmaras 23.78 40.09 S S3 Town in the peninsula of Chalkidiki (Neos) 1, 2, 3, Methoni 21.7 36.81 M, S 2, 4, 5, S3 Village known for its earthquake of 2008 S1, S4 Monemvasia 23.05 36.68 M 2 Village in the Peloponnese Mornos 22.09 38.53 S3 River and river delta Mountain range and village near the 1, 2, 3, 4, epicentre of the June 8, 2008 earthquake. Movri Mnt. 21.45 38.11 M, S 2, 5, S3, S8 6, S4 Used to designate the corresponding fault and fault zone Naoussa 22.07 40.62 M, S 2 4, S1 Town in northern Greece Naxos 25.37 37.1 M, S 2, S3 S1 Island in the Cyclades Northern 22.6 41.6 M 1, 2, S3 Country Macedonia Ohrid 20.8 41.12 S S3 S1 City and its lake in Northern Macedonia Mountain and mountain range in northern Olympus Mnt. 22.34 40.08 M 2 4 Greece Palliki Peninsula 20.39 38.23 M 2 Peninsula of the Cephalonia island Paramythia Mnt. 20.51 39.45 M 2 6 Mountain in Epirus Parga 20.39 39.28 M Town in western Greece Parnitha Mnt. 23.71 38.17 S Mountain north of Athens Patras 21.73 38.24 M 2 4, 5, 6, S1 City of western Greece Paxoi Island 20.16 39.21 M, S 2 S1 Island in the Ionian Sea Peloponnese 22.37 37.54 M, S 2, 9, S3, S30 4, 5, S5 Peninsula in southern Greece Permet 20.35 40.23 M, S 10 S1 Town in Albania Pilion 23.04 39.43 M 2 Mountain in northern Greece Mountain range located in north-western Pindos Mnt. 20.75 40.00 M 2, 10 4, 5 Greece and southern Albania Preveza 20.75 38.95 M,S 2, S3 1, 2, S4 City known for its earthquake of 2019 Psiloritis Mnt. 24.77 35.22 M 2 Mountain range in central Crete Pylos 21.69 36.91 M, S 2, S3; S30 S1 Town in the south-western Peloponnese Rio 21.78 38.3 M, S 2, S3 Town north of Patras Rodini 21.9 38.3 S S3 S1 Village in the Peloponnese Rodos 28.21 36.43 M, S 2, S3, S26, S27 S1 Island in the Dodecanese archipelago 1, 3, 4, Samos 26.97 37.75 M 2, 4, 5, S1, S3 Island of the eastern Aegean Sea S1, S4 1, 2, 3, 4, Island of the northern Aegean Sea, known Samothraki 25.59 40.44 M, S 2, S2, S3 S4 for its earthquake of 2014 2, 7, S3, S26, Sandanski 23.28 41.56 M, S S1 Town in southern Bulgaria S27 1, 2, S3, S12, Santorini 25.44 36.4 M, S S5 Island in Cyclades, known for its volcano S13 Selinountas 22.14 38.22 S S3 River and river delta Sitia 26.1 35.2 M, S 2 S1 Town in Crete Skopelos 23.67 39.12 M, S 2 S1 Island of the north Aegean Sea

45 1, 2, 3, Island in the Sporades archipelago, Skyros 24.59 38.81 M, S 2, S3 S1, S4 known for its earthquake of 2001 Sparta 22.42 37.07 M 2 City in Peloponnese Island in the entrance of the Gulf of Spetsopoula 23.17 27.22 S S3, S30 Argos, eastern Peloponnese Village in Aetolia-Acarnania (see Stamna 21.28 38.51 M 2 4 Katouna) Strofades 21 37.25 M, S 2 S1 Islands in the Ionian Sea Strymon 23.39 41.34 M 2 4 River in the Balkan Peninsula Thessalia 22.04 39.61 S S3 Region in central Greece Second largest city of Greece, located in Thessaloniki 22.94 40.64 M, S 2, S3 S1 central Macedonia Thiva 23.32 38.32 M, S 2 S1 City in central Greece Thrace 24.88 41.12 M 2 5 Region in north-eastern Greece Tilos 27.37 36.42 M 2 Island in Dodecanese archipelago Tinos 25.16 37.57 M 2 4 Island in Cyclades Archipelago 2, 7, S3, S28, Tirana 19.81 41.32 M, S Capital city of Albania S29 Trichonis Lake 21.55 38.56 M 2 4, 6 Lake in Aetolia-Acarnania Tripoli 22.37 37.51 M, S 2, S3, S17, S18 City in the central Peloponnese Trizonia Isl. 22.06 38.37 M, S 2, S3 S1 Islet in the Gulf of Corinth Turkey 28.88 38.96 M, S 1, 2, S3 Country Vlora 19.49 40.46 M, S 2 S1 City in Albania Volos 22.95 39.36 M, S 2 S1 City in Thessalia Volvi Lake 23.53 40.65 S S3 Lake in the root of Chalkidiki Peninsula Vrasna 23.65 40.7 S S3 S1 Village near the Volvi lake Zakros 26.24 35.09 M, S 2 S1 Ancient site of eastern Crete 1, 2, 4, 5, S3, 1, 2, 3, Island in Ionian Sea, known for its Zakynthos 20.84 37.78 M, S S10, S11, S14, S1, S4 earthquakes of 2006 and 2018 S24 923

924

925 Table S3. Scatter of velocities for 110 pairs of stations separated by less than 15 km. The formal

926 uncertainties provided by the calculations and the linear regression of the time series are scaled to

927 bet fit with the scatters of velocities. See Fig. 3 for graphic representation of the east, north and

928 vertical scatters.

929

Sta-1 Sta-2 Dist diff horizontal diff_E diff_N diff_U σ_vE σ_vN σ_vU

velocity km mm yr-1 DION DYNG 0.02 0.14 0.1 0.1 0.1 0.58 0.81 1.97 AGRI AGRU 0.13 0.00 0 0 0.6 0.42 0.50 1.08 THI1 THIV 0.13 0.85 0.3 0.8 0.6 0.45 0.50 1.22 MOLA MOLU 0.18 0.28 0.2 0.2 0.3 0.50 0.50 1.35 TRIP TRIU 0.26 1.63 0.3 1.6 1.6 0.50 0.57 1.70 SKOU SKP1 0.28 0.54 0.5 0.2 1.5 0.64 0.71 1.42 MARM MARU 0.31 1.20 1.2 0 1.4 0.64 0.64 1.84 NEAB NEAU 0.31 1.28 0.8 1 0.4 0.72 0.50 1.80 TEI1 TEIS 0.33 0.78 0.6 0.5 2.8 0.85 1.00 2.26

46 ZAKU ZAKY 0.35 0.22 0.1 0.2 0.3 0.50 0.50 1.28 THIU THIV 0.35 0.30 0.3 0 2 0.45 0.50 1.30 AGPA PRKV 0.37 0.73 0.2 0.7 0 0.50 0.50 1.22 SKOP SKOU 0.40 2.51 2.2 1.2 0.1 0.50 0.54 1.14 ITEA ITEU 0.41 0.32 0.3 0.1 0.5 0.50 0.58 1.42 KALM KLAM 0.44 0.00 0 0 0.7 0.58 0.58 1.22 AFY0 AFYT 0.45 1.28 0.8 1 0.3 0.57 0.57 1.27 NAXO NAXU 0.55 0.22 0.1 0.2 1.4 0.50 0.50 1.14 ASTU ASTY 0.56 0.20 0 0.2 1.2 0.50 0.58 1.39 CHI1 CHIU 0.58 1.14 0.3 1.1 7.2 0.94 0.81 1.72 LEUK LFKD 0.60 1.52 1.4 0.6 0.7 0.64 0.57 1.35 GREU GREV 0.62 0.85 0.6 0.6 1.1 0.45 0.58 1.53 ZAK1 ZAKY 0.68 0.72 0.6 0.4 5.2 0.58 0.58 1.44 HALK XALK 0.68 1.21 1.1 0.5 0.7 0.73 0.58 1.99 KAL1 KAL2 0.72 0.85 0.8 0.3 1.3 0.42 0.50 1.50 ARTA ARTU 0.74 1.84 1.6 0.9 0.9 0.50 0.50 1.35 SYR1 SYRO 0.75 0.36 0.2 0.3 1.4 0.57 0.64 1.41 CHIO CHIU 0.79 0.71 0.1 0.7 0.7 0.58 0.50 1.22 IOAN IOAU 1.06 0.85 0.3 0.8 0.9 0.58 0.50 1.22 KOM1 KOMO 1.13 0.71 0.7 0.1 0.2 0.50 0.72 1.56 LEFK LEUK 1.19 1.08 0.6 0.9 5.5 0.78 0.81 1.80 KATE KATU 1.20 0.50 0.5 0 2.7 0.45 0.50 1.30 SAND SNDA 1.41 1.10 0 1.1 0.6 0.50 0.58 1.63 AIGI EGIO 1.42 1.87 0.5 1.8 4.3 0.67 0.99 2.03 NOMI SANT 1.43 1.80 1.8 0.1 1.3 0.42 0.42 0.92 THIR WNRY 1.45 0.80 0 0.8 0.5 0.71 0.71 1.63 KALU KALY 1.49 1.14 0.7 0.9 0.8 0.58 0.58 1.42 ALE1 ALEX 1.67 0.92 0.2 0.9 0.3 0.63 0.85 1.80 PAT3 PATR 1.75 1.33 1.3 0.3 5.5 0.58 0.67 1.66 HGOU IGOU 1.97 0.61 0.1 0.6 1.8 0.57 0.50 1.35 MKMN PKMN 2.03 2.68 2.4 1.2 3.6 0.50 0.50 1.39 MOZI SANT 2.44 2.38 0.9 2.2 0.2 0.58 0.58 1.08 LARI LARM 2.72 1.00 0.8 0.6 0.4 0.45 0.50 1.14 LAMA LAMI 2.73 0.51 0.5 0.1 1.3 0.45 0.58 1.30 PAT0 PAT2 2.81 2.21 2.2 0.2 0.3 0.58 0.67 1.34 ALEU ALE1 2.97 1.55 0.4 1.5 0.4 0.72 0.89 2.03 DSLN MOZI 3.02 2.28 1.4 1.8 0.8 0.64 0.71 1.42 ARSA PAT0 3.26 2.42 0.3 2.4 0.1 0.50 0.58 1.08 KORO KORU 4.05 0.51 0.5 0.1 2 0.72 0.81 1.75 KERA PKMN 4.07 1.12 1 0.5 0.1 0.50 0.42 1.06 RIBA KERA 4.14 1.34 0.6 1.2 0.4 0.57 0.50 1.20 GEYB NAFP 4.20 1.70 1.7 0.1 2.9 0.42 0.64 1.22 KOUN AIGI 4.31 0.71 0.1 0.7 4.5 0.36 0.45 1.30 PYL1 PYLO 4.64 0.89 0.8 0.4 4 0.85 0.58 1.22 KITH KTHA 4.71 1.28 0.8 1 5.9 0.99 1.08 2.55 EYPA XILI 4.83 2.12 0.3 2.1 4.8 0.36 0.45 1.17 SNTR PKMN 4.94 0.80 0.8 0 1.3 0.42 0.42 1.06 IKAR IKAU 5.03 1.03 0.9 0.5 3.2 0.92 1.06 2.62 ART1 ARTA 5.09 2.20 1.4 1.7 0.2 0.67 0.67 1.58 FLOR FLOU 5.31 0.76 0.3 0.7 0.2 0.50 0.57 1.64 CHA1 TUC2 5.38 0.81 0.7 0.4 1.2 0.45 0.54 1.30 VALI AIGI 5.52 0.63 0.2 0.6 1.2 0.42 0.64 1.42 ARGO KEFA 5.60 2.75 0.9 2.6 0.2 0.57 0.57 1.41 EXAN SPAN 5.66 2.04 2 0.4 0.1 0.42 0.36 1.62 DUTH XANT 5.68 1.43 0.3 1.4 0.7 0.45 0.58 1.17 SPAR SPRT 5.78 1.89 1 1.6 5 0.72 0.72 2.46

47 DRAM DRAU 5.86 0.50 0.4 0.3 1.3 0.42 0.72 1.20 KARU KRPS 6.29 1.14 1.1 0.3 1.1 0.85 0.64 1.84 028A ZAK1 6.44 1.36 1.3 0.4 3.4 1.03 0.94 2.25 ABEL AMFI 6.60 0.85 0.6 0.6 5 0.57 0.64 3.27 NOA1 DYNG 6.95 0.36 0.3 0.2 0.1 0.36 0.50 1.06 GALA ITEA 7.03 2.00 0.1 2 0.5 0.50 0.58 1.42 LAMB ROD3 7.21 2.72 0.3 2.7 1.5 0.42 0.57 1.35 MOIR SIVA 7.23 1.27 0.9 0.9 1.4 0.58 0.67 1.48 AUT1 THS1 7.71 0.72 0.6 0.4 1.2 0.28 0.36 0.86 MILO MLOS 7.87 0.40 0.4 0 1.1 0.36 0.50 1.25 KIPO KEFA 7.90 2.43 0.4 2.4 2.2 0.64 0.64 1.41 LEMN MOUD 8.00 1.24 0.3 1.2 0.2 0.45 0.54 1.12 MESO RETS 8.01 1.61 0.2 1.6 3.2 0.50 0.64 1.49 NEOS TEIS 8.12 0.67 0.3 0.6 3.5 0.72 0.94 1.89 KRP1 KRPS 8.29 0.20 0 0.2 0.1 0.78 0.64 1.72 METH PYL1 8.45 2.50 2.5 0.1 2.8 0.85 0.58 1.35 VLSM ARGO 8.47 0.67 0.6 0.3 0.3 0.45 0.50 1.17 AKYR MOIR 8.72 1.01 0.1 1 2.2 0.50 0.58 1.39 MET0 NOA1 9.16 0.64 0.4 0.5 0.1 0.28 0.36 0.92 ATRS LEON 9.34 0.45 0.4 0.2 0.6 0.50 0.57 1.56 THES AUT1 9.40 1.27 0.9 0.9 2.2 0.45 0.45 1.03 RAFI DION 9.43 0.85 0.3 0.8 1.6 0.71 0.86 2.22 ASSO FISK 9.60 1.36 0.4 1.3 1.9 0.64 0.64 1.56 PVOG AGRI 10.02 1.10 0 1.1 0.4 0.50 0.50 1.00 ANAV KERT 10.43 0.50 0.4 0.3 1.2 0.57 0.64 1.49 SISS VLSM 10.46 2.94 0.9 2.8 3.9 0.54 0.58 1.62 LIMN LEMN 10.54 1.57 0.7 1.4 2.4 0.45 0.36 0.94 EPID LYGO 10.90 0.40 0.4 0 0.3 0.50 0.57 1.35 ARKI ATAL 11.71 0.14 0.1 0.1 0.1 0.63 0.58 1.57 ORES OREU 11.83 0.94 0.5 0.8 0.3 0.36 0.50 1.08 LESU LESV 11.93 1.03 0.9 0.5 0.4 0.50 0.58 1.35 KARD KLOK 12.07 0.76 0.7 0.3 0.9 0.45 0.54 1.25 KERY PTKG 12.42 1.65 1.6 0.4 0.3 0.99 0.72 1.84 KIAT KRYO 12.46 0.64 0.4 0.5 2.3 0.50 0.64 1.56 AGNI NEAP 12.51 0.85 0.6 0.6 1.3 0.81 0.64 1.56 PALC XRSO 12.60 0.32 0.3 0.1 0.1 0.72 0.72 1.49 ROZH SMLN 13.62 0.45 0.2 0.4 0 0.54 0.63 1.43 PROV TEI1 13.83 1.14 0.3 1.1 0.8 0.67 0.67 1.75 LEPE AGRU 14.03 0.45 0.2 0.4 3.3 0.67 0.81 1.75 SITI ZKRO 14.22 1.33 0.3 1.3 5.2 0.67 0.57 1.50 KASI KERU 14.23 1.87 1.8 0.5 0.3 0.50 0.50 1.08 ANOA RLSO 14.24 0.63 0.2 0.6 0.9 0.54 0.54 1.25 KERK KERU 14.38 1.02 1 0.2 0.3 0.50 0.50 1.42 AMAL 030A 14.67 2.30 1.9 1.3 2.3 0.76 0.81 1.79 ANIS PALC 14.91 0.73 0.2 0.7 1.1 0.64 0.78 1.80 930

931

932 Table S4. Observed and modelled coseismic and postseismic displacements for the earthquakes

933 listed in Tab. 1. The faults parameters are in Tab. 2. We provide those numbers only for the events

934 modelled in this paper. For the event modelled in previous works, the observed and modelled

48 935 displacements can be found in Ilieva et al. (2016), Serpetsidaki et al. (2014), Elias (2013), Briole et

936 al. (2015), Ganas et al. (2016, 2018b, 2019, 2020). For the Methoni 2008 earthquake the vertical

937 displacement cannot be assessed for the instantaneous coseismic but it can be assessed for the first

938 ten days.

939

Earthquake Phase Site Long. Lat. Observation (mm) Model (mm)

° ° East North Up East North Up Skyros NSKR 24.5430 38.8870 -58 14 -57.9 13.8 coseismic CG38 24.1090 38.6610 -16 -17 -7.6 -20.2 2001 DION 23.9330 38.0790 -1 -3 -0.9 -3.5 Lefkada coseismic see Ilieva et al. (2016) 2003 Zakynthos coseismic KERI 20.8080 37.6550 -10 -2 8 -8.9 -2 11 2006 ATRS 22.8348 37.2444 -3.9 -4.1 -2 -2.8 KERY 22.3836 36.4933 -17.7 -9.1 -13.8 -11.2 KITH 23.0151 36.2745 -2.9 -0.1 -1.3 -0.6 coseismic METH 21.7045 36.8255 -1.8 -19.7 1.3 -20.3 PYL1 21.7420 36.8955 -1.8 -14.9 0.3 -16.7 Methoni SPRT 22.3721 37.0460 -5.8 -8.9 -3.7 -7.2 VASS 21.8993 37.4304 -0.1 -4.9 -0.5 -4.6 2008 KERY 22.3836 36.4933 -22.5 -12.4 -7.3 -18.3 -14.9 -9.9 coseismic METH 21.7045 36.8255 -1.0 -25.4 6.0 2.1 -27.1 -9.9 PYL1 21.7420 36.8955 0.0 -23.5 0.0 0.6 -22.3 -7.2 + 10 days VASS 21.8993 37.4304 1.2 -2.5 5.8 -0.6 -6.1 -0.9 KERY 22.3836 36.4933 -43.1 -17.4 -5.8 -27.6 -22.2 -13.3 whole METH 21.7045 36.8255 -12.8 -42.3 -14.1 -16.9 -60.6 -21.6 PYL1 21.7420 36.8955 -22.5 -47.3 -11.3 -15.5 -50.3 -15.7 sequence VASS 21.8993 37.4304 -0.6 -12.7 8.9 -5.7 -16.5 -2.2 Movri coseismic see Serpetsidaki et al. (2014) 2008 Efpalion coseismic see Elias (2013) 2010 West Crete coseismic TUC2 24.0706 35.5332 -5.2 -4.2 -0.2 -5.4 -2.9 -1.7 2013 Cephalonia

2014 (26/1 coseismic see Briole et al. (2015)

& 3/2) Samo- coseismic LEMN 25.1805 39.8972 -21.4 -58.4 -6.9 -20.7 -56.1 -9.3 PRKV 26.2650 39.2457 -4.2 -0.4 -3.8 -1.1 thraki 2014 ALEX 25.8534 40.8491 13.5 18.2 -5.7 11.5 20.7 -5.1 ORES 26.5458 41.4142 3 7 2.8 7.1

49 KRUM 25.6519 41.4735 2.9 4.2 2.3 5.0 STRA 23.7908 40.5202 2.1 0.9 5.5 -1.9 THS1 22.9593 40.6275 1.6 1.4 2.5 -0.9 SKYR 24.5646 38.9045 -2.4 -4.5 -1.8 -4.5 CANA 26.4143 40.1112 -31.6 15.3 3.6 -34.0 12.8 8.0 IPSA 26.3798 40.9176 10.5 18.5 7.9 19.4 YENC 27.2417 39.9359 -7.3 4 -8.0 2.9 TEKR 27.4964 40.9583 0.8 2.8 0.0 2.5 LEMN 25.1805 39.8972 -37 -85.4 -13.2 -38.1 -85.8 -7.3 PRKV 26.2650 39.2457 -1.6 -6.5 -6.5 -1.2 ALEX 25.8534 40.8491 16.7 37 2.7 17.6 32.8 -6.7 whole ORES 26.5458 41.4142 3.4 11.5 4.5 11.4 KRUM 25.6519 41.4735 3.4 5.7 3.8 8.4 sequence THS1 22.9593 40.6275 2.3 0.5 4.5 -1.5 YENC 27.2417 39.9359 -11.4 4.2 -12.7 4.6 TEKR 27.4964 40.9583 1.1 7.1 0.1 3.9 CANA 26.4143 40.1112 -54 14.9 3.6 -51.9 19.9 9.1 coseismic see Ganas et al. (2016) PONT 20.5852 38.6190 -13 -18 -10 -22 DRAN 20.5746 38.6839 -5 -5 -8 -9 postseismi EXAN 20.6186 38.7540 -8 -12 -6 -5 SPAN 20.6736 38.7813 -5 -4 -6 -5 c 138 days FISK 20.5771 38.4597 0 -14 7 -16 ASSO 20.5417 38.3779 0 -13 2 -9 ARTA 20.9918 39.1588 0.9 -2 -1 -1.6 ARTU 20.9999 39.1568 0.1 -2.3 -1 -1.6 KTCH 21.2470 38.4116 0.7 3 -0.6 0.1 Lefkada ABEL 21.2150 38.9231 -1.1 2.8 -3.1 -1.1 AMFI 21.1498 38.8924 -1.1 2.7 -3.9 -1.3 2015 AGRU 21.4078 38.6234 -1.4 1.5 -1.3 0.1 postseismi VLSM 20.5886 38.1768 -1.8 -5.6 -0.2 -1.1 AGRI 21.4091 38.6240 -2.4 1.2 -1.3 0.1 c 471 days SPAN 20.6736 38.7813 -8.7 -4.5 -3.7 -3.5 PONT 20.5852 38.6190 -12.3 -0.1 -4.6 -3.3 FISK 20.5771 38.4597 -1.9 -3.2 -2.6 -4.6 ASSO 20.5417 38.3779 -2.5 -3.5 -1.4 -3 LFKD 20.7121 38.8345 -2.6 -4.8 -3 -3.3 SISS 20.6594 38.1009 -4.3 -5.6 0 -0.8 KEFA 20.4381 38.1957 -2.7 -3.3 -0.3 -0.9 Gulpinar coseismic see Ganas et al. (2018b) – this swarm was modelled with InSAR only 2017 AGPA 26.2690 39.2455 0.8 13.1 3.8 0.2 12.8 0.5 coseismic LESV 26.5538 39.1001 10.5 18.6 4.4 11 18.5 1.5 LESU 26.4491 39.0306 20.9 40.1 13 20.4 40.3 8.1 + 4 days Lesvos CHIO 26.1272 38.3679 -1.3 -5.6 2.8 -1.8 -7 -0.4 AGPA 26.2690 39.2455 -0.3 17.6 1.4 0.3 15.9 0.3 2017 coseismic LESV 26.5538 39.1001 11.6 23.4 3.8 14.0 24 1.1 LESU 26.4491 39.0306 25.6 57.2 3 27.2 55.5 8.5 + 140 days CHIO 26.1272 38.3679 -3.4 -10.5 -1.4 -2.2 -8.4 -0.6 Kos 2017 coseismic see Ganas et al. (2019) Zakynthos coseismic see Ganas et al. (2020) postseismi 2018 see Ganas et al. (2020) c 21 days postseismi AMAL 21.3554 37.7959 -8.9 -6.8 -4.1 -8.8 -10.3 -2.8 ZAKU 20.8891 37.7792 -8.5 0.8 1.2 -4.7 -3.1 1.8

50 c 181 days PYRG 21.4622 37.6788 -8 -5.8 1.3 -7.6 -6.8 -0.9 ZAKY 20.8851 37.7792 -7.2 -1.3 4.4 -4.7 -3 1.8 KALA 22.1023 38.0306 -6.1 4 0.2 -2.7 -1.5 -0.7 TROP 21.9595 37.7301 -5.3 -0.7 -0.5 -3.4 -0.7 -1.2 KOPA 21.8183 37.2895 -3.9 -1.6 0.4 -0.3 -0.4 -0.2 TRIU 22.3668 37.5093 -1 1 -0.1 -0.8 -0.1 -0.2 VLSM 20.5886 38.1768 -0.8 -1.2 -0.5 -0.1 -0.3 -0.2 TRIP 22.3694 37.5082 -0.1 0 -0.4 -0.8 -0.1 -0.2 STRF 21.0156 37.2454 1.5 -16.2 -0.9 0.4 -13 10.7 RLSO 21.4647 38.0558 -0.1 -8.5 -9.6 -4.4 -8.5 -5.5 GARG 21.6398 37.0662 -4.4 -1.5 -1.6 0.4 -0.7 0 MESO 21.4749 38.3663 -3.6 -3.5 -5 -1.2 -3.5 -1 LEUK 20.708 38.830 2.6 -2.4 2.4 1.6 -1.9 1.4 Preveza LFKD 20.712 38.835 2.7 -3.1 1.1 1.7 -1.9 1.4 coseismic SPAN 20.674 38.781 -2.6 -4.5 0.7 0.8 -1.4 0.7 2019 LEFK 20.695 38.830 1.4 -0.3 1.5 -2 1.5 Eratini VALI 22.1348 38.2338 3 0.4 0.6 0.1 coseismic PSAR 22.1843 38.3218 3 1 -3 2.4 -1 -3 2019 GALA 22.3915 38.3752 0 0 0 0.5 0.7 0.4 ATHI 23.6983 37.9268 1.8 -2.8 -0.7 0.1 -0.4 0.1 Magoula AIGU 23.4398 37.7335 4.4 -1.5 4.7 -0.1 -0.4 0 coseismic AGTH 23.1230 37.9386 0 0 0 -0.1 -0.2 0 2019 MET0 23.7623 38.0653 0 0 0 -0.3 0.1 0 Kanallaki coseismic GARD 20.5710 39.3506 -2.4 2.8 -1.1 1.7 1.1 2020 940

941

942 Table S5. Observed and modelled displacements corresponding for some of the transients presented

943 in the section S4.

944

Event Type of Site Long. Lat. Observation (mm) Model (mm)

event ° ° East North Up East North Up AMAL 21.3554 37.7959 1.4 -3 -6 1.5 -3.3 -4.5 West Slow-slip PYRG 21.4622 37.6788 1.6 -1 -2 1.0 -1.1 -1.2 TROP 21.9595 37.7301 1.5 0 -7 1.3 0.3 0.6 Peloponnese earthquake KOPA 21.8183 37.2895 0.9 -1.0 0.7 ZAKY 20.8851 37.7792 1.0 -0.3 -0.5 2018 ? STRF 21.0156 37.2453 0.0 0.0 0.0 GARG 21.6398 37.0662 0.1 -0.4 0.2 Santorini un- volcanic Observation Model

rest unrest Anomalous velocity Anomalous velocity (mm

12/2/ 2011 (mm yr-1) yr-1) SNTR 25.3571 36.3576 -11.2 -37.4 14.3 -7.9 -27.9 13 - PKMN 25.3846 36.3962 -11.4 -66 60.9 0.2 -68.2 71.3 MKMN 25.4012 36.4085 41.1 -50 50 41.1 -60 105.7

51 KERA 25.3481 36.4179 -83.1 -12.6 54.6 -63.2 -24 69.4 RIBA 25.3435 36.4550 -57.2 33.1 47.8 -44.9 31.9 41.8 DSLN 25.3892 36.4608 7.5 75.2 -7.7 1.8 67.2 67.9 MOZI 25.4222 36.4555 60.1 61.3 51.1 41.3 42.9 49.6 3/5/2012 SANT 25.4226 36.4336 67.5 16.9 61.7 68 16.2 77.9 NOMI 25.4286 36.4217 74.3 3.1 62.5 64.7 -9 62.5 WNRY 25.4367 36.3874 42.2 -38.5 16.8 30.4 -27.7 23.8 SNTR 25.3571 36.3576 0.4 1.9 0.7 0.1 0.4 -0.2 PKMN 25.3846 36.3962 1.2 1.9 -0.6 0 1.1 -1.2 MKMN 25.4012 36.4085 -1.2 0.7 -4.2 -0.7 1 -1.7 Santorini KERA 25.3481 36.4179 2.2 1.4 -0.7 1 0.4 -1.1 relaxation RIBA 25.3435 36.4550 1.6 0.2 -0.3 0.8 -0.5 -0.7 post-unrest DSLN 25.3892 36.4608 -1.4 -4.3 -1.4 0 -1.1 -1.1 post-unrest MOZI 25.4222 36.4555 -2.8 -2.5 -0.6 -0.7 -0.7 -0.8 period SANT 25.4226 36.4336 -1.9 -0.3 -0.5 -1.1 -0.3 -1.3 NOMI 25.4286 36.4217 -0.1 -0.2 0.9 -1 0.1 -1 WNRY 25.4367 36.3874 0.8 -1.6 3.6 -0.5 0.5 -0.4 THIR 25.4525 36.3849 0.8 -0.8 3.1 -0.5 0.4 -0.3 945

946

947 Table S6. Comparison of solutions (velocities and uncertainties) for the stations AUT1, DYNG,

948 PAT0, TUC2 (blue dots in Fig. 4). The solutions ROB-EUREF, LTK-EUREF and UGA-CNRS are

949 published by the European Plate Observing System (EPOS, https://gnssproducts.epos.ubi.pt). The

950 solution NGL is published by the Nevada Geodetic Laboratory (http://geodesy.unr.edu). Our initial

951 solution is the one before correcting the velocity biases induced by the coseismic and postseismic

952 displacements.

953

Solution vE vN vU σE σN σU mm yr-1 mm yr-1 Station TUC2 ROB-EUREF 7.91 -11.7 0.06 0.03 0.04 0.13 LTK-EUREF 7.83 -12.03 -0.19 0.05 0.07 0.07 UGA-CNRS 8.01 -12.12 -0.2 0.2 0.17 0.5 NGL 8.29 -11.65 -0.5 0.21 0.22 0.7

Average 8.01 -11.88 -0.21 r.m.s. deviation 0.25 0.23 0.23

Our - initial 8 -12.1 -0.8 Our - final 8.4 -11.5 -0.7 0.2 0.2 0.5 Station AUT1 ROB-EUREF 25.07 6.82 -1.41 0.01 0.01 0.05 LTK-EUREF 25.21 6.84 -1.56 0.02 0.03 0.05 UGA-CNRS 24.83 6.4 -1.58 0.13 0.15 0.46

52 NGL 25.36 7.19 -1.88 0.18 0.2 0.49

Average 25.12 6.81 -1.61 r.m.s. deviation 0.23 0.32 0.20

Our - initial 25.1 7.3 -1.00 Our – final 24.9 7.5 -1.2 0.2 0.2 0.5 Station DYNG ROB-EUREF 7.5 -13.04 2.3 0.13 0.14 0.53 LTK-EUREF 7.7 -12.37 1.34 0.1 0.14 0.15 UGA-CNRS 8.17 -11.7 1.1 0.34 0.44 1.03 NGL 8.57 -11.35 0.3 0.3

Average 7.99 -12.12 1.58 r.m.s. deviation 0.48 0.75 0.63

Our - initial 8.6 -11.2 0.2 Our – final 8.9 -10.7 0.3 0.3 0.4 0.8 Station PAT0 ROB-EUREF 7.6 -6.76 0.61 0.05 0.05 0.06 LTK-EUREF 7.49 -7.06 -0.54 0.06 0.1 0.1 UGA-CNRS 7.19 -7.13 0.84 0.25 0.3 0.7 NGL 8.25 -6.04 -0.29 0.24 0.24 0.8

Average 7.63 -6.75 0.16 r.m.s. deviation 0.45 0.50 0.67

Our - not corrected 7.4 -6.7 -0.3 Our – final 8.1 -5.8 -0.3 0.3 0.3 0.6 954

955

956 Table S7. ITRF2014 coordinates of the 329 points at the epoch 2020.0 and component of the

957 velocity due to the coseismic and postseismic displacements (plotted in Fig. 4).

958

Code X - ITRF2014 Y - ITRF2014 Z - ITRF2014 vE vN vU m m m mm yr-1 mm yr-1 mm yr-1 ABEL 4632182.9222 1798100.2154 3985882.4918 -1.1 -0.6 -0.1 AFY0 4482594.5277 1942982.0770 4086649.7545 0.1 -0.2 0.1 AFYT 4482813.2962 1943168.9184 4086310.3652 1.2 -0.6 0.5 AGNI 4701329.4692 2264416.1353 3655273.4283 -0.1 -0.4 -0.1 AGPA 4435359.3335 2189114.1461 4013633.6070 -0.5 4.3 0.6 AGRI 4645296.2043 1821316.3987 3959869.7497 -1.8 -0.9 0.4 AGRU 4645369.2138 1821233.6866 3959822.8667 -1.7 -1.1 0.2 AGTH 4632100.8797 1977957.7805 3900159.3453 -0.6 -0.4 0.1 AIGI 4648353.4739 1884923.1388 3926681.0337 -1.5 -1.2 0.2 AIGU 4633842.1748 2009072.0931 3882131.3356 0.1 -0.6 0.7 AKYR 4745189.2171 2203912.8223 3636387.2261 0 0 0

53 ALE1 4347527.6822 2108491.8087 4149545.7095 1.8 4.1 -0.6 ALEU 4345768.3947 2110782.7232 4150203.1076 0.4 1.5 -0.1 ALEX 4348024.5431 2106924.0086 4149806.1308 2.7 4.4 -1.0 ALXA 4480452.1054 1850696.4317 4131207.8992 0.1 -0.1 0.1 AMAL 4699931.7960 1837661.5956 3887632.4013 -3.7 -4.6 -2.5 AMFI 4636041.2909 1793524.2584 3983067.4408 -1.6 -0.9 -0.3 ANAV 4617420.9542 2046555.3683 3882113.9672 -0.6 -0.4 0 ANDR 4577719.5436 2109107.1865 3895528.1143 -0.3 -0.5 -0.2 ANIK 4729934.2300 1938158.2445 3801984.5358 0 -0.6 -0.2 ANIS 4768156.9470 2099021.1241 3668286.1936 0 -0.2 -0.1 ANKY 4752742.6501 2046958.0638 3716342.6346 0.2 -0.3 -0.1 ANOA 4670647.6677 1848260.6093 3917540.6792 -2.1 -2.8 -0.4 ANOC 4631275.8476 1864118.7115 3957893.8269 -1 -0.8 0.2 ANOP 4762454.0124 2129948.9600 3658026.0587 0 0 0 ANTL 4393226.5518 2605026.6410 3807552.5083 0 0.1 0 ARGO 4703171.5177 1757741.3892 3920241.9594 -1 0.6 -0.8 ARGU 4668377.6121 1952120.2026 3869886.3704 -0.6 -0.5 0 ARKI 4583365.5696 1948697.2309 3971175.0016 0 0 0 ARKL 4722575.2536 2229220.4316 3650304.8176 0 -0.3 -0.1 ARSA 4652935.0023 1862616.7519 3931816.2442 -1.1 -1.7 0.3 ART1 4625771.0735 1769494.2308 4005529.7086 0 0 0 ARTA 4623628.7732 1774079.9941 4006075.3868 -1.6 -1.2 0.1 ARTU 4623474.5638 1774784.6212 4005880.1201 -0.3 -0.2 0 ASPR 4456276.1885 1857791.4466 4154022.0800 0.1 -0.1 0.1 ASSO 4688202.2939 1756738.9890 3938431.4482 -0.6 -1 0.1 ASTU 4597082.4486 2276913.0624 3777345.3013 -0.2 -0.6 0.2 ASTY 4597137.0136 2277372.4222 3777033.3316 -0.2 -0.5 0.1 ATAL 4591113.8439 1948751.3047 3962396.6750 -0.4 -0.4 0 ATHI 4612201.5459 2024479.6727 3899754.8002 -0.4 -0.3 0.1 ATRS 4685505.7513 1972962.7698 3839262.5782 -0.7 -1.1 -0.2 AUT1 4466283.3342 1896167.2420 4126096.9751 0.2 -0.2 0.2 AYD1 4459716.7455 2355111.7916 3891560.4024 -0.1 0.1 0 AYVL 4415163.2970 2219261.4904 4019162.6612 -0.8 0.2 -0.3 BAJR 4433855.9736 1620286.5265 4275336.9737 0.1 0 0 BALK 4347061.4874 2301032.0377 4047360.3381 -0.6 0.3 0.2 BAND 4299255.5772 2285638.3346 4106154.1327 0 0 0 BELE 4545267.9253 1941847.2508 4017743.2261 -0.6 -0.6 0.1 BERA 4551491.4265 1652679.3814 4137533.8221 0.1 0 -0.1 BURS 4265347.8176 2365803.3118 4096299.4211 -0.3 0.1 0.1 CANA 4374891.3160 2173079.5382 4087527.8082 0.7 0.9 2.1 CESM 4490049.3519 2226201.3361 3931995.5779 -0.2 -0.2 0 CHA1 4747721.2526 2115577.5157 3684206.0725 -0.1 -0.3 0 CHA2 4748703.0940 2117759.3000 3681704.6577 0.1 -0.1 0 CHIO 4495609.8641 2205020.5340 3937602.2080 -0.8 -1.9 0.1 CHI1 4495217.3546 2206279.9503 3937295.4856 -0.1 -0.1 0 CHIU 4495334.0146 2205747.5585 3937460.3361 -1.1 -2.6 0 DATC 4533056.9617 2379084.3017 3791550.2344 0.1 0.0 -0.1 DEIR 4354177.7763 2378749.0799 3995884.2367 -0.2 0.1 0 DIDI 4511089.0010 2325222.8889 3850343.3982 -0.1 0 0 DION 4595216.4667 2039453.0250 3912626.7627 -0.3 -0.2 0 DNZL 4412026.2812 2454909.1837 3884891.1659 -0.1 0.1 0 DRAM 4390543.5174 1961124.3211 4176408.0993 0.2 -0.1 0.1 DRAN 4667626.1903 1752081.3787 3965215.0995 -3.2 -2.8 0.4 DRAU 4388885.3808 1966684.2662 4175552.4658 0.2 -0.1 0.1 DSLN 4639882.2578 2202109.2771 3769559.5164 -0.1 -0.3 -0.1

54 DUTH 4362689.6393 2026648.4463 4174234.5235 1 0 0.2 DYNG 4595221.4853 2039434.8334 3912627.1696 -0.3 -0.5 -0.1 EDES 4481402.7923 1815268.9531 4146235.2775 0.2 -0.1 -2.9 EGIO 4647361.4592 1884782.3286 3927705.4833 -1.5 -1.3 0.1 ELAS 4537158.5769 1852148.5606 4069049.6817 -0.3 -0.4 0.1 EPID 4649652.6003 1988698.9244 3873750.2695 -0.6 -0.3 0 EXAN 4661896.0610 1754021.1575 3971442.6914 1.9 1.5 5.5 EYPA 4641299.7931 1868463.2881 3942785.9626 0.1 -0.9 -3.9 FARS 4570303.9823 1882285.9274 4017895.4863 -0.3 -0.2 0.1 FETH 4476942.1067 2494267.9366 3784205.8169 -0.1 0 0 FINI 4450189.3326 2584503.1466 3755291.2518 0 0 0 FISK 4681830.8291 1757645.0081 3945556.9200 -0.8 -1.7 1 FLOR 4503588.9483 1765092.8538 4144411.3643 0 0 0 FLOU 4500256.4067 1768719.9518 4146401.6751 0 0 0.1 FREN 4338874.6626 2916762.4650 3641321.3105 0 0 0 GALA 4629221.0900 1907238.2141 3938236.7831 -0.8 -0.6 0.2 GARD 4623897.8135 1735337.7251 4022543.7624 -0.4 0.4 -0.3 GARG 4736664.9453 1879186.5214 3823449.8703 0.3 -1.2 -0.2 GEYB 4652384.5643 1857653.8554 3934628.8549 -0.9 -1.0 0.1 GODE 4375326.7073 1923754.3024 4210078.6788 0.2 -0.2 0.2 GOUM 4459119.2071 1842234.2348 4158151.0064 0 -0.1 0.1 GREU 4549208.4689 1785635.3245 4085560.8390 -0.1 -0.1 0 GREV 4548969.9361 1786215.4270 4085549.9555 -0.1 -0.1 0.1 GYTH 4724277.2787 1963208.4566 3796180.0974 0.1 -0.4 -0.1 HALK 4582206.1384 2002158.6984 3946082.1171 -0.4 -0.5 0 HARC 4293416.2732 2394865.9577 4050966.5544 -0.2 0.1 0.1 HERA 4716223.0674 2213407.7007 3667317.4873 0 -0.3 -0.1 HGOU 4623593.6288 1707374.8465 4034710.5568 -0.3 0.5 -0.2 HRA1 4714991.1721 2217766.8280 3666396.8893 -0.1 -0.4 0 IDI0 4728461.1492 2193915.6758 3664532.1953 0 -0.2 0 IERA 4706217.0083 2274793.7888 3642764.4500 -0.1 -0.2 0 IERI 4447711.3458 1968986.6293 4111977.9278 0.3 -0.3 0.2 IGOU 4622591.7533 1706526.6340 4036186.5731 -0.5 1 -0.3 IKAR 4537348.5736 2235036.6037 3873017.6217 -0.1 -0.1 0.0 IKAU 4536702.6128 2239552.3036 3870917.6576 0.5 -0.3 0.7 IOAN 4595018.1559 1750120.7588 4049693.4410 -0.2 0.1 0 IOAU 4595524.1121 1750642.7706 4048912.9005 -0.2 0.1 -0.1 IPSA 4324027.5111 2144564.6780 4155539.4334 0 0.1 -0.1 ISTA 4208830.0003 2334850.5752 4171267.3826 0.7 -0.3 0.1 ISTI 4566494.3500 1952665.9365 3988574.8681 -0.5 -0.5 0 ISTN 4223660.7099 2325015.2573 4161716.3341 -0.1 0.1 0 ITEA 4624377.5354 1908884.1992 3943058.1908 -0.6 -0.5 0.1 ITEU 4624299.0006 1908571.5269 3943305.2627 -0.8 -0.5 0.3 IZMI 4456584.2957 2278765.5673 3939930.8419 -0.3 0.1 0.1 KAL1 4466060.3442 1882577.9495 4132364.4825 0.2 -0.1 0.3 KAL2 4466275.5813 1881896.1850 4132436.8650 0.1 -0.1 0.2 KALA 4661224.5304 1892941.0929 3908575.1256 -1.3 -1.1 0.1 KALM 4722865.8546 1918848.9528 3820418.1381 -0.1 -0.5 -0.2 KALU 4547987.8155 2313483.6373 3814119.2364 0.1 0.7 1.4 KALV 4660885.1198 1893799.3384 3908621.5695 0 0 0 KALY 4547764.9373 2314818.1821 3813512.8381 0.3 0.5 1.0 KARD 4568777.0391 1853558.3929 4032787.5407 -0.5 -0.5 0.1 KARP 4614696.6560 1846019.2489 3985316.5863 -1 -0.7 0.3 KARU 4622674.0063 2376852.8424 3683909.0837 -0.1 -0.9 -0.1 KARY 4581342.3905 2079978.3270 3906876.2782 -0.4 -0.4 0.1

55 KASI 4616572.4548 1674415.7337 4056441.4463 -0.1 0.2 -0.1 KAST 4525141.0294 1760929.8097 4122756.4540 0.1 0 0.1 KATC 4573403.4024 2409322.2238 3723878.4494 0.0 -0.8 0.0 KATE 4501779.5917 1866655.4639 4100967.0953 0.1 -0.2 0.1 KATU 4502079.1544 1865508.9518 4101174.0614 0 -0.2 0.1 KAVA 4394295.6752 1993277.9271 4157207.5462 0.4 -0.2 0.2 KEFA 4703108.3646 1752637.2274 3922567.7774 -0.5 -0.4 -0.3 KERA 4644059.3266 2200009.9107 3765747.3565 -0.1 -0.1 -0.1 KERK 4635244.0547 1675501.2652 4034866.3800 -0.1 0.2 -0.1 KERT 4610361.4718 2050628.8640 3888621.2786 -0.3 -0.3 0.1 KERU 4625602.7214 1675829.7527 4045532.1787 -0.2 0.4 -0.2 KERY 4746993.4849 1954971.8403 3772452.4343 0 0 0 KIAT 4640093.1778 1945712.3707 3906677.2832 -1 -0.5 0.2 KIKA 4389217.2131 2301675.1916 4001606.8927 0 0 0 KILK 4441924.1710 1874330.3188 4162250.1989 0.1 -0.1 0.1 KIPO 4705430.9673 1745107.5790 3923271.0275 -0.8 12.5 7.2 KIRL 4238828.4981 2180141.2539 4224130.5524 0 0 0 KISM 4761965.4177 2085815.4355 3682780.5406 0 0 0 KITH 4738508.9677 2012858.0408 3752992.1296 0 0 0 KLAM 4723123.9979 1918891.3749 3820063.7887 -0.1 -0.6 -0.1 KLOK 4564746.8948 1845610.9558 4040935.1667 -0.2 -0.3 0.1 KOM1 4346766.9895 2063591.2869 4172738.8095 0.2 0.4 0 KOMO 4346497.8852 2064660.5092 4172502.5725 0.3 0.5 -0.1 KOPA 4716817.3248 1888333.5416 3843136.5209 -1.1 -0.8 -0.3 KORI 4638456.3597 1962313.1952 3900363.8579 -0.3 -0.4 0 KORO 4742622.1149 1911936.9793 3799439.0772 0.9 -1.6 -0.6 KORU 4743405.2909 1908128.7613 3800588.8022 0.4 -0.8 -0.3 KOUN 4651613.2934 1883702.4852 3924106.2758 -0.8 -1.0 -0.1 KOZA 4523932.7689 1808094.7813 4103851.2459 -0.1 -0.1 0 KRDI 4580688.8367 1843528.2724 4023954.6984 -0.5 -0.6 0.2 KRIN 4655852.2402 1877291.6582 3922476.7978 -1 -1.1 0.1 KRP1 4628005.7263 2370338.0861 3681437.2862 -0.1 -0.9 -0.1 KRPS 4622868.7282 2371861.8970 3687744.3651 -0.1 -0.1 -0.1 KRUM 4314415.9685 2071931.2977 4202147.3804 0.3 1 -0.3 KRYO 4647807.3478 1936397.6871 3903591.3797 -0.7 -0.8 0 KTCH 4664075.2869 1813472.1490 3941385.5656 -1 -1.7 0 KTHA 4737986.6255 2017258.8610 3751411.6697 0 0 0 KTIM 4589632.7325 1865522.6492 4004201.7969 0 -0.1 0 KUKE 4443700.6345 1654086.7734 4252246.1469 0.1 0.0 0.0 KUST 4359748.8462 1824895.9882 4269354.2357 0 0 0 KYMI 4553884.2305 2037546.5212 3960901.9210 -0.4 -0.4 0.1 LAMA 4594430.3564 1896710.8415 3983618.1512 -0.6 -0.5 0.1 LAMB 4646524.2961 1874785.1310 3933433.5602 -0.6 -0.6 0.1 LAMI 4595888.6010 1894529.2056 3982860.6631 -0.6 -0.4 0.1 LARI 4547498.5850 1874350.9735 4047090.4741 -0.2 -0.3 0.1 LARM 4549397.1290 1874003.6492 4045167.9618 -0.2 -0.4 0.2 LEFK 4654193.1103 1758184.1658 3977667.2140 -0.8 -0.4 -0.1 LEMN 4434466.0611 2084864.4793 4069305.5210 -4.4 -9.8 -0.7 LEON 4689385.8337 1977142.7428 3831865.3529 -0.3 -0.3 0 LEPE 4643596.7482 1809529.1775 3967365.3381 -8.2 -3.2 1.6 LERO 4541890.9353 2299708.5316 3829500.7593 -0.2 0.4 0.8 LESU 4442031.0578 2209792.1946 3995066.5492 8.0 17.5 2.4 LESV 4433553.6763 2215687.6151 4000985.6645 3.1 6.9 0.8 LEUK 4653778.5321 1759300.8946 3977650.4440 -4.7 -5 2 LFKD 4653359.5632 1759483.0896 3978044.6946 -6.5 -7.5 1

56 LIBR 4508692.8290 1669107.8875 4177626.5209 0.1 0 0 LIDO 4626106.9712 1887975.4180 3951917.6800 -0.4 -0.6 0.1 LIMN 4440204.0751 2076246.0788 4067382.6631 -1.7 -3.8 0.2 LIVA 4609275.0812 1943760.3011 3943990.2562 -0.1 -0.1 0 LOVE 4234958.2049 1948619.1487 4338767.3956 0.1 0.1 0 LYGO 4655629.0699 1979810.7523 3871701.4074 -0.4 -0.3 0.2 MANT 4565396.7625 1983180.0766 3974825.9909 -0.5 -0.4 0.2 MAOR 4673991.9909 1830738.8543 3921748.5701 -0.5 0.3 0.5 MARM 4471213.1470 1970221.4997 4085990.7879 0 -0.3 0.2 MARU 4470998.5943 1970344.9829 4086179.3189 0.2 -0.2 0.1 MEGI 4483542.6960 2546274.4502 3741787.3841 -0.1 0 0 MENA 4646460.0785 2009359.9455 3866971.5710 0 0 0 MESA 4648874.2787 1839340.1046 3948013.2317 -1.3 -1.6 0.2 MESO 4659693.8174 1833149.0852 3937428.5676 -1.1 -1.6 0.0 MET0 4601878.2879 2026055.7098 3911285.3007 -0.4 -1.1 -0.3 METH 4749280.5141 1890407.1585 3801955.8268 0 -0.2 -0.1 MILO 4659071.5573 2116515.9551 3794601.0229 -0.1 -0.4 -0.1 MKMN 4642463.0390 2204528.6284 3764821.6853 0 -0.2 0 MLOS 4655589.1743 2123551.4587 3795046.2009 -0.1 0 0 MOIR 4742066.4292 2198956.2562 3642845.7649 -0.1 -0.3 0 MOLA 4711863.5316 1985602.6309 3800289.0137 -0.1 -0.3 -0.1 MOLU 4711870.0054 1985755.6331 3800182.0938 0 -0.3 0 MOUD 4432736.8130 2092390.7217 4067260.0140 -1.8 -1.7 0.1 MOZI 4639049.3902 2204980.2823 3769175.6722 -0.1 -0.2 0 MTSV 4578206.1540 1774174.1165 4059209.9276 -0.2 -0.2 -0.1 MUGL 4475458.9781 2416283.2852 3836973.7806 -0.1 0 0 MURG 4568077.9482 1804409.6053 4055914.6283 -0.1 -0.2 0 MURT 4751742.5362 2155790.0690 3656050.9084 0 -0.3 0 MYKN 4582963.4126 2169202.2630 3856468.9216 -0.3 -0.4 0 NAFP 4649596.4424 1857938.4714 3937758.5284 -1.1 -1.1 0.1 NAOU 4492475.6234 1821587.4821 4131598.8169 0 -0.1 0 NAXO 4601830.5469 2183256.6483 3826123.5606 -0.1 -0.4 -0.1 NAXU 4601796.7078 2182796.3023 3826435.6500 -0.2 -0.4 0.1 NEAB 4722484.8086 2010445.0131 3773780.0096 0 -0.3 -0.1 NEAP 4701748.2040 2253746.8138 3661764.1595 0 0 0 NEAU 4722475.2059 2010176.3278 3773940.1067 0 -0.3 -0.3 NEOS 4409468.4173 1930722.4602 4170777.0748 0.2 -0.1 0.1 NOA1 4599643.3342 2034828.0650 3909891.6831 -0.4 -0.4 -0.1 NOMI 4640783.4944 2206447.6642 3766127.0722 -0.1 -0.1 0 NVRK 4386262.7079 1940952.7242 4190907.0682 -1 1.9 5.1 ORES 4285376.0133 2140887.4609 4197120.1610 0.6 1.7 -0.5 OREU 4279034.1881 2136144.3445 4205870.7478 0.2 0.6 -0.1 ORID 4498451.4972 1708267.2983 4173592.0468 0 -0.1 0 OROP 4710648.7787 2245605.4826 3656293.7411 0 -0.1 0 PALC 4777538.1716 2091297.7819 3659632.0486 0 -0.2 0 PARO 4611210.6885 2165372.6877 3825002.3367 -0.2 -0.4 0.1 PAT0 4655002.7778 1860619.7634 3930281.8108 -0.7 -0.9 0 PAT1 4655020.7070 1860614.1048 3930261.8161 -3.9 -0.7 -2.1 PAT2 4656879.0354 1859254.0459 3928704.4499 -2.6 -2.6 0.1 PAT3 4660614.3390 1856776.2265 3925338.3093 -3.9 -3.9 -0.1 PATR 4659433.5436 1857319.0759 3926511.7439 -0.8 -1 0 PAXO 4645392.9561 1705853.4545 4010584.8614 -0.9 1.8 -0.6 PAZA 4309560.6625 1949464.6431 4264810.4596 0.2 0 0.1 PERI 4565086.6631 2045605.4033 3943763.0177 -0.4 -0.3 0 PERM 4571819.3914 1695983.9408 4097965.4691 0 0.1 0

57 PESH 4470863.1675 1665038.9431 4220035.2811 0.1 0 0 PIRG 4522824.7069 1881638.1910 4071243.7531 -0.2 -0.2 0.1 PKMN 4643869.0824 2203539.6573 3763744.1540 -0.1 -0.1 -0.1 POLI 4464666.1858 1935542.0644 4110273.8257 0.2 -0.2 0.2 POLY 4464686.1571 1936157.0847 4110050.1285 -0.1 0 0 PONT 4671272.5416 1754437.2551 3959389.5175 -22.2 -39.3 -4.2 PRKV 4435581.3201 2188830.5574 4013585.9556 -0.6 1.2 0.5 PROV 4420020.1595 1911521.1193 4168199.7087 0.5 -0.3 0.3 PSAR 4639531.0752 1891874.9223 3933580.1287 -0.6 -0.8 0.0 PSAT 4649321.8454 1866320.7902 3934160.1386 -1.4 -1.3 0.1 PTKG 4747739.8347 1965562.8675 3766013.1174 0.3 -0.6 -0.3 PTOL 4513036.6158 1793961.1610 4121755.8525 0.2 0.1 0.1 PVOG 4642285.9678 1830835.4296 3958994.5083 -1.8 -1.1 0.3 PYL1 4743818.9926 1891821.1150 3808253.3346 -3.4 -9.7 -2.9 PYLO 4744074.1990 1887446.3062 3809805.6410 0.2 -0.9 -0.2 PYRG 4703872.4104 1849321.6304 3877320.6698 -1.8 -2.2 -0.8 RAFI 4595861.4992 2047113.4224 3907152.5176 -0.6 -0.4 0.1 RETH 4732718.0313 2168114.4918 3673072.1244 0 -0.2 -0.1 RETS 4657391.3354 1827963.2164 3943088.8583 -1.4 -1.7 0.1 RIBA 4641874.6593 2198523.6984 3768942.3816 0 0 0 RGNI 4653408.5679 1808559.1044 3956194.0301 -2.1 -1.5 0.2 RLSO 4679939.0291 1840151.2799 3910407.6730 -0.8 -0.8 0.6 ROD3 4650274.3400 1868670.8252 3932631.1813 -1 -1 0.1 RODO 4537597.5093 2429130.5797 3754440.1724 0.3 -0.6 0.1 ROZH 4332980.7709 1996660.8797 4221598.9318 0.2 0 0 SALH 4409101.5292 2356570.0460 3947657.3691 -0.2 0.1 0.1 SAMO 4508287.4143 2267958.4859 3887335.3144 -0.5 0.1 -0.3 SAMU 4499749.7797 2290112.1315 3884235.7864 -1.3 1.2 -2.4 SAND 4391640.6319 1888416.4485 4208542.0087 0 0 0.1 SANT 4640364.4259 2205651.9684 3767241.4413 -0.1 -0.2 0 SARY 4231111.1589 2241812.5601 4199557.9863 0.1 0.2 -0.1 SERR 4417216.4162 1878405.9247 4186617.7370 0.5 -0.3 0.2 SIFN 4634604.5571 2135191.0982 3814041.6158 -0.2 -0.2 0 SISS 4702472.6167 1773110.2163 3914367.8584 1.7 -4.7 -2.6 SITI 4685118.3183 2295830.5066 3656471.0251 -0.1 -0.4 0 SIVA 4746676.8380 2194481.6266 3639548.4335 0.1 -0.1 0 SKOP 4536006.1238 1993792.8655 4002852.4821 -0.3 -0.8 -0.1 SKOU 4536091.3314 1994089.8806 4002602.0292 -0.4 -0.4 0.0 SKP1 4536100.6365 1994334.5626 4002472.0404 -0.4 -0.4 0 SKYR 4520263.9622 2066159.8002 3984144.0817 -0.7 -1.2 -0.2 SMLN 4341528.0878 1997939.1938 4211061.7300 0 0 0.1 SNDA 4390644.8737 1888820.2655 4209453.7998 0 -0.1 0.1 SNTR 4647218.3861 2202396.7792 3760306.7976 0 -0.1 0 SOFI 4319371.8301 1868688.0569 4292064.0639 0 0 0.1 SPAN 4658312.1002 1757780.8360 3973702.6921 -8.9 -7.7 -0.3 SPAR 4709762.4183 1943511.2108 3824424.4883 -0.2 -0.4 0 SPET 4672716.1990 1998423.5033 3841084.6270 -0.1 -0.3 -0.1 SPRT 4713881.7457 1940236.8166 3821960.6752 0 0 0 STEF 4547444.9237 1906215.0808 4032293.9464 -0.2 -0.4 0.1 STRA 4442959.4853 1958723.6538 4122226.4888 1.2 -0.5 0.5 STRF 4745440.9416 1823089.1702 3839131.6786 1.7 -1.4 -1.8 SVI1 4274925.4826 2104754.6097 4225699.2551 0 0.2 0 SYR1 4597678.3367 2138015.8530 3856344.1041 -0.3 -0.3 0 SYRO 4597169.1470 2138292.8297 3856806.4933 -0.3 -0.3 0 TEI1 4414135.2642 1924030.5409 4168706.6841 -0.2 -0.3 0

58 TEIS 4414007.1274 1924331.8077 4168705.4259 0 0 0 TEKR 4278776.5676 2227057.9747 4158962.1136 0 0.4 -0.1 THAS 4404239.0615 2010414.6202 4138456.2937 0.7 -0.2 0.2 THES 4472857.7900 1892082.3186 4120761.2749 0.1 -0.2 0.1 THI1 4601697.9587 1983418.3720 3933201.2546 -0.8 -0.9 0 THIR 4642084.2031 2209438.0915 3762870.5860 -0.1 -0.4 -0.2 THIU 4601397.7384 1983528.5635 3933489.4022 -0.5 -0.4 0.1 THIV 4601629.8249 1983525.5307 3933234.5061 -0.4 -0.5 0 THS1 4463674.5051 1890981.9396 4131173.8594 0.3 -0.3 0.3 TILO 4564714.3614 2365481.9481 3762366.5774 -0.1 -0.1 0 TINU 4583046.9485 2153543.0674 3865062.8813 0 0 0 TRIK 4572319.3270 1828602.9351 4040134.0081 -0.2 -0.4 0.1 TRIP 4685169.2148 1928160.3127 3862714.9723 -0.8 -1.3 2.5 TRIU 4685194.9932 1927920.6227 3862822.4076 -2.3 -2.6 7.8 TRIZ 4640396.8967 1881698.2903 3937360.0884 -0.7 -0.6 0.1 TROP 4684940.7715 1888987.4318 3882268.3726 -2.2 -0.5 -1.8 TUC2 4744543.8434 2119412.0506 3686258.7205 -0.4 -0.6 -0.1 VALI 4646733.3842 1890136.2542 3925882.0935 -0.9 -0.9 -2 VALY 4719043.7730 1905038.1714 3831960.3917 -0.2 -0.7 -0.1 VAM0 4747348.2515 2133418.5625 3674795.7775 0 -0.2 -0.1 VASS 4706041.5069 1891746.0853 3856151.1910 0 0 0 VERI 4495880.4961 1838353.3913 4120155.2145 0.2 -0.2 0.1 VLSM 4699991.5545 1765547.8674 3921162.2947 -2.9 -4.4 -0.9 VOLO 4549875.3166 1920748.9819 4022789.4949 -0.2 -0.4 0.1 VOLU 4565221.5349 1915377.9723 4007961.3457 -0.3 -0.3 0.1 VOUR 4604202.5048 1784424.3871 4024774.8621 -0.6 -0.4 0.1 VRAS 4433123.2002 1946359.8975 4138235.7969 0.2 -0.2 0.1 WNRY 4642530.8627 2208076.5946 3763086.9613 0 -0.2 0 XALK 4581712.0101 2002387.0708 3946497.0255 -0.6 -0.4 0.1 XANT 4366480.1686 2022971.9081 4172140.9631 0.2 0.1 0.1 XILI 4644380.1014 1868139.2762 3939073.1413 -1.0 -1.0 0.1 XRSO 4777260.3404 2080454.6045 3666053.2612 0 0 0 YENC 4354314.9762 2241822.1827 4072748.6380 -1.2 0.5 0.4 YSST 4225148.1215 2341157.9537 4151211.7118 -0.2 0.2 0.1 YUND 4338116.8010 1918205.9923 4251887.8045 0 -0.1 0 ZAGO 4536866.3880 1935307.4955 4031029.8680 -0.6 -0.5 0.1 ZAKU 4715769.1725 1799745.5741 3886137.1326 -12 -13.8 0.1 ZAK1 4715987.8731 1798769.5154 3886307.2524 0 -0.1 -0.1 ZAKY 4715888.8271 1799412.8017 3886130.4757 -5 -8.2 -1.4 ZKRO 4685994.9539 2307519.3947 3648460.3418 -0.1 -0.4 -0.1 028Α 4720317.0119 1797293.0719 3881778.3162 0.9 -2.1 -0.6 029Α 4710031.0674 1869466.2358 3860256.7488 0.4 -0.8 -0.3 030Α 4695504.5027 1827746.7423 3897484.6561 1 -1.9 -0.6 959

960

961 FIGURES 962

59 963 964

965 Figure S1. Greece and surroundings with the location of the twenty one earthquakes listed in Tab.

966 1, and the seismicity from the NOA catalogue for the period 2000-2020, magnitude ≥ 4 and depth ≤

967 35 km. The beachballs show the focal mechanism of the studied events, they are centred at the

968 location of the epicentres.

969

60 970

971 972

973 Figure S2. Time series of the north coordinate of the GPS station LEMN. The large offset

974 corresponds to the nearby Samothraki earthquake of May 24, 2014 (Tab. 1, Tab. 2). After the

975 earthquake the velocity of the station has changed. The transient can be modelled with a sum of two

976 exponential terms of time constant 65 and 1295 days respectively.

977 978

61 979 980 981 Figure S3. Additional toponymy for features not mentioned in the main text but in the 982 supplementary material (coordinates and description of the sites in Tab. S2). 983 984

62 985 986 987 Figure S4. Duration of the time series of GPS data for the 329 stations. The average duration is 6.5 988 years. 989 990

991 992 993 Figure S5. Completeness of the time series acquired at the 329 GPS stations. The average for the 994 whole set of stations is 85%. 995

63 996 997 998 Figure S6. Comparison of our solution with the solution of the Nevada Geodetic Laboratory 999 (NGL), here for the time series of the station LEMN. The curves are de-trended from the modelled 1000 effect of the earthquakes and the average secular velocity. The difference between the slopes is 1001 below 0.05 mm yr-1. See the native time series of the north component of LEMN in Fig. S2. 1002

1003 1004 1005 Figure S7. Histogram of east and up velocity uncertainties at the 329 GPS stations.

64 1006

1007 Figure S8. Time series of the RLSO station showing the coseismic offset of the Mw = 6.4 June 8, 1008 2008 Movri earthquake of amplitude -3.4 mm in east and 9.8 mm in north, and the postseimic 1009 deformation with relaxation time of 103 days and amplitudes -1.1 and 1.4 mm in east and north. 1010

1011 1012 Figure S9. Time series of the stations MET0, AGTH, AIGU and ATHI. The vertical line

1013 corresponds to the Mw = 5.3 Magoula earthquake of July 19, 2019. At AIGU and ATHI we measure 1014 an offset at the date of the earthquake with amplitudes that are reported in Tab. S4.

65 1015

1016 Figure S10. Residual time series of various stations corrected for coseismic, postseismic and 1017 secular terms (top) and common mode noise correction (bottom) estimated using the eight top

1018 stations (MET0 to PAT0). The solid vertical line is on October 25, 2018, day of the M w = 6.8 1019 Zakynthos earthquake. The dashed line corresponds to a transient visible at the stations AMAL and 1020 PYRG on May 10, 2018 (see S4.1 and see Fig. S11 for the north and up components).

66 1021

1022 Figure S11. North and up components of the same stations whose east component is plotted in Fig.

1023 S10. The solid vertical line is on October 25, 2018, day of the M w = 6.8 Zakynthos earthquake. The

1024 dashed line corresponds to the transient observed in Fig. S10 at the stations AMAL and PYRG on

1025 May 10, 2018. Here, a subsidence is observed at TROP, PYRG and AMAL around May 10, 2018.

67 1026

1027 Figure S12. Time series of the east and up components of the station NOMI used to assess the

1028 linear character of the deformation during unrest and the starting and ending dates of the volcanic

1029 unrest. We see also a large annual term of the vertical component that was modelled and removed.

68 1030

1031 Figure S13. Time series of the east component of the GPS stations of Santorini. (a) Original time

1032 series, (b) time series corrected from a linear signal starting on February 12, 2011 and ending on

1033 May 3, 2012 (446 days). See Tab. S5 for the amplitude of the slope at the stations.

69 1034

1036 Figure S14. Time series of the vertical coordinate of the GPS stations VALI and XILI. The secular

1037 trends are 2.4 mm yr-1 for VALI and -6.2 mm yr-1 for XILI. The Zakynthos earthquake of October

1038 2018 (vertical line) produced a subsidence of the two stations that we model with a decreasing

1039 exponential starting the day of the earthquake with time constant 36 days for VALI and 16 days for

1040 XILI and amplitude 10.2 and 9.6 mm respectively.

1041

1042

70 1043 1044 1045 Figure S15. Sentinel-1 interferogram of the western Gulf of Corinth made with the scenes acquired 1046 on October 20 and November 13, 2018 on the descending track 80, using DIAPASON InSAR 1047 Sentinel-1 TOPSAR version 1.1.8 on the GEP-TEP (http://geohazards-tep.eu). One cycle of colour 1048 corresponds to 28 mm of displacement along the line of sight. The geo-information layer is fom 1049 OpenTopoMap (http://opentopomap.org). 1050 1051

71 1052 1053 1054 Figure S16. InSAR ground velocities in the region of Larissa (Thessaly) in the line of sight of the 1055 Sentinel-1 descending track 7, measured for the period October 16, 2019 to May 31, 2010, using the 1056 CNR-IREA P-SBAS processor on the GEP-TEP (Berardino et al. 2002; Casu et al. 2014). Red dots 1057 are permanent GPS stations. Blue colour on the map indicates range increase from -5 cm yr-1 (uplift) 1058 to 5 cm yr -1 (subsidence). 1059 1060

72 1061

1062

1063 Figure S17. Time series of the station TRIP. On January 12, 2019 (vertical line) a transient

1064 deformation starts with a sudden motion of -1.9 mm in east, 0.2 mm in north and 0.0 mm in up,

1065 followed by a progressive deformation that we model with an exponential of time constant 17 days

1066 and amplitude -3.3 mm in east, -9.7 mm in north and 29.7 mm in up.

1067

1068

73 1070 Figure S18. Sentinel-1 interferogram of the Tripoli area made with the scenes acquired on 1071 December 31, 2018 and January 24, 2019 on the ascending track 175. The two GPS stations TRIP 1072 and TRIU (time series in Fig. S17) are located in the centre of Tripoli. 1073 1074

74 1075

1076

1077 Figure S19. Time series of the horizontal coordinates of the GPS stations KALU and KALY,

1078 Kalymnos Island. The vertical line indicates the date of the Lesvos 2017 earthquake. The amplitude

1079 of the offset is -1.7 mm and 3 mm in east and north for KALU, and 0.2 mm and 3.5 mm for KALY.

1080 The relaxation time of the signal is 16 days. See also Fig. S20.

1081

1082

75 1083 1084 1085 Figure S20. Sentinel-1 interferogram of the Kalymnos Island made with the scenes acquired on 1086 May 25 and June 30, 2017 on the descending track 36. The vectors correspond to the offset plotted 1087 in Fig. S19. 1088 1089

76 1090 1091 1092 Figure S21. Time series of the three components of the vector KRIN (located on a landslide) - 1093 KOUN (the closest permanent station, considered here as a local reference). See also Fig. S22 and 1094 Fig. S23. 1095 1096

77 1097 1098 1099 Figure S22. Ground velocities in the line of sight of the Sentinel-1 ascending track 175 over the 1100 area of KRIN calculated for the period January 1, 2017 to May 15, 2020, using the TRE- 1101 ALTAMIRA FASTVEL processor on the GEP-TEP. The amplitude of the velocity of the station 1102 KRIN is 101.5, -21.2 and -18.7 mm yr-1 in east, north and up respectively. The approximate 1103 boundary of the landslide (according also to the interferogram made with the CNR-IREA P-SBAS 1104 processor, see Fig. S23) is plotted as a light red contour. Red indicates motion away from the 1105 satellite. 1106 1107

78 1108 1109 1110 Figure S23. InSAR ground velocities around the station KRIN, represented in the line of sight of 1111 the Sentinel-1 ascending track 175 for the period January 1, 2017 to May 15, 2020, using the CNR- 1112 IREA P-SBAS processor on the GEP-TEP. Red indicates motion away from the satellite. 1113 1114

79 1115 1116 1117 Figure S24. Time series of the horizontal coordinates of the EUREF (https://www.epncb.oma.be/)

1118 station PAT0. Top is corrected for the coseismic and postseismic effects of the various earthquakes,

1119 bottom is not. For example, the Zakynthos 2018 earthquake is clear in the uncorrected east

80 1120 component. The difference in slope between the two solutions is -0.67 mm yr-1 in east and -0.85 mm

1121 yr-1 in north. East and north ITRF2014 velocity of PAT0: EUREF: 8.35 mm yr-1, -6.59 mm yr-1; our

1122 uncorrected velocity: 7.43 mm yr-1, -6.65 mm yr-1; our corrected velocity: 8.1 mm yr-1, -5.8 mm yr-1.

1123 1124

1125 1126 1127 Figure S25. Istanbul profile. Velocities plotted with X axis along the line Istanbul - pole AT/EU and

1128 Y profile perpendicular. The solid and dashed line represents the best fitting linear increase of the

1129 velocity (using the five points in red) as a function of the distance to the pole. The corresponding

1130 rate is 27.3 10-9 rad yr-1 counter clockwise, i.e. 1.56° Ma-1.

1131

1132

81 82 1134 Figure S26. Velocities along a profile Sandanski (Bulgaria) - Rodos. We plotted of a sub-set of 1135 stations located at less than 55 km from the profile. The axis is aligned with the azimuth of the 1136 Euler pole of rotation of Anatolia with respect to Eurasia, located at 31.96°E and 32.02°N according 1137 to Le Pichon & Kreemer (2010). 1138 1139

1140 1141 1142 Figure S27. Sandanski profile (location of the points and vectors in Fig. S26). Velocities of the

1143 stations located in the vicinity (± 55 km) of the line Sandanski-Rodos aligned with the Euler pole of

1144 rotation of Anatolia with respect to Eurasia according to Le Pichon & Kreemer (2010). The oblique

1145 solid line shows the tangential velocities that are expected for a rigid rotation of Anatolia around the

1146 pole plotted in Fig. 7. The grey dots show the values at LIMN, LEMN and MOUD after correcting

1147 for the expected component of strain accumulation due to the vicinity of the North Anatolian fault.

1148 1149 1150

83 1151 Figure S28. Velocities along a profile Tirana-Karpathos. We plotted of a sub-set of stations located 1152 at less than 55 km from the axis connecting Tirana and Karpathos aligned with the Euler pole of 1153 rotation of Anatolia with respect to Eurasia located at 31.96°E and 32.02°N according to Le Pichon 1154 & Kreemer (2010).

84 1155 1156 Figure S29. Tirana-Karpathos profile. Velocities of the stations located in the vicinity (± 55 km) of

1157 the line. See location of the points and line and the vectors in Fig. S28.

1158

1159 1160 Figure S30. Section of the south Peloponnese Pylos-Spetsopoula

85 1161

1162

1163

1164 Figure S31. East and north velocities in Crete with respect to average.

1165

1166

86 1167

1168

1169 Figure S32. North-south profile across central Macedonia (see points in Fig. 10).

1170 1171 1172

87 1173 1174 1175 Figure S33. Cumulative distribution of the whole set of 322 vertical velocities (grey) and the subset 1176 of 177 (black) with velocity uncertainties lower than 1 mm yr-1 (see Tab. S1). Best Gaussian fits are 1177 calculated for 90% of the population, which corresponds to velocities comprises between -2.7 and 1178 1.4 mm yr-1 in the case of the 177-stations distribution. The remaining 6% on the negative velocities 1179 side and 4% on the positive velocities side (drawn with + symbol) do not fit with a Gaussian 1180 distribution and were not used to calculate the best fit. The median µ of the distribution is -0.38 mm 1181 yr-1 with the best 177 stations and -0.27 mm yr-1 for the whole set of 322 stations. The 1182 corresponding standard deviation σ is 0.79 and .95 mm yr-1, respectively. Two stations (GEYB and 1183 PSAT) present anomalous subsidence discussed in the text, GEYB (mentioned in S4.7) is built on a 1184 recent thick (>10 m) embankment built during the construction of the bridge Rio-Antirio, PSAT is 1185 located near (~1 km) the fault plane of the active Psathopirgos fault, on its footwall. Three stations 1186 present anomalous uplift, ZAKY and AMAL presumably affected local uplift (insufficiently 1187 modelled by our post-seismic model, see S3.14) during the post-seismic relaxation following the 1188 2018 Zakynthos earthquake, VALI affected by the uplift of the Selinoutas river delta (see S4.3) as 1189 shown in Fig. S14, induced by the activity of the offshore fault north of Aigion. 1190

88 1191 1192

1193 Figure S34. The boundaries of the ten blocks of our model. The North Aegean block is the

1194 combination of two sub-blocks one for the north (Lemnos and Lesvos), one for the south (Chios,

1195 Samos). Here we use one single block (despite the large north-south extension and the localization

1196 zone between Lesvos and Chios) because we do not have yet, at the moment, enough GPS stations

1197 in the area and long enough time series to divide the zone in two. The focal mechanisms are for

1198 events located between 0 and 35 km of depth, published by the National and Kapodistrian

1199 University of Athens (http://www.geophysics.geol.uoa.gr/) for the period 2003-2018.

89