scholarly journals Neo-Tethys geodynamics and mantle convection: from extension to compression in Africa and a conceptual model for obduction

2016 ◽  
Vol 53 (11) ◽  
pp. 1190-1204 ◽  
Author(s):  
Laurent Jolivet ◽  
Claudio Faccenna ◽  
Philippe Agard ◽  
Dominique Frizon de Lamotte ◽  
Armel Menant ◽  
...  

Since the Mesozoic, Africa has been under extension with shorter periods of compression associated with obduction of ophiolites on its northern margin. Less frequent than “normal” subduction, obduction is a first order process that remains enigmatic. The closure of the Neo-Tethys Ocean, by the Upper Cretaceous, is characterized by a major obduction event, from the Mediterranean region to the Himalayas, best represented around the Arabian Plate, from Cyprus to Oman. These ophiolites were all emplaced in a short time window in the Late Cretaceous, from ∼100 to 75 Ma, on the northern margin of Africa, in a context of compression over large parts of Africa and Europe, across the convergence zone. The scale of this process requires an explanation at the scale of several thousands of kilometres along strike, thus probably involving a large part of the convecting mantle. We suggest that alternating extension and compression in Africa could be explained by switching convection regimes. The extensional situation would correspond to steady-state whole-mantle convection, Africa being carried northward by a large-scale conveyor belt, while compression and obduction would occur when the African slab penetrates the upper–lower mantle transition zone and the African plate accelerates due to increasing plume activity, until full penetration of the Tethys slab in the lower mantle across the 660 km transition zone during a 25 Myr long period. The long-term geological archives on which such scenarios are founded can provide independent time constraints for testing numerical models of mantle convection and slab–plume interactions.

2021 ◽  
Author(s):  
Matteo Desiderio ◽  
Anna J. P. Gülcher ◽  
Maxim D. Ballmer

<p>According to geochemical and geophysical observations, Earth's lower mantle appears to be strikingly heterogeneous in composition. An accurate interpretation of these findings is critical to constrain Earth's bulk composition and long-term evolution. To this end, two main models have gained traction, each reflecting a different style of chemical heterogeneity preservation: the 'marble cake' and 'plum pudding' mantle. In the former, heterogeneity is preserved in the form of narrow streaks of recycled oceanic lithosphere, stretched and stirred throughout the mantle by convection. In the latter, domains of intrinsically strong, primordial material (enriched in the lower-mantle mineral bridgmanite) may resist convective entrainment and survive as coherent blobs in the mid mantle. Microscopic scale processes certainly affect macroscopic properties of mantle materials and thus reverberate on large-scale mantle dynamics. A cross-disciplinary effort is therefore needed to constrain present-day Earth structure, yet countless variables remain to be explored. Among previous geodynamic studies, for instance, only few have attempted to address how the viscosity and density of recycled and primordial materials affect their mutual mixing and interaction in the mantle.</p><p>Here, we apply the finite-volume code <strong>STAGYY</strong> to model thermochemical convection of the mantle in a 2D spherical-annulus geometry. All models are initialized with a lower, primordial layer and an upper, pyrolitic layer (i.e., a mechanical mixture of basalt and harzburgite), as is motivated by magma-ocean solidification studies. We explore the effects of material properties on the style of mantle convection and heterogeneity preservation. These parameters include (i) the intrinsic strength of basalt (viscosity), (ii) the intrinsic density of basalt, and (iii) the intrinsic strength of the primordial material.</p><p>Our preliminary models predict a range of different mantle mixing styles. A 'marble cake'-like regime is observed for low-viscosity primordial material (~30 times weaker than the ambient mantle), with recycled oceanic lithosphere preserved as streaks and thermochemical piles accumulating near the core-mantle boundary. Conversely, 'plum pudding' primordial blobs are also preserved when the primordial material is relatively strong, in addition to the 'marble cake' heterogeneities mentioned above. Most notably, however, the rheology and the density anomaly of basalt affect the appearance of both recycled and primordial heterogeneities. In particular, they control the stability, size and geometry of thermochemical piles, the enhancement of basaltic streaks in the mantle transition zone, and they influence the style of primordial material preservation. These results indicate the important control that the physical properties of mantle constituents exert on the style of mantle convection and mixing over geologic time. Our numerical models offer fresh insights into these processes and may advance our understanding of the composition and structure of Earth's lower mantle.</p>


2021 ◽  
Author(s):  
Marcel Thielmann ◽  
Gregor Golabek ◽  
Hauke Marquardt

