Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.     

Spatial patch structure and adaptive strategy for desert shrub of Reaumuria soongorica in arid ecosystem of the Heihe River Basin

In many arid ecosystems, vegetation frequently occurs in high-cover patches interspersed in a matrix of low plant cover. However, theoretical explanations for shrub patch pattern dynamics along climate gradients remain unclear on a large scale. This context aimed to assess the variance of the Reaumuria soongorica patch structure along the precipitation gradient and the factors that affect patch structure formation in the middle and lower Heihe River Basin (HRB). Field investigations on vegetation patterns and heterogeneity in soil properties were conducted during 2014 and 2015. The results showed that patch height, size and plant-to-patch distance were smaller in high precipitation habitats than in low precipitation sites. Climate, soil and vegetation explained 82.5% of the variance in patch structure. Spatially, R. soongorica shifted from a clumped to a random pattern on the landscape towards the MAP gradient, and heterogeneity in the surface soil properties (the ratio of biological soil crust (BSC) to bare gravels (BG)) determined the R. soongorica population distribution pattern in the middle and lower HRB. A conceptual model, which integrated water availability and plant facilitation and competition effects, was revealed that R. soongorica changed from a flexible water use strategy in high precipitation regions to a consistent water use strategy in low precipitation areas. Our study provides a comprehensive quantification of the variance in shrub patch structure along a precipitation gradient and may improve our understanding of vegetation pattern dynamics in the Gobi Desert under future climate change.

     

Melting of the subcontinental lithospheric mantle by the Emeishan mantle plume; evidence from the basal alkaline basalts in Dongchuan, Yunnan, Southwestern China

The Emeishan continental flood basalt (ECFB) sequence in Dongchuan, SW China comprises a basal tephrite unit overlain by an upper tholeiitic basalt unit. The upper basalts have high TiO₂ contents (3.2–5.2 wt.%), relatively high rare-earth element (REE) concentrations (40 to 60 ppm La, 12.5 to 16.5 ppm Sm, and 3 to 4 ppm Yb), moderate Zr/Nb and Nb/La ratios (9.3–10.2 and 0.6–0.9, respectively) and relatively high _{Nd (t)} values, ranging from − 0.94 to 2.3, and are comparable to the high-Ti ECFB elsewhere. The tephrites have relatively high P₂O₅ (1.3–2.0 wt.%), low REE concentrations (e.g., 17 to 23 ppm La, 4 to 5.3 ppm Sm, and 2 to 3 ppm Yb), high Nb/La (2.0–3.9) ratios, low Zr/Nb ratios (2.3–4.2), and extremely low _{Nd (t)} values (mostly ranging from − 10.6 to − 11.1). The distinct compositional differences between the tephrites and the overlying tholeiitic basalts cannot be explained by either fractional crystallization or crustal contamination of a common parental magma. The tholeiitic basalts formed by partial melting of the Emeishan plume head at a depth where garnet was stable, perhaps > 80 km. We propose that the tephrites were derived from magmas formed when the base of the previously metasomatized, volatile-mineral bearing subcontinental lithospheric mantle was heated by the upwelling mantle plume.     

He–Ar and Nd–Sr isotopic compositions of late Pleistocene felsic plutonic back arc basin rocks from Ulleungdo volcanic island, South Korea: Implications for the genesis of young plutonic rocks in a back arc basin

We report analyses of noble gases and Nd–Sr isotopes in mineral separates and whole rocks of late Pleistocene (< 0.2 Ma) monzonites from Ulleungdo, South Korea, a volcanic island within the back arc basin of the Japan island arc. A Rb–Sr mineral isochron age for the monzonites is 0.12 ± 0.01 Ma. K–Ar biotite ages from the same samples gave relatively concordant ages of 0.19 ± 0.01and 0.22 ± 0.01 Ma. ⁴⁰Ar/³⁹Ar yields a similar age of 0.29 ± 0.09 Ma. Geochemical characteristics of the felsic plutonic rocks, which are silica oversaturated alkali felsic rocks (av., 12.5 wt% in K₂O + Na₂O), are similar to those of 30 alkali volcanics from Ulleungdo in terms of concentrations of major, trace and REE elements. The initial Nd–Sr isotopic ratios of the monzonites (⁸⁷Sr/⁸⁶Sr = 0.70454–0.71264, ¹⁴³Nd/¹⁴⁴Nd = 0.512528–0.512577) are comparable with those of the alkali volcanics (⁸⁷Sr/⁸⁶Sr = 0.70466–0.70892, ¹⁴³Nd/¹⁴⁴Nd = 0.512521–0.512615) erupted in Stage 3 of Ulleungdo volcanism (0.24–0.47 Ma). The high initial ⁸⁷Sr/⁸⁶Sr values of the monzonites imply that seawater and crustally contaminated pre-existing trachytes may have been melted or assimilated during differentiation of the alkali basaltic magma.A mantle helium component (³He/⁴He ratio of up to 6.5 RA) associated with excess argon was found in the monzonites. Feldspar and biotite have preferentially lost helium during slow cooling at depth and/or during their transportation to the surface in a hot host magma. The source magma noble gas isotopic features are well preserved in fluid inclusions in hornblende, and indicate that the magma may be directly derived from subcontinental lithospheric mantle metasomatized by an ancient subduction process, or may have formed as a mixture of MORB-like mantle and crustal components. The radiometric ages, geochemical and Nd–Sr isotopic signatures of the Ulleungdo monzonites as well as the presence of mantle-derived helium and argon, suggests that these felsic plutonic rocks evolved from alkali basaltic magma that formed by partial melting of subcontinental lithospheric mantle beneath the back arc basin located along the active continental margin of the southeastern part of the Eurasian plate.     

