Progressive weakening within the overriding plate during dual inward dipping subduction

Dual inward dipping subduction often produces complex deformation patterns in the overriding plate. However, the geodynamic process of how dual inward dipping subduction relates to this deformation is still poorly understood. Here we apply a composite viscosity, dependent on multiple parameters, e.g., temperature, pressure, strain rate etc., in 2-D thermo-mechanical numerical modelling to investigate how dual inward dipping subduction modifies the rheological structure of the overriding plate. Three variables are investigated to understand what controls the maximum degree of weakening. We find that the initial length and thickness of the overriding plate are negatively correlated with the magnitude of viscosity reduction. While the initial thickness of the subducting plate positively relates to the magnitude of viscosity reduction. The progressive weakening can result in a variety of stretching states ranging from 1) little or no lithosphere thinning and extension, to 2) limited thermal lithosphere thinning, and 3) localised rifting followed by spreading extension. Compared with single sided subduction, dual inward dipping subduction further reduces the magnitude of viscosity of the overriding plate. It does this by creating a dynamic fixed boundary condition for the overriding plate and forming a stronger upwelling mantle flow underlying the overriding plate. Three types of feedback weakening cycles are recognised, among which the strain rate weakening mechanism plays the dominant role in lowering the viscosity of the overriding plate throughout the simulation. Strain rate weakening is also a precondition for initiating thermal weakening, strain localisation and lithosphere thinning.

High 3He/4He in central Panama reveals a distal connection to the Galápagos plume

It is well established that mantle plumes are the main conduits for upwelling geochemically enriched material from Earth's deep interior. The fashion and extent to which lateral flow processes at shallow depths may disperse enriched mantle material far (>1,000 km) from vertical plume conduits, however, remain poorly constrained. Here, we report He and C isotope data from 65 hydrothermal fluids from the southern Central America Margin (CAM) which reveal strikingly high 3He/4He (up to 8.9RA) in low-temperature (≤50 °C) geothermal springs of central Panama that are not associated with active volcanism. Following radiogenic correction, these data imply a mantle source 3He/4He >10.3RA (and potentially up to 26RA, similar to Galápagos hotspot lavas) markedly greater than the upper mantle range (8 ± 1RA). Lava geochemistry (Pb isotopes, Nb/U, and Ce/Pb) and geophysical constraints show that high 3He/4He values in central Panama are likely derived from the infiltration of a Galápagos plume–like mantle through a slab window that opened ∼8 Mya. Two potential transport mechanisms can explain the connection between the Galápagos plume and the slab window: 1) sublithospheric transport of Galápagos plume material channeled by lithosphere thinning along the Panama Fracture Zone or 2) active upwelling of Galápagos plume material blown by a “mantle wind” toward the CAM. We present a model of global mantle flow that supports the second mechanism, whereby most of the eastward transport of Galápagos plume material occurs in the shallow asthenosphere. These findings underscore the potential for lateral mantle flow to transport mantle geochemical heterogeneities thousands of kilometers away from plume conduits.

Geochronology and Isotope Geochemistry of the Yingfang Pb-Zn-Ag Deposit: Implications for Large-Scale Metallogeny along the Northern Flank of the North China Craton

The northern flank of the North China Craton (NCC) hosts a linear zone of gold, molybdenum, silver, lead, and zinc polymetallic ore deposits. Among these, the Yingfang Pb-Zn-Ag deposit is located in the central part of the Yanshan–Liaoxi metallogenic belt (YLMB) which extends for approximately 1000 km and forms part of the major mineralized zone. In this study, we characterize the mineralization and trace the ore genesis based on new sulfur and lead isotopic geochemistry and evaluate the timing of mineralization from Rb-Sr isotope dating of sulfides. The pyrite δ34S values range from +3.2‰ to +5.8‰ with a mean at +4.07‰, close to the values of mantle and meteorite sulfur. The 206Pb/204Pb values range from 16.833 to 18.956, 207Pb/204Pb from 15.374 to 15.522, and 208Pb/204Pb from 37.448 to 37.928. Five samples of sulfide, from the Yingfang deposit, yield a Rb-Sr isochron age of 135.7 ± 4.1 Ma. This age is close to the age of the adjacent Niujuan Ag-Au deposit and the associated Er’daogou granite, suggesting a close relationship between magmatism and metallogeny in this region. The S and Pb isotopes of the regional silver polymetallic deposits show similar sources of ore-forming materials. According to a compilation of the available age data on the Mesozoic ore deposits in the northern flank of the NCC, we divide the mineralization into the following four periods: 240–205 Ma, 190–160 Ma, 155–135 Ma, and 135–100 Ma. Mesozoic magmatism and mineralization in the Yingfang deposit mainly took place at 245 Ma and 145–135 Ma. We correlate the Pb-Zn-Ag mineralization to metallogeny associated with large-scale inhomogeneous lithosphere thinning beneath the NCC.

