Geographical and vertical variation of the earth's seismicity

Information concerning a total number of 13700 instrumentally recorded earthquakes is used to study the geographical and the vertical distribution of the Earth's seismicity. From these earthquakes, which form four complete samples of data (M 7.0, 1894–1992; M 6.5, 1930–1992; M 6.0, 1953–1992; M 5.5, 1966–1992), 11511 are shallow (h 60 km), 2085 are of intermediate focal depth (61 h 300 km) and 564 are deep focus earthquakes (301 h 720 km). The parameters a and b of the frequency-magnitude relationship were calculated in a grid of equally spaced points at 1° by using the data of earthquakes located inside circles centered at each point. The radius of the circles increased from 30 km with a step of 10 km until the information for the earthquakes located inside the circle fulfil three criteria which concern the size of the sample used to compute these parameters at each point of the grid. The results are given in a qualitative way (epicenter maps) as well as in a quantitative way (mean return periods).     

Seismotectonics of the northern Aegean area

Data concerning the focal mechanism and the spatial distribution of earthquakes have been used to investigate the active tectonics of the northern Aegean and the surrounding area.A thrust region, which includes the northernmost part of the Aegean and at least part of the Marmara Sea, has been defined. An amphitheatrical Benioff zone dipping towards the thrust region from south, east and probably from west, at a mean angle of about 30°, has been detected.The thrust region is surrounded by a region of normal faulting. An eastward progression of the seismic activity in this normal faulting region between 1954 and 1971 has been observed.A correspondence between the earthquake occurrence in the thrust and normal faulting regions has also been observed. Each large shock produced by tensional mechanism in the region of normal faulting is preceded or followed by one or more shocks of compressional mechanism in the thrust region.The focal mechanism, the distribution of the earthquake foci with intermediate focal depth, as well as some magnetic and gravimetric observations can be interpreted by assuming that dense oceanic crust sinks in the northern part of this area and that the adjacent lithosphere moves by segmentation to fill the void with the consequence of producing tensile stresses associated with normal faulting. Such a mechanism of lithospheric interaction suggests that accretion probably takes place in this area.     

A lithospheric model to interpret focal properties of intermediate and shallow shocks in central greece

Summary An attempt has been made to interpret the striking difference, in focal properties, between the intermediate and shallow earthquakes in Central Greece and an observed time sequence of these shocks by a lithospheric model. This model consists of a lithospheric slab descending from the Ionian to the Aegean and a back-arc expanding Aegean lithosphere. Thrust faulting near the top surface of the slab, caused by the sinking of the slab, triggers spreading and normal faulting in the back-arc Aegean region.     

Spatial Distribution of Time-independent Seismicity in China

—The seismicity in the territory of China, a seismotectonically complicated region, has been examined by using three complete samples of earthquakes which occurred during the last two centuries (1800–1995). The b value of the Gutenberg-Richter relation was estimated by using this data sample. Taking into account the fact that the b value is spatially more stable than the a value, the b values were calculated at the nodes of a normal grid superposing on the entire area studied, and their distribution was examined. The results showed that the b value increases smoothly from 0.4 in inner-Mongolia to 0.8 in the east, south and southwest of China with higher values (b>0.8) in the Taiwan region. Furthermore, keeping fixed the obtained b values, the a value distribution was also examined. In order to display more detailed information about the seismicity, smaller cell surface (10,000 km²) for the calculation of the a values was chosen. The mean return periods for different cutoff magnitudes were also calculated for each of these small cells. It was observed that the mean return periods are the shortest ones in China, which are 10 and 50 years for the magnitude larger than or equal to 6.0 and 7.0, respectively.     

An algorithm for fast 3D inversion of surface electrical resistivity tomography data: application on imaging buried antiquities

In this work a new algorithm for the fast and efficient 3D inversion of conventional 2D surface electrical resistivity tomography lines is presented. The proposed approach lies on the assumption that for every surface measurement there is a large number of 3D parameters with very small absolute Jacobian matrix values, which can be excluded in advance from the Jacobian matrix calculation, as they do not contribute significant information in the inversion procedure. A sensitivity analysis for both homogeneous and inhomogeneous earth models showed that each measurement has a specific region of influence, which can be limited to parameters in a critical rectangular prism volume. Application of the proposed algorithm accelerated almost three times the Jacobian (sensitivity) matrix calculation for the data sets tested in this work. Moreover, application of the least squares regression iterative inversion technique, resulted in a new 3D resistivity inversion algorithm more than 2.7 times faster and with computer memory requirements less than half compared to the original algorithm. The efficiency and accuracy of the algorithm was verified using synthetic models representing typical archaeological structures, as well as field data collected from two archaeological sites in Greece, employing different electrode configurations. The applicability of the presented approach is demonstrated for archaeological investigations and the basic idea of the proposed algorithm can be easily extended for the inversion of other geophysical data.     

