Geochronological and petrological constraints from the evolution in the Saxon Granulite Massif,Germany, on the Variscan continental collision orogeny

Controversy over the plate tectonic affinity and evolution of the Saxon granulites in a two‐ or multi‐plate setting during inter‐ or intracontinental collision makes the Saxon Granulite Massif a key area for the understanding of the Palaeozoic Variscan orogeny. The massif is a large dome structure in which tectonic slivers of metapelite and metaophiolite units occur along a shear zone separating a diapir‐like body of high‐P granulite below from low‐P metasedimentary rocks above. Each of the upper structural units records a different metamorphic evolution until its assembly with the exhuming granulite body. New age and petrologic data suggest that the metaophiolites developed from early Cambrian protoliths during high‐P amphibolite facies metamorphism in the mid‐ to late‐Devonian and thermal overprinting by the exhuming hot granulite body in the early Carboniferous. A correlation of new Ar–Ar biotite ages with published P–T–t data for the granulites implies that exhumation and cooling of the granulite body occurred at average rates of ~8 mm/year and ~80°C/Ma, with a drop in exhumation rate from ~20 to ~2.5 mm/year and a slight rise in cooling rate between early and late stages of exhumation. A time lag of c. 2 Ma between cooling through the closure temperatures for argon diffusion in hornblende and biotite indicates a cooling rate of 90°C/Ma when all units had assembled into the massif. A two‐plate model of the Variscan orogeny in which the above evolution is related to a short‐lived intra‐Gondwana subduction zone conflicts with the oceanic affinity of the metaophiolites and the timescale of c. 50 Ma for the metamorphism. Alternative models focusing on the internal Variscan belt assume distinctly different material paths through the lower or upper crust for strikingly similar granulite massifs. An earlier proposed model of bilateral subduction below the internal Variscan belt may solve this problem.     

Carboniferous eclogite and garnet–omphacite granulite from northeastern Hainan Island,South China: Implications for the evolution of the eastern Palaeo‐Tethys

High‐P (HP) eclogite and associated garnet–omphacite granulite have recently been discovered in the Mulantou area, northeastern Hainan Island, South China. These rocks consist mainly of garnet, omphacite, hornblende, quartz and rutile/ilmenite, with or without zoisite and plagioclase. Textural relationships, mineral compositions and thermobarometric calculations demonstrate that the eclogite and garnet–omphacite granulite share the same three‐stage metamorphic evolution, with prograde, peak and retrograde P?T conditions of 620–680°C and 8.7–11.1 kbar, 820–860°C and 17.0–18.2 kbar, and 700–730°C and 7.1–8.5 kbar respectively. Sensitive high‐resolution ion microprobe U–Pb zircon dating, coupled with the identification of mineral inclusions in zircon, reveals the formation of mafic protoliths before 355 Ma, prograde metamorphism at c. 340–330 Ma, peak to retrograde metamorphism at c. 310–300 Ma, and subsequent pegmatite intrusion at 295 Ma. Trace element geochemistry shows that most of the rocks have a MORB affinity, with initial ε_Nd values of +2.4 to +6.7. As with similar transitional eclogite–HP granulite facies rocks in the thickened root in the European Variscan orogen, the occurrence of relatively high P?T metamorphic rocks of oceanic origin in northeastern Hainan Island suggests Carboniferous oceanic subduction leading to collision of the Hainan continental block, or at least part of it, with the South China Block in the eastern Palaeo‐Tethyan tectonic domain.     

Generation of continuous surface soil moisture dataset using combined optical and thermal infrared images

Surface soil moisture (SSM) is a critical variable for understanding water and energy flux between the atmosphere and the Earth's surface. An easy to apply algorithm for deriving SSM time series that primarily uses temporal parameters derived from simulated and in situ datasets has recently been reported. This algorithm must be assessed for different biophysical and atmospheric conditions by using actual geostationary satellite images. In this study, two currently available coarse‐scale SSM datasets (microwave and reanalysis product) and aggregated in situ SSM measurements were implemented to calibrate the time‐invariable coefficients of the SSM retrieval algorithm for conditions in which conventional observations are rare. These coefficients were subsequently used to obtain SSM time series directly from Meteosat Second Generation (MSG) images over the study area of a well‐organized soil moisture network named REMEDHUS in Spain. The results show a high degree of consistency between the estimated and actual SSM time series values when using the three SSM dataset‐calibrated time‐invariable coefficients to retrieve SSM, with coefficients of determination (R²) varying from 0.304 to 0.534 and root mean square errors ranging from 0.020 m³/m³ to 0.029 m³/m³. Further evaluation with different land use types results in acceptable debiased root mean square errors between 0.021 m³/m³ and 0.048 m³/m³ when comparing the estimated MSG pixel‐scale SSM with in situ measurements. These results indicate that the investigated method is practical for deriving time‐invariable coefficients when using publicly accessed coarse‐scale SSM datasets, which is beneficial for generating continuous SSM dataset at the MSG pixel scale.     

Concentration–discharge relationships in the coal mined region of the New River basin and Indian Fork sub‐basin,Tennessee, USA

