Weak Shear Zones Research Articles

Structural Evolution, Exhumation Rates, and Rheology of the European Crust During Alpine Collision: Constraints From the Rotondo Granite—Gotthard Nappe

AbstractThe rheology of crystalline units controls the large‐scale deformation geometry and dynamics of collisional orogens. Defining a time‐constrained rheological evolution of such units may help unravel the details of collisional dynamics. Here, we integrate field analysis, pseudosection calculations and in situ garnet U–Pb and mica Rb–Sr geochronology to define the structural and rheological evolution of the Rotondo granite (Gotthard nappe, Central Alps). We identify a sequence of four (D1–D4) deformation stages. Pre‐collisional D1 brittle faults developed before Alpine peak metamorphism, which occurred at 34–20 Ma (U–Pb garnet ages) at 590 ± 25°C and 0.9 ± 0.1 GPa. The reactivation of D1 structures controlled the rheological evolution, from D2 reverse mylonitic shearing at amphibolite facies (520 ± 40°C and 0.8 ± 0.1 GPa) at 18–20 Ma (white mica Rb–Sr ages), to strike‐slip, brittle‐ductile shearing at greenschist‐facies D3 (395 ± 25°C and 0.4 ± 0.1 GPa) at 14–15 Ma (white mica and biotite Rb–Sr ages), and then to D4 strike‐slip faulting at shallow conditions. Although highly misoriented for the Alpine collisional stress orientation, D1 brittle structures controlled the localization of D2 ductile mylonites accommodating fast (∼3 mm/yr) exhumation rates due to their weak shear strength (<10 MPa). This structural and rheological evolution is common across External Crystalline Massifs (e.g., Aar, Mont Blanc), suggesting that the European upper crust was extremely weak during Alpine collision, its strength controlled by weak ductile shear zones localized on pre‐collisional deformation structures, that in turn controlled localized exhumation at the scale of the orogen.

A case study on the stability of a big underground powerhouse cavern cut by an interlayer shear zone in the China Baihetan hydropower plant

AbstractThe big underground powerhouse cavern of the China Baihetan hydropower plant is 438 m long, 34 m wide, and 88.7 m high. It is cut by a weak interlayer shear zone and its high sidewall poses a huge stability problem. This paper reports our successful solution of this problem through numerical simulations and a replacement‐tunnel scheme in the detailed design stage and close site monitoring in the excavation stage. Particularly, in the detail design stage, mechanical parameters of the shear zone were carefully determined through laboratory experiments and site tests. Then, deformation of the surrounding rocks and the shear zone under high in situ stress conditions was predicted using 3 Dimensional Distinct Element Code (3DEC). Subsequently, a replacement‐tunnel scheme was proposed for the treatment on the shear zone to prevent severe unloading relaxation of surrounding rocks. In the construction period, excavation responses were closely monitored on deformations of surrounding rocks and the shear zone. The effect of local cracking in the replacement tunnels on sidewall stability was evaluated using the strength reduction method. These monitoring results were compared with the predicted numerical simulation in the detailed design stage. It is found that the shear zone greatly modified the deformation mode of the cavern surrounding rocks. Without any treatment, rock mass deformation on the downstream sidewall was larger than 125 mm and the shearing deformation of the shear zone was 60–70 mm. These preset replacement tunnels can reduce not only the unloading and relaxation of rock masses but also the maximum shearing deformation of the shear zone by 10–20 mm. The predictions by numerical simulation were in good agreement with the monitoring results. The proposed tunnel‐replacement scheme can not only restrain the shear zone deformation but also enhance the safety of surrounding rocks and concrete tunnels. This design procedure offers a good reference for interaction between a big underground cavern and a weak layer zone in the future.

Weakness of the Indian Lower Crust Beneath the Himalaya Inferred From Postseismic Deformation of the 2015 Mw 7.8 Gorkha Earthquake

