A microearthquake study of the Mygdonian graben (northern Greece)

During March and April 1984, a temporary network of 29 portable stations was operated in the region of the Mygdonian graben near Thessaloniki (northern Greece), where a destructive earthquake (M_s = 6.5) had occurred in the Summer of 1978. During a period of six weeks we recorded 540 earthquakes with magnitudes ranging from −0.2 to 3.0. From this set of data, 254 events are selected which according to us have a precision in epicenter and depth better than 1.5 km. A total of 54 single-event focal mechanisms have been determined.The seismicity and focal mechanisms show a rather complex pattern. There are no clear individual faults, but the E-W and NW-SE striking zones show N-S extension. Zones striking NNE-SSW show dextral strike-slip motion but NW-SE zones with sinistral strike-slip are also observed.In the center of the graben where the 1978 earthquake was located, we observe several thrust mechanisms distributed in two groups showing either NW-SE or E-W compression; these earthquakes seem to be located 2 km above the earthquakes showing normal mechanisms.The mean direction of the T-axes, found from the focal mechanisms, trends N15° and dips sub-horizontal.We propose a model for the formation and evolution of a complex graben system comprising several stages. In the initial stage the deformation occurs along pre-existing NW-SE or NNE-SSW faults, with normal or strike-slip movements. In the second stage, a new, E-W trending group of normal faults is formed over the ancient fault network. These new faults have a direction perpendicular to the mean T-axis and accommodate better the actual state of stress. At this stage the initial faults adjust to the deformation produced by the E-W trending new faults, and may constitute geometric barriers to the evolution of the new normal faults.     

On-chip earthquake simulation model using potentials

A two-dimensional (2-D) Cellular Automata (CA) dynamic system constituted of cells-charges has been proposed for the simulation of the earthquake process. The CA model has been calibrated with the use of real data. The calibration incorporates major seismic characteristics of the area under test. The simulation results are found in good quantitative and qualitative agreement with the recorded Gutenberg–Richter (GR) scaling relations. The model is enriched with a powerful multi-parameter interface that enables the user to load real data from different regions. This paper examines the on-chip realisation of the model and its instrumentation. The CA model hardware implementation is based on Field Programmable Gate Array (FPGA) logic. It utilises an array of 32 × 32 cells. Parameters that construct the local CA rule constitute the input data. The initial seed, which to some extent corresponds to the seismic features of the area under test, is loaded in a semi-parallel way and the process is completed in a certain number of time steps. The automatic response of the processor provides the corresponding GR scaling law of the area under study. The hardware implementation of the CA-based earthquake simulation model is advantageous in terms of low-cost, high-speed, compactness and portability features. It can operate as a preliminary data-acquisition filter that accelerates the evaluation of recorded data as far as its origin time, spatial and magnitude completeness and quality are concerned. Software that performs reliable automatic phase picking, as well as data elaboration, can be assembled next to the earthquake recording instruments (the whole network) output to assure a quick and reliable iteration of the focal parameters of a recorded earthquake (epicentre coordinates, focal depth and magnitude). The dedicated processor can be accommodated right after this stage (before any manual elaboration) focusing on the near real-time development of a reliable qualitative dynamical seismic record and a mapping of the seismic characteristics of the area.     

Decelerating precursory seismicity in Vrancea

Preshock seismic excitation followed by seismic quiescence has been observed in the seismogenic region of strong shallow mainshocks. The strain released by such preshocks is decelerating with the time to the mainshock and is fitted by a power-law with a power value larger than unit. This model is tested in the present work for the intermediate-depth earthquakes of the Vrancea region, generated in an isolated seismogenic zone proper for such testing. A backward application of this “decelerating preshock strain” model for the case of 4.3.1977 (M = 7.5) earthquake, for which reliable data are available, shows a good fit of the power-law pattern to the seismic activity preceding the main shock. The occurrence rate of recent intermediate-depth shocks in Vrancea indicates that this region is currently in a state of decelerating seismic deformation, which may lead to the generation of a strong intermediate-depth mainshock there at about the beginning of the third decade of the present century. The respective uncertainties are unknown due to lack of previous relative studies.     

Small-scale spatial variation of the stress field in the back-arc Aegean area: Results from the seismotectonic study of the broader area of Mygdonia basin (N. Greece)