Travel times of Pn -waves in the Aegean and surrounding area

Summary. Based on accurately located 23 very shallow earthquakes ( h = 1–14 km) in northern and central Greece by portable networks of seismic stations and by the joint epicentre method, the travel times of the P_n -waves from the foci of these earthquakes to the sites of 54 permanent stations in the Balkan region have been determined. The travel times of P_n -waves in the central and eastern part of the area (eastern Greece, south-eastern Yugoslavia, the Aegean Sea, Bulgaria, southern Romania, western Turkey) fit a straight line very well with the P_n velocity equal to 7.9 ± 0.1 km s^-1. On the contrary, the travel times of P_n -waves to stations in the western part of the area (Albania, western Greece) do not fit this curve because the P_n -waves travelling to these stations are delayed by more than 1 s due to the thicker crust under the Dinarides–Hellenides mountain range. Time delays for P_n -waves have been calculated for each permanent station in the Balkan area with respect to the mean travel-time curve of these waves in the central and eastern part of the area. Corrections of the travel times for these delays contribute very much to the improvement of the accuracy in the location of the shallow earthquakes in the Aegean and surrounding area.     

Tsunamis in and near Greece and their relation to the earthquake focal mechanisms

The major earthquake-induced tsunamis reliable known to have occurred in and near Greece since antiquity are considered in the light of the recently obtained reliable data on the mechanisms and focal depths of the earthquakes occurring here. (The earthquake data concern the major shocks of the period 1962–1986.) First, concise information is given on the most devastating tsunamis. Then the relation between the (estimated) maximum tsunami intensity and the earthquake parameters (mechanism and focal depth) is examined. It is revealed that the most devastating tsunamis took place in areas (such as the western part of the Corinthiakos Gulf, the Maliakos Gulf, and the southern Aegean Sea) where earthquakes are due to shallow normal faulting. Other major tsunamis were nucleated along the convex side of the Hellenic arc, characterized by shallow thrust earthquakes. It is probably somewhere there (most likely south of Crete) that the region's largest known tsunami occurred in AD 365, claiming many lives and causing extensive devastation in the entire eastern Mediterranean. Such big tsunamis seem to have a return period of well over 1000 years and can be generated by large shallow earthquakes associated with thrust faulting beneath the Hellenic trench, where the African plate subduces under the Euroasian plate. Lesser tsunamis are known in the northernmost part of the Aegean Sea and in the Sea of Marmara, where strike-slip faulting is observed. Finally, an attempt is made to combine the tsunami and earthquake data into a map of the region's main tsunamigenic zones (areas of the sea bed believed responsible for past tsunamis and expected to nucleate tsunamis in the future).     

Long-term earthquake prediction in the circum-pacific convergent belt

Investigation of the time-dependent seismicity in 274 seismogenic regions of the entire continental fracture system indicates that strong shallow earthquakes in each region exhibit short as well as intermediate term time clustering (duration extending to several years) which follow a power-law time distribution. Mainshocks, however (interevent times of the order of decades), show a quasiperiodic behaviour and follow the ‘regional time and magnitude predictable seismicity model’. This model is expressed by the following formulas $$\begin{gathered} \log T_t = 0.19 M_{\min } + 0.33 M_p - 0.39 \log m_0 + q \hfill \\ M_f = 0.73 M_{\min } - 0.28 M_p + 0.40 \log m_0 + m \hfill \\ \end{gathered} $$ which relate the interevent time,T _t (in years), and the surface wave magnitude,M _f , of the following mainshock: with the magnitude,M _min, of the smallest mainshock considered, the magnitude,M _p , of the preceded mainshock and the moment rate,m ₀ (in dyn.cm.yr^?1), in a seismogenic region. The values of the parametersq andm vary from area to area. The basic properties of this model are described and problems related to its physical significance are discussed. The first of these relations, in combination with the hypothesis that the ratioT/T _t , whereT is the observed interevent time, follows a lognormal distribution, has been used to calculate the probability for the occurrence of the next very large mainshock (M _s ≥7.0) during the decade 1993–2002 in each of the 141 seismogenic regions in which the circum-Pacific convergent belt has been separated. The second of these relations has been used to estimate the magnitude of the expected mainshock in each of the regions.     

Surface fault traces and fault plane solutions of the May–June 1978 major shocks in the Thessaloniki area, Greece.

It is shown that the foci of the recent earthquakes in the Thessaloniki area of northern Greece are located in an arcuate seismic zone which is associated with the Serbomacedonian geologic zone. Three main lines of fracture have been observed in the epicentral area after the May–June 1978 earthquakes. Field and macroseismic observations as well as fault plane solutions for the main shock and for the largest foreshock show that both earthquakes are due to a strike slip sinistral motion with a small reverse component on a steeply dipping and trending southeast-northwest fault.