<p>The rheology of the Earth’s lower mantle is poorly constrained due to a lack of knowledge of the rheological behaviour of its constituent minerals. In addition, the lower mantle does not consist of only a single, but of multiple mineral phases with differing deformation behaviour. The rheology of Earth’s lower mantle is thus not only controlled by the rheology of its individual constituents (bridgmanite and ferropericlase), but also by their interplay during deformation. This is particularly important when the viscosity contrast between the different minerals is large. Experimental studies have shown that ferropericlase may be significantly weaker than bridgmanite and may thus exert a strong control on lower mantle rheology.</p><p>Here, we thus explore the impact of phase morphology on the rheology of a ferropericlase-bridgmanite mixture using numerical models. We find that elongated ferropericlase structures within the bridgmanite matrix significantly lower the effective viscosity, even in cases where no interconnected network of weak ferropericlase layers has been formed. In addition to the weakening, elongated ferropericlase layers result in a strong viscous anisotropy. Both of these effects may have a strong impact on lower mantle dynamics, which makes is necessary to develop upscaling methods to include them in large-scale mantle convection models. We develop a numerical-statistial approach to link the statistical properties of a ferropericlase-bridgmanite mixture to its effective viscosity tensor. With this approach, both effects are captured by analytical approximations that have been derived to describe the evolution of the effective viscosity (and its anisotropy) of a two-phase medium with aligned elliptical inclusions, thus allowing to include these microscale processes in large-scale mantle convection models.</p>


2007 ◽  
Vol 20 (2) ◽  
pp. 233-254 ◽  
Author(s):  
Paul R. Field ◽  
Robert Wood

Abstract Composite mean fields and probability distribution functions (PDFs) of rain rate, cloud type and cover, cloud-top temperature, surface wind velocity, and water vapor path (WVP) are constructed using satellite observations of midlatitude cyclones from four oceanic regions (i.e., the North Pacific, South Pacific, North Atlantic, and South Atlantic). Reanalysis surface pressure fields are used to ascertain the locations of the cyclone centers, onto which the satellite fields are interpolated to give a database of ∼1500 cyclones from a two-year period (2003–04). Cyclones are categorized by their strength, defined here using surface wind speed, and by their WVP, and it is found that these two measures can explain a considerable amount of the intercyclone variability of other key variables. Composite cyclones from each of the four ocean basins exhibit similar spatial structure for a given strength and WVP. A set of nine composites is constructed from the database using three strength and three WVP ranges and is used to demonstrate that the mean column relative humidity of these systems varies only slightly (0.58–0.62) for a doubling in WVP (or equivalently a 7-K rise in sea surface temperature) and a 50% increase in cyclone strength. However, cyclone-mean rain rate increases markedly with both cyclone strength and WVP, behavior that is explained with a simple warm conveyor belt model. Systemwide high cloud fraction (tops above 440 hPa) increases from 0.23 to 0.31 as cyclone strength increases by 50%, but does not vary systematically with WVP. It is suggested that the composite fields constitute useful diagnostics for evaluating the behavior of large-scale numerical models, and may provide insight into how precipitation and clouds in midlatitude cyclones respond under a changed climate.


2021 ◽  
Author(s):  
Anna J. P. Gülcher ◽  
Maxim D. Ballmer ◽  
Paul J. Tackley

<p>The nature of compositional heterogeneity in Earth’s lower mantle is a long-standing puzzle that can inform about the thermochemical evolution and dynamics of our planet. On relatively small scales (<1km), streaks of recycled oceanic crust (ROC) and lithosphere are distributed and stirred throughout the mantle, creating a “marble cake” mantle. On larger scales (10s-100s of km), compositional heterogeneity may be preserved by delayed mixing of this marble cake with either intrinsically-dense or -strong materials of e.g. primordial origin. Intrinsically-dense materials may accumulate as piles at the core-mantle boundary, while intrinsically viscous (e.g., enhanced in the strong mineral MgSiO<sub>3 </sub>bridgmanite) may survive as blobs in the mid-mantle for large timescales (i.e., as plums in the mantle “plum pudding”). So far, only few, if any, studies have quantified mantle dynamics in the presence of different types of heterogeneity with distinct physical properties.<br><br>Here, we use 2D numerical models of global-scale mantle convection to investigate the coupled evolution and mixing of (intrinsically-dense) recycled and (intrinsically-strong) primordial material. We explore the effects of ancient compositional layering of the mantle, as motivated by magma-ocean solidification studies, and the physical parameters of the primordial material. Over a wide parameter range, primordial and recycled heterogeneity is predicted to coexist with each other. Primordial material usually survives as mid-to-large scale blobs in the mid-mantle, and this preservation is largely independent on the initial primordial-material volume. In turn, recycled oceanic crust (ROC) persists as piles at the base of the mantle and as small streaks everywhere else. The robust coexistence between recycled and primordial materials in the models indicate that the modern mantle may be in a hybrid state between the “marble cake” and “plum pudding” styles.<br><br>Finally, we put our model predictions in context with geochemical studies on early Earth dynamics as well as seismic discoveries of present-day lower-mantle heterogeneity. For the latter, we calculate synthetic seismic velocities from output model fields, and compare these synthetics to tomography models, taking into account the limited resolution of seismic tomography. Because of the competing effects of compositional and thermal anomalies on S-wave velocities, it is difficult to identify mid-mantle bridgmanitic domains in seismic tomography images. This result suggests that, if present, bridgmanitic domains in the mid-mantle may be “hidden” from seismic tomographic studies, and other approaches are needed to establish the presence/absence of these domains in the present-day deep Earth.</p>


