Subduction-related mafic to felsic magmatism in the Malayer–Boroujerd plutonic complex,western Iran

The Malayer–Boroujerd plutonic complex (MBPC) in western Iran, consists of a portion of a magmatic arc built by the northeast verging subduction of the Neo-Tethys plate beneath the Central Iranian Microcontinent (CIMC). Middle Jurassic-aged felsic magmatic activity in MBPC is manifested by I-type and S-type granites. The mafic rocks include gabbroic intrusions and dykes and intermediate rocks are dioritic dykes and minor intrusions, as well as mafic microgranular enclaves (MMEs). MBPC Jurassic-aged rocks exhibit arc-like geochemical signatures, as they are LILE- and LREE-enriched and HFSE- and HREE-depleted and display negative Nb–Ta anomalies. The gabbro dykes and intrusions originated from metasomatically enriched garnet-spinel lherzolite [Degree of melting (f_mel) ~ 15%] and exhibit negative Nd and positive to slightly negative ε_Hf(T) (+ 3.0 to ? 1.6). The data reveal that evolution of Middle Jurassic magmatism occurred in two stages: (1) deep mantle-crust interplay zone and (2) the shallow level upper crustal magma chamber. The geochemical and isotopic data, as well as trace element modeling, indicate the parent magma for the MBPC S-type granites are products of upper crustal greywacke (f_mel: 0.2), while I-type granites formed by partial melting of amphibolitic lower crust (f_mel: 0.25) and mixing with upper crustal greywacke melt in a shallow level magma chamber [Degree of mixing (f_mix): 0.3]. Mixing between andesitic melt leaving behind a refractory dense cumulates during partial crystallization of mantle-derived magma and lower crustal partial melt most likely produced MMEs (f_mix: 0.2). However, enriched and moderately variable ε_Nd(T) (? 3.21 to ? 4.33) and high (⁸⁷Sr/⁸⁶Sr)_i (0.7085–0.7092) in dioritic intrusions indicate that these magmas are likely experienced assimilation of upper crustal materials. The interpretations of magmatic activity in the MBPC is consistent with the role considered for mantle-derived magma as heat and mass supplier for initiation and evolution of magmatism in continental arc setting, elsewhere.     

Flooding a landscape: impact of Holocene transgression on coastal sedimentology and underwater archaeology in Kiladha Bay (Greece)

Franchthi Cave, bordering Kiladha Bay, in Greece, is a key archaeological site, due to its long occupation time, from?~?40,000 to?~?5000 year BP. To date, no clear evidence of Neolithic human dwellings in the cave was found, supporting the assumption that Neolithic people may have built a village where there is now Kiladha Bay. During the Neolithic period/Early Holocene, wide areas of the bay were indeed emerged above sea level. Bathymetric and seismic data identified a terrace incised by a valley in?~?1 to 2 m sediment depth. Eight sediment cores, up to 6.3-m-long, were retrieved and analysed using petrophysical, sedimentological, geochemical, and chronostratigraphic methods. The longest core extends into the exposure surface, consisting of a layer of carbonate rubble in a finer matrix, representing weathering processes. Dated organic remains place this unit at?~?8500 cal year BP. It is overlain by stiff silty mud representing an estuarine environment. This mud is capped by reduced sediments with roots marking an exposure surface. A shell-layer, dated to?~?6300 cal year BP, overlies this terrestrial sequence, reflecting the marine transgression. This layer occurs at 10.8 mbsl, 7.7 m deeper than the global sea level at that time, suggesting tectonic subsidence in the area. It is overlain by finer-grained marine carbonate-rich sediments. The top of the core shows traces of eutrophication, pebbles and marine shells, all likely a result of modern anthropogenic processes. These results are interpreted in the context of human occupation: the exposed surface contains pottery sherds, one dating to the Early to Middle Neolithic period, indicating that Neolithic people were present in this dynamic landscape interacting with a migrating coastline. Even if the artefacts are isolated, future investigations of the submerged landscape off Franchthi Cave might lead to the discovery of a Neolithic village, which eventually became buried under marine sediments.     

