Fluid Infiltration Through Oceanic Lower Crust in Response to Reaction‐Induced Fracturing: Insights From Serpentinized Troctolite and Numerical Models

2020 ◽  
Vol 125 (11) ◽  
Author(s):  
Kazuki Yoshida ◽  
Atsushi Okamoto ◽  
Hiroyuki Shimizu ◽  
Ryosuke Oyanagi ◽  
Noriyoshi Tsuchiya ◽  
...  
2021 ◽  
Vol 9 ◽  
Author(s):  
Lin Chen

The continental lower crust is an important composition- and strength-jump layer in the lithosphere. Laboratory studies show its strength varies greatly due to a wide variety of composition. How the lower crust rheology influences the collisional orogeny remains poorly understood. Here I investigate the role of the lower crust rheology in the evolution of an orogen subject to horizontal shortening using 2D numerical models. A range of lower crustal flow laws from laboratory studies are tested to examine their effects on the styles of the accommodation of convergence. Three distinct styles are observed: 1) downwelling and subsequent delamination of orogen lithosphere mantle as a coherent slab; 2) localized thickening of orogen lithosphere; and 3) underthrusting of peripheral strong lithospheres below the orogen. Delamination occurs only if the orogen lower crust rheology is represented by the weak end-member of flow laws. The delamination is followed by partial melting of the lower crust and punctuated surface uplift confined to the orogen central region. For a moderately or extremely strong orogen lower crust, topography highs only develop on both sides of the orogen. In the Tibetan plateau, the crust has been doubly thickened but the underlying mantle lithosphere is highly heterogeneous. I suggest that the subvertical high-velocity mantle structures, as observed in southern and western Tibet, may exemplify localized delamination of the mantle lithosphere due to rheological weakening of the Tibetan lower crust.


2020 ◽  
Author(s):  
Ryan Stoner ◽  
Mark Behn ◽  
Bradley Hacker

<p>Geochronological and thermobarometric data from a lower crustal xenolith suite in the Pamir offer a unique record of the transport of lower crust to mantle depths after an episode of slab breakoff. We compare petrologically constrained pressure-temperature-time paths from the xenoliths to pressure-temperature-time (P-T-t) paths of tracked markers in 2-D numerical geodynamic models of density foundering with thermodynamically calculated densities. We investigate whether gravitational “drip” instabilities or the peeling back of a dense layer—delamination—can reproduce the P-T-t paths seen in the xenoliths, with the ancillary goal of capturing the positive feedback between mechanical thickening and densification of the lower crust. Key thermobarometric observations from the xenoliths we try to match in our numerical study are: (1) initial heating at near-constant pressure followed by (2) a sharp increase in pressure with continued heating. We find that thick crustal sections develop P-T-t paths in numerical models of delamination that match the observations from xenoliths: the lower crust initially heats due to return flow from upwelling asthenosphere, and then foundering mantle lithosphere and crust show a marked increase in pressure with additional heating. Initial gravitational drip instabilities founder with relatively little heating yet may thin the mantle lithosphere sufficiently to allow for subsequent delamination or asymmetric drips to nucleate in the region of hotter, thinner mantle lithosphere. Such subsequent asymmetric drips or delamination entrain crust that closely follows the P-T-t path from xenoliths. This suggests that the xenoliths were not derived from an initial drip instability, but instead from later instabilities or delamination enabled by thinning of the lithosphere. In all models where density foundering occurs, the positive feedback between contraction and densification of the lower crust leads to the loss of initially positively buoyant lower crust. The combination of geological and numerical methods constrains the geometry and triggers of lower crustal foundering during collision. Contraction alone does not match the record of foundering; the lithosphere must have also been asymmetrically thinned.</p><p> </p>


2020 ◽  
Author(s):  
Pauline Chenin ◽  
Gianreto Manatschal ◽  
Stefan M. Schmalholz ◽  
Thibault Duretz

