Using geodynamic modeling to test plate tectonic scenarios of the Mediterranean-Alpine area

Author(s):  
Boris Kaus ◽  
Eline Le Breton ◽  
Georg Reuber ◽  
Christian Schuler

<p>Using geological and geophysical data, it is possible to reconstruct the past motion of tectonic plates involved in the Alpine orogeny and propose possible scenarios for their geological evolution. However, those scenarios have not yet been tested for geodynamic consistency.</p><p>Here, we perform 3D thermomechanical geodynamic simulations of the Mediterranean-Alpine area starting with a plate tectonic reconstruction at 20 Ma based on the work of Le Breton et al. (2017). The models include viscoelastoplastic rheologies and a free surface, and thus simulate the spontaneous occurrence of shear zones as well as the development of topography. Whereas some aspects of the tectonic reconstruction are well constrained (i.e. past position of the plates and subduction-collision fronts), many details such as the dip and length of the subducted plates, their thermal structure as well as their rheology, are unknown. The models are run forward in time to see to which extent they are consistent with the kinematic reconstructions. Perhaps unsurprisingly, our initial modelling attempts show a wide variety of behavior, including slab break off events and slab rollbacks in the wrong directions. Yet, in all cases tested so far, the model evolution does not reproduce the present-day geological setting, with Adria frequently moving too far towards the east and breaking apart internally, frequently no Alpine chain forming and in some cases new subduction zones developing within the Western Mediterranean that swallow Sardinia and Corsica.</p><p>Reproducing geological scenarios with thermomechanical geodynamic modelling thus requires substantial additional work, both from the modelling side (testing the effect of uncertain parameters on the behaviour of plates and subduction zones), as well as from the plate reconstruction side (assessing which parameters are well constrained and need to be reproduced).  Nevertheless, interesting insights can already be obtained from our models, and in our presentation, we will highlight some of the links between interacting subducting plates and plate motion.</p><p>Le Breton E, Handy MR, Molli G, Ustaszewski K (2017) Post-20 Ma Motion of the Adriatic Plate: New Constraints From Surrounding Orogens and Implications for Crust-Mantle Decoupling. Tectonics 36:3135–3154. doi: 10.1002/2016TC004443</p><div> <div> <div> </div> </div> </div>

2020 ◽  
Author(s):  
Derya Gürer ◽  
Roi Granot ◽  
Douwe J.J. van Hinsbergen

<p>The relative motions of the tectonic plates show remarkable variation throughout Earth’s history. Major changes in relative motion between the tectonic plates are traditionally viewed as spatially and temporally isolated events linked to forces acting on plate boundaries (i.e., formation of same-dip double subduction zones, changes in the strength of the boundary), or thought to be associated with mantle dynamics. A Cretaceous global plate reorganization event has been postulated to have affected all major plates. The Cretaceous ‘swing’ in Africa-Eurasia relative plate motion provides an ideal test-bed for assessing the temporal and spatial evolution of both relative plate motions and surrounding geological markers. Here we show a novel plate kinematic model for the closure of the Tethys Ocean by implementing intra-Cretaceous Quiet Zone time markers and combine the results with the geological constraints found along the convergent plate boundary. Our results allow to assess the order, causes and consequences of geological events and unravel a chain of tectonic events that set off with the onset of horizontally-forced double subduction ~105 Myr ago, followed by a 40 Myr long period of acceleration of the Africa relative to Eurasia that peaked at 80 Myr ago (at rates four times as high as previously predicted). This acceleration, which was likely caused by the pull of two same-dip subduction zones was followed by a sharp decrease in plate velocity, when double subduction terminated with ophiolite obduction onto the African margin. These tectonic forces acted on the eastern half of the Africa-Eurasia plate boundary, which led to counterclockwise rotation of Africa and sparked new subduction zones in the western Mediterranean region. Our analysis identifies the Cretaceous double subduction episode between Africa and Eurasia as a link in the global plate tectonic chain reaction and provides a dynamic view on plate reorganizations.</p>


2021 ◽  
Author(s):  
Rosalia Lo Bue ◽  
Manuele Faccenda ◽  
Jianfeng Yang

