George Veis

Global Positioning System constraints on plate kinematics and dynamics in the eastern Mediterranean and Caucasus

We present and interpret Global Positioning System (GPS) measurements of crustal motions for the period 1988–1997 at 189 sites extending east-west from the Caucasus mountains to the Adriatic Sea and north-south from the southern edge of the Eurasian plate to the northern edge of the African plate. Sites on the northern Arabian platform move 18±2 mm/yr at N25°±5°W relative to Eurasia, less than the NUVEL-1A circuit closure rate (25±1 mm/yr at N21°±7°W). Preliminary motion estimates (1994–1997) for stations located in Egypt on the northeastern part of Africa show northward motion at 5–6±2 mm/yr, also slower than NUVEL-IA estimates (10±1 mm/yr at N2°±4°E). Eastern Turkey is characterized by distributed deformation, while central Turkey is characterized by coherent plate motion (internal deformation of <2 mm/yr) involving westward displacement and counterclockwise rotation of the Anatolian plate. The Anatolian plate is de-coupled from Eurasia along the right-lateral, strike-slip North Anatolian fault (NAF). We derive a best fitting Euler vector for Anatolia-Eurasia motion of 30.7°± 0.8°N, 32.6°± 0.4°E, 1.2°±0.1°/Myr. The Euler vector gives an upper bound for NAF slip rate of 24±1 mm/yr. We determine a preliminary GPS Arabia-Anatolia Euler vector of 32.9°±1.2°N, 40.3°±1.1°E, 0.8°±0.2°/Myr and an upper bound on left-lateral slip on the East Anatolian fault (EAF) of 9±1 mm/yr. The central and southern Aegean is characterized by coherent motion (internal deformation of <2 mm/yr) toward the SW at 30±1 mm/yr relative to Eurasia. Stations in the SE Aegean deviate significantly from the overall motion of the southern Aegean, showing increasing velocities toward the trench and reaching 10±1 mm/yr relative to the southern Aegean as a whole.

Active deformation of the Corinth rift, Greece : Results from repeated Global Positioning System surveys between 1990 and 1995

Between 1990 and 1995, we carried out seven Global Positioning System (GPS) campaigns in the Corinth rift area in order to constrain the spatial and temporal crustal deformation of this active zone. The network, 193 points over ∼10,000 km2, samples most of the active faults. In order to estimate the deformation over a longer period, 159 of those points are also Greek triangulation pillars previously measured between 1966 and 1972. Two earthquakes of magnitude 6.2 and 5.9 have occurred in the network since it was installed. The extension rate deduced from the analysis of the different GPS data sets is 14±2 mm/yr oriented N9° in the west, 13±3 mm/yr oriented S-N in the center, and 10±4 mm/yr oriented N19°W in the east of the gulf. The comparison between GPS and triangulation gives higher rates and less angular divergence (25±7 mm/yr, N4°E; 22±7 mm/yr, S-N; 20±7 mm/yr, N15°W, respectively). Both sets of data indicate that the deforming zone is very narrow (10–15 km) in the west, might be wider in the center (15–20 km), and is more diffuse in the east. The analysis of the displacements observed after the Ms = 6.2, June 15, 1995, and the Ms = S.9, November 18, 1992, earthquakes, both located in the west of the gulf, together with seismological and tectonic observations shows that these two earthquakes occurred on low-angle (≤35°) north dipping normal faults located between 4.5 and 10 km depth in the inner part of the rift. Assuming that the deformation is concentrated in relatively narrow deforming zones, we use a simple model of a dislocation in an elastic half-space to study the implication of the localization. Using the geometry of the known seismogenic faults, our observations imply continuous aseismic deformation in the uppermost crust of the inner rift. This model predicts geodetic strain rates close to seismic strain rates in opposition to previous estimates. This is because our model takes into account the activity on low-angle normal faults in the inner rift and an effective seismogenic layer of 6–7 km, about half that usually assumed.

GPS‐derived strain rate field within the boundary zones of the Eurasian, African, and Arabian Plates