2020 ◽  
Author(s):  
Gregor J. Golabek ◽  
Anna J. P. Gülcher ◽  
Marcel Thielmann ◽  
Paul J. Tackley ◽  
Maxim D. Ballmer

<p>Rocks in the Earth’s interior are not homogeneous but consist of different mineralogical phases with different rheological properties. Deformation of heterogeneous rocks is thus also heterogeneous, and strongly depends on the rheological contrasts and spatial distribution of the mineral phases. In Earth’s lower mantle, the main rock constituents are bridgmanite (Br) and smaller amounts of ferropericlase (Fp). Bridgmanite is substantially stronger than ferropericlase [1]. Recent studies propose that lower mantle rheology is highly dependent on the relative mineral abundances and distribution of these two phases [1,2]. It has been suggested that for bridgmanite-depleted compositions, the viscosity decreases with accumulating strain due to the interconnection of the weaker ferropericlase. This implies that deformation may localize in the lower mantle, potentially aiding the formation and preservation of compositionally distinct and “hidden” reservoirs away from these regions of localized deformation [3]. Therefore, understanding the rheological nature of Br-Fp aggregates on a small-scale is crucial for assessing the dynamics of global mantle convection. Here, we address this objective with multi-scale numerical approaches.  </p><p>Using a numerical-statistical approach, a connection between ferropericlase morphology and effective rheology of Earth’s lower mantle has recently been established [4]. Results show that bulk-rock weakening depends on the topology of the weak phase as well on its rheology, but also that significant rheological weakening can already be achieved when ferropericlase does not (yet) form an interconnected weak layer.</p><p>In a second suite of models, we implement a macro-scale description of strain-weakening based on the micro-scale solutions found in [4] in a global mantle convection model to test the first-order effect of strain weakening on convection dynamics in the lower mantle. We present 2D numerical models of thermochemical convection in a spherical annulus geometry [5] that include a new implementation of tracking the strain ellipse at each tracer through time. We further allow lower mantle materials to rheologically weaken once a certain strain threshold has been reached. Preliminary results indicate that strain localizes along both up- and downwellings in the lower mantle and that rheological weakening has a stabilizing effect on these conduits. </p><p>This multi-scale approach is essential for addressing lower-mantle rheological behavior and our results form an important step towards addressing the feasibility of isolated, long-lived geochemical reservoirs in Earth’s lower mantle.</p><p>[1] Yamazaki and Karato (2001), Am. Mineral. 86, 385-391. [2] Girard et al. (2016), Science 351, 144-147. [3] Ballmer et al. (2017), Nat. Geosci. 10, 236-240. [4] Thielmann et al. (2020), Geochem. Geophys. Geosyst., doi:10.1029/2019GC008688. [5] Hernlund and Tackley (2008), Phys. Earth Planet. Int. 171, 48–54.</p>


2019 ◽  
Vol 218 (3) ◽  
pp. 2089-2105 ◽  
Author(s):  
Tianyang Lyu ◽  
Zhiyuan Zhu ◽  
Benjun Wu

