Origin of the Yinshan epithermal-porphyry Cu–Au–Pb–Zn–Ag deposit,southeastern China: insights from geochemistry,Sr–Nd and zircon U–Pb–Hf–O isotopes

The Yinshan deposit is a large epithermal-porphyry polymetallic deposit, and the timing and petrogenesis of ore-hosting porphyries have been hotly debated. We present new results from geochemical, whole-rock Sr–Nd and zircon U–Pb–Hf–O isotopic investigations. Zircon U–Pb data demonstrate that the quartz porphyry, dacitic porphyry, and quartz dioritic porphyry formed at ?172.2 ± 0.4 Ma, ?171.7 ± 0.5 Ma, and ?170.9 ± 0.3 Ma, respectively. Inherited zircon cores show significant age spreads from ?730 to ?1390 Ma. Geochemically, they are high-K calc-alkaline or shoshonitic rocks with arc-like trace element patterns. They have similar whole-rock Nd and zircon Hf isotopic compositions, yet an increasing trend in ?_Nd(t) and ?_Hf(t) values typifies the suite. Older (inherited) zircons of the three porphyries display Hf compositions comparable to those of the Jiangnan Orogen basement rocks. In situ zircon oxygen isotopic analyses reveal that they have similar oxygen isotopic compositions, which are close to those of mantle zircons. Moreover, a decreasing trend of δ¹⁸O values is present. We propose that the ore-related porphyries of the Yinshan deposit were emplaced contemporaneously and derived from partial melting of Neoproterozoic arc-derived mafic (or ultra-mafic) rocks. Modelling suggests that the quartz porphyries, dacitic porphyries, and quartz dioritic porphyries experienced ?25%, ?10%, and ?10% crustal contaminations by Shuangqiaoshan rocks. Our study provides important constraints on mantle–crust interaction in the genesis of polymetallic mineralization associated with Mesozoic magmatism in southeastern China.     

Deciphering the geochronology of a large granitoid pluton (Karkonosze Granite,SW Poland): an assessment of U–Pb zircon SIMS and Rb–Sr whole-rock dates relative to U–Pb zircon CA-ID-TIMS

Granitoid plutons are often difficult to radiometrically date precisely due to the possible effects of protracted and complex magmatic evolution, crustal inheritance, and/or partial re-setting of radiogenic clocks. However, apart from natural/geological issues, methodological and analytical problems may also contribute to blurring geochronological data. This may be exemplified by the Variscan Karkonosze Pluton (SW Poland). High-precision chemical abrasion (CA) ID-TIMS zircon data indicate that the two main rock types, porphyritic and equigranular, of this igneous body were both emplaced at ca. 312 Ma, while field evidence points to a younger age for the latter. This is in contrast to the earlier reported SIMS (SHRIMP) zircon dates that scattered mainly between ca. 322 and 302 Ma. In an attempt to overcome this dispersion, at least in part caused by radiogenic lead loss, the CA technique was used before SHRIMP analysis. The ²⁰⁶Pb/²³⁸U age obtained in this way from a sample of porphyritic granite is 322 ± 3 Ma, ~16 Ma older than the untreated zircons; another porphyritic sample yielded a mean age of 319 ± 3 Ma, and the mean age was 318 ± 4 Ma for an equigranular granite sample – all three somewhat older than the age obtained by ID-TIMS. Older SIMS dates of ca. 318–322 Ma might indicate either faint inheritance or that zircon domains crystallized during earlier stages of Karkonosze igneous evolution. The ID-TIMS results have been used to re-assess the whole-rock Rb–Sr data. Excluding a porphyritic granite with excess radiogenic ⁸⁷Sr, it appears that isotopic homogeneity was achieved for most samples during the 312 Ma event, as shown by a pooled 21-point isochron with an age of 311 ± 3 Ma and an initial ⁸⁶Sr/⁸⁶Sr of 0.7067 ± 4. Local crustal contamination by stopping of metapelitic material might account for the more radiogenic Sr isotope signature observed in biotite-rich schlieren. A critical re-evaluation of all available SHRIMP data using the ID-TIMS age of 312 Ma as a benchmark suggests that the observed scatter may be partly attributed to analytical and methodological problems, in particular failing to distinguish subtly discordant spots from truly concordant ones, which is a serious limitation of the microbeam analytical approach. Other likely pitfalls contributing to geochronological scatter are identified in the published Re–Os ages on molybdenite and the ⁴⁰Ar/³⁹Ar data on micas. A scenario postulating a 15–20 milliion year evolution of the Karkonosze Pluton cannot be established on the basis of available geochronological data, which rather supports a brief igneous event, although a more protracted pre-emplacement evolution is possible. A short timescale for crystallization of large igneous bodies, as suggested by the ID-TIMS data from the Karkonosze Granite, is in line with models of transport of granitic magmas through dikes to form large plutons.     