We use the GPS velocity field (1988–1998) for eastern Mediterranean and Asia Minor to determine the crustal deformation strain rate field in an area bounded by 35°N and 43°N, and 20°E and 48°E. We calculate the normal and shear strain rate components associated with the major faults and compare these qualitatively with seismological data. Uncertainties in the calculation of the strain rates reach 50 nstrain yr−1 in sparsely observed parts of Anatolia, whereas we estimate errors <20 nstrain yr−1 in the Aegean and Marmara regions. The largest compressional strain rate components in the eastern part of the study area occur along the Greater Caucasus mountain front reaching 70 nstrain yr−1. (1 nstrain yr−1 = 0.0317 × 10−15 s−1). The North Anatolian Fault Zone is the clearest feature in the shear strain rate field. It is expressed as a pronounced dextral strike-slip fault zone, reaching rates of up to 170 nstrain yr−1. This holds true also for the Izmit area, where the August 17, 1999, earthquake occurred. Central Anatolia is almost strain-free, whereas extension prevails in western Anatolia. The principal axes of extension vary around the N-S direction with strain rates of up to 85 nstrain yr−1. These extensional areas coincide with graben features and normal faulting earthquakes. The central and southwestern Aegean Sea is strain-free with values far below 40 nstrain yr−1. The seismic cluster around the Dodekanissa islands, southeastern Aegean Sea, coincides with NW-SE oriented extension, attaining strain rates of up to 90 nstrain yr−1. This area of extension also exhibits recent active volcanism. The entire Hellenic arc shows compressional strain rates perpendicular to the arc. The Pliny-Strabo troughs along the eastern segment of the arc show left-lateral shear strain rates reaching 80 nstrain yr−1. Significant extension is found in central Greece, with a NNE-SSW oriented maximum of 120 nstrain yr−1 centered around the Gulf of Corinthos. The Kephalonia Fault Zone in NW Greece is a distinct dextral fault zone, separating Apulia from the rapidly moving Aegean microplate. Right-lateral shear strain rates reach 150 nstrain yr−1.

Geodetic strain of Greece in the interval 1892–1992

A first-order triangulation of Greece was carried out in the 1890s. Reoccupation, using Global Positioning System receivers, of 46 of the 93 original markers yielded estimates of the deformation of the region over the intervening interval. Broad regions have similar geodetic strain over the 100-year time span. Strain north of the Gulf of Korinthos is predominantly north-south extension, though with a significant east-west component. The central Peloponnisos is relatively stable, whereas the gulfs of the southern Peloponnisos are all characterized by uniaxial east-west extension. The seismic expression of strain for the entire region, calculated from the seismic moment tensors of earthquakes of M S ≥ 5.8 during the past 100 years, accounts for only 20-50% of the geodetically determined strain. At a scale of 50-100 km, the fraction of the strain that is expressed seismically varies much more than this range. In particular, whereas seismic strain in the eastern Gulf of Korinthos over the past 100 years is commensurate with the geodetic strain, there is rapid extension across the western Gulf of Korinthos (∼0.3 μstrain yr -1 ), with negligible seismic strain for the 100 year period prior to 1992. The Egion earthquake of June 1995 in the western Gulf of Korinthos released only a small proportion (≤20%) of the elastic strain that had accumulated in that region. The observed distribution of displacements can be explained by the relative rotation of two plates with a broad accommodation zone between them, but it is equally consistent with the deformation that would be expected of a sheet of fluid moving toward a low-pressure boundary at the Hellenic Trench. A simple calculation implies that if the region does behave as a fluid, then its effective viscosity is ∼10 22 ∼10 23 Pa s. Such viscosities are consistent with the deformation of a lithosphere obeying a rheological law similar to that obtained for olivine in the laboratory.

The strain rate field in the eastern Mediterranean region, estimated by repeated GPS measurements

We use the combined GPS velocity field of the eastern Mediterranean for the period 1988 to 1996 to determine crustal deformation strain rates in a region comprising the Hellenic arc, the Aegean Sea, and western Anatolia. We interpret the velocity field and determine the strain rate tensor by the spatial derivatives of the collocated motion vectors. The region following the line Marmara Sea, North Aegean Trough, northern central Greece, and the central Ionian islands is associated with strong right-lateral shear motion, with maximum shear strain rates of 180 nano-strain/a (180×10−9/a). In the central Aegean Sea, N–S-oriented extensional processes prevail, reaching 100 nano-strain/a. The southern Aegean is characterized by relatively small strain rates. Maximum extensional components of the strain rate tensor, reaching 150 nano-strain/a in a N–S direction, are found in central Greece. The Hellenic arc is associated with moderate arc-parallel extension and strong compression perpendicular to it. Projections of the strain rates parallel to the major fault zones reveal that the northern Aegean is governed by the westward continuation of the North Anatolian Fault Zone which is associated with strong dextral shearing (maximum 220 nano-strain/a), accompanied by numerous large earthquakes in this century.

Geodetic estimate of seismic hazard in the Gulf of Korinthos

The recent 15 June 1995, M0 = 6.0 × 1018 N m, Aigion earthquake in the western Gulf of Korinthos has focussed attention on the seismic hazard of the region. Although there have been few large earthquakes in the region during this century, the historical record suggests that there may have been many large earthquakes there in the interval 1750–1900. We present geodetic data that give estimates of the rate of extension of the Gulf of Korinthos during this century and which suggest that less than half of the elastic strain in the central and western Gulf of Korinthos has been released by earthquakes during this century. In contrast, the seismic and geodetic strains in the eastern Gulf of Korinthos are in agreement with each other. If the discrepancy between seismic and geodetic strains in the western Gulf of Korinthos that has accumulated during this century is removed in earthquakes, the moment release will be equivalent to several Ms > 6.5 earthquakes.