<p>In the geodynamic context of the slow Africa-Europe plates convergence, the Central-Western Mediterranean region has been involved in a complex subduction process, which in the last 30 Myr was characterized by the rapid retreat of the Ionian slab, the opening of back-arc extensional basins (i.e., Liguro-Provençal, Algerian, Alboran, and Tyrrhenian basins) and episodes of slab lateral tearing, segmentation and break-off.  A proper modelling of 3-D mantle flow evolution beneath the Mediterranean could provide important clarifications about the complex mantle dynamics of this region and help us understanding the interaction between surface tectono-magmatic processes and mantle convection patterns. </p><p>The mantle flow and its relations with plate horizontal and vertical motions can be determined by measuring seismic anisotropy generated by strain-induced lattice/crystal preferred orientation (LPO/CPO) of intrinsically anisotropic minerals. Seismic anisotropy is widespread in the Mediterranean and it shows an intricate pattern that likely has some relations with the recent (20-30 Myr) behavior of subducting slabs. The extrapolation of the mantle flow from seismic anisotropy is neither simple nor always warranted, especially at subduction zones where complex and non-steady-state 3D flow patterns may establish.  A promising approach, which helps reducing the number of plausible models that can explain a given anisotropy dataset, is to compare seismic measurements with predictions of numerical and experimental flow models (Long et al.,2007). Recently, Faccenda and Capitanio (2013) and Faccenda (2014) have extended this methodology to account for the non-steady state evolution typical of many subduction zones, yielding mantle fabrics that are physically consistent with the deformation history.</p><p>In this study, we apply a similar modelling approach to the complex Central-Western Mediterranean convergent margin. We use the wealth of observations from the Mediterranean region available in the literature to design and calibrate 3D thermo-mechanical subduction modelling. We test different initial configurations defined at 30 Ma according to the paleogeographic and tectonic reconstructions derived from (Lucente and Speranza, 2001; Carminati et al., 2012; van Hinsbergen et al., 2014) in order to optimize the fit between predicted and observed slabs position and obtain a final model configuration resembling the present-day surface and deeper structures.</p><p>In particular, here we want to evaluate the influence on rollback rates, trench shape and the occurrence and timing of slab tears (Mason et al., 2010) of structural heterogeneities within the Adria plate as proposed by (Lucente and Speranza, 2001). In all models, subduction migrates south-eastward driven by the subducting oceanic lithosphere, and slab lateral tearing or break-off occurs when a continental margin enters the trench. More importantly, we show that the presence of a stiffer continental promontory in central Adria together with a thinned continental margin in the Umbria-Marche region plays a fundamental role on (i) the development of a slab window below the Central Apennines, (ii) the retreat of the Northern Apenninic trench till the Adriatic Sea, and (iii) the retreat of the Ionian slab till the present-day position.</p>


Geosphere ◽  
2020 ◽  
Vol 16 (6) ◽  
pp. 1408-1424 ◽  
Author(s):  
Geoffrey A. Abers ◽  
Peter E. van Keken ◽  
Cian R. Wilson

Abstract The plate interface undergoes two transitions between seismogenic depths and subarc depths. A brittle-ductile transition at 20–50 km depth is followed by a transition to full viscous coupling to the overlying mantle wedge at ∼80 km depth. We review evidence for both transitions, focusing on heat-flow and seismic-attenuation constraints on the deeper transition. The intervening ductile shear zone likely weakens considerably as temperature increases, such that its rheology exerts a stronger control on subduction-zone thermal structure than does frictional shear heating. We evaluate its role through analytic approximations and two-dimensional finite-element models for both idealized subduction geometries and those resembling real subduction zones. We show that a temperature-buffering process exists in the shear zone that results in temperatures being tightly controlled by the rheological strength of that shear zone’s material for a wide range of shear-heating behaviors of the shallower brittle region. Higher temperatures result in weaker shear zones and hence less heat generation, so temperatures stop increasing and shear zones stop weakening. The net result for many rheologies are temperatures limited to ≤350–420 °C along the plate interface below the cold forearc of most subduction zones until the hot coupled mantle is approached. Very young incoming plates are the exception. This rheological buffering desensitizes subduction-zone thermal structure to many parameters and may help explain the global constancy of the 80 km coupling limit. We recalculate water fluxes to the forearc wedge and deep mantle and find that shear heating has little effect on global water circulation.


2021 ◽  
Author(s):  
Erik van der Wiel ◽  
Cedric Thieulot ◽  
Wim Spakman ◽  
Douwe van Hinsbergen

