Abstract. The Alboran system reflects the interplay of slow Nubia–Iberia convergence, inherited structures, and ongoing lithospheric attenuation. Earthquake occurrence mainly tracks active crustal fault networks that partition the oblique plate motion into strike-slip and extension. The 2021–2022 seismic swarm in the western Alboran Sea represents an exceptional episode characterized by complex temporal-spatial evolution, primarily influenced by fluid-driven processes interacting with inherited fault systems. We analyzed arrivals of approximately 7,000 seismic events recorded by Spanish (IGN) and Moroccan (CNRST) seismic networks. Both bulletins have been individually and jointly processed using the double-difference algorithm (HypoDD) and a regionally optimized velocity model. This approach significantly improved hypocentral precision, reducing event scatter and delineating a clearly defined, near-vertical seismic conduit-oriented NW–SE. Spatiotemporal analyses revealed distinct episodes of seismic migration, consistent with episodic fluid overpressure pulses, confirmed by diffusivity values (1.2–13.9 m2/s) characteristic of fluid-controlled swarms. Focal mechanisms predominantly indicated strike-slip motion, aligning with the regional transtensional tectonics and pinpointing the unrecognized Ras Tarf Fault as the primary seismogenic structure likely linked to the Al-Idrissi Fault System (AIFS). Integration with vertical and horizontal shear-wave velocity models (VSH and VSV) highlighted velocity anomalies at depths of 30–70 km, suggesting the presence of partially serpentinized mantle wedges above a delaminating slab segment, further supporting fluid involvement. Our results emphasize the critical interplay between deep lithospheric fluids, inherited fault structures, and regional tectonic stress, providing a comprehensive framework for understanding the 2021–2022 swarm dynamics which could improve the seismic hazard assessment in the region.

The 2020-2021 seismic sequence at the Western Gulf of Corinth, Central Greece, has been thoroughly analyzed using seismological and geodetic data processing. We present a high-resolution dataset of over 4000 relocated earthquakes between June 2020 and February 2021, delineating the activated structures. We examine the evolution of various clusters that were triggered during a three-stage sequence which began on December 23, 2020, with an Mw = 4.6 event near Marathias, migrated eastwards, following an Mw = 5.0 event on January 12, 2021 near Trizonia Island, and culminated on February 17 with an Mw = 5.3 offshore event north of Psatho-pyrgos. Focal mechanisms for the 20 stronger events, determined by moment tensor inversion, and 36 weaker ones, obtained through first motion polarities measurements, revealed dominant normal faulting. The fault plane of the February 17, 2021 major event is modeled through the inversion of geodetic data. Results suggest a very shallow geodetic centroid at 1.5 km depth, consistent with the seismic centroid (3.5 km), but different from the relocated hypocenter (7.5 km). Spatiotemporal analysis reveals seismic migration, following a diffusion law with D values in the range 0.06-0.25 m2s  1, consistent with triggering due to pore-pressure diffusion by fluids intrusion. This enables failure on non-optimally oriented faults, which explains the significant proportion of the observed strike-slip and oblique-normal faulting. Coulomb stress transfer due to the major events shows only weak stress-loading on the fault of the January 12, 2021 earthquake and a higher level of positive stress transfer to the south-dipping than to the north-dipping nodal plane of the 17 February event. The latter could have facilitated the nucleation of the rupture at 7.5 km and its further propagation to shallower depths along a south-dipping plane, where the bulk of its seismic and geodetic moment was released.

The North Anatolian Fault is the similar to 1200-km-long active continental transform boundary between Anatolia and Eurasia. This strike-slip system initiated around 10-12 Ma and experienced diachronous episodes of strain localization along its strike. The structural evolution of the similar to 350-km-long fault segments crossing the North Aegean Sea remains to be accurately investigated. There, the moder n Nor th Anatolian Fault is localized along two main branches: the nor ther n branch ends at the Nor th Aegean Trough and the souther n branch ends at the Edremit-Skyros Trough. The Evia Basin is located in the North Aegean Domain between the North Anatolian Fault and the Corinth Rift. This study presents seismic reflection lines crossing the aforementioned structures of the North Aegean Domain, which document their subsurface structure and the sedimentary record of their activity since the Messinian. The seismic-reflection data set is tied to regional-scale stratigraphic markers, which constrains the age of main tectonic events related to the formation of the North Anatolian Fault. The seismic-reflection lines show that the two main branches of the North Anatolian Fault became localized structures at 1.3-2 Ma, coe v all y with the formation of the Evia Basin. Since 2 Ma, the North Aegean Troughs developed as a series of horsetail basins propagating westwards at the termination of the branches of the North Anatolian Fault. On a regional scale, the wide and diffuse North Anatolian transtensive shear zone active from Serravalian to Late Pliocene turned into a narrower shear zone at the two main branches of the North Anatolian Fault since the Early Pleistocene. This abrupt episode of strain localization occurred in the frame of the major Early Pleistocene change in stress regime from NE-SW to N-S extension, which has been observed throughout the Aegean Sea.