SUMMARY Lithospheric plates on the Earth's surface interact with each other, producing distinctive structures comprising two descending slabs. Double-slab subduction with inward-dipping directions represents an important multiplate system that is not yet well understood. This paper presents 2-D numerical models that investigate the dynamic process of double-slab subduction with inward dipping, focussing on slab geometry and mantle transition zone upwelling flow. This unique double-slab configuration limits trench motion and causes steep downward slab movement, thus forming fold piles in the lower mantle and driving upward mantle flow between the slabs. The model results show the effects of lithospheric plate properties and lower-mantle viscosity on subducting plate kinematics, overriding plate stress and upward mantle flow beneath the overriding plate. Appropriate lower-mantle strength (such as an upper–lower mantle viscosity increase with a factor of 200) allows slabs to penetrate into the lower mantle with periodical buckling. While varying the length and thickness of a long overriding plate (≥2500) does not have a substantial effect on slab geometry, its viscosity has a marked impact on slab evolution and mantle flow pattern. When the overriding plate is strong, slabs exhibit an overturned geometry and hesitate to fold. Mantle transition zone upwelling velocity depends on the speed of descending slabs. The downward velocity of slabs with a large negative buoyancy (caused by thickness or density) is very fast, inducing a significant transition zone upwelling flow. A stiff slab slowly descends into the deep mantle, causing a small upward flow in the transition zone. In addition, the temporal variation of mantle upwelling velocity shows strong correlation with the evolution of slab folding geometry. In the double subduction system with inward-dipping directions, the mantle transition zone upwelling exhibits oscillatory rise with time. During the backward-folding stage, upwelling velocity reaches its local maximum. Our results provide new insights into the deep mantle source of intraplate volcanism in a three-plate interaction system such as the Southeast Asia region.


2018 ◽  
Vol 31 (7) ◽  
pp. 2969-2975 ◽  
Author(s):  
Catherine M. Naud ◽  
Derek J. Posselt ◽  
Susan C. van den Heever

In Naud et al., a compositing method was utilized with CloudSat– CALIPSO observations to obtain mean transects of cloud vertical distribution and surface precipitation across cold fronts, and to examine their sensitivity to the large-scale properties of the parent extratropical cyclone. This reply demonstrates the value of compositing for evaluating numerical models, and presents additional results that address the issue of the sensitivity of the initial results to the frontal detection methodology and the potential misclassification of occlusions as cold fronts. Here a sensitivity study of the cold front composite transects of cloud cover to the input datasets or the method utilized to locate the cold fronts demonstrates that these composite transects are robust and only marginally sensitive to cold front location methods. The same conclusion is reached for the robustness of the contrast between Northern and Southern Hemisphere cloud transects. While occlusions cannot directly be flagged within the database at this point, comparisons of transects obtained for subsets of cyclones of different age indicate that the misclassification of occluded fronts as cold fronts does not explain the predominance of cloud and precipitation on the warm side of the cold fronts. The strong signal on the warm side might be better explained by a predominance of forward sloping cold fronts, or the presence of the warm conveyor belt.


2020 ◽  
Vol 191 ◽  
pp. 37
Author(s):  
Adrien Romagny ◽  
Laurent Jolivet ◽  
Armel Menant ◽  
Eloïse Bessière ◽  
Agnès Maillard ◽  
...  

Slab retreat, slab tearing and interactions of slabs are first-order drivers of the deformation of the overriding lithosphere. An independent description of the tectonic evolution of the back-arc and peripheral regions is a pre-requisite to test the proposed conceptual, analogue and numerical models of these complex dynamics in 3-D. We propose here a new series of detailed kinematics and tectonic reconstructions from 35 Ma to the Present shedding light on the driving mechanisms of back-arc rifting in the Mediterranean where several back-arc basins all started to form in the Oligocene. The step-by-step backward reconstructions lead to an initial situation 35 Ma ago with two subduction zones with opposite direction, below the AlKaPeCa block (i.e. belonging to the Alboran, Kabylies, Peloritani, Calabrian internal zones). Extension directions are quite variable and extension rates in these basins are high compared to the Africa-Eurasia convergence velocity. The highest rates are found in the Western Mediterranean, the Liguro-Provençal, Alboran and Tyrrhenian basins. These reconstructions are based on shortening rates in the peripheral mountain belts, extension rates in the basins, paleomagnetic rotations, pressure-temperature-time paths of metamorphic complexes within the internal zones of orogens, and kinematics of the large bounding plates. Results allow visualizing the interactions between the Alps, Apennines, Pyrenean-Cantabrian belt, Betic Cordillera and Rif, as well as back-arc basins. These back-arc basins formed at the emplacement of mountain belts with superimposed volcanic arcs, thus with thick, hot and weak crusts explaining the formation of metamorphic core complexes and the exhumation of large portions of lower crustal domains during rifting. They emphasize the role of transfer faults zones accommodating differential rates of retreat above slab tears and their relations with magmatism. Several transfer zones are identified, separating four different kinematic domains, the largest one being the Catalan-Balearic-Sicily Transfer Zone. Their integration in the wider Mediterranean realm and a comparison of motion paths calculated in several kinematic frameworks with mantle fabric shows that fast slab retreat was the main driver of back-arc extension in this region and that large-scale convection was a subsidiary driver for the pre-8 Ma period, though it became dominant afterward. Slab retreat and back-arc extension was mostly NW-SE until ∼ 20 Ma and the docking of the AlKaPeCa continental blocks along the northern margin of Africa induced a slab detachment that propagated eastward and westward, thus inducing a change in the direction of extension from NW-SE to E-W. Fast slab retreat between 32 and 8 Ma and induced asthenospheric flow have prevented the transmission of the horizontal compression due to Africa-Eurasia convergence from Africa to Eurasia and favored instead upper-plate extension driven by slab retreat. Once slab retreat had slowed down in the Late Miocene, this N-S compression was felt and recorded again from the High Atlas to the Paris Basin.