In the present work a detailed seismotectonic study of the broader area of the Mygdonia basin (N. Greece) is performed. Digital data for earthquakes which occurred in the broader Mygdonia basin and were recorded by the permanent telemetric network of the Geophysical Laboratory of the Aristotle University of Thessaloniki during the period 1989–1999 were collected and fault plane solutions for 50 earthquakes which occurred in the study area were calculated with a modified first motions approach which incorporates amplitude and radiation pattern information. Fault plane solutions for the 3 main shocks of Volvi (23/05/78, M_W = 5.8 and 20/06/78, M_W = 6.5) and Arnaia (04/05/95, M_W = 5.8) events and the 1978 aftershock sequence were additionally used. Moreover, data from two local networks established in the Mygdonia basin were also incorporated in the final dataset.Determination of the stress field was realized by the use of the method of Gephart and Forsyth [Gephart, J.W., Forsyth, D.W., 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: application to the San Fernando earthquake sequence: Jour. Geophys. Res., v.89, no. B11, p. 9305–9320] for the stress tensor inversion and the results were compared with independent estimates based on the calculation of the average moment tensor [Papazachos, C.B.,Kiratzi, A.A., 1992. A formulation for reliable estimation of active crustal deformation and its application to central Greece. Geophys. J. Int. 111, 424–432]. The obtained stress results show a relatively good agreement between the two approaches, with differences in the azimuth of the dominant extension axis of the order of 10°. Furthermore, comparison with independent information for the mean stress axes provided by the study of kinematics on neotectonic faults [Mountrakis, D., Kilias, A., Tranos, M., Thomaidou, E., Papazachos, C., Karakaisis, G., Scordilis, E., Chatzidimitriou, P., Papadimitriou, E., Vargemezis, G., Aidona, E., Karagianni, E., Vamvakaris, D. Skarlatoudis, A. 2003. Determination of the settings and the seismotectonic behavior of the main seismic-active faults of Northern Greece area using neotectonic and seismological data. Earthquake Planning and Protection Organisation (OASP) (in Greek)] shows a similar agreement with typical misfit of the order 10°. The stress inversion method was modified in order to select one or both nodal planes of the focal mechanism which corresponds to the “true” fault plane of the occurred earthquakes and was able to select a single fault plane in the majority of examined cases. Using this approach, the obtained fault plane rose diagrams are in agreement with results from various neotectonic studies. Moreover, several secondary active fault branches were identified, which are still not clearly observed in the field.     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

Surface fault traces,fault plane solution and spatial distribution of the aftershocks of the September 13, 1986 earthquake of Kalamata (Southern Greece)

A shallow earthquake ofM _S=6.2 occurred in the southern part of the Peloponnesus, 12 km north of the port of the city of Kalamata, which caused considerable damage. The fault plane solution of the main shock, geological data and field observations, as well as the distribution of foci of aftershocks, indicate that the seismic fault is a listric normal one trending NNE-SSW and dipping to WNW. The surface ruptures caused by the earthquake coincide with the trace of a neotectonic fault, which is located 2–3 km east of the city of Kalamata and which is related to the formation of Messiniakos gulf during the Pliocene-Quaternary tectonics. Field observations indicate that the earthquake is due to the reactivation of the same fault.A three-days aftershock study in the area, with portable seismographs, recorded many aftershocks of which 39 withM _S1.7 were very well located. The distribution of aftershocks forms two clusters, one near the epicenter of the main shock in the northern part of the seismogenic volume, and the other near the epicenter of the largest aftershock (M _S=5.4) in the southern part of this volume. The central part of the area lacks aftershocks, which probably indicates that this is the part of the fault which slipped smoothly during the earthquake.     

Evidence for transform faulting in the Ionian sea: The Cephalonia island earthquake sequence of 1983

Accurate locations of aftershocks of the January 17, 1983 (M _s=7.0) main shock in the Ionian islands have been determined, as well as fault plane solutions for this main shock and its largest aftershock, which are interpreted as a right-lateral, strike-slip motion with a thrust component, on a fault striking in about a NE-SW direction.This is considered as a transform fault in the northwesternmost part of the Hellenic arc.     

Crustal and upper mantle structure of the Kozani-Grevena area obtained by non-linear inversion of P and S travel times

The present study focuses on the P and S crustal and uppermost mantle velocity structure in the broader Kozani-Grevena area. The velocity structure is derived from the inversion of travel times of local events. The main data source is the travel times from the aftershock sequence of the large event of 13 of May 1995 (M_w = 6.6) which occurred in the study area. An appropriate preconditioning of the final linearized system is used to reduce ray density effects on the results. An attempt is made to interpret the features and details of the crustal structure in terms of the geotectonic setting of the area. The observed features of the deeper crustal and uppermost mantle structure are in very good agreement with previous results. Specifically, a crustal thickening is observed along a ENE-WSW direction, perpendicular to the well-known Dinaric trend (NNW-SSE) of the geological formations of the area, in accordance with the theoretical expectation of a thicker crust under the accretion prism which starts at the SW edge of the study area.     

Recent reliable observations and improved tests on synthetic catalogs with spatiotemporal clustering verify precursory decelerating–accelerating seismicity

We examined the seismic activity which preceded six strong mainshocks that occurred in the Aegean (M?=?6.4–6.9, 33–43° N, 19–28° E) and two strong mainshocks that occurred in California (M?=?6.5–7.1, 32–41° N, 115–125° W) during 1995–2010. We find that each of these eight mainshocks has been preceded by a pronounced decelerating and an equally easily identifiable accelerating seismic sequence with the time to the mainshock. The two preshock sequences of each mainshock occurred in separate space, time, and magnitude windows. In all eight cases, very low decelerating seismicity, as well as very low accelerating seismicity, is observed around the actual epicenter of the ensuing mainshock. Statistical tests on the observed measures of decelerating, q _d, and accelerating, q _a, seismicity against similar measures calculated using synthetic catalogs with spatiotemporal clustering based on the ETAS model show that there is an almost zero probability for each one of the two preshock sequences which preceded each of the eight mainshocks to be random. These results support the notion that every strong shallow mainshock is preceded by a decelerating and an accelerating seismic sequence with predictive properties for the ensuing mainshock.