The Aetolia-Akarnanian region, in Western Greece, is considered to be part of a micro-plate in formation, named the Ionian Island-Akarnanian Block (IAB), in the larger-scale Central Mediterranean tectonic context. The IAB accommodates the deformations between the surrounding tectonic structures that are the Corinth Gulf, the Hellenic subduction, the Kefalonia Transform Fault and the Apulian collision. This work presents the first results of a dense temporary seismic survey in the Aetolia-Akarnanian region (from the Amvrakikos Gulf to the Patras Gulf). Our local dense network has been designed in order to avoid gaps and to allow the recording of a major part of the Akarnania seismicity. With a semi-automatic events detection and picking program, we detected more than 15000 events from October 2015 to December 2018. With this important data set we constrained a 1D local velocity model. The comparison with the previous published models shows a possible significant velocity variation inside the region and especially at the Trichonis lake graben. Thanks to our data set and our velocity model, we precisely located 12723 seismic events with magnitude 0 < ML < 4.6, and a magnitude of completeness Mc = 1.0, that represents actually the most important catalogue for the Aetolia-Akarnania. Seismicity highlights specific seismic structures as clusters and a seismic plane below the West of Corinth Gulf that are briefly discussed.

SUMMARY An unusual seismic activity has recently occurred in the Gripp valley, located in the central part of the French Pyrenees. Since spring 2020, two new swarms appeared, clearly outside the usual location of the seismicity in this area. On 20 September 2020, almost concomitantly with the activation of the second seismic swarm, a hole suddenly opened in the bed of a local river, the Adour de Payolle. This hole drained the water from the river, which dried up over 500 m. We follow and study the spatial and temporal evolution of these clusters, using four temporary stations deployed a few days or months after the beginning of the crisis to complete the regional network. These additional data lead to the construction of a comprehensive catalogue of more than 4900 earthquakes, using both a template matching approach and a deep-learning based phase picking method to complete and improve the initial catalogue available from the French seismological agency. This allows highlighting a slow and clear migration of the seismicity during 1 yr. Precise absolute and relative event locations reveal a dipping faulting structure, confirmed by the focal mechanism estimated for the highest magnitude event of the sequence (ML 3). We propose to explain the observed migration of the seismicity by deep fluids going up through a newly discovered faulting structure.

inversion, we explore the ranges of geometric parameters that can explain the data. We obtain an average final model and its standard deviation, with small slip amplitude at the surface, consistent with the field observations, and with slip as large as 2.5 m at depth. This model is compared with those previously published. We conclude that an antithetic fault is not required to explain the SAR data.