GeoArabia ◽  
2015 ◽  
Vol 20 (4) ◽  
pp. 67-122
Author(s):  
Philipp Wolpert ◽  
Martin Bartenbach ◽  
Peter Suess ◽  
Randolf Rausch ◽  
Thomas Aigner ◽  
...  

ABSTRACT Uppermost Jurassic–Lower Cretaceous carbonates of the Sulaiy Formation are well exposed at the type locality Dahal Hit, and along the entire natural escarpment near Ar Riyad, the capital of the Kingdom of Saudi Arabia. This study provides a facies and sequence-stratigraphic analysis based on detailed sedimentological and gamma-ray logging of 12 outcrop sections. The sections represent the Sulaiy Formation along a 60 km-long outcrop belt, including the Hith-Sulaiy transition in a large solution cavity named Dahal Hit, situated south of Ar Riyad. The latter section is studied in detail because it is the only locality in Saudi Arabia where the Hith Anhyrite (Hith Formation in this study) to the Sulaiy Formation transition crops out. Ten lithofacies types were identified for the Sulaiy Formation including potential reservoirs such as oolitic cross-bedded grainstones, biostromal boundstones, and bioclast-rich, graded pack-to-grainstones. Lithofacies types are grouped into five facies associations: (1) offshoal, (2) transition zone, foreshoal, (4) shoal margin, and (5) shoal, distributed along a carbonate ramp. Their vertical stacking pattern revealed a systematic hierarchy of cyclicity consisting of small-scale cycles, medium-scale cycle sets and two large-scale sequences for the Sulaiy Formation. Four cycle motifs, with an average thickness of 2–4 m, are present: (1) offshoal to transition zone cycle motif, (2) offshoal to foreshoal cycle motif, (3) transition zone to shoal margin cycle motif, and foreshoal to shoal margin cycle motif. A total of 15 cycle sets, ranging between 8 and 12 m in thickness each, were interpreted. They were correlated, where possible, across the study area. Three types of medium-scale cycle sets are observed: (1) offshoal to shoal cycle set motif, (2) offshoal to foreshoal cycle set motif, and (3) shoal margin to offshoal cycle set motif. The Lower Sulaiy Sequence consists of twelve cycle sets and is interpreted to contain two Arabian Plate maximum flooding surfaces (MFS): (1) Upper Tithonian MFS J110 (147 Ma) in its lowermost part, which is interpreted to be the time-equivalent of the Manifa reservoir in subsurface Arabia. (2) Lower Berriasian MFS K10 (144 Ma) in the seventh-up cycle set. The Upper Sulaiy Sequence is only represented in the Wadi Nisah Section and is believed to be incomplete because the Sulaiy/Yamama Formation boundary was not included in our study. It is presumed to contain Upper Berriasian MFS K20 (141 Ma).


I wish to suggest a mechanism for ocean-floor spreading, other than deep mantle convection cells. If the geotherm intersects the mantle solidus at the base of the upper mantle, a partly melted layer will result. It will be less dense than the overlying unmelted garnet peridotite lithosphere and the situation will be gravitationally unstable. It will be more voluminous than the overlying unmelted mantle and will distend the lithosphere which will break up into plates. It will have low rigidity and decouple those lithospheric plates from the underlying transition zone and lower mantle. Melted material, perhaps with crystals, will escape in two ways.


Sign in / Sign up

Export Citation Format

Share Document