AbstractThe 2015 Mw 7.8 Gorkha (Nepal) earthquake induced prolonged postseismic deformation extending northward beyond the Yarlung Zangbo Suture, which provides unique opportunities to better understand the lithospheric rheology in Himalaya and southern Tibet. Here, we used the first 5‐year Global Positioning System observations to study the main postseismic processes following this event, including viscoelastic relaxation and afterslip, based on a three‐dimensional finite element model. We considered a realistic geometry of the underthrusting Indian plate according to various geophysical images. We found that the models with a uniform elastic Indian lower crust fail to fit the vertical displacements. A heterogeneous Indian lower crust with the transition from elastic (high‐viscosity) to low‐viscosity approximately under the Main Central Thrust is required to reproduce the observed postseismic uplift between China‐Nepal border and Peiku Lake, indicating the weakness of the Indian lower crust from the Lesser to High Himalaya. The afterslip simulated using a weak shear zone takes place in the adjacent area downdip of the rupture zone. The preferred model suggested that viscosities of the Tibetan lower crust, weakened Indian lower crust, and shear zone are 3 × 1018, 1019, and 4 × 1018 Pa s, respectively. The viscosity of the underthrusting Indian upper mantle was roughly estimated to be greater than 1021 Pa s. The model results imply that the near‐field deformation is dominated by both afterslip and viscoelastic relaxation of the weakened Indian lower crust, not only afterslip as suggested by previous studies.

Characterizing large-scale weak interlayer shear zones using conditional random field theory

The shear behavior of large-scale weak intercalation shear zones (WISZs) often governs the stability of foundations, rock slopes, and underground structures. However, due to their wide distribution, undulating morphology, complex fabrics, and varying degrees of contact states, characterizing the shear behavior of natural and complex large-scale WISZs precisely is challenging. This study proposes an analytical method to address this issue, based on geological fieldwork and relevant experimental results. The analytical method utilizes the random field theory and Kriging interpolation technique to simplify the spatial uncertainties of the structural and fabric features for WISZs into the spatial correlation and variability of their mechanical parameters. The Kriging conditional random field of the friction angle of WISZs is embedded in the discrete element software 3DEC, enabling activation analysis of WISZ C2 in the underground caverns of the Baihetan hydropower station. The results indicate that the activation scope of WISZ C2 induced by the excavation of underground caverns is approximately 0.5–1 times the main powerhouse span, showing local activation. Furthermore, the overall safety factor of WISZ C2 follows a normal distribution with an average value of 3.697.

Mantle exhumation at magma-poor rifted margins controlled by frictional shear zones

The transition zone from continental crust to the mature mid-ocean ridge spreading center of the Iberia-Newfoundland magma-poor rifted margins is mostly composed of exhumed mantle characterized by highs and domes with varying elevation, spacing and shape. The mechanism controlling strain localization and fault migration explaining the geometry of these peridotite ridges is poorly understood. Here we show using forward geodynamic models that multiple out-of-sequence detachments with recurring dip reversal form during magma-poor rifting and mantle exhumation as a consequence of the strength competition between weak frictional-plastic shear zones and the thermally weakened necking domain beneath the exhuming footwall explaining geometry of these peridotite ridges. Model behaviour also shows that fault types and detachment styles vary with spreading rate and fault strength and confirm that these results can be compared to other magma poor passive margins such as along Antarctica-Australia and to ultra-slow mid-ocean spreading systems as the South-West Indian Ridge.

Technological innovations in construction of underground caverns in basaltic rocks at Baihetan Hydropower Station on Jinsha River