Analysis of lithospheric deformation is key to understanding current tectonics and other active deformation processes. The Alceste project, conducted in the framework of the Résif Seismicity Transverse Action, aims at proposing an updated seismic hazard model in metropolitan France built from the most recent data and academic consensus. One of the contributions will come from geodetic observations through strain rate integration in seismotectonic zoning and seismic hazard models. Most of Western Europe in general, and metropolitan France in particular, is located within the Eurasian Plate, which has very low deformation and seismicity rates. These GNSS-derived secular velocity field could be related to the combination of different deformation processes, with a minor contribution from plate tectonics (relative plate motions, mantle convection, etc), while most of the measured velocities could be explained by non-tectonic long-term or transient processes (gravitational motions, Glacial Isostatic Adjustment, erosion, anthropogenic deformation, etc). Some of these physical processes causing surface deformation also reflect stress changes at depth that may be associated with loading on active faults and seismicity. Properly mapping this deformation is therefore a key to better assess seismic hazard in slow straining areas. In order (i) to assess the variability due to the diversity of the strain-rate calculation methods used in the scientific community and (ii) to test their capacity to resolve low-amplitude consistent surface deformation, we conduct a benchmark exercise. We build sets of synthetic velocity fields sampled at the existing GNSS permanent stations from the RENAG (REseau NAtional GNSS Permanent), RGP (Réseau GPS Permanent) and other permanent and non-permanent benchmarks. Our synthetic velocity fields have the same characteristics (noise, uncertainties) as the observed velocities in metropolitan France but they contain surface deformation signals from known physical processes (block rotations, fault elastic loading, large scale flexure, etc). We compare the strain rate invariants derived independently by nine different RENAG research teams (using different software) to the expected strain rate patterns and discuss drawbacks and advantages of each approach. In a second step, we analyzed the strain rate tensors derived from the synthetic velocity fields to discuss potential regional style of the deformation in metropolitan France. Previous studies have shown that the computation of strain rate tensors is impacted by the user-defined parameters and the algorithm specificity used. Exploring these different biases in the strain rate solutions represent the opportunity to improve the understanding of the conventional problem of the standard interpolation.

Moderate‐to‐large earthquakes in rifts may occur on leading boundary faults or inner antithetic faults. Here we show a rare case of the 2020–2021 seismic sequence in the Corinth rift, that culminated in the shallow rupture of the antithetic fault, neither preceded nor followed by the leading fault rupture. The hypocenter of the largest shock (Mw 5.3 of 17 February 2021) was located at ∼8 km depth. However, seismic waveform data, supported by satellite‐geodetic and tide gauge measurements, pointed to rupture at shallow depth (∼3 km), where no earthquakes were previously observed. We show that the earthquake most probably ruptured two orthogonal, conjugate fault segments: a weak nucleation phase occurred in the microseismically highly active sub‐horizontal detachment layer, followed – a few seconds later – by a larger, shallow moment release on a high‐angle, south‐dipping normal fault. The latter is the Mornos offshore fault, antithetic to the leading, north‐dipping Psathopyrgos fault. Our study presents the first instrumental/observational evidence of a very shallow Mw 5+ event in this rift – and one of the few reported worldwide. The depth limit of the main shallow slip patch coincides with the expected crossing of the Mornos fault with the Psathopyrgos fault, stressing the importance of fault segmentation and rooting inherited from the rift history. This unusual shallow slip in a depth range with little background seismicity and few aftershocks needs to be further investigated by dynamic modeling as a possible prototype of hazardous events in rift environments.

Following the installation of a temporary seismological network in western Greece north of the Gulf of Patras, we determined the quality of the sites of each of the 10 stations in the network. For this, we used the horizontal-to-vertical spectral ratio (HVSR) method and calculated an average curve over randomly selected days between 0 and 10 Hz. The daily HVSR curve is determined by the HVSR 12-hr calculation (1 hr every two) without distinction between seismic ambient noise and earthquake signal. The HVSR curves obtained can be classified in three categories: flat curves without amplification, curves with a amplification peaks covering a large frequency range, and curves with one or more narrow peaks. In this third category C3, one station has one peak, two have two and one has three. On the contrary of what it is commonly assumed, the amplitudes and the resonance frequencies of these narrow peaks are not stable over time in C3. We determined the maximum of the amplitude of each peak with the corresponding central frequency for each day during 2.5 yr. Except for the station with three peaks, which finally appears stable within the uncertainties, the principal peak exhibits a seasonal variation, with a maximum in winter and a minimum in summer, the observations being more dispersed during winter. The second peak, when it exists, varies in the same way except at one station where it varies oppositely. These variations are clearly correlated with the loading and unloading cycle of the underlying aquifers as shown by the comparison with water level and yield measurements from wells located close to the stations. Moreover, they are also correlated with the vertical surface displacements observed at continuously recording GPS stations. The dispersion of the observed maximum amplitude in winter is probably related to the rainfall and the soil moisture modifying the S-wave velocity as revealed by other studies. From this study, we would like to emphasize that the use the HVSR method to constrain the S-wave velocity and the thickness of the sediment layer over the bedrock in the basin, has to be done with caution. Upon further confirmation of its robustness, the HVSR methodology presented here could be a good and easy-to-use tool for a qualitative survey of the aquifer backdrop and its seasonal behaviour, and of the soil moisture conditions.