<p>Long-lived, Mesozoic-Cenozoic subduction zones such as the Pacific slab under the Americas and the Tethyan slab under Eurasia consumed thousands of kms of lithosphere of which remnants are detected in today’s mantle by seismic tomography. Major differences, however, in subduction zone evolution occurred between these systems which include strong variations in subduction rate, slab morphological evolution, and trench motion, which all appear mostly to be accommodated in the upper 1000 km of the mantle (van der Meer et al. 2018). Furthermore, sinking rates of slabs below this zone tend to be similar for different subduction systems and an order of magnitude smaller than their plate/subduction velocities. Working from the premise that the mantle rheology that accommodated these subduction systems is basically similar, although still poorly constrained, we test the hypothesis that the contrasting evolution of these subduction systems is primarily tied in with the global plate tectonic forcing of subduction.</p><p>It is generally accepted that plate motion is primarily driven by slab pull with contributions from ridge push, rather than the drag of the underlying mantle. If correct, numerical subduction models should be able to obtain upper as well as lower mantle subduction velocities and sinking rates similar to those reconstructed from geological records. We are at the start of this investigation and will present the numerical model setup, modeling strategy, and preliminary results of a 2-D subduction modelling experiment. We implement a 2D-cylindrical model setup for solving the conservation of momentum, mass and energy with the open-source geodynamics code ASPECT (Kronbichler et al. 2012) using a nonlinear visco-plastic rheology and including the major phase changes. Our focus is on the possible role of the absolute motion of the subducting and overriding plates in concert with slab pull variation reconstructed from plate tectonic evolution models, while in both subduction cases the same (partly nonlinear) mantle rheological processes are required to accommodate slab morphology change and slab sinking. Kinematic modelling constraints are derived from global plate tectonic evolution models, while the tomographically inferred present-day stage provides the end-stage geometry of slabs.</p><p>van der Meer, D. G., Van Hinsbergen, D. J., & Spakman, W. (2018). Atlas of the underworld: Slab remnants in the mantle, their sinking history, and a new outlook on lower mantle viscosity. Tectonophysics, 723, 309-448.</p><p>Kronbichler, M., Heister, T., & Bangerth, W. (2012). High accuracy mantle convection simulation through modern numerical methods. Geophysical Journal International, 191(1), 12-29.</p>


2016 ◽  
Author(s):  
Simin Gao ◽  
◽  
Margarete Jadamec

2017 ◽  
Vol 14 (23) ◽  
pp. 5425-5439 ◽  
Author(s):  
Wenchao Cao ◽  
Sabin Zahirovic ◽  
Nicolas Flament ◽  
Simon Williams ◽  
Jan Golonka ◽  
...  

Abstract. Paleogeographic reconstructions are important to understand Earth's tectonic evolution, past eustatic and regional sea level change, paleoclimate and ocean circulation, deep Earth resources and to constrain and interpret the dynamic topography predicted by mantle convection models. Global paleogeographic maps have been compiled and published, but they are generally presented as static maps with varying map projections, different time intervals represented by the maps and different plate motion models that underlie the paleogeographic reconstructions. This makes it difficult to convert the maps into a digital form and link them to alternative digital plate tectonic reconstructions. To address this limitation, we develop a workflow to restore global paleogeographic maps to their present-day coordinates and enable them to be linked to a different tectonic reconstruction. We use marine fossil collections from the Paleobiology Database to identify inconsistencies between their indicative paleoenvironments and published paleogeographic maps, and revise the locations of inferred paleo-coastlines that represent the estimated maximum transgression surfaces by resolving these inconsistencies. As a result, the consistency ratio between the paleogeography and the paleoenvironments indicated by the marine fossil collections is increased from an average of 75 % to nearly full consistency (100 %). The paleogeography in the main regions of North America, South America, Europe and Africa is significantly revised, especially in the Late Carboniferous, Middle Permian, Triassic, Jurassic, Late Cretaceous and most of the Cenozoic. The global flooded continental areas since the Early Devonian calculated from the revised paleogeography in this study are generally consistent with results derived from other paleoenvironment and paleo-lithofacies data and with the strontium isotope record in marine carbonates. We also estimate the terrestrial areal change over time associated with transferring reconstruction, filling gaps and modifying the paleogeographic geometries based on the paleobiology test. This indicates that the variation of the underlying plate reconstruction is the main factor that contributes to the terrestrial areal change, and the effect of revising paleogeographic geometries based on paleobiology is secondary.


Author(s):  
J.A. Reina-Hervás ◽  
J.E. García Raso ◽  
M.E. Manjón-Cabeza

The capture of a specimen of Sphoeroides spengleri (Osteichthyes: Tetraodontidae), 17 December 2000 and 29·7 mm total length, from the Málaga coast (Alborán Sea, western Mediterranean) represents the first record of a new alien species for Mediterranean waters.


2021 ◽  
Author(s):  
Nalan Lom ◽  
Abdul Qayyum ◽  
Derya Gürer ◽  
Douwe G. van der Meer ◽  
Wim Spakman ◽  
...  