The two underground powerhouses on the left and right banks of Baihetan Hydropower Station (BHT-HS) are symmetrically arranged inside basaltic rocks in the alpine gorge reach of the lower Jinsha River. Each of them is equipped with 8 sets of 1,000 MW hydraulic turbine generators, with the largest unit capacity, total installed capacity, span of main powerhouse (i.e. 34 m), diameter of dome-shaped tail-regulating chamber (max diameter of 48 m) and scale of underground caverns around the world. Due to complex layout, interaction between caverns and high hollowing rate, the construction process of the underground cavern group is faced many challenges. For example, the cavern excavation is located in the high ground stress area dominated by horizontal tectonic stress; the underlying basalt fractures are well developed, and the rock is hard and brittle, with low crack initiation strength; columnar jointed basalts develop locally; a number of large and weak gently-sloped shear zones develop diagonally cutting the main powerhouse. Following the construction ideology of “understanding, utilization, protection, monitoring, and feedback analysis on the surrounding rock masses during the construction period”, and with the engineering construction procedure of “the excavation, analysis, prediction and inspection of each rock layer”, this paper proposes a set of innovation construction technologies including “advanced pre-control, thin vertical layer, small planar partition, short excavation progress, fine blasting, strong and quick support, short interval monitoring, quick feedback, precise dynamic optimization” to control the internal cracking of the hard and brittle basaltic rocks, discontinuous deformation of staggered zones and spatio-temporal deformation of underground caverns in high geostress. It has promoted the safe and orderly construction of BHT-HS, the technical progress of the construction of hydropower station projects in China, and may also be applicable to the design and construction of similar projects worldwide.

Modelovanje zatvaranja uskih okeanskih basena numeričkim simulacijama

Numerical simulations represent an established approach in investigation of geodynamic processes in the Earth's interior and it is applied to a vast variety of geological problems. The closure of a narrow oceanic realm is modelled in this paper. Firstly, a simplified numerical model of a narrow ocean is established. The thermomechanical system is then simulated by solving Navier-Stokes and temperature equations. By imposing boundary conditions for velocity the convergent regime is simulated which leads to a subduction process along a predefined weak shear zone mantle. The subduction subsequently results in development of magmatism, the final closure of the oceanic realm and the onset of continental collision. Considerable amount of ophiolites were successfully obducted onto the overriding continent during the course of the simulation.

Viscoelastic relaxation of the upper mantle and afterslip following the 2014 MW8.1 Iquique earthquake

An improved understanding of postseismic crustal deformation following large subduction earthquakes may help to better understand the rheological properties of upper mantle and the slip behavior of subduction interface. Here we construct a three-dimensional viscoelastic finite element model to study the postseismic deformation of the 2014 MW8.1 Iquique, Chile earthquake. Elastic units in the model include the subducting slab, continental and oceanic lithospheres. Rheological units include the mantle wedge, the oceanic asthenosphere and upper mantle. We use a 2 ​km thick weak shear zone attached to the subduction fault to simulate the time-dependent stress-driven afterslip. The viscoelastic relaxation in the rheological units is represented by the Burgers rheology. We carry out grid-searches on the shear zone viscosity, thickness and viscosity of the asthenosphere, and they are determined to be 1017 ​Pa ​s, 110 ​km and 2 ​× ​1018 ​Pa ​s, respectively. The stress-driven afterlsip within the first two years is up to ~47 ​cm and becomes negligible after two years (no more than 5 ​cm/yr). Our results suggest that a thin, low-viscosity oceanic asthenosphere together with a weak shear zone attached to the fault are required to better reproduce the observed postseismic deformation.

Decadal Viscoelastic Postseismic Deformation of the 1964 Mw9.2 Alaska Earthquake

AbstractViscoelastic postseismic deformation after the 1964 Mw9.2 Alaska earthquake extends thousands of kilometers from the rupture region and lasts for decades, providing unique opportunities to better understand the three‐dimensional rheological properties of the Alaska subduction zone. We have optimized a three‐dimensional viscoelastic finite element model to study processes that control the postseismic deformation of the 1964 event. The model includes an elastic continental plate and an elastic oceanic plate, a two‐layered viscoelastic oceanic upper mantle, and a uniform viscoelastic mantle wedge. Stress‐driven afterslip is simulated by a thin weak shear zone. The viscoelastic relaxation of the upper mantle and shear zone is represented by the bi‐viscous Burgers rheology. The model has determined the viscosities of the mantle wedge and shear zone to be 3 × 1019 Pa s and 8 × 1016 Pa s, respectively. The afterslip takes place mostly within the first 5 years after the earthquake and is up to 4 m, equivalent to a modeled earthquake of Mw8.5. Model results reveal a spatial and temporal correlation between the afterslip distribution and later slow slip events. The model predicts that the surface deformation about 200 years after the earthquake will be controlled mostly by the relocking of the fault. Further tests on the impact of lateral variation in the mantle wedge viscosity indicate that the viscosity in the continental upper mantle thousands of kilometers from the rupture area may be about an order of magnitude higher than that of the mantle wedge in the subduction zone.