Seismic monitoring of southwestern France began in the 1960s, and homogeneous coverage by observation networks has been in place since the 1990s. The accumulation of data now allows a refined understanding of regional seismicity, not only on its spatial aspects, but also on the regularity of the earthquake distribution over time. This paper is both a review of the work carried out on the subject, and a progress report on the current knowledge of the regional seismicity in its seismotectonic context. With the support of maps, the available catalogs are exploited at different nested scales, from the region as a whole to the numerous clusters that characterize the seismicity of southwestern France, and more specifically that of the Pyrenees. An exhaustive study of these Pyrenean clusters and their temporal behavior is proposed, allowing in particular a better description of the prominent seismicity stripe to the northwest of the range.

We investigate a seismic crisis that occurred in the western Gulf of Corinth (Greece) between December 2020 and February 2021. This area is the main focus of the Corinth Rift Laboratory (CRL) network, and has been closely monitored with local seismological and geodetic networks for 20 yr. The 2020–2021 seismic crisis evolved in three stages: It started with an Mw 4.6 event near the northern shore of the Gulf, opposite of Aigion, then migrated eastward toward Trizonia Island after an Mw 5.0 event, and eventually culminated with an Mw 5.3 event, ∼3 km northeast of the Psathopyrgos fault. Aftershocks gradually migrated westward, triggering another cluster near the junction with the Rion–Patras fault. Moment tensor inversion revealed mainly normal faulting; however, some strike-slip mechanisms also exist, composing a complex tectonic regime in this region dominated by east–west normal faults. We employ seismic and geodetic observations to constrain the geometry and kinematics of the structures that hosted the major events. We discuss possible triggering mechanisms of the second and third stages of the sequence, including fluids migration and aseismic creep, and propose potential implications of the Mw 5.3 mainshock for the seismic hazard of the region.

Characterized for the first time in 2013, the Island Akarnanian Block (IAB) is a micro-plate located in the western Greece. This micro-plate accommodates the deformation in between larger scale tectonic structures as the Gulf of Corinth (South-East), the Hellenic subduction (South) and the Apulian Collison (West). We started a micro-seismic survey (MADAM) at the end of 2015 with a dense seismological network over the area, between the Gulf of Patras and the Gulf of Amvrakikos. In order to obtain precise locations of the recorded events, we better constrained the local velocity model. In fact, several velocity models (local or regional) have been proposed for this area. However, the velocity model generally used by the scientific community remains the Hasslinger 98 velocity model. This model, nevertheless, raises some questions about its physical meaning, mainly due to a low velocity layer between 4 and 7 km-depth. Thanks to our seismological network and permanent networks of the Corinth Rift Laboratory and the Hellenic Unified Seismic Network, we collected and analysed a huge quantity of data acquired between October 2015 and December 2017. Those analyses of more than 10,000 events allowed us to develop a new and robust local velocity model, which is consistent with the seismic data and the geophysical observations. The observed seismic activity is characterized by the presence of numerous clusters. The clusters are analysed in detail by relative relocations in order to appraise their physical processes and their possible implications in the fault activity to finally have a better understanding of the deformation mode(s) of the IAB micro-plate.

The North-Eastern zone of the Gulf of Corinth in Greece is characterized by the rotation of a micro-plate in formation. The Island Akarnanian Block (IAB) have been progressively individualized since the Pleistocene (less than ~ 1.5 My ago). This micro-plate is the result of a larger-scale tectonic context with, on one side the N-S extension of the Gulf of Corinth to the East, and on the other side the Hellenic subduction to the South and the Apulian collision to the West. To the Northeast, the IAB micro-plate is bounded by a large North-South sinistral strike-slip fault system, the Katouna-Stamna Fault (KSF) and by several normal faults. To the North, normal faults reach the limit between Apulian and Eurasian plates and to the East, they form the East-West graben of Trichonis lake.Although the structures and dynamics behind the Gulf of Corinth extension are today relatively known, nevertheless, the set of faults linking the Gulf of Corinth to the Western subduction structures remain poorly studied. The seismicity recorded by the Greek national network shows discrepancies regarding to the faults mapped on the surface.At the end of 2015, a new micro-seismicity campaign started with the deployment of a temporary seismological network in an area ranging from the Gulf of Patras to the Amvrakikos Gulf toward the North. This network includes 17 seismic stations, recording continuously, added to the permanent stations of the Corinth Rift Laboratory (CRL) and of the Hellenic Unified Seismic Network (HUSN).The analysis of the seismological records is still in process for the 2016 and 2017 years. Our study consists first in picking the P- and S- waves, and then to precisely localize the seismic events recorded by our temporary seismological network combined with the permanent ones. We will present here the event location map obtained for the 2016-2017 period, a new seismic velocity model, and focal mechanisms. The seismic activity including thousands of events, is characterized by the presence of numerous clusters of few days to few weeks duration. The clusters are analysed in detail by relative relocations in order to appraise their physical processes and their implications in the fault activity. We will discuss the deformation mode of the region and build a seismotectonic model consistent with the regional geodynamics and observations.

