The dynamics of the eastern Mediterranean and eastern Turkey

• Published in: Geophysical journal international 2010-12, Vol.183 (3), p.1165-1184
• Main Authors: Özeren, M. Sinan, Holt, William E.
• Format: Article
• Language: English
• Subjects: Boundary conditions; Dynamics; Dynamics of lithosphere and mantle; Estimates; Fault zone rheology; Global Positioning System; Lithosphere; Mathematical models; Space geodetic surveys; Strain rate; Stresses
• ISSN: 0956-540X; 1365-246X
• DOI: 10.1111/j.1365-246X.2010.04819.x

In this study we investigate the dynamics of the region that includes Greece, the Aegean Sea and Asia Minor. In a least-squares inversion we solve for a continuous strain rate field, and corresponding velocity field, that satisfies 872 GPS data. The estimate of the geodetic strain rate field provides constraints for our dynamic analysis. Next, we separately solve the depth integrated 3-D force balance equations for depth-integrated deviatoric stresses within the lithosphere, in which body force input comes from differences in vertically integrated vertical stress, or differences in gravitational potential energy per unit area (GPE). These GPE estimates calibrate the absolute magnitudes of deviatoric stresses that are acting within the lithosphere. Further, we investigate the sensitivity of our stress field solutions by using two different crustal structure models: one from compiled crustal structure estimates obtained primarily from relatively recent seismic observations and the other from the Crust 2.0 model. In an iterative least-squares inversion we then solve for stress field boundary conditions that, when added to the contribution of deviatoric stresses associated with GPE differences, provides a best fit to the directions of principal axes and relative magnitudes of the principal axes of the rates of strain obtained in the kinematic analysis. Robust features that arise from the boundary condition solution are NNE forcing along the southern boundary east of about 33° E (0.5–1.2 × 1012 N m−1), with a rapid anticlockwise rotation of forces to the west of this, along with an outward pulling force (∼0.4 × 1012 N m−1) directed SSW along the entire Hellenic Arc segment. This force system along the Hellenic Arc can be interpreted as a result of slab rollback. The total depth integrated 3-D deviatoric stresses in the final dynamic solution provides an excellent match to the deformation indicators throughout the region, with vertically integrated stress magnitudes of order 0.5–2.5 × 1012 N m−1. We use constraints from derived stress magnitudes, together with GPS-defined scalar values of strain rate magnitude, to define bulk effective viscosities of the lithosphere. Depth-averaged effective viscosities for the entire lithosphere are high within the Black Sea, of order 0.7–3 ×1023 Pa-s, relative to surrounding continental lithosphere. North Anatolian shear zone, northern Aegean Sea and Gulf of Corinth are characterized by low depth averaged viscosities of order