Evolution of subcontinental lithospheric mantle beneath eastern China: Re–Os isotopic evidence from mantle xenoliths in Paleozoic kimberlites and Mesozoic basalts

The ages of subcontinental lithospheric mantle beneath the North China and South China cratons are less well-constrained than the overlying crust. We report Re–Os isotope systematics of mantle xenoliths entrained in Paleozoic kimberlites and Mesozoic basalts from eastern China. Peridotite xenoliths from the Fuxian and Mengyin Paleozoic diamondiferous kimberlites in the North China Craton give Archean Re depletion ages of 2.6–3.2 Ga and melt depletion ages of 2.9–3.4 Ga. No obvious differences in Re and Os abundances, Os isotopic ratios and model ages are observed between spinel-facies and garnet-facies peridotites from both kimberlite localities. The Re–Os isotopic data, together with the PGE concentrations, demonstrate that beneath the Archean continental crust of the eastern North China Craton, Archean lithospheric mantle of spinel- to diamond-facies existed without apparent compositional stratification during the Paleozoic. The Mesozoic and Cenozoic basalt-borne peridotite and pyroxenite xenoliths, on the other hand, show geochemical features indicating metasomatic enrichment, along with a large range of the Re–Os isotopic model ages from Proterozoic to Phanerozoic. These features indicate that lithospheric transformation or refertilization through melt-peridotite interaction could be the primary mechanism for compositional changes during the Phanerozoic, rather than delamination or thermal-mechanical erosion, despite the potential of these latter processes to play an important role for the loss of garnet-facies mantle. A fresh garnet lherzolite xenolith from the Yangtze Block has a Re depletion age of ∼1.04 Ga, much younger than overlying Archean crustal rocks but the same Re depletion ages as spinel lherzolite xenoliths from adjacent Mesozoic basalts, indicating Neoproterozoic resetting of the Re–Os system in the South China Craton.     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

In situ strength measurements on natural upper-mantle minerals

Using in situ strength measurements at pressures up to 10 GPa and at room temperature, 400, 600, and 700°C, we examined rheological properties of olivine, orthopyroxene, and chromian-spinel contained in a mantle-derived xenolith. Mineral strengths were estimated using widths of X-ray diffraction peaks as a function of pressure, temperature, and time. Differential stresses of all minerals increase with increasing pressure, but they decrease with increasing temperature because of elastic strain on compression and stress relaxation during heating. During compression at room temperature, all minerals deform plastically at differential stress of 4–6 GPa. During subsequent heating, thermally induced yielding is observed in olivine at 600°C. Neither orthopyroxene nor spinel shows complete stress relaxation, but both retain some stress even at 700°C. The strength of the minerals decreases in the order of chromian-spinel ≈ orthopyroxene > olivine for these conditions. This order of strength is consistent with the residual pressure of fluid inclusions in mantle xenoliths.     

青藏高原东南木里地区二叠纪苦橄岩的Os-Sr-Nd同位素地球化学研究

本文报道了在青藏高原东南木里地区发现的二叠纪苦橄岩和与其共生玄武岩的主微量元素地球化学特征以及Os-Sr-Nd同位素组成。苦橄岩和与其共生玄武岩受地壳混染作用影响较小。根据苦橄岩的Ti/Y比值和初始的Os同位素组成,将木里苦橄岩分为两类:高Ti/Y型苦橄岩和低Ti/Y型苦橄岩,其中高Ti/Y型苦橄岩具有高的γ_(Os)= 5.3～ 10.7和ε_(Nd)= 5.9～ 6.4,与全球典型洋岛玄武岩的Os和Nd同位素组成接近,代表了地幔柱源区的同位素特征;而低Ti/Y型苦橄岩具有低的γ_(Os)=-4.1～ 1.2和ε_(Nd)= 3.2～ 5.0,可能表明受到了SCLM(大陆岩石圈地幔)源区物质的混染。与其共生的玄武岩具有低的γ_(Os)=-3.5～-1.6和ε_(Nd)=-0.6～ 0.7,表明其来自于不同于低Ti/Y型苦橄岩也有异于高Ti/Y型苦橄岩的地幔源区,但是也可能受到了SCLM物质的混染。基于Nd-Os同位素的地幔柱与SCLM的二端元混合模型显示:低Ti/Y型苦橄岩可能是SCLM物质组分与地幔柱起源的苦橄质原始岩浆混合形成的;与苦橄岩共生的玄武岩可能是由地幔柱来源的玄武质岩浆与SCLM小比例熔融的熔体混合形成的。