The seismicity in the Pyrenees is continuous and well surveyed since more than 20 years. We use the catalogue of seismicity between 1997 and 2013 to explore the spatial variations of the b-value, which corresponds to the slope of the frequency-magnitude distribution of the earthquakes. Especially, variations of the b-value characterise the state of stress of the crust, possibly highlighting a deficit of large earthquake occurrence. We estimate the differential crustal stress from the b-value using a relationship published by Scholz (Geophys Res Lett 42:1399–1402, 2015). We also estimate the stress drop variations by determining a power law which links the magnitude to the seismic source radius in the Pyrenees. We focus on the depth variations and we analyse vertical profiles of b-value, differential stress and stress drop, first in the whole Pyrenean belt, then in 10 subregions. The b-values are generally smaller than 1, except in the uppermost 3–5 km where the obtained high values could be linked to the presence of fluids. Downward, the b-values decrease slowly or remain constant until a depth of increase, which could correspond to the brittle-ductile limit of the crust. We propose that this depth and the regional and vertical variations of the b-values are related to the regional tectonic context and possibly to the density heterogeneities. We also suggest that stress drop and differential stress are linearly correlated and that the stress drop is at least 1.8‰ of the differential stress.

Normal faulting mechanisms observed in the northern foothills of the Central–Western Pyrenees are remarkable, since one expects thrust faults at a convergent plate boundary. To understand the mechanisms involved, we used numerical modeling and investigated the impact of the following processes: gravitational potential energy associated with topography and dense crustal blocks; isostatic compensation in response to denudation and/or sedimentation. To decipher the effects of each process, we designed three model geometries and added a pre-existing weak fault where most of the seismicity occurs. To evaluate our model results, we derived the fault slip rate from the focal mechanisms in the region where we have the fault in our model. We found a slip rate of ∼15m/Ma, which is in agreement with our modelling results. We conclude that flexural rebound induced by surface processes is able to explain the seismicity in Central–Western Pyrenees.

The Yellowstone–East Snake River Plain hotspot track has been intensely studied since several decades and is widely considered to result from the interaction of a mantle plume with the North American plate. An integrated conclusive geodynamic interpretation of this extensive data set is however presently still lacking, and our knowledge of the dynamical processes beneath Yellowstone is patchy. It ...

The Béarn range, located to the north of the Axial Zone in the Western Pyrenees, is affected by numerous small‐magnitude seismic events. These events overlap an area characterized by specific geological structures which are interpreted to have resulted from multistage extensional and compressional deformation. An analysis of surface geology draped over digital elevation model, together with field investigations, allow identification of two main shortening episodes with differing direction of contraction: D1 represents the inversion of the North Pyrenean Basin, whose Mesozoic infill was detached from highly extended crust and transported southward over the necking zone and the northern margin of the Iberian plate; D2 corresponds to the collision stage and is characterized in the study area by backfolding and backthrusting deformation coeval with uplift in the axial part of the chain due to thickening of the Iberian plate. The microseismicity appears to concentrate along the basal part of the inverted basin units (D1) where this initially low angle thrust has been tilted and steepened during collision (D2). We propose that local steepening of this ancient inversion structure, which should not be named “North Pyrenean Fault,” provided the suitable dip for extensional solicitation in association with the present uplift of the Axial Zone, whatever the driving mechanism of this uplift could be.

