M. J. R. Wortel

Bio: M. J. R. Wortel is an academic researcher from Utrecht University. The author has contributed to research in topics: Subduction & Lithosphere. The author has an hindex of 38, co-authored 57 publications receiving 6162 citations.

Subduction and Slab Detachment in the Mediterranean-Carpathian Region

Abstract: Seismic tomography models of the three-dimensional upper mantle velocity structure of the Mediterranean-Carpathian region provide a better understanding of the lithospheric processes governing its geodynamical evolution. Slab detachment, in particular lateral migration of this process along the plate boundary, is a key element in the lithospheric dynamics of the region during the last 20 to 30 million years. It strongly affects arc and trench migration, and causes along-strike variations in vertical motions, stress fields, and magmatism. In a terminal-stage subduction zone, involving collision and suturing, slab detachment is the natural last stage in the gravitational settling of subducted lithosphere.

The Hellenic Subduction Zone: A tomographic image and its geodynamic implications

Abstract: New tomographic images of the Hellenic subduction zone demonstrate slab penetration in the Aegean Upper Mantle to depths of at least 600 km. Beneath Greece the lower part of the slab appears to be detached at a depth of about 200 km whereas it still seems to be unruptured beneath the southern Aegean. Schematically we derive minimum time estimates for the duration of the Hellenic subduction zone that range from 26 to 40 Ma. This is considerably longer than earlier estimates which vary between 5 and about 13 Ma.

Subduction history of the Tethyan region derived from seismic tomography and tectonic reconstructions

Abstract: [1] In the mantle underneath the Tethyan suture zone, large volumes of positive velocity anomalies have been imaged by seismic tomography and interpreted as the present-day signature of subducted Tethyan lithosphere. We investigate the Mesozoic-Cenozoic subduction history of the region by integrating independent information from mantle tomography and tectonic reconstructions. Three different subduction scenarios for the Tethyan oceanic lithosphere, representative for the available tectonic reconstructions, are used to predict the present thermally anomalous volumes associated with the lithospheric surface subducted since the late Mesozoic. Next, these predicted thermal volumes and their expected positions are compared to the relevant anomalous volumes derived from seismic tomographic images. In this analysis we include, among others, the possible effects of ridge subduction and slab detachment after the Cenozoic continental collisions, absolute plate motion, and slab thickening in the mantle. Our preferred subduction model comprises the opening of large back-arc oceanic basins within the Eurasian margin. The model points to slab thickening by a factor of 3 in the mantle, in which case the estimated volumes allow for active oceanic spreading (∼1–2.5 cm/yr) in the Tethyan lithosphere during convergence. Our results further indicate the occurrence of early Oligocene slab detachment underneath the northern Zagros suture zone, followed by both westward and eastward propagation of the slab tear and diachronous Eocene to Miocene slab detachment below the eastern to western Himalayas. Free sinking rates of the detached material of ∼2 cm/yr in the lower mantle provide the best fit between the tomographic mantle structure and our Tethyan subduction model.

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

Abstract: 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.

First‐ and second‐order patterns of stress in the lithosphere: The World Stress Map Project