Progressive,episodic deformation in the Mexican Fold–Thrust Belt (central Mexico): evidence from isotopic dating of folds and faults

We used illite Ar/Ar dating to obtain absolute ages of folds and shear zones formed within the Mexican Fold–Thrust Belt (MFTB). The methodology takes advantage of illite dating in folded, clay-bearing layers and the ability to obtain accurate ages from small-size fractions of illite using encapsulated Ar analysis. We applied our approach to a cross-section that involves folded Aptian–Cenomanian shale-bentonitic layers interbedded with carbonates of the Zimapán (ZB) and Tampico–Misantla (TMB) Cretaceous basins in central-eastern Mexico. Basinal carbonates were buried by syn-tectonic turbidites and inverted during the formation of the MFTB in the Late Cretaceous. Results from folds and shear zones record different pulses of deformation within this thin-skinned orogenic wedge.

Mineralogical compositions, variations in illite polytypes, illite crystallite size (CS), and Ar/Ar ages were obtained from several size fractions in limbs and hinges of the folds and in the shear zones. 1M_d-illite polytype (with CS of 6–9 nm) dominates in two folds in the TMB while 2M₁-illlite (with CS of 14–30 nm) dominates in the third fold, in the ZB, and in the fold/shear zone. From west (higher grade) to east (lower grade): Ar retention ages indicate shearing occurred at ~84 Ma in the westernmost shear zone, folding at ~82 Ma in the ZB with subsequent localized shearing at ~77 Ma, and Ar total gas ages constrain the time of folding at ~64 Ma on the west side of the TMB and ~44 Ma on the eastern edge. These results are consistent with the age and distribution of syn-tectonic turbidites and indicate episodic progression of deformation from west to east.     

Oligocene–Miocene geodynamic evolution of the central part of Urumieh-Dokhtar Arc of Iran

Basic volcanic rocks from Tafresh, west Kashan, and west Nain volcanic successions in the central part of Urumieh-Dokhtar Magmatic Assemblage (UDMA) of Iran yield K–Ar ages ranging from 26.8 to 18.2 Ma. These ages indicate significant Late Oligocene–Early Miocene basic volcanism in the UDMA. These ages, combined with K–Ar ages of 26.0 and 14.1 Ma, respectively, for associated low-silica and high-silica adakites, help constrain reconstructions of the UDMA geodynamic evolution. Late Oligocene–Early Miocene slab roll-back associated with an asthenospheric mantle influx are suggested as the major processes responsible for concurrent volcanism showing Nb–Ta-depleted, Nb–Ta-enriched and low-silica adakite signatures. Slab roll-back, the likely consequence of a decrease in subduction velocity, led to partial melting of the subducted slab and produced Early–Middle Miocene high-silica (dacitic) adakites. Oligocene to Miocene volcanic rocks do not conform to the Oligocene continental collisional model for the UDMA, rather they suggest a decrease in the subduction rate that prompted the asthenospheric mantle influx.     

典型隧道断面形态对拟建滇西大理至瑞丽铁路 高黎贡山隧道稳定性影响的数值模拟分析