On the dynamics of coupled grain size evolution and shear heating in lithospheric shear zones

The formation of weak, localized shear zones in the lithosphere via tectonic forces is a requirement for the operation of plate tectonics on Earth. However, how shear zones form is not well understood. Shear heating and grain size reduction are two commonly cited mechanisms, and both processes are likely to be active in lithospheric shear zones. However, shear heating and grain size reduction can potentially act against each other, as elevated temperatures increase grain growth rates, and fast grain growth can eliminate weak shear zones formed via grain size reduction. Here, an effectively one-dimensional simple shear model that considers shear heating and grain size evolution is used to determine how grain size reduction and shear heating interact. The results indicate that for constant velocity boundary conditions at typical tectonic strain rates and mid- to lower-lithosphere temperatures, grain size reduction suppresses shear heating by weakening the shear zone. As a result shear heating at these conditions is small (∼10 K or less), unless grain size reduction via damage is significantly less effective than the models assume. For constant stress boundary conditions, grain size reduction enhances shear heating and lowers the stress threshold required to trigger thermal runaway instabilities, as previous studies have found. Finally, combining shear heating and grain size reduction allows weak shear zones to form over a wider range of conditions than would result from either mechanism acting in isolation; coupling between shear heating and grain size reduction does not impede weak shear zone formation.

Unexpected earthquake hazard revealed by Holocene rupture on the Kenchreai Fault (central Greece): Implications for weak sub-fault shear zones

High-resolution elevation models, palaeoseismic trenching, and Quaternary dating demonstrate that the Kenchreai Fault in the eastern Gulf of Corinth (Greece) has ruptured in the Holocene. Along with the adjacent Pisia and Heraion Faults (which ruptured in 1981), our results indicate the presence of closely-spaced and parallel normal faults that are simultaneously active, but at different rates. Such a configuration allows us to address one of the major questions in understanding the earthquake cycle, specifically what controls the distribution of interseismic strain accumulation? Our results imply that the interseismic loading and subsequent earthquakes on these faults are governed by weak shear zones in the underlying ductile crust. In addition, the identification of significant earthquake slip on a fault that does not dominate the late Quaternary geomorphology or vertical coastal motions in the region provides an important lesson in earthquake hazard assessment.

Collapse of passive margins by lithospheric damage and plunging grain size

The collapse of passive margins has been proposed as a possible mechanism for the spontaneous initiation of subduction. In order for a new trench to form at the junction between oceanic and continental plates, the cold and stiff oceanic lithosphere must be weakened sufficiently to deform at tectonic rates. Such rates are especially hard to attain in the cold ductile portion of the lithosphere, at which the mantle lithosphere reaches peak strength. The amount of weakening required for the lithosphere to deform in this tectonic setting is dictated by the available stress. Stress in a cooling passive margin increases with time (e.g., due to ridge push), and is augmented by stresses present in the lithosphere at the onset of rifting (e.g., due to drag from underlying mantle flow). Increasing stress has the potential to weaken the ductile portion of the lithosphere by dislocation creep, or by decreasing grain size in conjunction with a grain-size sensitive rheology like diffusion creep. While the increasing stress acts to weaken the lithosphere, the decreasing temperature acts to stiffen it, and the dominance of one effect or the other determines whether the margin might weaken and collapse. Here, we present a model of the thermal and mechanical evolution of a passive margin, wherein we predict formation of a weak shear zone that spans a significant depth-range of the ductile portion of the lithosphere. Stiffening due to cooling is offset by weakening due to grain size reduction, driven by the combination of imposed stresses and grain damage. Weakening via grain damage is modest when ridge push is the only source of stress in the lithosphere, making the collapse of a passive margin unlikely in this scenario. However, adding even a small stress-contribution from mantle drag results in damage and weakening of a significantly larger portion of the lithosphere. We posit that rapid grain size reduction in the ductile portion of the lithosphere can enable, or at least significantly facilitate, the collapse of a passive margin and initiate a new subduction zone. We use this model to estimate the conditions for passive margin collapse for modern and ancient Earth, as well as for Venus.

