Three strong earthquakes ruptured the northwest Thessaly area, Central Greece, on the 3, 4 and 12 March 2021. Since the area did not rupture by strong earthquakes in the instrumental period of seismicity, it is of great interest to understand the seismotectonics and source properties of these earthquakes. We combined relocated hypocenters, inversions of teleseismic P-waveforms and of InSAR data, and moment tensor solutions to produce three fault models. The first shock (M_w = 6.3) occurred in a fault segment of strike 314° and dip NE41°. It caused surface subsidence −40 cm and seismic slip 1.2–1.5 m at depth ~10 km. The second earthquake (M_w = 6.2) occurred to the NW on an antithetic subparallel fault segment (strike 123°, dip SW44°). Seismic slip of 1.2 m occurred at depth of ~7 km, while surface subsidence −10 cm was determined. Possibly the same fault was ruptured further to the NW on 12 March (M_w = 5.7, strike 112°, dip SSW42°) that caused ground subsidence −5 cm and seismic slip of 1.0 m at depth ~10 km. We concluded that three blind, unknown and unmapped so far normal fault segments were activated, the entire system of which forms a graben-like structure in the area of northwest Thessaly.

On 30 October 2020, an Mw = 7.0 earthquake struck the eastern Aegean Sea. It triggered earthquake environmental effects (EEEs) on Samos Island detected by field surveys, relevant questionnaires, and Interferometric Synthetic Aperture Radar (InSAR) analysis. The primary EEEs detected in the field comprise coseismic uplift imprinted on rocky coasts and port facilities around Samos and coseismic surface ruptures in northern Samos. The secondary EEEs were mainly observed in northern Samos and include slope failures, liquefaction, hydrological anomalies, and ground cracks. With the contribution of the InSAR, subsidence was detected and slope movements were also identified in inaccessible areas. Moreover, the type of the surface deformation detected by InSAR is qualitatively identical to field observations. As regards the EEE distribution, effects were generated in all fault blocks. By applying the Environmental Seismic Intensity (ESI-07) scale, the maximum intensities were observed in northern Samos. Based on the results from the applied methods, it is suggested that the northern and northwestern parts of Samos constitute an almost 30-km-long coseismic deformation zone characterized by extensive primary and secondary EEEs. The surface projection of the causative offshore northern Samos fault points to this zone, indicating a depth–surface connection and revealing a significant role in the rupture propagation.

Planned decommissioning of coal-fired plants in Europe requires innovative technical and economic strategies to support coal regions on their path towards a climate-resilient future. The repurposing of open pit mines into hybrid pumped hydro power storage (HPHS) of excess energy from the electric grid, and renewable sources will contribute to the EU Green Deal, increase the economic value, stabilize the regional job market and contribute to the EU energy supply security. This study aims to present a preliminary phase of a geospatial workflow used to evaluate land suitability by implementing a multi-criteria decision making (MCDM) technique with an advanced geographic information system (GIS) in the context of an interdisciplinary feasibility study on HPHS in the Kardia lignite open pit mine (Western Macedonia, Greece). The introduced geospatial analysis is based on the utilization of the constraints and ranking criteria within the boundaries of the abandoned mine regarding specific topographic and proximity criteria. The applied criteria were selected from the literature, while for their weights, the experts’ judgement was introduced by implementing the analytic hierarchy process (AHP), in the framework of the ATLANTIS research program. According to the results, seven regions were recognized as suitable, with a potential energy storage capacity from 1.09 to 5.16 GWh. Particularly, the present study’s results reveal that 9.27% (212,884 m²) of the area had a very low suitability, 15.83% (363,599 m²) had a low suitability, 23.99% (550,998 m²) had a moderate suitability, 24.99% (573,813 m²) had a high suitability, and 25.92% (595,125 m²) had a very high suitability for the construction of the upper reservoir. The proposed semi-automatic geospatial workflow introduces an innovative tool that can be applied to open pit mines globally to identify the optimum design for an HPHS system depending on the existing lower reservoir.

