According to previous studies (e.g.,
[3][6]), the Hellenic Trench has undergone a retreat (approx. 580 km) since the Eocene–Miocene
[7][8]. This geological phenomenon has resulted in a tectonic regime, which is characterized by N–S and NE–SW extensional occurrences in the Aegean and Anatolian regions
[9][10]. The southwest Anatolian microplate motion is considered to be related to the exertion of pulling and pushing forces resulting from the HSZ in its western part
[5]. This movement is affecting the lithospheric thickness across Anatolia and the Aegean
[11][12]. The Hellenic Trench (HT) retreat results in the extensional stress field generation, causing the crustal thinning of central and northern Greece and western Anatolia (
[13][14] and references therein).
Since it is clear to the scientific community that satellite geodesy can provide detailed information about the motions of plate tectonics and the deformation of the crust, it is widely applied. Particularly, the Aegean and the eastern Mediterranean constitute an ideal laboratory site for geodetic analyses, and especially for the determination of strain
[15][16][17][18][19][20][21]. Therefore, various studies have been conducted in the Aegean region by implementing measurements recorded by either permanent or campaign GPS/GNSS stations
[22][23][24][25][26][27]. While preceding seismological studies only implied the gross geodynamic setting of the East Mediterranean (the relation between the Arabian, Nubian and Eurasian plates), the geodetic studies confirmed and measured plate boundaries and motions. Furthermore, they resulted in the separating of the region into several blocks, proposing potential geodynamic models (e.g.,
[28][29][30][31]), which are complicated. Focusing on the Aegean, two distinct velocity clusters are observed throughout the area, distributed on both sides of the NAT (and its continuation) and considered as the boundary between the North and South Aegean
[32][33]. The most recent studies predominantly focus on developing crustal block modelling (e.g.,
[31][34][35][36]), attempting to interpret their geodynamic activity, as well as their relationship to seismic events (e.g.,
[37][38][39][40]). Notably, the block modelling considers that each block partially receives the brittle tectonics effect
[31][41][42], caused by tectonically active structures, while it is imprinted on the upper-crust strain. This North Aegean region strain is documented on the seismic slip and the gravitational collapse of the Aegean crustal wedge in response to slab roll-back at the Eurasian–Nubian convergence boundary
[43][44][45], while in tandem, strike-slip faulting dominates in this region
[46][47][48].
2. Tectonic Geodesy Synthesis of the North Aegean Region
2.1. Synthesis of Strain Results
2.1.1. Maximum (MaHE) and Minimum (MiHE) Horizontal Extension
The MaHE and MiHE for the North Aegean Sea region are characterized by a wide range of strain field values. In particular, these parameters demonstrate a heterogeneous strain distribution within the North Aegean Sea region, as depicted in
Figure 2a. The MaHE values range from approximately −50 to 1000 nano-strain, with the highest value observed in the Thrace region. Similarly, the MiHE values range approximately from −1250 to 70 nano-strain
[49]. The North Aegean Sea is predominantly characterized by shear, with comparable values of MaHE and MiHE. Moreover, the MaHE and MiHE exhibit an increase in the northwestern region of Turkey, showing a slight NE–SW (MaHE orientation) extension dominance. These results indicate that the crustal dynamics of the North Aegean Trough (NAT) are more complex ranging from 100 to 700 nano-strain, particularly under a shear setting. At the same time, its extension is oriented in an NNE–SSW to N–S direction
[49]. The shear occurrence can also be observed in the North Aegean Basin (NAB), characterized by moderate magnitudes and a gradual transition in an approximately NW–SE extension direction. The NAB exhibits a prevalence of local extension dominance over compaction. The area exhibits a distinct extensional trend in contrast to the adjacent regions, which follows an ENE–WSW direction
[49]. The strain regime prevalent in the northern Aegean region, except for the wider Skyros Island area, is the transition from pure to quasi-pure extension.
Figure 2. Maximum (MaHE; blue arrows) and Minimum (MiHE; red arrows) Horizontal Extension in the regions of: (a) North Aegean Sea; (b) Strymon basin; and (c) Thessalian basin.
A broad spectrum of strain field values is exhibited in the Strymon basin region. In particular, the MaHE values show a range of variability from approximately −60 to 2650 × 10
−9, whereas the corresponding MiHE values demonstrate a range of variability from approximately −650 to 150 × 10
−9, as depicted in
Figure 2b. The primary analysis of the distribution of the values
[46] indicates that the absolute maximum values of both MaHE and MiHE are primarily concentrated along the eastern boundary and the northernmost part of the western boundary of the Strymon basin. Moreover, the analysis of MaHE and MiHE suggests that extension is the dominant process in the study area, resulting in a normal faulting activity. In particular, the Strymon basin exhibits a notable concentration of maximum values at its eastern and northwestern boundaries. At the same time, the high MaHE values of these regions, coupled with exceedingly low absolute values of corresponding MiHE, suggest the prevalence of dip-slip normal faulting activity
[46]. The documented data reveals a gradual decrease in both MaHe and MiHE values on either side of the area. Simultaneously, the difference in values between MaHE and MiHE decreases, with MaHE values remaining greater than the corresponding MiHE values. These findings suggest oblique-slip normal faulting activity in the Strymon and the adjacent (e.g., Drama) basins. Additionally, the northwestern part of the Strymon basin exhibits positive values for both MaHE and MiHE, indicating a prevalent, pure extensional setting characterized by normal faulting. Finally, it is worth mentioning that the direction of the MaHE vectors is nearly perpendicular to the strike of the faults and fault zones, thereby providing evidence for the existence of these tectonic structures and their normal faulting characteristics.
Focusing on the wider Thessalian basin region, the values show a range of −40 to 380 nano-strain, with the extensional regime prevailing as the dominant one (
Figure 2a,b)
[50]. The MaHE vectors exhibit an approximate N–S orientation in the southern Thessalian basin. This region also displays the most significant values, exceeding 200 nano-strain. A notable alteration in direction is discernible towards the westernmost Thessalian basin, while similar noticeable changes are documented in the MaHE direction of the northernmost part
[50]. Remarkable values exceeding 150 nano-strain are observed within the Larissa plain and the Thessalian basin southern boundary. The findings of the MaHE study indicate a progressive rise of tectonic activity in the central and eastern Thessalian basin.
2.1.2. Maximum Shear Strain (MaxSS)
The MaxSS parameter exhibits a wide distribution of values within the North Aegean Sea region, as depicted in
Figure 3. These values range from 10 to 1650 nano-strain, with both peak values observed in the Samothraki Plateau. The regions that show the highest values are documented in the eastern part of the Aegean Sea (western Turkey), the CAT (Central Aegean Trough) and the NAT, as shown in
Figure 3a. The MaxSS is characterized by slightly reduced values within the vicinity of the NAB, in the transitional zone that links the NAB and the CAT, and the broader region of northwestern Turkey. The regions showing great values are primarily associated with shear faulting. The extension may exhibit a greater magnitude in regions with possible oblique-to-normal dip-slip faulting. An example of such a scenario can be observed in the eastern Aegean Sea region, suggesting the occurrence of normal and oblique-slip faulting. Local kinematic variations (from pure normal dip-slip to strike-slip faulting) can be observed in regions with moderate values.
Figure 3. Maximum Shear Strain (MaxSS) in the wider North Aegean region: (a) MaxSS values (red circles) of all triangles, constructed within the three study areas; (b) MaxSS values grid, based on the interpolation (kriging) geostatistical analysis. Dark lines in (b) are faults.
Concerning the Strymon basin, the MaxSS value analysis shows a range of approximately 1 to 2750 × 10
−9 nano-strain (
Figure 3a). The highest MaxSS values are recorded in the eastern region of the Strymon basin, with comparable values noted in the northwestern margin
[46], respectively, as depicted in
Figure 3a. Furthermore, high MaxSS values are documented in the southern region of the western boundary of the Strymon basin. The results of the geostatistical analysis, as presented in
Figure 3b, provide the distribution of MaxSS. In particular, the analysis reveals that higher values of MaxSS are concentrated along the WNW–ESE axis of the Strymon basin
[46]. Notably, the adjacent basins (such as the Drama basin) exhibit exceedingly low, nearly negligible values. Consequently, the significant shear values highlight the tectonic structures that are actively distributed along the boundaries of the Strymon basin, with particular emphasis on the eastern one.
Finally, the Thessaly area exhibits a range of MSS values from 5 to 730 nano-strain. According to the analysis of MSS values (
Figure 3a,b), the regions with the highest values, exceeding 400 nano-strain, are adjacent to the western boundary of the Thessalian basin
[50]; additionally, high MSS values are observed in the southern part of this region. Moreover, intermediate values (approximately 200–400 nano-strain) have been detected all along the western margin, while the northern part is characterized by low MSS values, measuring less than 100 nano-strain.
2.1.3. Area Strain (AS)
Two distinct AS types are distinguished: (i) dilatation, which is directly related to extension, and (ii) compaction, which is directly related to compression. Compaction is the predominant form of strain in the North Aegean Sea region, while the nano-strain values for AS range from −1180 (compaction) to 970 (dilatation) (
Figure 4a). The dilatation encompasses a larger region (
Figure 4b); the most significant values are documented in the eastern and northeastern parts of the North Aegean Sea. Samothraki Plateau and western NAB have smaller but notable dilatation values. Most compaction occurs in the southerly area of the North Aegean Sea (NAB and NAT). In particular, the CAT exhibits remarkably low values. Compaction is limited to the north of the NAB and NAT, showing a slight presence in the northeast of the North Aegean Sea and along the eastern coast of Thessaly. In a case of a vast contrast between dilatation and compaction, extensional or compressional tectonics occur, respectively. However, shear tectonics are documented in a case of similar dilatation and compaction values. For instance, the broader NAT region is characterized by high dilatation and compaction values, defining a more complex setting, as the two strain types do not extend uniformly or equally.
Figure 4. Area Strain (AS) in the wider North Aegean area: (a) Dilatation (blue circles) and compaction (red circles) calculated values; (b) dilatation values grid, based on the interpolation (kriging) geostatistical methodology; and (c) compaction values grid, based on the interpolation (kriging) geostatistical methodology. Dark lines in (b,c) are faults.
Accordingly, the AS results of the Strymon basin reveal a dilatation dominance over compaction (
Figure 4a), with values ranging from approximately −600 (compaction) to 2750 × 10
−9 (dilatation). Even though the dilatation type dominates in the Strymon basin (
Figure 4b), it partially coexists with compaction, indicating shear occurrence
[46]. Specifically, the greatest concentrations are only documented within the Strymon basin area, especially along the eastern margin and the northwestern part of the Strymon basin. Concerning the high-value concentrations of compaction (
Figure 4b), they are predominantly documented at the eastern Strymon basin boundary
[46]. The difference between dilatation and compaction defines extension (dilatation > compaction) or compression (compaction > dilatation), whereas comparable values indicate shear tectonics
[49]. The Strymon basin is predominantly characterized by dilatation over compaction; nevertheless, the compaction occurrence is noteworthy, indicating that despite the dilatation dominance, a remarkable shear component is documented, leading to oblique-slip, normal faulting. A typical example is the southernmost part of the Strymon basin, which is adjacent to the North Aegean Sea and potentially affected by the corresponding strike-slip tectonic regime.
Regarding the Thessalian basin area, it exhibits AS values that vary from −420 to 360 nano-strain
[50], implying occurrences of both dilatation and compaction
(Figure 6a,b). The phenomenon of compaction is primarily spotted at the western boundary of Thessalian basin, where the absolute strain (AS) exceeds 100 nano-strain. However, there are also some areas in the central part, where the absolute strain ranges between 50 and 150 nano-strain. The dilatation shows a greater dispersion when higher values are considered, specifically with values exceeding 200 nano-strain in the southeastern Thessalian basin
[50]; the maximum dilatation is observed to the south of the Thessalian basin. Dilatation values ranging from 100 to 200 nano-strain have been observed in the northwestern part of the Thessalian basin, while at its northern part, small values (approximately 50 nano-strain) for both compaction and dilatation are documented.