Friction of Shear-Fracture Zones.

A shear fracture of brittle solids under compression undergoes a substantial evolution from the initial microcracking to a fully formed powder-filled shear zone. Experiments covering the entire process are relatively easy to conduct, but they are very difficult to investigate in detail. Numerically, the large strain limit has remained a challenge. An efficient simulation model and a custom-made experimental device are employed to test to what extent a shear fracture alone is sufficient to drive material to spontaneous self-lubrication. A "weak shear zone" is an important concept in geology, and a large number of explanations, specific for tectonic conditions, have been proposed. We demonstrate here that weak shear zones are far more general, and that their emergence only demands that a microscopic, i.e., fragment-scale, stress relaxation mechanism develops during the fracture process.

Mineral fabrics in high-level intrusions recording crustal strain and volcano–tectonic interactions: the Shellenbarger pluton, Sierra Nevada, California

The Minarets caldera is a volcano–plutonic complex in the Sierra Nevada, California, that exemplifies complex interactions between volcanism and tectonic deformation in continental-margin arcs. Caldera evolution commenced with emplacement of pre-collapse rhyolitic ash-flow tuff, followed by collapse and deposition of volcanic breccia and rhyodacitic ash-flow tuff. Subsequently, the volcanic rocks were deformed along the regional Bench Canyon shear zone. The caldera centre was then intruded by the resurgent c. 100 Ma steep-sided Shellenbarger granite pluton, which steepened the shear zone foliation. The pluton was overprinted by syn- to post-magmatic ∼NNE–SSW horizontal shortening; the same shortening was documented in several other Late Cretaceous syntectonic plutons in the Sierra Nevada and interpreted to record dextral transpression during convergence of the Farallon and North American plates. To explain the unusual tectonic fabric in the shallow-level Shellenbarger pluton, we develop a general model for strain partitioning in syntectonic magma bodies emplaced at various crustal levels. We propose that shallow intrusions, isolated within stiff crust, may tend to accommodate minor pure shear strain whereas simple shear dominates along weak faults and shear zones. By contrast, a rheological reversal is crossed deeper in the crust and magma bodies become the weakest, simple shear-dominated parts of the system. Supplementary material: Analytical methods and anisotropy of magnetic susceptibility and U–Th–Pb isotopic data are available at https://doi.org/10.6084/m9.figshare.c.3582749

Reaction-induced rheological weakening enables oceanic plate subduction.

Earth is the only terrestrial planet in our solar system where an oceanic plate subducts beneath an overriding plate. Although the initiation of plate subduction requires extremely weak boundaries between strong plates, the way in which oceanic mantle rheologically weakens remains unknown. Here we show that shear-enhanced hydration reactions contribute to the generation and maintenance of weak mantle shear zones at mid-lithospheric depths. High-pressure friction experiments on peridotite gouge reveal that in the presence of hydrothermal water, increasing strain and reactions lead to an order-of-magnitude reduction in strength. The rate of deformation is controlled by pressure-solution-accommodated frictional sliding on weak hydrous phyllosilicate (talc), providing a mechanism for the ‘cutoff' of the high peak strength at the brittle-plastic transition. Our findings suggest that infiltration of seawater into transform faults with long lengths and low slip rates is an important controlling factor on the initiation of plate tectonics on terrestrial planets.

Stress‐driven relaxation of heterogeneous upper mantle and time‐dependent afterslip following the 2011 Tohoku earthquake