Abstract: To date, more than 7300 in situ stress orientations have been compiled as part of the World Stress Map project. Of these, over 4400 are considered reliable tectonic stress indicators, recording horizontal stress orientations to within <±25°. Remarkably good correlation is observed between stress orientations deduced from in situ stress measurements and geologic observations made in the upper 1–2 km, well bore breakouts extending to 4–5 km depth and earthquake focal mechanisms to depths of ∼20 km. Regionally uniform stress orientations and relative magnitudes permit definition of broad-scale regional stress patterns often extending 20–200 times the approximately 20–25 km thickness of the upper brittle lithosphere. The “first-order” midplate stress fields are believed to be largely the result of compressional forces applied at plate boundaries, primarily ridge push and continental collision. The orientation of the intraplate stress field is thus largely controlled by the geometry of the plate boundaries. There is no evidence of large lateral stress gradients (as evidenced by lateral variations in stress regime) which would be expected across large plates if simple resistive or driving basal drag tractions (parallel or antiparallel to absolute motion) controlled the intraplate stress field. Intraplate areas of active extension are generally associated with regions of high topography: western U.S. Cordillera, high Andes, Tibetan plateau, western Indian Ocean plateau. Buoyancy stresses related to crustal thickening and/or lithospheric thinning in these regions dominate the intraplate compressional stress field due to plate-driving forces. These buoyancy forces are just one of several categories of “second-order” stresses, or local perturbations, that can be identified once the first-order stress patterns are recognized. These second-order stress fields can often be associated with specific geologic or tectonic features, for example, lithospheric flexure, lateral strength contrasts, as well as the lateral density contrasts which give rise to buoyancy forces. These second-order stress patterns typically have wavelengths ranging from 5 to 10+ times the thickness of the brittle upper lithosphere. A two-dimensional analysis of the amount of rotation of regional horizontal stress orientations due to a superimposed local stress constrains the ratio of the magnitude of the horizontal regional stress differences to the local uniaxial stress. For a detectable rotation of 15°, the local horizontal uniaxial stress must be at least twice the magnitude of the regional horizontal stress differences. Examples of local rotations of SHmax orientations include a 75°–85° rotation on the northeastern Canadian continental shelf possibly related to margin-normal extension derived from sediment-loading flexural stresses, a 50°–60° rotation within the East African rift relative to western Africa due to extensional buoyancy forces caused by lithospheric thinning, and an approximately 90° rotation along the northern margin of the Paleozoic Amazonas rift in central Brazil. In this final example, this rotation is hypothesized as being due to deviatoric compression oriented normal to the rift axis resulting from local lithospheric support of a dense mass in the lower crust beneath the rift (“rift pillow”). Estimates of the magnitudes of first-order (plate boundary force-derived) regional stress differences computed from modeling the source of observed local stress rotations magnitudes can be compared with regional stress differences based on the frictional strength of the crust (i.e., “Byerlee's law”) assuming hydrostatic pore pressure. The examples given here are too few to provide a definitive evaluation of the direct applicability of Byerlee's law to the upper brittle part of the lithosphere, particularly in view of uncertainties such as pore pressure and relative magnitude of the intermediate principal stresses. Nonetheless, the observed rotations all indicate that the magnitude of the local horizontal uniaxial stresses must be 1–2.5+ times the magnitude of the regional first-order horizontal stress differences and suggest that careful evaluation of such local rotations may be a powerful technique for constraining the in situ magnitude stress differences in the upper, brittle part of the lithosphere.

GPS constraints on continental deformation in the Africa-Arabia-Eurasia continental collision zone and implications for the dynamics of plate interactions

Abstract: [1] The GPS-derived velocity field (1988–2005) for the zone of interaction of the Arabian, African (Nubian, Somalian), and Eurasian plates indicates counterclockwise rotation of a broad area of the Earth's surface including the Arabian plate, adjacent parts of the Zagros and central Iran, Turkey, and the Aegean/Peloponnesus relative to Eurasia at rates in the range of 20–30 mm/yr. This relatively rapid motion occurs within the framework of the slow-moving (∼5 mm/yr relative motions) Eurasian, Nubian, and Somalian plates. The circulatory pattern of motion increases in rate toward the Hellenic trench system. We develop an elastic block model to constrain present-day plate motions (relative Euler vectors), regional deformation within the interplate zone, and slip rates for major faults. Substantial areas of continental lithosphere within the region of plate interaction show coherent motion with internal deformations below ∼1–2 mm/yr, including central and eastern Anatolia (Turkey), the southwestern Aegean/Peloponnesus, the Lesser Caucasus, and Central Iran. Geodetic slip rates for major block-bounding structures are mostly comparable to geologic rates estimated for the most recent geological period (∼3–5 Myr). We find that the convergence of Arabia with Eurasia is accommodated in large part by lateral transport within the interior part of the collision zone and lithospheric shortening along the Caucasus and Zagros mountain belts around the periphery of the collision zone. In addition, we find that the principal boundary between the westerly moving Anatolian plate and Arabia (East Anatolian fault) is presently characterized by pure left-lateral strike slip with no fault-normal convergence. This implies that “extrusion” is not presently inducing westward motion of Anatolia. On the basis of the observed kinematics, we hypothesize that deformation in the Africa-Arabia-Eurasia collision zone is driven in large part by rollback of the subducting African lithosphere beneath the Hellenic and Cyprus trenches aided by slab pull on the southeastern side of the subducting Arabian plate along the Makran subduction zone. We further suggest that the separation of Arabia from Africa is a response to plate motions induced by active subduction.