拟建的大理—瑞丽铁路穿越横断山区南段的滇西南地区,地形、地貌和地质条件都极为复杂,其中高黎贡山深埋超长隧道的工程稳定性问题一直是困扰铁路选线、设计和施工的重大工程难点。针对不同形态隧道断面可能对高黎贡山超长隧道工程稳定性产生的影响问题,在充分综合该区野外地质调查、地应力测量、岩石力学实验等成果与资料的基础上,利用ANSYS有限元应力分析软件对不同形态隧道断面的应力分布进行了有限元计算,给出了应力分布图像,分析了断面应力分布的特点和断面形状对应力分布的影响。同时,计算了围岩的应力与强度比,对2种不同断面的围岩稳定性进行了分析对比,最后根据分析结果对不同隧道断面形态下的隧道稳定性进行了综合评价,并据此提出了铁路隧道断面设计与施工中应重点关注的问题。     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

二重钢管防屈曲耗能支撑的有限元分析

设计了12组41种二重钢管防屈曲耗能支撑试件,应用有限元软件ABAQUS对其进行了有限元分析,研究了径厚比、约束比、长细比和边界条件对支撑性能的影响.分析结果表明:二重钢管防屈曲耗能支撑的径厚比取值不大于24、约束比的取值不小于3时,支撑的滞回曲线饱满、稳定;当支撑的长细比增加时,可以通过提高约束比来避免支撑发生整体屈...     

Pseudo-dynamic test and numerical simulation of high-strength concrete frame structure reinforced with high-strength rebars

This paper describes an investigation of a high-strength concrete frame reinforced with high-strength rebars that was tested in the structure engineering laboratory at Shenyang Jianzhu University. The frame specimen was pseudo- dynamically loaded to indicate three earthquake ground motions of different hazard levels, after which the test specimen was subjected to a pseudo-static loading. This paper focuses on the design, construction and experiment of the test frame and validation of the simulation models. Research shows that a high-strength concrete frame reinforced with high-strength rebars is more eff icient and economical than a traditional reinforced concrete frame structure. In addition to the economies achieved by effective use of materials, research shows that the frame can provide enough strength to exceed conventional reinforced concrete frames and provide acceptable ductility. The test study provides evidence to validate the performance of a high- strength concrete frame designed according to current seismic code provisions. Based on previous test research, a nonlinear FEM analysis is completed by using OpenSees software. The dynamic responses of the frame structure are numerically analyzed. The results of the numerical simulation show that the model can calculate the seismic responses of the frame by OpenSees. At the same time, the test provides additional opportunities to validate the performance of the simulation models.     

泥岩渗流-应力耦合蠕变损伤模型研究(Ⅱ):数值仿真和参数反演

进一步分析了第Ⅰ部分[1]提出的泥岩渗流-应力耦合蠕变损伤模型。在连续损伤力学理论和比奥(Biot)理论的基础上,导出了考虑渗流-应力-损伤耦合的蠕变损伤有限元格式,建立了弹性预测、塑性修正、损伤修正-渗透系数修正的数值分析框架,编制了非线性有限元分析程序。根据监测的衬砌长期变形数据,采用优化反分析法获得了蠕变损伤模型中的待定参数,并应用于比利时核废料库施工过程中泥岩巷道围岩渗流-应力耦合过程、损伤演化以及长期稳定性分析,研究结果表明,泥岩开挖后渗透性明显增大,约为原岩的120倍,蠕变效应导致泥岩裂隙和渗透性自愈合,约3.5年后渗透性基本恢复到原岩的数量级,围岩中部的蠕变明显大于顶部和底部。研究成果对软岩隧洞长期稳定性的预测与预报具有一定的参考意义。     

A Kinematic Thermal Model for Descending Slabs with Velocity Boundary Layers: A Case Study for the Tonga Subducting Slab

For the purpose of investigating the influence of metastable olivine (MO) phase transformations on both deep seismicity and stagnation of slabs, we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50o dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers, which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800oC and 1200oC isotherm, respectively. The elevation of olivine transformation near ~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive, zonal, negative density anomalies, respectively. These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab: (i) slab pull force exerts above a depth of ~230 km; (ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km; (iii) as the subducting materials accumulated over 660 km, deepest earthquakes occur due to MO transformation; (iv) a flattened ‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.