Heat Flow Data in an Area of the Eastern Southern Basin and Range in Arizona Contribute to an Analysis of Neogene Lithosphere Thinning Greater than 100 Km

Heat flow data and thermochronologic derived paleotemperature gradient data are examined to calculate heat flow ~25 Ma and, at present, for a southern Basin and Range location north of Tucson, Arizona. An increase in the surface heat flow is estimated from ~25 Ma to the present; changing from ~47 to ~83 mW m-2. Steady-state conduction temperature vs. depth profiles provide estimates of lithosphere thicknesses both for the present and for ~25 Ma. Different heat transfer models for present heat flow predict present LAB depth that agrees with seismic studies. From these temperature profiles, lithosphere thinning from ~184 km to ~70 km is suggested during the Neogene. Mantle lithosphere thinning caused by thermal phenomena is likely a fundamental driving force for southern Basin and Range extension. Because the mantle lithosphere has likely thinned much more than the crust, it is shown that additional vertical advection, such as an asthenosphere plume, delaminating part of the mantle lithosphere, convection cells, and rising magmas along conduits, add to the vertical advection component of upper mantle lithosphere extension. Interestingly, values of heat flow 25 Ma, lithosphere thicknesses 25 Ma, and Neogene lithosphere thinning are somewhat similar for the Four Corners area of the Colorado Plateau and the southern Basin and Range, even though Neogene tectonic development was quite different, i.e., no Neogene extension in the Colorado Plateau vs. ~57% in the southern Basin and Range. Neogene lithosphere thinning phenomena are likely different in the two regions.

Clinopyroxene and Garnet Mantle Cargo in Kimberlites as Probes of Dharwar Craton Architecture and Geotherms, with Implications for Post-1·1 Ga Lithosphere Thinning Events Beneath Southern India

The Wajrakarur Kimberlite Field (WKF) on the Eastern Dharwar Craton in southern India hosts several occurrences of Mesoproterozoic kimberlites, lamproites and ultramafic lamprophyres, for which mantle-derived xenoliths are rare and only poorly preserved. The general paucity of mantle cargo has hampered the investigation of the nature and evolution of the continental lithospheric mantle (CLM) beneath cratonic southern India. We present a comprehensive study of the major and trace element compositions of clinopyroxene and garnet xenocrysts recovered from heavy mineral concentrates for three c.1·1 Ga old WKF kimberlite pipes (P7, P9, P10), with the goal to improve our understanding of the cratonic mantle architecture and its evolution beneath southern India. The pressure-temperature conditions recorded by peridotitic clinopyroxene xenocrysts, estimated using single-pyroxene thermobarometry, suggest a relatively moderate cratonic mantle geotherm of 40 mW/m2 at 1·1 Ga. Reconstruction of the vertical distribution of clinopyroxene and garnet xenocrysts, combined with some rare mantle xenoliths data, reveals a compositionally layered CLM structure. Two main lithological horizons are identified and denoted as layer A (∼80–145 km depth) and layer B (∼160–190 km depth). Layer A is dominated by depleted lherzolite with subordinate amounts of pyroxenite, whereas layer B comprises mainly refertilised and Ti-metasomatized peridotite. Harzburgite occurs as a minor lithology in both layers. Eclogite stringers occur within the lower portion of layer A and at the bottom of layer B near the lithosphere–asthenosphere boundary at 1·1 Ga. Refertilisation of layer B is marked by garnet compositions with enrichment in Ca, Ti, Fe, Zr and LREE, although Y is depleted compared to garnet in layer A. Garnet trace element systematics such as Zr/Hf and Ti/Eu indicate that both kimberlitic and carbonatitic melts have interacted with and compositionally overprinted layer B. Progressive changes in the REE systematics of garnet grains with depth record an upward percolation of a continuously evolving metasomatic agent. The intervening zone between layers A and B at ∼145–160 km depth is characterized by a general paucity of garnet. This ‘garnet-paucity’ zone and an overlying type II clinopyroxene-bearing zone (∼115–145 km) appear to be rich in hydrous mineral assemblages of the MARID- or PIC kind. The composite horizon between ∼115–160 km depth may represent the product of intensive melt/rock interaction by which former garnet was largely reacted out and new metasomatic phases such as type II clinopyroxene and phlogopite plus amphibole were introduced. By analogy with better-studied cratons, this ‘metasomatic horizon’ may be a petrological manifestation of a former mid-lithospheric discontinuity at 1·1 Ga. Importantly, the depth interval of the present-day lithosphere–asthenosphere boundary beneath Peninsular India as detected in seismic surveys coincides with this heavily overprinted metasomatic horizon, which suggests that post-1·1 Ga delamination of cratonic mantle lithosphere progressed all the way to mid-lithospheric depth. This finding implies that strongly overprinted metasomatic layers, such as the ‘garnet-paucity’ zone beneath the Dharwar craton, present structural zones of weakness that aid lithosphere detachment and foundering in response to plate tectonic stresses.