Coastal environments are highly recognized for their spectacular morphological features and economic activities, such as agriculture, maritime traffic, fishing, and tourism. In the context of climate change and the evolution of physical processes, the occurrence of intense natural phenomena adjacent to populated coastal areas may result in natural hazards, causing human and/or structural losses. As an outcome, scientific interest in researching and assessing multi-hazard susceptibility techniques has increased rapidly in an effort to better understand spatial patterns that are threatening coastal exposed elements, with or without temporal coincidence. The islands of Milos and Thira (Santorini Island) in Greece are prone to natural hazards due to their unique volcano-tectonic setting, the high number of tourist visits annually, and the unplanned expansion of urban fabric within the boundaries of the low-lying coastal zone. The main goal of this research is to analyze the onshore coastal terrain’s susceptibility to natural hazards, identifying regions that are vulnerable to soil erosion, torrential flooding, landslides and tsunamis. Therefore, the objective of this work is the development of a multi-hazard approach to the South Aegean Volcanic Arc (SAVA) islands, integrating them into a superimposed susceptibility map utilizing Multi-Criteria Decision-Making (MCDM) analysis. The illustrated geospatial workflow introduces a promising multi-hazard tool that can be implemented in low-lying coastal regions globally, regardless of their morphometric and manmade characteristics. Consequently, findings indicated that more than 30% of built-up areas, 20% of the transportation network, and 50% of seaports are within the high and very high susceptible zones, in terms of the Extended Low Elevation Coastal Zone (ELECZ). Coastal managers and decision-makers must develop a strategic plan in order to minimize potential economic and natural losses, private property damage, and tourism infrastructure degradation from potential inundation and erosion occurrences, which are likely to increase in the foreseeable future.

The 27 September 2021 damaging mainshock (M_w6.0) is the first known strong earthquake that ruptured the Arkalochori area, Crete Isl., Greece, during the entire historical period, making it an unexpected event in the long-term sense. The area is characterized by the presence of the normal active Kastelli Fault (KF) striking NNE-SSW and dipping towards ~WNW. The KF, of surface exposure only ~6 km, at its southern tip is truncated by the nearly perpendicular active Nipiditos fault. The main shock was preceded by foreshock activity lasting for ~3.9 months, thus the mainshock turned out to be an expected event in the short-term sense. Maximum ground subsidence of ~20 cm was estimated from InSAR images, but this also incorporates deformation that may have been caused by the largest aftershock (M_w5.1) of 28 September 2021. The fault model produced from the inversion of InSAR observations indicated strike 216°, dip towards ~NW at angle 53°, rake −95°, and is consistent with fault-plane solutions obtained from routine moment tensor analysis. The geodetic seismic moment calculated from the Okada’s formalism is 1.14 × 10¹⁸ N·m (M_w6.0), while a maximum slip of 1.03 m was found at depths from 3.5 km to 5 km. The entire aftershock epicenters cloud strikes in a ~SW-NE direction but is distributed in two clusters, the southern and the northern ones. The foreshock cloud, the main slip patch, the deformation area, and the strongest aftershocks all fall within the southern cluster. The foreshocks concentration at the deepest edge of the main slip patch was a foreshadow of the mainshock nucleation area. The northern cluster, which is very likely due to the gradual expansion of aftershocks, is situated in the KF hanging wall block. To interpret the main seismic slip in the southern cluster area we propose the existence of a buried KF segment at the SSW-wards prolongation of the emerged at the surface segment. Assuming a rectangular seismic fault stress drop Δσ~7 bars was found. However, for a circular fault area, which in this case is more realistic, we get Δσ = 55 bars. This is a relatively large value for Greek earthquakes but is explainable by increased fault rigidity as a result of the long repeat time of strong earthquakes in KF.

We investigate an earthquake sequence involving an M_w = 4.6 mainshock on 2 December 2020, followed by a seismic swarm in July–October 2021 near Thiva, Central Greece, to identify the activated structures and understand its triggering mechanisms. For this purpose, we employ double-difference relocation to construct a high-resolution earthquake catalogue and examine in detail the distribution of hypocenters and the spatiotemporal evolution of the sequence. Furthermore, we apply instrumental and imaging geodesy to map the local deformation and identify long-term trends or anomalies that could have contributed to stress loading. The 2021 seismic swarm was hosted on a system of conjugate normal faults, including the eastward extension of the Yliki fault, with the main activated structures trending WNW–ESE and dipping south. No pre- or coseismic deformation could be associated with the 2021 swarm, while Coulomb stress transfer due to the M_w = 4.6 mainshock of December 2020 was found to be insufficient to trigger its nucleation. However, the evolution of the swarm is related to stress triggering by its major events and facilitated by pore-fluid pressure diffusion. The re-evaluated seismic history of the area reveals its potential to generate destructive M_w = 6.0 earthquakes; therefore, the continued monitoring of its microseismicity is considered important.