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.     

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.     

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.     

Re-Os同位素对峨眉山大火成岩省成因制约的探讨

峨眉山大火成岩省（ELIP）主要由玄武岩、玄武质火山碎屑岩及少量的苦橄岩（包括越南的科马提岩）、长英质岩石以及层状岩体和岩墙组成,其物质来源直接关系到其成因是否与地幔柱活动有关。Re-Os同位素体系是地核、地幔和地壳物质的最佳示踪剂。前人对ELIP内的Re-Os同位素研究表明,低Ti玄武岩的Os含量为0．006&#215;10^-9-0．40010^-9,^187Os／^188Os初始值为0．1371～1．403,并提出其与地幔柱活动有关;而高Ti玄武岩的Os含量为0．00410^-9～0．56010^-9,^187Os／^188Os初始值为0．1271～5．19,认为起源于大陆岩石圈地幔或地幔柱上升过程中受到大量岩石圈地幔“混染”（xu JF et al．,2007）;科马提岩的0s含量为1．2410^-9～7．0010^-9,^187Os／^188Os初始值为0．1251～0．1261,苦橄岩的Os含量为0．3210^-9～2．32910^-9,^187Os／^188Os初始值为0．1233～0．1266,指示苦橄岩和科马提岩均来自亏损地幔源区（Hanski et al．,2004;陈雷等,2007）。本文利用Os含量最低、^187Os／^188Os最高的高Ti玄武岩作为地壳端员,用铁质陨石、原始上地幔（PUM）和亏损地幔（DMM）作为地核和各种地幔端员,分别做二元混合计算,结果显示绝大多数玄武岩和所有苦橄岩及科马提岩均落在地壳和DMM混合曲线附近,并且邻区特提斯洋地幔岩与DMM具有相近的Os含量和^187Os／^188Os组成,据此推测峨眉山火成岩的形成与特提斯洋的活动有关,主要受控于地壳和亏损地幔的相互作用。     

Seasonal use of coastal resources by otters: Comparing sandy and rocky stretches

Seasonal patterns of use of coastal resources by otters have been studied, simultaneously performing diet analyses and spraint counts (expressed as marking intensity, MI). In spite of the small size of the study area, two different sectors could be clearly identified, a sandy and a rocky one. The origin (freshwater or coastal) of most prey types could also be determined. Surveys were performed bimonthly and classified as autumn–winter (A–W) or spring–summer (S–S).In the sandy sector, eels, flatfish and grey mullets were more frequently consumed in A–W, while crayfish was more frequently consumed in S–S. None of the identified prey types in the rocky sectors showed different frequency of occurrence values between the two periods. In this sector otter fed mainly on coastal prey throughout the year, while in the sandy one there was a clear shift from coastal prey to freshwater ones in S–S. The decrease of coastal prey consumption in S–S in the sandy sector coincided with the lowest MI values, showing that in this period otters use coastal areas less intensely. Though shifts in otter diet and use of space in the sandy sector are probably influenced by the population dynamics of estuarine fish species and crayfish, the increased tourist disturbance during spring and summer in the sandy sector could be an important determinant of the observed patterns. Finally, the results show that MI can be used as a reliable measure of habitat use by otters, at least in small areas.     

Rare species, habitat diversity and functional redundancy in marine benthos

Macro-ecological theories relating species richness, abundance, range size, biological traits and environmental tolerance have rarely been tested in marine soft-sediments, despite the spatial extent of these habitats and the inherent richness of resident communities. This study examines the contribution of rare species to marine soft-sediment communities from New Zealand, focussing on the relationships of range size with abundance, environment, habitat diversity and life history traits. 54% of the 351 species sampled exhibited restricted ranges (found at ≤ 2 sites). In contrast to many terrestrial systems, we observed only a weak positive relationship between abundance and frequency of occurrence. Restricted-range species were not randomly distributed, with their distribution related to habitat characteristics, suggesting an important link between habitat diversity and rarity. They exhibited a similar range of traits to the total observed species pool, suggesting that they are not only important to biodiversity but could play a role in stability. Restricted range species were generally not small and this, together with the number of different biological traits represented, suggests that rare species are important to the functioning of marine systems. Thus, our results highlight the importance of considering rare species in habitat-based approaches to conservation.