Wim Spakman

Bio: Wim Spakman is an academic researcher from Utrecht University. The author has contributed to research in topics: Subduction & Lithosphere. The author has an hindex of 67, co-authored 159 publications receiving 18726 citations. Previous affiliations of Wim Spakman include École Normale Supérieure & University of Oslo.

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.

Closing the gap between regional and global travel time tomography

Abstract: Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al. [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [Spakman and Bijwaard, 1998]. Second, we exploit the global data set of Engdahl et al. [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used (P, pP, and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Zagros orogeny: a subduction-dominated process

Abstract: This paper presents a synthetic view of the geodynamic evolution of the Zagros orogen within the frame of the Arabia-Eurasia collision. The Zagros orogen and the Iranian plateau preserve a record of the long-standing convergence history between Eurasia and Arabia across the Neo-Tethys, from subduction/obduction processes to present-day collision (from ~ 150 to 0 Ma). We herein combine the results obtained on several geodynamic issues, namely the location of the oceanic suture zone, the age of oceanic closure and collision, the magmatic and geochemical evolution of the Eurasian upper plate during convergence (as testified by the successive Sanandaj-Sirjan, Kermanshah and Urumieh-Dokhtar magmatic arcs), the P-T-t history of the few Zagros blueschists, the convergence characteristics across the Neo-Tethys (kinematic velocities, tomographic constraints, subduction zones and obduction processes), together with a survey of recent results gathered by others. We provide lithospheric-scale reconstructions of the Zagros orogen from ~ 150 to 0 Ma across two SW-NE transects. The evolution of the Zagros orogen is also compared to those of the nearby Turkish and Himalayan orogens. In our geotectonic scenario for the Zagros convergence, we outline three main periods/regimes: (1) the Mid to Late Cretaceous (115-85 Ma) corresponds to a distinctive period of perturbation of subduction processes and interplate mechanical coupling marked by blueschist exhumation and upper-plate fragmentation, (2) the Paleocene-Eocene (60-40 Ma) witnesses slab break-off, major shifts in arc magmatism and distributed extension within the upper plate, and (3) from the Oligocene onwards (~ 30-0 Ma), collision develops with a progressive SW migration of deformation and topographic build-up (Sanandaj-Sirjan Zone: 20-15 Ma, High Zagros: ~12-8 Ma; Simply Folded Belt: 5-0 Ma) and with partial slab tear at depths (~10 Ma to present). Our reconstructions underline the key role played by subduction throughout the whole convergence history. We finally stress that such a long-lasting subduction system with changing boundary conditions also makes the Zagros orogen an ideal natural laboratory for subduction processes.

Geodynamics of flat subduction: Seismicity and tomographic constraints from the Andean margin

Abstract: The cause and geodynamic impact of fiat subduction are investigated. First, the 1500 km long Peru fiat slab segment is examined. Earthquake hypocenter data image two morphologic highs in the subducting Nazca Plate which correlate with the posi- tions of subducted oceanic plateaus. Travel time tomo- graphic images confirm the three-dimensional slab ge- ometry and suggest a lithospheric tear may bound the NW edge of the fiat slab segment, with possible slab de- tachment occurring down dip as well. Other fiat slab re- gions worldwide are discussed: central Chile, Ecuador, NW Colombia, Costa Rica, Mexico, southern Alaska, SW Japan, and western New Guinea. Flat subduction is shown to be a widespread phenomenon, occuring in 10% of modern convergent margins. in nearly all these cases, as a spatial and temporal correlation is observed between subducting oceanic plateaus and fiat subduc- tion, we conclude that fiat subduction is caused pri- marily by (1) the buoyancy of thickened oceanic crust of moderate to young age and (2) a delay in the basalt to eclogite transition due to the cool thermal structure of two overlapping lithospheres. A statistical analysis of seismicity along the entire length of the Andes demon- strates that seismic energy release in the upper plate at a distance of 250-800 km from the trench is on aver- age 3-5 times greater above fiat slab segments than for adjacent steep slab segments. We propose this is due to higher interplate coupling and the cold, strong rhe- ology of the overriding lithosphere which thus enables stress and deformation to be transmitted hundreds of kilometers into the heart of the upper plate.

Greater India Basin hypothesis and a two-stage Cenozoic collision between India and Asia

Abstract: Cenozoic convergence between the Indian and Asian plates pro- duced the archetypical continental collision zone comprising the Himalaya mountain belt and the Tibetan Plateau. How and where India-Asia convergence was accommodated after collision at or be- fore 52 Ma remains a long-standing controversy. Since 52 Ma, the two plates have converged up to 3,600 35 km, yet the upper crustal shortening documented from the geological record of Asia and the Himalaya is up to approximately 2,350-km less. Here we show that the discrepancy between the convergence and the shortening can be explained by subduction of highly extended continental and oceanic Indian lithosphere within the Himalaya be- tween approximately 50 and 25 Ma. Paleomagnetic data show that this extended continental and oceanic "Greater India" promontory resulted from 2,675 700 km of North-South extension between 120 and 70 Ma, accommodated between the Tibetan Himalaya and cratonic India. We suggest that the approximately 50 Ma "India"- Asia collision was a collision of a Tibetan-Himalayan microconti- nent with Asia, followed by subduction of the largely oceanic Greater India Basin along a subduction zone at the location of the Greater Himalaya. The "hard" India-Asia collision with thicker and contiguous Indian continental lithosphere occurred around 25-20 Ma. This hard collision is coincident with far-field deforma- tion in central Asia and rapid exhumation of Greater Himalaya crys- talline rocks, and may be linked to intensification of the Asian monsoon system. This two-stage collision between India and Asia is also reflected in the deep mantle remnants of subduction imaged with seismic tomography.

Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book

Abstract: This book is a tutorial written by researchers and developers behind the FEniCS Project and explores an advanced, expressive approach to the development of mathematical software. The presentation spans mathematical background, software design and the use of FEniCS in applications. Theoretical aspects are complemented with computer code which is available as free/open source software. The book begins with a special introductory tutorial for beginners. Followingare chapters in Part I addressing fundamental aspects of the approach to automating the creation of finite element solvers. Chapters in Part II address the design and implementation of the FEnicS software. Chapters in Part III present the application of FEniCS to a wide range of applications, including fluid flow, solid mechanics, electromagnetics and geophysics.

Geologically current plate motions

Abstract: SUMMARY We describe best-fitting angular velocities and MORVEL, a new closure-enforced set of angular velocities for the geologically current motions of 25 tectonic plates that collectively occupy 97 per cent of Earth's surface. Seafloor spreading rates and fault azimuths are used to determine the motions of 19 plates bordered by mid-ocean ridges, including all the major plates. Six smaller plates with little or no connection to the mid-ocean ridges are linked to MORVEL with GPS station velocities and azimuthal data. By design, almost no kinematic information is exchanged between the geologically determined and geodetically constrained subsets of the global circuit—MORVEL thus averages motion over geological intervals for all the major plates. Plate geometry changes relative to NUVEL-1A include the incorporation of Nubia, Lwandle and Somalia plates for the former Africa plate, Capricorn, Australia and Macquarie plates for the former Australia plate, and Sur and South America plates for the former South America plate. MORVEL also includes Amur, Philippine Sea, Sundaland and Yangtze plates, making it more useful than NUVEL-1A for studies of deformation in Asia and the western Pacific. Seafloor spreading rates are estimated over the past 0.78 Myr for intermediate and fast spreading centres and since 3.16 Ma for slow and ultraslow spreading centres. Rates are adjusted downward by 0.6–2.6 mm yr−1 to compensate for the several kilometre width of magnetic reversal zones. Nearly all the NUVEL-1A angular velocities differ significantly from the MORVEL angular velocities. The many new data, revised plate geometries, and correction for outward displacement thus significantly modify our knowledge of geologically current plate motions. MORVEL indicates significantly slower 0.78-Myr-average motion across the Nazca–Antarctic and Nazca–Pacific boundaries than does NUVEL-1A, consistent with a progressive slowdown in the eastward component of Nazca plate motion since 3.16 Ma. It also indicates that motions across the Caribbean–North America and Caribbean–South America plate boundaries are twice as fast as given by NUVEL-1A. Summed, least-squares differences between angular velocities estimated from GPS and those for MORVEL, NUVEL-1 and NUVEL-1A are, respectively, 260 per cent larger for NUVEL-1 and 50 per cent larger for NUVEL-1A than for MORVEL, suggesting that MORVEL more accurately describes historically current plate motions. Significant differences between geological and GPS estimates of Nazca plate motion and Arabia–Eurasia and India–Eurasia motion are reduced but not eliminated when using MORVEL instead of NUVEL-1A, possibly indicating that changes have occurred in those plate motions since 3.16 Ma. The MORVEL and GPS estimates of Pacific–North America plate motion in western North America differ by only 2.6 ± 1.7 mm yr−1, ≈25 per cent smaller than for NUVEL-1A. The remaining difference for this plate pair, assuming there are no unrecognized systematic errors and no measurable change in Pacific–North America motion over the past 1–3 Myr, indicates deformation of one or more plates in the global circuit. Tests for closure of six three-plate circuits indicate that two, Pacific–Cocos–Nazca and Sur–Nubia–Antarctic, fail closure, with respective linear velocities of non-closure of 14 ± 5 and 3 ± 1 mm yr−1 (95 per cent confidence limits) at their triple junctions. We conclude that the rigid plate approximation continues to be tremendously useful, but—absent any unrecognized systematic errors—the plates deform measurably, possibly by thermal contraction and wide plate boundaries with deformation rates near or beneath the level of noise in plate kinematic data.