<p>Iran is a mosaic of continental blocks that are surrounded by Tethyan oceanic relics. Remnants of these oceanic rock assemblages are exposed around the Central Iranian Microcontinent (CIM), discretely along the Sanandaj-Sirjan Zone and in Jaz-Murian. The ophiolite belts surrounding the CIM are mainly assumed to represent narrow back-arc basins that opened in Cretaceous and closed before the Eocene. Although these ophiolites are exposed as small pieces on continental crust today, they represent oceans wide enough to form supra-subduction ophiolites and arc-related magmatic rocks which suggest that their palaeogeographic width was at least some hundreds of kilometers. Current models for the palaeogeographic dimension, opening and closure of these basins are highly schematic. They usually seem plausible in two-dimensional reconstructions, however a single three-dimensional model explaining whole Iran and its surrounding regions has not been fully accomplished.  This is mostly because while the geological record provides constraints on the origin and ages of the subducted ocean floor, it provides limited information about onset and cessation of the subduction and almost no constraints on the dimension of these oceans and the subduction zones that consumed them.</p><p>In this study, we follow a novel approach in estimating the dimension and evolution of these back-arc basin by using seismic tomography. Seismic tomography has revealed that we can image and trace subducted lithosphere relics. Imaged mantle structure is now being used to link sinking slabs with sutures and to define shape of a slab. Systematic comparison of regions where the timing of subduction is reasonably well constrained by geological data showed that slabs sink gradually through the mantle at rates more or less the same. This perspective enabled us to study slab shape as a function of absolute trench motion. While mantle stationary trenches tend to create steep slabs or slab walls, the flat-lying segments are formed where the overlying trenches are mobile relative to the mantle, normal facing during roll-back, overturned during slab advance.  Under the assumption of vertical sinking after break-off, it is also possible to locate the palaeo-trenches.  When combined with absolute plate motion reconstructions, tomographically determined volume and size of the subducted lithosphere can also be used to estimate the size/width of the prehistoric oceans. To this end, we build on and further develop concepts that relate absolute trench motion during subduction to modern slab geometry to evaluate the possible range of dimensions associated with opening and closure of the Iranian back-arc basins.</p>


2018 ◽  
Vol 18 (10) ◽  
pp. 7287-7312 ◽  
Author(s):  
Arineh Cholakian ◽  
Matthias Beekmann ◽  
Augustin Colette ◽  
Isabelle Coll ◽  
Guillaume Siour ◽  
...  

Abstract. The simulation of fine organic aerosols with CTMs (chemistry–transport models) in the western Mediterranean basin has not been studied until recently. The ChArMEx (the Chemistry-Aerosol Mediterranean Experiment) SOP 1b (Special Observation Period 1b) intensive field campaign in summer of 2013 gathered a large and comprehensive data set of observations, allowing the study of different aspects of the Mediterranean atmosphere including the formation of organic aerosols (OAs) in 3-D models. In this study, we used the CHIMERE CTM to perform simulations for the duration of the SAFMED (Secondary Aerosol Formation in the MEDiterranean) period (July to August 2013) of this campaign. In particular, we evaluated four schemes for the simulation of OA, including the CHIMERE standard scheme, the VBS (volatility basis set) standard scheme with two parameterizations including aging of biogenic secondary OA, and a modified version of the VBS scheme which includes fragmentation and formation of nonvolatile OA. The results from these four schemes are compared to observations at two stations in the western Mediterranean basin, located on Ersa, Cap Corse (Corsica, France), and at Cap Es Pinar (Mallorca, Spain). These observations include OA mass concentration, PMF (positive matrix factorization) results of different OA fractions, and 14C observations showing the fossil or nonfossil origins of carbonaceous particles. Because of the complex orography of the Ersa site, an original method for calculating an orographic representativeness error (ORE) has been developed. It is concluded that the modified VBS scheme is close to observations in all three aspects mentioned above; the standard VBS scheme without BSOA (biogenic secondary organic aerosol) aging also has a satisfactory performance in simulating the mass concentration of OA, but not for the source origin analysis comparisons. In addition, the OA sources over the western Mediterranean basin are explored. OA shows a major biogenic origin, especially at several hundred meters height from the surface; however over the Gulf of Genoa near the surface, the anthropogenic origin is of similar importance. A general assessment of other species was performed to evaluate the robustness of the simulations for this particular domain before evaluating OA simulation schemes. It is also shown that the Cap Corse site presents important orographic complexity, which makes comparison between model simulations and observations difficult. A method was designed to estimate an orographic representativeness error for species measured at Ersa and yields an uncertainty of between 50 and 85 % for primary pollutants, and around 2–10 % for secondary species.


Sign in / Sign up

Export Citation Format

Share Document