The Pyrenean mountain range is a slowly deforming belt with continuous and moderate seismic activity. To quantify its deformation field, we present the velocity field estimated from a GPS survey of the Pyrenees spanning 18 yr. The PotSis and ResPyr networks, including a total of 85 GPS sites, were installed and first measured in 1992 and 1995-1997, respectively, and remeasured in 2008 and 2010. We obtain a deformation field with velocities less than 1 mm yr(-1) across the range. The estimated velocities for individual stations do not differ significantly from zero with 95 per cent confidence. Even so, we estimate a maximum extensional horizontal strain rate of 2.0 +/- 1.7 nanostrain per year in a N-S direction in the western part of the range. We do not interpret the vertical displacements due to their large uncertainties. In order to compare the horizontal strain rates with the seismic activity, we analyse a set of 194 focal mechanisms using three methods: (i) the 'r' factor relating their P and T axes, (ii) the stress tensors obtained by fault slip inversion and (iii) the strain-rate tensors. Stress and strain-rate tensors are estimated for: (i) the whole data set, (ii) the eastern and western parts of the range separately, and (iii) eight zones, which are defined based on the seismicity and the tectonic patterns of the Pyrenees. Each of these analyses reveals a lateral variation of the deformation style from compression and extension in the east to extension and strike-slip in the west of the range. Although the horizontal components of the strain-rate tensors estimated from the seismic data are slightly smaller in magnitude than those computed from the GPS velocity field, they are consistent within the 2 sigma uncertainties. Furthermore, the orientations of their principal axes agree with the mapped active faults.

The Pyrenees results from the convergence of the Iberian and Eurasian plates since Cenozoic, but the present stress regime of the range is debated, as most of the recent earthquakes exhibit normal solutions. We analyse the seismicity in central-western Pyrenees, which is the most active part of the range with the largest events. Seismicity maps obtained at different depths reveal quasi-periodic features in focus distribution, and a preferential occurrence of large earthquakes at the base of the upper crust. The superimposition of the seismicity to the gravity anomalies shows that earthquakes are mostly located on the southern border of positive Bouguer anomalies, which correspond to dense lower crust blocks trapped in the upper crust during the Pyrenean convergence. We propose that the seismicity results from the subsidence of these blocks previously exhumed inside pull-apart basins. This scenario explains all together the geographic distribution of the seismicity, its magnitude distribution and the predominance of normal focal solutions. It shows that the normal mechanisms do not necessarily imply a general north–south extension of the range, but may be compatible with a weak compressive regime, and that the stress field may not be uniform along the range.

Space geodetic techniques such as interferometric synthetic aperture radar (InSAR) and global positioning systems (GPS) have demonstrated to be useful in mapping the displacement fields of large earthquakes (moment magnitude (M-w) approximate to 6 or higher). However, the displacement fields of smaller earthquakes (M-w < 5.5), such as those that typically result from the collision of the European and African plates, are less often analysed by space geodetic techniques, and their characterization, in terms of slip along the fault plane at depth and focal depth location, often challenges current seismological techniques. This letter presents the results of InSAR analysis of the 11 May 2011, M-w 5.1, Lorca earthquake. The Lorca earthquake occurred close to an area undergoing rapid subsidence due to sediment compaction related to water pumping. Therefore, it is challenging to separate the InSAR signals due to the earthquake from those due to human activity. We used four sets of SAR data acquired from the European C-band Advanced SAR (ASAR) sensor on board the Environmental Satellite (ENVISAT) to map the surface-displacement field in the Lorca region. Then, we use a simple elastic dislocation model to characterize the fault plane geometry and the fault slip at depth. We find that the InSAR signals can be explained by approximate to 21 cm reverse slip with a approximate to 6 cm left-lateral component on a 3 kmx3 km segment centred at 4.2 km depth dipping 45 degrees NW and striking N65 degrees E, consistent with the rupture of a segment of the Alhama de Murcia fault and consistent with recent published analyses. Interestingly, the InSAR signal can also be explained by approximate to 21 cm reverse slip with a approximate to 6 cm left-lateral component on a 3 kmx3 km segment centred at approximate to 4.2 km depth dipping 50 degrees SE and striking N230 degrees E, consistent with preliminary focal plane solutions indicating a rupture on a previously unmapped blind structure. We conclude that the second model cannot be rejected on the base of the InSAR results, the complex surface-displacement pattern (containing both seismic and non-seismic displacement), the different preliminary moment tensor solutions and the published locations of aftershocks at depth.