<p>Although so-called "magma-poor" rifted margins display a large variability on a local scale, they are characterized by a number of common primary features worldwide such as their first-order architecture (proximal, necking, hyperextended, exhumation and oceanic domains), their lithological evolution along dip and the deformation processes associated with their different rifting stages. In this contribution, we first emphasize the primary morphological and lithological architecture of magma-poor rifted margins and how they relate to specific deformation modes (pure shear thinning, mechanical necking, frictional extensional wedge, detachment faulting and seafloor spreading). Second, we focus on the necking stage of rifting, which corresponds to the first major thinning event (when the crust is thinned from its initial thickness to ~ 10 km). We display the range of possible topographic and thermal evolutions of "magma-poor" and "sedimentary starved" rift systems depending on their lithosphere rheology. Our two-dimensional thermo-mechanical numerical models suggest that extension of lithospheres where the crust and the mantle are mechanically decoupled by a weak lower crust results in a complex morphotectonic evolution of rift systems, with formation of temporary restricted sub-basins framed by uplifted parts of the future distal margin. Mechanical decoupling between the crust and the mantle controls also largely the thermal evolution of rift systems during the necking phase since for equivalent extension rates and initial geotherms: (i) weak/decoupled lithospheres have a higher geothermal gradient at the end of the necking phase than strong/coupled lithospheres; and (ii) weak/decoupled lithospheres show intense heating of the lower crust at the rift center and intense cooling of the crust on either side of the rift center, unlike strong/coupled lithospheres. These behaviors contrast with the continuous subsidence and cooling predicted by the commonly used depth-uniform thinning model.</p>


2021 ◽  
Author(s):  
L B Harris ◽  
P Adiban ◽  
E Gloaguen