AbstractUnderstanding of postseismic deformation following great subduction zone earthquakes is complicated by the combined effects of viscoelastic relaxation of earthquake‐induced stresses in the upper mantle and time‐dependent afterslip on the megathrust. We integrate geodetic observations and constraints on afterslip from small repeating earthquakes on the megathrust to better distinguish contributions from these two postseismic processes. We have developed a three‐dimensional, spherical viscoelastic finite element model to study the postseismic deformation of the 2011 Mw9.0 Tohoku earthquake that has been recorded at unprecedented high resolution in space and time. We model stress‐driven afterslip in a 2 km thick weak shear zone away from historic rupture zones on the megathrust. We model both the viscoelastic relaxation of the upper mantle and shear zone deformation with a transient Burgers body rheology. The transient Kelvin viscosity is assumed to be one order of magnitude lower than that of the Maxwell viscosity. Viscoelastic relaxation of the mantle wedge alone causes postseismic uplift and seaward motion in the upper plate, opposite to the pattern from relaxation of just the oceanic upper mantle. Afterslip on the fault produces uplift updip of the afterslip zone and subsidence over its downdip edge and mostly seaward motion above the afterslip zone. The best fit Maxwell viscosity of the shear zone at depths ≤50 km is 1017 Pa s, constrained by afterslip estimates from repeating earthquakes. The optimal viscosities of the deep weak shear zone, continental mantle wedge, and oceanic upper mantle are determined to be 5 × 1017 Pa s, 3 × 1019 Pa s, and 5 × 1019 Pa s, respectively. The stress‐driven afterslip in the shear zone is up to ~3.5 m in the first 2 years after the earthquake, equivalent to an Mw8.4. Our model reproduces the first‐order pattern of the GPS observations both in horizontal and in vertical directions. Seafloor geodetic observations of subsidence and landward motions near the high‐slip zone of the earthquake provide evidence for a low‐viscosity asthenosphere below the oceanic lithosphere.

Macro- and microstructural analysis of the North Tea Lake Mylonite Zone: an extensional shear zone in the Central Gneiss Belt, Grenville Province, Ontario

The North Tea Lake Mylonite Zone is a late extensional ductile fault that is concordant with and has reworked fabrics of the North Tea Lake Shear Zone, the frontal thrust shear zone of the upper amphibolite–granulite facies Kiosk domain within the interior of the Central Gneiss Belt. North Tea Lake Mylonite Zone fabric is an anomalously fine-grained mylonite compared to Central Gneiss Belt gneisses, and consists of three microstructural domains that display progressive recrystallization and grain size refinement of the protolith granitoid. On the basis of petrography and electron backscatter diffraction, these microdomains are inferred to represent a transition from dominantly dislocation creep to diffusion creep and diffusion-accommodated grain boundary sliding at elevated stress (>100 MPa), low fluid activity, and temperatures ∼500 °C. The North Tea Lake Mylonite Zone is interpreted to mark a step in the progressive transition in deformation mode during late- to post-Ottawan extension and cooling of the Grenville orogen from weak, wide, wet, and warm shear zones to Rigolet-phase cooler, narrow, ultrafine, high-stress shear zones reworking dry protoliths.

Modes of continental extension in a crustal wedge

We ran numerical experiments of the extension of a crustal wedge as an approximation to extension in an orogenic belt or a continental margin. We study the effects of the strength of the lower crust and of a weak mid-crustal shear zone on the resulting extension styles. A weak mid-crustal shear zone effectively decouples upper crustal extension from lower crustal flow. Without the mid-crustal shear zone, the degree of coupling between the upper and the lower crust increases and extension of the whole crust tends to focus on the thickest part of the wedge. We identify three distinct modes of extension determined by the strength of the lower crust, which are characterized by 1) localized, asymmetric crustal exhumation in a single massif when the lower crust is weak, 2) the formation of rolling-hinge normal faults and the exhumation of lower crust in multiple core complexes with an intermediate strength lower crust, and 3) distributed domino faulting over the weak mid-crustal shear zone when the lower crust is strong. A frictionally stronger mid-crustal shear zone does not change the overall model behaviors but extension occurred over multiple rolling-hinges. The 3 modes of extension share characteristics similar to geological models proposed to explain the formation of metamorphic core complexes: 1) the crustal flow model for the weak lower crust, 2) the rolling-hinge and crustal flow models when the lower crust is intermediate and 3) the flexural uplift model when the lower crust is strong. Finally we show that the intensity of decoupling between the far field extension and lower crustal flow driven by the regional pressure gradient in the wedge control the overall style of extension in the models.

