Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Fast genetic mapping in barley: case studies of cuticle mutants using RNA-sequencing

Barley(Hordeum vulgare L.) is one of the earliest domesticated crop species and ranked as the fourth largest cereal production worldwide. Forward genetic studies in barley have greatly advanced plant genetics during the last century; however, most genes are identified by the conventional mapping method. Array genotyping and exome-capture sequencing have also been successfully used to target the causal mutation in barley populations, but these techniques are not widely adopted because of associated costs and partly due to the huge genome size of barley. This review summarizes three mapping cases of barley cuticle mutants in our laboratory with the help of RNA-sequencing. The causal mutations have been successfully identified for two of them and the target genes are located in the pericentromeric regions. Detailed information on the mapping-by-sequencing, mapping-and-sequencing, and RNA-sequencing assisted linkage mapping are presented and some limitations and challenges on the mapping assisted by RNA sequencing are also discussed. The alternative and elegant methods presented in this review may greatly accelerate forward genetics of barley mapping, especially for laboratories without large funding.     

Metamorphic phase relations in orthopyroxene-bearing granitoids: implication for high-pressure metamorphism and prograde melting in the continental crust

In this work, the factors controlling the formation and preservation of high-pressure mineral assemblages in the metamorphosed orthopyroxene-bearing metagranitoids of the Sandmata Complex, Aravalli-Delhi Mobile Belt (ADMB), northwestern India have been modelled. The rocks range in composition from farsundite through quartz mangerite to opdalite, and with varying K₂O, Ca/(Ca + Na)_rock and FeO_tot + MgO contents. A two stage metamorphic evolution has been recorded in these rocks.
An early hydration event stabilized biotite with or without epidote at the expense of magmatic orthopyroxene and plagioclase. Subsequent high-pressure granulite facies metamorphism (∼15 kbar, ∼800 °C) of these hydrated rocks produced two rock types with contrasting mineralogy and textures. In the non-migmatitic metagranitoids, spectacular garnet ± K-feldspar ± quartz corona was formed around reacting biotite, plagioclase, quartz and/or pyroxene. In contrast, biotite ± epidote melting produced migmatites, containing porphyroblastic garnet incongruent solids and leucosomes.
Applying NCKFMASHTO T–M (H₂O) and P–T pseudosection modelling techniques, it is demonstrated that the differential response of these magmatic rocks to high-pressure metamorphism is primarily controlled by the scale of initial hydration. Rocks, which were pervasively hydrated, produced garnetiferous migmatites, while for limited hydration, the same metamorphism formed sub-solidus garnet-bearing coronae. Based on the sequence of mineral assemblage evolution and the mineral compositional zoning features in the two metagranitoids, a clockwise metamorphic P–T path is constrained for the high-pressure metamorphic event. The finding has major implications in formulating geodynamic model of crustal amalgamation in the ADMB.     

Effect of metamorphism and partial melting of host rocks on zircons

Abstract Zircons have been studied from different layers of migmatites (from Arvika, western Sweden and Nelaug, southern Norway) and from a paragneiss (from Arvika) associated with one of the migmatites. The main purpose of the investigation is to establish whether or not information about zircons can help in the elucidation of the parentage and rock-forming processes of migmatites.
The elongation ratio of zircons from all layers is small and characteristic of sedimentary zircons. Further, the absence of characteristic colours and the growth trends of the zircons (indicated by the reduced major axes) observed in the various samples both support a sedimentary parentage for these rocks. The zircons of all layers exhibit secondary growth (overgrowth, outgrowth and multiple growth) due to metamorphism. Compared with the zircons from the paragneiss, those of the migmatite layers are more clouded and less rounded, some of them becoming opaque or even skeletal; this is especially true of the zircons from the leucosomes. These observations indicate an alteration of the original sedimentary zircons in the migmatite, especially in the leucosomes, in response to the migmatization process, previously interpreted as partial melting.     

Petrogenesis of Adakitic Intrusive Rocksin Xu-Huai Area, Eastern China

1IntroductionAdakite has specially sense in the resuming structure background of rock formed,because of its distinctive geochemistry feature and implying partic-ular petrogenesis and ore genesis mechanism (De-fant and Drummond,1990),as a result,it arouses numerous researchers widely concern. After adakite conception was drawn by Chinese (Wang etal.,2000), it was widely concerned by domestic re-searchers (Wang etal.,2000,2001a,2001b;Xu etal.,2000;Pan etal.,2001;Qian,2001;Xu etal.,2001;X…     

The Expansion of the Distance Between Two Arbitrary Gravitating Points as Applied to Specific Models of Celestial Bodies

For the complete system of biharmonic functions a suitable representation in spheroidal coordinates is found. It is used for expanding the distance between two arbitrary gravitating points and its inverse quantity in appropriate series. Such expansions can be of interest and useful in a number of problems of celestial mechanics and stellar dynamics.     

Application of stress release models to historical earthquakes from North China

Univariate and multivariate stress release models are fitted to historical earthquake data from North China. It is shown that a better fit is obtained by treating separately the Eastern part of the region, including the North China Plain and Bohai Sea, and the Western part of the region, including the Ordos Plateau and its Eastern boundary. Further improvement is obtained by fitting the large events (M7.6) and smaller events in the Western region by different stress release models. The comparisons are made by computing the likelihoods of the fitted models and discounting the number of parameters used by Akaike's AIC criterion. The models are used to develop long-term risk scenarios for the East and West regions.     

Trace element distribution among rock-forming minerals from metamorphosed to partially molten basic igneous rocks in a contact aureole (Fuerteventura, Canaries)

Low pressure partial melting of basanitic and ankaramitic dykes gave rise to unusual, zebra-like migmatites, in the contact aureole of a layered pyroxenite–gabbro intrusion, in the root zone of an ocean island (Basal Complex, Fuerteventura, Canary Islands). These migmatites are characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations. Their mineralogy consists of plagioclase (An_32–36), diopside, biotite, oxides (magnetite, ilmenite), +/− amphibole, dominated by plagioclase in the leucosome and diopside in the melanosome. The melanosome is almost completely recrystallised, with the preservation of large, relict igneous diopside phenocrysts in dyke centres. Comparison of whole-rock and mineral major- and trace-element data allowed us to assess the redistribution of elements between different mineral phases and generations during contact metamorphism and partial melting.

Dykes within and outside the thermal aureole behaved like closed chemical systems. Nevertheless, Zr, Hf, Y and REEs were internally redistributed, as deduced by comparing the trace element contents of the various diopside generations. Neocrystallised diopside – in the melanosome, leucosome and as epitaxial phenocryst rims – from the migmatite zone, are all enriched in Zr, Hf, Y and REEs compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of enriched primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post-solidus reequilibration on the cooling of the migmatite terrain i.e. mineral–melt equilibria were reset to mineral–mineral equilibria.     

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.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.