Aeromagnetic and ground gravity data for the Canadian Superior Province, filtered to extract long wavelength components and converted to pseudo-gravity, highlight deep, N-S trending regional-scale, rectilinear faults and margins to discrete, competent mafic or felsic granulite blocks (i.e. at high angles to most regional mapped structures and sub-province boundaries) with little to no surface expression that are spatially associated with lode ('orogenic') Au and Ni-Cu-PGE-Cr occurrences. Statistical and machine learning analysis of the Red Lake-Stormy Lake region in the W Superior Province confirms visual inspection for a greater correlation between Au deposits and these deep N-S structures than with mapped surface to upper crustal, generally E-W trending, faults and shear zones. Porphyry Au, Ni, Mo and U-Th showings are also located above these deep transverse faults. Several well defined concentric circular to elliptical structures identified in the Oxford Stull and Island Lake domains along the S boundary of the N Superior proto-craton, intersected by N- to NNW striking extensional fractures and/or faults that transect the W Superior Province, again with little to no direct surface or upper crustal expression, are spatially associated with magmatic Ni-Cu-PGE-Cr and related mineralization and Au occurrences. The McFaulds Lake greenstone belt, aka. 'Ring of Fire', constitutes only a small, crescent-shaped belt within one of these concentric features above which 2736-2733 Ma mafic-ultramafic intrusions bodies were intruded. The Big Trout Lake igneous complex that hosts Cr-Pt-Pd-Rh mineralization west of the Ring of Fire lies within a smaller concentrically ringed feature at depth and, near the Ontario-Manitoba border, the Lingman Lake Au deposit, numerous Au occurrences and minor Ni showings, are similarly located on concentric structures. Preliminary magnetotelluric (MT) interpretations suggest that these concentric structures appear to also have an expression in the subcontinental lithospheric mantle (SCLM) and that lithospheric mantle resistivity features trend N-S as well as E-W. With diameters between ca. 90 km to 185 km, elliptical structures are similar in size and internal geometry to coronae on Venus which geomorphological, radar, and gravity interpretations suggest formed above mantle upwellings. Emplacement of mafic-ultramafic bodies hosting Ni-Cr-PGE mineralization along these ringlike structures at their intersection with coeval deep transverse, ca. N-S faults (viz. phi structures), along with their location along the margin to the N Superior proto-craton, are consistent with secondary mantle upwellings portrayed in numerical models of a mantle plume beneath a craton with a deep lithospheric keel within a regional N-S compressional regime. Early, regional ca. N-S faults in the W Superior were reactivated as dilatational antithetic (secondary Riedel/R') sinistral shears during dextral transpression and as extensional fractures and/or normal faults during N-S shortening. The Kapuskasing structural zone or uplift likely represents Proterozoic reactivation of a similar deep transverse structure. Preservation of discrete faults in the deep crust beneath zones of distributed Neoarchean dextral transcurrent to transpressional shear zones in the present-day upper crust suggests a 'millefeuille' lithospheric strength profile, with competent SCLM, mid- to deep, and upper crustal layers. Mechanically strong deep crustal felsic and mafic granulite layers are attributed to dehydration and melt extraction. Intra-crustal decoupling along a ductile décollement in the W Superior led to the preservation of early-formed deep structures that acted as conduits for magma transport into the overlying crust and focussed hydrothermal fluid flow during regional deformation. Increase in the thickness of semi-brittle layers in the lower crust during regional metamorphism would result in an increase in fracturing and faulting in the lower crust, facilitating hydrothermal and carbonic fluid flow in pathways linking SCLM to the upper crust, a factor explaining the late timing for most orogenic Au. Results provide an important new dataset for regional prospectively mapping, especially with machine learning, and exploration targeting for Au and Ni-Cr-Cu-PGE mineralization. Results also furnish evidence for parautochthonous development of the S Superior Province during plume-related rifting and cannot be explained by conventional subduction and arc-accretion models.


Geology ◽  
2020 ◽  
Vol 48 (7) ◽  
pp. 673-677 ◽  
Author(s):  
Benjamin Malvoisin ◽  
Håkon Austrheim ◽  
György Hetényi ◽  
Julien Reynes ◽  
Jörg Hermann ◽  
...  

Abstract The densification of the lower crust in collision and subduction zones plays a key role in shaping the Earth by modifying the buoyancy forces acting at convergent boundaries. It takes place through mineralogical reactions, which are kinetically favored by the presence of fluids. Earthquakes may generate faults serving as fluid pathways, but the influence of reactions on the generation of seismicity at depth is still poorly constrained. Here we present new petrological data and numerical models to show that in the presence of fluids, densification reactions can occur very fast, on the order of weeks, and consume fluids injected during an earthquake, which leads to porosity formation and fluid pressure drop by several hundreds of megapascals. This generates a mechanically highly unstable system subject to collapse and further seismic-wave emission during aftershocks. This mechanism creates new pathways for subsequently arriving fluids, and thus provides a route for self-sustained densification of the lower crust.


2020 ◽  
Vol 82 ◽  
pp. 149-160
Author(s):  
N Kargapolova

Numerical models of the heat index time series and spatio-temporal fields can be used for a variety of purposes, from the study of the dynamics of heat waves to projections of the influence of future climate on humans. To conduct these studies one must have efficient numerical models that successfully reproduce key features of the real weather processes. In this study, 2 numerical stochastic models of the spatio-temporal non-Gaussian field of the average daily heat index (ADHI) are considered. The field is simulated on an irregular grid determined by the location of weather stations. The first model is based on the method of the inverse distribution function. The second model is constructed using the normalization method. Real data collected at weather stations located in southern Russia are used to both determine the input parameters and to verify the proposed models. It is shown that the first model reproduces the properties of the real field of the ADHI more precisely compared to the second one, but the numerical implementation of the first model is significantly more time consuming. In the future, it is intended to transform the models presented to a numerical model of the conditional spatio-temporal field of the ADHI defined on a dense spatio-temporal grid and to use the model constructed for the stochastic forecasting of the heat index.


2003 ◽  
Vol 59 (3-4) ◽  
pp. 10
Author(s):  
D. Yu. Kulik ◽  
S. L. Senkevich ◽  
Victor Ivanovich Tkachenko
Keyword(s):  

Sign in / Sign up

Export Citation Format

Share Document