Plagioclase peridotite or olivine-plagioclase assemblage in orogenic peridotites: its implications on high-temperature decompression of the subcontinental lithosphere-asthenosphere boundary zone

<p>Orogenic peridotites are expected to provide direct information with high spatial resolution for a better understanding of the processes taking place in the lithosphere and asthenosphere boundary zones (LABZ), where the transfer mechanisms of heat, material, and momentum from the Earth’s interior to the surface drastically change. Plagioclase peridotite or olivine-plagioclase assemblage <em>sensu lato</em> has been reported from some orogenic peridotites. The olivine-plagioclase assemblage in fertile systems is in principle not stable even at the depth of the upper most subcontinental lithospheric mantle (SCLM) because (1) the common crustal thickness in normal non-cratonic SCLM is ~35km, (2) the Moho temperature for the mean steady-state continental geotherm is much lower than 600°C, (3) the upper stability limit of plagioclase (plagioclase to spinel facies transition) becomes shallower with decrease in temperature, and (4) kinetic barrier for subsolidus reactions in the peridotite system becomes enormous at temperatures below 600°C. The occurrence of olivine-plagioclase assemblage in some orogenic peridotite bodies, therefore, implies transient and dynamic high-temperature (>800°C) processing at depth shallower than 20km (plagioclase-spinel facies boundary at ~800°C), i.e., high-temperature decompression of LABZ up to the depth closer to the Moho. Adiabatic decompression of high-temperature LABZ leading to decompressional melting with inefficient melt segregation may give rise to plagioclase peridotite. Decompression along moderately high temperature adiabatic path or heating to allow subsolidus reactions leading to transformation of either spinel peridotites or garnet peridotites may give rise to plagioclase peridotite. However, decompression of LABZ associated with efficient cooling does not produce any olivine-plagioclase assemblage. Plagioclase peridotites thus could provide precious information on the dynamics of shallowing LABZ and underlying asthenosphere.</p><p>We have examined several orogenic peridotite complexes, Ronda, Pyrenees, Lanzo, and Horoman, to clarify the extent of shallow thermal processing based on olivine-plagioclase assemblage. The key approach of this study is searching olivine-plagioclase assemblage not only in various lithologies but also in microstructures, whose scale and mode of occurrence provide extent and strength of thermal processing in the shallow upper mantle. The wide-spread occurrence of plagioclase peridotites and localized partial melting in Lanzo suggest exhumation along high temperature adiabatic paths from the thermally structured <span>LABZ in the </span>Seiland subfacies; the predominance of plagioclase peridotites and its localized partial melting in Horoman <span>suggest </span> exhumation along variously heated paths from the garnet stability field; the moderate development of plagioclase peridotites without partial melting in Ronda suggest exhumation along variously but weekly heated paths from the spinel-garnet stability field, and the occurrence of minor plagioclase peridotites in Pyrenees suggests exhumation along cold path from the garnet-spinel facies boundaries. We propose that the extent of shallower thermal processing decreases, and thus lithosphere thinning becomes less extensive in this order.</p>