THE STRAIN FIELD IN NORTHWESTERN GREECE AND THE IONIAN ISLANDS : RESULTS INFERRED FROM GPS MEASUREMENTS

Abstract Recent crustal movements detected by the analysis of repeated satellite geodetic measurements reflect the ongoing geodynamic processes in the Alpine-Mediterranean area. Superimposed on the large-scale counterclockwise rotation of the African plate, complex dynamic processes are affecting the lithospheric fragments between the African and Eurasian plates. Key features to better understand the driving forces and associated seismic activity in the Africa/Eurasia collision zone are the Calabrian and Hellenic arcs. In this paper geodynamic investigations along the West Hellenic arc are discussed. They are based on two epochs (1989 and 1993) of satellite geodetic measurements carried out using the US Global Positioning System (GPS). The results are presented in terms of relative displacements and strain rates. Within the time span of 4 years southwestern Greece has moved to the southwest relative to southeastern Italy by an average of 120 mm, increasing from 80 mm at Lefkada, in the center of the Ionian Islands, to 160 mm at the Peloponnesus. The maximum strain rate is 0.18 μstrain/a located in the vicinity of Lefkada, where anomalously high earthquake activity is observed. The data provide strong evidence for dextral strike-slip motion on the order of 25 mm/a along the Kephalonia Fault Zone (KFZ). The deformation field of the KFZ is interpreted as a transition zone between the kinematics of the Apulian platform and the West Hellenic fold and thrust belts.

Geodetic investigation of the 13 May 1995 Kozani‐Grevena (Greece) Earthquake

The Ms=6.6 13 May 1995 Kozani-Grevena earthquake struck a region of low historical seismic activity which includes a 10-year-old triangulation network in northern Greece. After the earthquake, monuments from this network were occupied with GPS to measure co-seismic displacements. Inversion of the co-seismic displacement field to yield a source mechanism is achieved by use of a hybrid simplex-Monte-Carlo method which requires no a priori constraints. The model focal mechanism agrees well with the global CMT solution and locally observed aftershocks, but implies a significantly higher scalar moment than do seismological or SAR interferometry studies, and has a longer fault length than the model based on SAR interferometry.

Sea level in the Mediterranean: a first step towards separating crustal movements and absolute sea-level variations

Abstract The SELF (SEa Level Fluctuations: geophysical interpretation and environmental impact) project has been developed and realized in the framework of the Environment Programme designed by the Commission of the European Communities. The SELF project was aimed at providing a reliable base for the determination, in the Mediterranean area, of sea-level variations which could then be used as a possible indicator of climate changes and to study the interactions taking place among the ocean, the atmosphere: and the solid Earth. The project has made it possible to define a consistent network of well-established tide gauges encompassing the Mediterranean Basin as far as the Black Sea and to determine to centimeter accuracy the tide gauge benchmark heights in a global well-defined reference system such as the one provided by the SLR/VLBI space techniques. The SELF network constitutes, for the Mediterranean, the necessary prerequisite towards achieving the actual capability to separate vertical crustal movements from true sea-level variations. This has been accomplished through the use of space techniques namely SLR, VLBI and GPS in conjunction with Water Vapor Radiometer observations and absolute gravity measurements. The analysis of the available tide gauge records has shown a high spatial coherence of the annual to multidecadal sea-level variability. Sea-level fluctuations at periods longer than two months were found to be strongly correlated with air pressure. The seasonal cycle was found to be variable in time. Relative sea-level trends determined from records longer than 30 years are less than 1.5 mm/yr. Crustal movement rates as determined from the tide gauge records are in general of the order ± 1.0 mm/yr. The geological observations have shed light on the fact that a marked variability of crustal movements occurs on both the temporal and spatial scale, and it represents a major contribution to relative sea-level fluctuations. This fact has been verified for the selected sector which belongs to one of the more geodynamically active areas of the Central Mediterranean (Aeolian Archipelago). However, this work has shown that, at least at the tide gauges included in the present study, crustal movements are small compared to the decadal to multidecadal sea-level variability but of the same order as the long-term trend in sea level, thus necessitating a careful monitoring if crustal movement is to be separated from the oceanographic contribution to relative sea-level changes.

Displacement field and fault model for the September 7, 1999 Athens earthquake inferred from ERS2 satellite radar interferometry

On September 7, 1999, a moderate (Mw=5.9) normal faulting earthquake occurred in the northwest of Athens (Hellas) causing heavy damages and casualties. Using interferometric combinations of ERS2 SAR images, we analyzed the coseismic deformation field. Two fringes are observed south of the Fill mountain, up to the coastline of the Elefsis gulf. They correspond to 56 mm increase in slant range. Modeling the earthquake as a dislocation in an elastic half-space, we inverted the interferometric data to assess the fault location and geometry and the amplitude of the coseismic slip. The model suggests -300 mm slip on an 18 km long blind fault composed of two pieces. The intersection of the fault plane with the Earth surface is located in the Fill mountain with a -N 120 o orientation.