High-pressure phase behavior of MnTiO<Subscript>3</Subscript>: decomposition of perovskite into MnO and MnTi<Subscript>2</Subscript>O<Subscript>5</Subscript>

The phase relations and compression behavior of MnTiO₃ perovskite were examined using a laser-heated diamond-anvil cell, X-ray diffraction, and analytical transmission electron microscopy. The results show that MnTiO₃ perovskite becomes unstable and decomposes into MnO and orthorhombic MnTi₂O₅ phases at above 38 GPa and high temperature. This is the first example of ABO₃ perovskite decomposing into AO + AB₂O₅ phases at high pressure. The compression behavior of volume, axes, and the tilting angle of TiO₆ octahedron of MnTiO₃ perovskite are consistent with those of other A²⁺B⁴⁺O₃ perovskites, although no such decomposition was observed in other perovskites. FeTiO₃ is also known to decompose into two phases, instead of transforming into the CaIrO₃-type post-perovskite phase and we argue that one of the reasons for the peculiar behavior of titanate is the weak covalency of the Ti–O chemical bonds.     

Role of Crystal Interlocking on the Strength of Brittle Rocks

Results of an experimental programme on heterogeneous rock-like specimens of dental plaster confirm the pronounced role of tensile microcracks on brittle failure. Microbuckling of very small rock-columns formed amid closely located tensile cracks was observed as the key incident connecting stable phenomenon of tensile cracking to unstable phenomenon of shearing and subsequent macroscopic failure. Using the classical beam and buckling theories and considering geometry of the problem a new failure criterion is proposed. As a novel attempt, this new failure criterion relates the compressive strength of rock to three basic microstructural properties, i.e. degree of crystal interlocking, average Young modulus and average tensile strength of rock forming minerals.     

Hydraulic Conductivity Tests for Evaluating Compatibility of Lateritic Soil—Fly Ash Mixtures with Municipal Waste Leachate

Long term competent performance of liner systems is a critical issue in the design and construction of waste repositories due to adverse interactions associated with leachate generated by wastes. This study was conducted to verify the efficacy of fly ash stabilization in enhancing compatibility between lateritic soil and municipal waste leachate. Applications investigated include soil mixtures containing 0, 5, 10, 15, and 20% fly ash compacted at approximately 2% wet of optimum moisture content with modified proctor energy. Baseline hydraulic conductivity was first established at every level of fly ash content by permeating soil mixtures with tap water before permeation with leachate in a compaction mould permeameter using the falling head test method. Results show that the trend in hydraulic conductivity of specimen containing 0% fly ash was characterized by a gradual but erratic decrease which may suggests partial entry of the leachate cations into the double layer. Conversely, specimens containing fly ash showed a general trend consisting of an initial drop in k (up to an order of magnitude) that was followed by slight decrease sustained until k stabilized and later terminated. Above 10% fly ash content, the relatively high values of k observed was not connected with the reactivity of the soil mixtures with leachate, rather it may be attributed to excessive fly ash content that altered their textural and hydraulic properties. The result of this study is potentially significant in the assessment of fly ash as a compatibility enhancing agent which can be admixed in barrier materials that are susceptible to adverse reactions with the liquid to be contained.     

Helium and argon isotope geochemistry of natural gases in China’s petroliferous basins