Future accreted terranes: a compilation of island arcs, oceanic plateaus, submarine ridges, seamounts, and continental fragments

Abstract. Allochthonous accreted terranes are exotic geologic units that originated from anomalous crustal regions on a subducting oceanic plate and were transferred to the overriding plate by accretionary processes during subduction. The geographical regions that eventually become accreted allochthonous terranes include island arcs, oceanic plateaus, submarine ridges, seamounts, continental fragments, and microcontinents. These future allochthonous terranes (FATs) contribute to continental crustal growth, subduction dynamics, and crustal recycling in the mantle. We present a review of modern FATs and their accreted counterparts based on available geological, seismic, and gravity studies and discuss their crustal structure, geological origin, and bulk crustal density. Island arcs have an average crustal thickness of 26 km, average bulk crustal density of 2.79 g cm−3, and three distinct crustal units overlying a crust–mantle transition zone. Oceanic plateaus and submarine ridges have an average crustal thickness of 21 km and average bulk crustal density of 2.84 g cm−3. Continental fragments presently on the ocean floor have an average crustal thickness of 25 km and bulk crustal density of 2.81 g cm−3. Accreted allochthonous terranes can be compared to these crustal compilations to better understand which units of crust are accreted or subducted. In general, most accreted terranes are thin crustal units sheared off of FATs and added onto the accretionary prism, with thicknesses on the order of hundreds of meters to a few kilometers. However, many island arcs, oceanic plateaus, and submarine ridges were sheared off in the subduction interface and underplated onto the overlying continent. Other times we find evidence of terrane–continent collision leaving behind accreted terranes 25–40 km thick. We posit that rheologically weak crustal layers or shear zones that were formed when the FATs were produced can be activated as detachments during subduction, allowing parts of the FAT crust to accrete and others to subduct. In many modern FATs on the ocean floor, a sub-crustal layer of high seismic velocities, interpreted as ultramafic material, could serve as a detachment or delaminate during subduction.

Assessment of rock slope stability and structurally controlled failures along Samma escarpment road, Asir Region (Saudi Arabia)

Samma escarpment road is located in Asir Region, Saudi Arabia. It is located NW of Abha city. This escarpment road represents a major corridor in the area which connects different cities and touristic resorts in the region. It is descended from Sudah plateau at about 2700 m above sea level (asl) toward Wadi al Aws at about 1500 m asl. The total length of the road section is about 8 km which is passing through a highly mountainous area characterized by a complex geological and structural elements. This road has been exposed to frequent slope failures from time to time due to various factors such as intense rain storms, different geological and structural elements (weak rocks, shear zones, and faults), road characteristics (different horizontal/vertical curvatures and narrow road section), and human activities (uncontrolled rock cuts). Many sliding events have been documented along this escarpment road, particularly during and following rainstorms. The purpose of this paper is to evaluate the geology and structures by applying the rock mass rating (RMR) and slope stability (structurally controlled failures) along the Samma escarpment road. The stability analysis was performed using two Dips and RockPack III programs with the help of RocLab software. The Summa escarpment rock cuts were classified into 51 stations which were investigated in detail. Results indicated that most rock stations are poor quality, and 31, 32, and 41 stations are stable, 4, 6, and 7 stations are marginally stable, and 16, 13, and 3 stations are potentially unstable due to planar, wedge, and toppling failures, respectively. Finally, different recommendations and remediation methods were suggested as mitigation measures.