For many basins, identifying changes to water quality over time and understanding current hydrologic processes are hindered by fragmented and discontinuous water‐quality and hydrology data. In the coal mined region of the New River basin and Indian Fork sub‐basin, muted and pronounced changes, respectively, to concentration–discharge (C–Q) relationships were identified using linear regression on log‐transformed historical (1970s–1980s) and recent (2000s) water‐quality and streamflow data. Changes to C–Q relationships were related to coal mining histories and shifts in land use. Hysteresis plots of individual storms from 2007 (New River) and the fall of 2009 (Indian Fork) were used to understand current hydrologic processes in the basins. In the New River, storm magnitude was found to be closely related to the reversal of loop rotation in hysteresis plots; a peak‐flow threshold of 25 cubic meters per second (m³/s) segregates hysteresis patterns into clockwise and counterclockwise rotational groups. Small storms with peak flow less than 25 m³/s often resulted in dilution of constituent concentrations in headwater tributaries like Indian Fork and concentration of constituents downstream in the mainstem of the New River. Conceptual two or three component mixing models for the basins were used to infer the influence of water derived from spoil material on water quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Slope stability analysis using evolutionary optimization techniques

Slope stability optimization, in the presence of a band of a weak layer between two strong layers, is accounted for in complicated geotechnical problems. Classical optimization algorithms are not suitable for solving such problems as they need a proper preliminary solution to converge to a valid result. Therefore, it is necessary to find a proper algorithm which is capable of finding the best global solution. Recently a lot of metaheuristic algorithms have been proposed which are able to evade local minima effectively. In this study four evolutionary algorithms, including well‐known and recent ones, such as genetic algorithm, differential evolution, evolutionary strategy and biogeography‐based optimization (BBO), are applied in slope stability analysis and their efficiencies are explored by three benchmark case studies. Result show BBO is the most efficient among these evolutionary algorithms and other proposed algorithms applied to this problem. Copyright © 2016 John Wiley & Sons, Ltd.     

New P–T constraints on the Tamayen glaucophane‐bearing rocks,eastern Taiwan: Perple_X modelling results and geodynamic implications

New pseudosection modelling was applied to better constrain the P–T conditions and evolution of glaucophane‐bearing rocks in the Tamayen block of the Yuli belt, recognized as the world's youngest known blueschist complex. Based on the predominant clinoamphibole, textural relationships, and mineral compositions, these glaucophane‐bearing high‐P rocks can be divided into four types. We focused on the three containing garnet. The chief phase assemblages are (in decreasing mode): amphibole + quartz + epidote + garnet + chlorite + rutile/titanite (Type‐I), phengite + amphibole + quartz + garnet + chlorite + epidote + titanite + biotite + magnetite (Type‐II), and amphibole + quartz + albite + epidote + garnet + rutile + hematite + titanite (Type‐III). Amphibole exhibits compositional zoning from core to rim as follows: glaucophane → pargasitic amphibole → actinolite (Type‐I), barroisite → Mg‐katophorite/taramite → Fe‐glaucophane (Type‐II), glaucophane → winchite (Type‐III). Using petrographic data, mineral compositions and Perple_X modelling (pseudosections and superimposed isopleths), peak P–T conditions were determined as 13 ± 1 kbar and 550 ± 40 °C for Type‐I, 10.5 ± 0.5 kbar and 560 ± 30 °C for Type‐II (thermal peak) and 11 ± 1 kbar and 530 ± 30 °C for Type‐III. The calculations yield higher pressures and temperatures than previously thought; the difference is ~1–6 kbar and 50–200 °C. The three rock types record similar P–T retrograde paths with clockwise trajectories; all rocks followed trajectories with substantial pressure decrease under near‐isothermal conditions (Type‐I and Type‐III), with the probable exception of Type‐II where decompression followed colder geotherms. The P–T paths suggest a tectonic environment in which the rocks were exhumed from maximum depths of ~45 km within a subduction channel along a relative cold geothermal gradient of ~11–14 °C km^?1.     

A bounding surface elasto‐viscoplastic constitutive model for non‐isothermal cyclic analysis of asphaltic materials

An elasto‐viscoplastic constitutive model for asphaltic materials is presented within the context of bounding surface plasticity theory, taking into account the effects of the stress state, void binder degree of saturation, temperature and strain rate on the material behaviour. A stress state dependent non‐linear elasticity model is introduced to represent time‐independent recoverable portion of the deformation. The consistent visco‐plasticity framework is utilised to capture the rate‐dependent, non‐recoverable strain components. The material parameters introduced in the model are identified, and their determination from conventional laboratory tests is discussed. The capability of the model to reproduce experimentally observed response of asphaltic materials is demonstrated through numerical simulations of several laboratory test data from the literature. Copyright © 2016 John Wiley & Sons, Ltd.     

Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Modelling the flow of self‐compacting concrete

A Lagrangian particle‐based method, smooth particle hydrodynamics (SPH), is used in this paper to model the flow of self‐compacting concretes (SCC) with or without short steel fibres. An incompressible SPH method is presented to simulate the flow of such non‐Newtonian fluids whose behaviour is described by a Bingham‐type model, in which the kink in the shear stress vs shear strain rate diagram is first appropriately smoothed out. The viscosity of the SCC is predicted from the measured viscosity of the paste using micromechanical models in which the second phase aggregates are treated as rigid spheres and the short steel fibres as slender rigid bodies. The basic equations solved in the SPH are the incompressible mass conservation and Navier–Stokes equations. The solution procedure uses prediction–correction fractional steps with the temporal velocity field integrated forward in time without enforcing incompressibility in the prediction step. The resulting temporal velocity field is then implicitly projected on to a divergence‐free space to satisfy incompressibility through a pressure Poisson equation derived from an approximate pressure projection. The results of the numerical simulation are benchmarked against actual slump tests carried out in the laboratory. The numerical results are in excellent agreement with test results, thus demonstrating the capability of SPH and a proper rheological model to predict SCC flow and mould‐filling behaviour. Copyright © 2010 John Wiley & Sons, Ltd.