The age-accumulation effect of 40Ar in hydrocarbon source rocks was discussed in accordance with the decay law of radioactive elements. In terms of the mean values of 40Ar/36Ar, the old Sinian gas reservoirs (mean values of 40Ar/36Ar: 7009) were definitely distinguished from the Permian gas reservoirs (mean values of 40Ar/36Ar: 1017) in Weiyuan, Sichuan Province, and the gas source of the Permian gas reservoir (mean values of 40Ar/36Ar: 5222) in well Wei-7 with the Weiyuan structure is defined as the Sinian system. Based on the values of 40Ar/36Ar, the coal-type gases (The source rocks are of the C-P system; mean values of 40Ar/36Ar: 1125) are definitely distinguished from the oil-type gases (The source rocks are of the Tertiary system; mean values of 40Ar/36Ar: 590) in the Tertiary reservoirs of the Zhongyuan Oilfield. Besides, 40Ar/36Ar values also have a positive effect on the oil-source correlation of oil reservoirs in ancient hidden mountains. According to the crust-mantle interchange information reflected by 3He/4He values, petroliferous provinces in China can be divided into three major tectonic regions. (1) The eastern active region: The crust-mantle volatile matter exchanges actively, and the 3He/4He values are mainly around 10-6, partly around 10-7. (2) The central stable region: The 3He/4He values are all around 10-8. (3) The western sub-stable region: The 3He/4He values are mainly around 10-8, and around 10-7 on the edges of the basins. Helium contents of some gas wells in China’s eastern petroliferous region reach the industrial abundance (He≈0.05%–0.1%), the 3He/4He values reach 10-6, and the equivalent values for the mantle-source components in helium gas can reach 30%–50%. As viewed from this, a new type of crust-mantle composite helium resources has been proposed. Geneses of some CO2 gas reservoirs in the east of China and some issues concerning mantle-source methane were discussed in the light of the helium and carbon isotopes of CO2 and CH4 in natural gases. In the discussion on helium isotopic characteristics of inclusions in the reservoirs, it was discovered that the 3He/4He values are close to those in natural gases. That is to say, this phenomenon is related to regional tectonism. The 3He/4He, CO2/3He and CH4/3He data were used to discuss the tectonic activities of fault zones in a certain number of regions in China.     

Zircon U-Pb ages and Hf isotopic characteristics of the Dehe biotite monzonitic gneiss pluton in the North Qinling orogen and their geological significance

The Dehe granitic pluton intruded the Xiahe Group which is in the core complex of the North Qinling Orogenic Belt(NQOB).It shows gneissic bedding and possesses typical S-type granite minerals such as muscovite and garnet.LA-ICP-MS U-Pb isotopic dating of the Dehe granite yielded a weighted average age of 925±23 Ma which represents the emplacement age of the pluton.Most of the εHf(t) values are negative,and the two-stage model ages are consistent with the age of the Qinling Group.The isotope data show that the Dehe granite was formed in the following geological setting:in the syn-collision setting of the NQOB in the Neoproterozoic,crustal thickening induced partial melting of materials derived from the Qinling complex,and then the maga upwelled and intruded into the Xiahe Group.The formation of the Dehe S-type granite implied the occurrence of a convergent event in the QOB during the Neoproterozoic.     

Beyond dual-porosity modeling for the simulation of complex flow mechanisms in shale reservoirs

The state of the art of modeling fluid flow in shale reservoirs is dominated by dual-porosity models which divide the reservoirs into matrix blocks that significantly contribute to fluid storage and fracture networks which principally control flow capacity. However, recent extensive microscopic studies reveal that there exist massive micro- and nano-pore systems in shale matrices. Because of this, the actual flow mechanisms in shale reservoirs are considerably more complex than can be simulated by the conventional dual-porosity models and Darcy’s law. Therefore, a model capturing multiple pore scales and flow can provide a better understanding of the complex flow mechanisms occurring in these reservoirs. This paper presents a micro-scale multiple-porosity model for fluid flow in shale reservoirs by capturing the dynamics occurring in three porosity systems: inorganic matter, organic matter (mainly kerogen), and natural fractures. Inorganic and organic portions of shale matrix are treated as sub-blocks with different attributes, such as wettability and pore structures. In kerogen, gas desorption and diffusion are the dominant physics. Since the flow regimes are sensitive to pore size, the effects of nano-pores and micro-pores in kerogen are incorporated into the simulator. The multiple-porosity model is built upon a unique tool for simulating general multiple-porosity systems in which several porosity systems may be tied to each other through arbitrary connectivities. This new model allows us to better understand complex flow mechanisms and eventually is extended into the reservoir scale through upscaling techniques. Sensitivity studies on the contributions of the different flow mechanisms and kerogen properties give some insight as to their importance. Results also include a comparison of the conventional dual-porosity treatment and show that significant differences in fluid distributions and dynamics are obtained with the improved multiple-porosity simulation.     

A semi-implicit method for incompressible three-phase flow in porous media

In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken.